Subversion Repositories f9daq

Rev

Go to most recent revision | Blame | Compare with Previous | Last modification | View Log | RSS feed

  1. //---------------------------------------------------
  2. int plot2d(char *opt="colz", float max=0){
  3. char hname[50];
  4.   char hn[50];
  5.   TCanvas *c;
  6.   TH1I *gICh = (TH1I * ) gDirectory->Get("mapping");  
  7.   for (int i=0;i<4;i++) {
  8.     sprintf(hname,"cz%d",i);
  9.     sprintf(hn,"Corr Pmt %d",i);
  10.     c =new TCanvas(hname,hn, 1000,1000);
  11.     c->Divide(4,4);
  12.     for (int j=0;j<16;j++){
  13.       c->cd(gICh->GetBinContent(j+1)+1);
  14.       sprintf(hn,"corrch%d",i*16+j);
  15.        TH1F*h=    (TH1F * ) gDirectory->Get(hn);
  16.       if (max>0) h->SetMaximum(max);
  17.       h->Draw(opt);
  18.     }
  19.     c->Modified();
  20.     c->Update();
  21.   }
  22. return 0;
  23. }
  24.  
  25. //---------------------------------------------------
  26. int plot( char *fname=NULL, float max=0){
  27.   if (fname!=NULL){
  28.     cout << "Rootname:" << fname << endl;
  29.     delete gROOT->GetListOfFiles()->FindObject(fname); // clear memory of file name
  30.     if( gSystem->AccessPathName(fname) ) {
  31.       cout << endl << "File: " << fname << " not there!!!" << endl << endl;
  32.       return 0;
  33.     }
  34.     TFile *rootfile = new TFile(fname);
  35.   }
  36.  //if (buf)  delete buf;
  37.   char hn[255];
  38.  
  39.   TCanvas *c1;
  40.   for (int i=0;i<4;i++){
  41.     sprintf(hn,"csumadc%d",i);
  42.     c1=new TCanvas(hn,hn,1200,1200);
  43.     c1->Divide(9,9,0,0);
  44.     sprintf(hn,"sumadc%d",i);
  45.     TH2F *h = ((TH2F * )gDirectory->Get(hn));
  46.     for (int j=0;j<81;j++){
  47.       c1->cd(1+j);
  48.       TH1D *h0 = h->ProjectionX(" ",1+j,1+j);
  49.       h0->DrawCopy();
  50.     }
  51.     c1->Modified();
  52.     c1->Update();
  53.   }  
  54.  
  55.  
  56.  
  57.   TCanvas *c2=new TCanvas("c2","c2",800,800);
  58.   c2->Divide(2,2);
  59.   for (int i=0;i<4;i++){
  60.     c2->cd(1+i);
  61.     sprintf(hn,"pmt1%d",i);
  62.      ((TH1F * ) gDirectory->Get(hn))->Draw("colz");
  63.   }
  64.   c2->Modified();
  65.   c2->Update();
  66.  
  67.   TCanvas *c6=new TCanvas("c6","c6",800,800);
  68.   c6->Divide(2,2);
  69.   for (int i=0;i<4;i++){
  70.     c6->cd(1+i);
  71.     sprintf(hn,"pmt3%d",i);
  72.      ((TH1F * ) gDirectory->Get(hn))->Draw("colz");
  73.   }
  74.   c6->Modified();
  75.   c6->Update();
  76.  
  77.   TCanvas *c5=new TCanvas("c5","c5",800,800);
  78.  
  79.  
  80.   ((TH1F * ) gDirectory->Get("globalxy"))->Draw("colz");
  81.  
  82.   c5->Modified();
  83.   c5->Update();
  84.  
  85.  
  86.   TCanvas *c3=new TCanvas("c3","c3",800,800);
  87.   c3->Divide(2,2);
  88.   for (int i=0;i<4;i++){
  89.     c3->cd(1+i);
  90.     sprintf(hn,"pmt2%d",i);
  91.      ((TH1F * ) gDirectory->Get(hn))->Draw("colz");
  92.   }
  93.  
  94.  
  95.   c3->Modified();
  96.   c3->Update();
  97.  
  98.   TCanvas *c4=new TCanvas("c4","c4",800,800);
  99.   c4->Divide(4,4);
  100.   for (int i=0;i<16;i++){
  101.     c4->cd(1+i)->SetLogy();
  102.     sprintf(hn,"ach%d",i);
  103.      ((TH1F * ) gDirectory->Get(hn))->Draw();
  104.   }
  105.   c4->Modified();
  106.   c4->Update();
  107.   plot2d("colz",max);
  108.   return 0;
  109. }
  110.  
  111. //---------------------------------------------------
  112. int plotadc(){
  113.   char hn[255];
  114.   TCanvas *c;
  115.   char hname[50];
  116.   char hn[50];
  117.   TH1I *gICh = (TH1I * ) gDirectory->Get("mapping");
  118.   for (int i=0;i<4;i++) {
  119.     sprintf(hname,"cx%d",i);
  120.     sprintf(hn,"Pmt %d",i);
  121.     c =new TCanvas(hname,hn, 1000,1000);
  122.     c->Divide(4,4);
  123.     for (int j=0;j<16;j++){
  124.       c->cd(gICh->GetBinContent(j+1)+1)->SetLogy();
  125.       sprintf(hn,"ach%d",i*16+j);  
  126.        ((TH1F * ) gDirectory->Get(hn))->Draw();
  127.     }
  128.     c->Modified();
  129.     c->Update();
  130.   }
  131.  
  132.   return 0;
  133. }
  134.  
  135. //---------------------------------------------------
  136. int plotsingles(char *opt="colz", float max=0){
  137.  char hname[50];
  138.  char hn[50];
  139.  TCanvas *c;
  140.   TH1I *gICh = (TH1I * ) gDirectory->Get("mapping");
  141.  for (int i=0;i<4;i++) {
  142.     sprintf(hname,"cy%d",i);
  143.     sprintf(hn,"Single Pmt %d",i);
  144.     c =new TCanvas(hname,hn, 1000,1000);
  145.     c->Divide(4,4);
  146.     for (int j=0;j<16;j++){
  147.        c->cd(gICh->GetBinContent(j+1)+1);
  148.       sprintf(hn,"singlech%d",i*16+j);
  149.       TH1F*h=    (TH1F * ) gDirectory->Get(hn);
  150.       if (max>0) gHistoSingle[i*16+j]->SetMaximum(max);
  151.       h->Draw(opt);
  152.     }
  153.     c->Modified();
  154.     c->Update();
  155.   }
  156.  
  157.   return 0;
  158. }
  159.  
  160.  
  161.  
  162.  
  163. //---------------------------------------------------
  164. int plotcut(char *opt="", float max=0){
  165. char hname[50];
  166.   char hn[50];
  167.   TCanvas *c;
  168.   TH1I *gICh = (TH1I * ) gDirectory->Get("mapping");
  169.   for (int i=0;i<4;i++) {
  170.     sprintf(hname,"cz%d",i);
  171.     sprintf(hn,"Cut Pmt %d",i);
  172.     c =new TCanvas(hname,hn, 1000,1000);
  173.     c->Divide(4,4);
  174.     for (int j=0;j<16;j++){
  175.        c->cd(gICh->GetBinContent(j+1)+1)->SetLogy();
  176.        sprintf(hn,"corrch%d",i);
  177.        TH1F*h=    (TH1F * ) gDirectory->Get(hn);
  178.       if (max>0) h->SetMaximum(max);
  179.       h->Draw(opt);
  180.     }
  181.     c->Modified();
  182.     c->Update();
  183.   }
  184. return 0;
  185. }
  186. //---------------------------------------------------
  187.