Subversion Repositories f9daq

Rev

Blame | Last modification | View Log | RSS feed

  1. #include "TROOT.h"
  2. #include "TFile.h"
  3. #include "TBenchmark.h"
  4. #include "TH1F.h"
  5. #include "TH2F.h"
  6. #include "TCanvas.h"
  7. #include "TStyle.h"
  8. #include "TPad.h"
  9. #include "TF1.h"
  10. #include "TGraph.h"
  11. #include "TSpectrum.h"
  12.  
  13. #include "RTUtil.h"
  14.  
  15. #define TDC_BIN (25./1000.) //1 TDC bin in ns
  16. #define MIKRO_BIN 0.49609/1000. //1 mikro step in mm
  17.  
  18. #define HFILL_COLOR 18
  19.  
  20. void mppc(char *fname, char *plopt="t", int nChannels = 1, double fitw=1.0, char *fitf="g")
  21. {
  22.         char fullname[256];
  23.        
  24.         //get ROOT file with histograms
  25.         char fnameroot[1024];
  26.         TFile * rootfile;
  27.         TDirectory *dir;
  28.        
  29.         sprintf(fnameroot, "root/%s.root", fname);
  30.         rootfile = (TFile *) gROOT->FindObject(fname);
  31.         if(rootfile==NULL) rootfile = new TFile(fnameroot);
  32.         if(rootfile==NULL) {
  33.           printf("Cannot open root file %s!!!\n",fnameroot);
  34.           return;
  35.         }
  36.         dir = (TDirectory*) rootfile;
  37.        
  38.        
  39.         // set draw style
  40.         gStyle->SetOptStat("ne");
  41.         gStyle->SetPalette(1, 0);
  42.        
  43.         gStyle->SetPaperSize(TStyle::kA4);
  44.         gStyle->SetStatBorderSize(1);
  45.         gStyle->SetFrameBorderMode(0);
  46.         gStyle->SetFrameFillColor(0);
  47.         gStyle->SetCanvasBorderMode(0);
  48.         gStyle->SetPadBorderMode(0);
  49.         gStyle->SetPadColor(0);
  50.         gStyle->SetCanvasColor(0);
  51.         gStyle->SetStatColor(0);
  52.         gStyle->SetOptFit(11);
  53.         gStyle->SetOptStat();
  54.         gStyle->SetPadRightMargin(0.15);
  55.         gStyle->SetPadLeftMargin(0.12);
  56.         //gStyle->SetTitleYOffset(1.4);
  57.        
  58.   TCanvas *c[16];
  59.   int cc=-1;
  60.   char hname[256];
  61.   TH1F *hp1d; TH2F *hp2d;
  62.  
  63.  
  64. // 2d Scan --------------------------------------------------------------------------------------------
  65.  
  66.   if( strchr(plopt, '2') != NULL ) {                   
  67.     c[++cc] = new TCanvas("2D Scan", fname, 650, 400, 600, 600);
  68.                        
  69.   int i=0;
  70.                         gPad->SetLogz();
  71.                                 sprintf(hname, "h2d%d",i); hp2d = (TH2F *) dir->Get(hname);
  72.                                
  73.                                 (hp2d->GetYaxis())->SetLabelOffset(0.01);
  74.                                 hp2d->SetTitle("; X [mm]; Y[mm]");
  75.                                 hp2d->DrawCopy("COLZ");
  76.                
  77.                 sprintf(fullname, "ps/%s_2d.eps", fname);
  78.                 c[cc]->SaveAs(fullname);
  79.   }
  80.  
  81. // 1d Scan Y ------------------------------------------------------------------------------------------
  82.  
  83.   if( strchr(plopt, 'y') != NULL ) {                   
  84.     c[++cc] = new TCanvas("Y Scan", fname, 250, 200, 600, 600);
  85.                        
  86.   int i=0;
  87.                         //gPad->SetLogy();
  88.                                 sprintf(hname, "hnhitsy%d",i); hp1d = (TH1F *) dir->Get(hname);
  89.                                
  90.                                 hp1d->SetTitle("Y Scan; Y; Count");
  91.                                 hp1d->DrawCopy();
  92.                
  93.                 sprintf(fullname, "ps/%s_y.eps", fname);
  94.                 c[cc]->SaveAs(fullname);
  95.   }  
  96.  
  97. // 1d Scan X ------------------------------------------------------------------------------------------
  98.  
  99.   if( strchr(plopt, 'x') != NULL ) {                   
  100.     c[++cc] = new TCanvas("X Scan", fname, 250, 200, 600, 600);
  101.                        
  102.   int i=0;
  103.                         //gPad->SetLogy();
  104.                                 sprintf(hname, "hnhitsx%d",i); hp1d = (TH1F *) dir->Get(hname);
  105.                                
  106.                                 hp1d->SetTitle("X Scan; X; Count");
  107.                                 hp1d->DrawCopy();
  108.                
  109.                 sprintf(fullname, "ps/%s_x.eps", fname);
  110.                 c[cc]->SaveAs(fullname);
  111.   }  
  112.  
  113. // Fits of corrected TDCs -----------------------------------------------------------------------------
  114.  
  115.   TF1 *fg = new TF1("fg", "gaus");
  116.         TF1 *fgg = new TF1("fgg", "gaus(0)+gaus(3)");
  117.         fgg->SetParNames("Constant","Mean","Sigma","Constant2","Mean2","Sigma2");
  118.   if( strchr(plopt, 'f') != NULL ) {
  119.         c[++cc] = new TCanvas("Corrected TDCs", fname, 50*(cc+1), 0, 600, 850);
  120.                 //c[cc]->Divide(2,4);
  121.                 c[cc]->Divide(nChannels);
  122.  
  123.                 for(int i=0;i<nChannels;i++) {
  124.                   c[cc]->cd(i+1); gPad->SetLogy();
  125.                         sprintf(hname, "hctdc%d",i); hp1d = (TH1F *) dir->Get(hname);
  126.                         hp1d->SetTitle("cTDC; TDC [ns]; N");
  127.                         (hp1d->GetXaxis())->SetRangeUser(-40*TDC_BIN,160*TDC_BIN);
  128.                         if( hp1d->GetMaximum() < 50 ) continue;
  129.                         //(hp1d->GetXaxis())->SetRangeUser(-0.25,0.75);
  130.                        
  131.                         if( strcmp(fitf, "g")==0 ) {
  132.                         fg->SetParameters(hp1d->GetMaximum(), hp1d->GetBinCenter(hp1d->GetMaximumBin()), 0.05);
  133.                         hp1d->Fit(fg, "0ql", "", fg->GetParameter(1)-0.1*fitw, fg->GetParameter(1)+0.05*fitw);
  134.                         hp1d->Fit(fg, "ql", "", fg->GetParameter(1)-0.1*fitw, fg->GetParameter(1)+0.05*fitw);
  135.                        
  136.                         printf("Ch[%d] Sigma =%6.1lfps\n",i,fg->GetParameter(2)*1000.);
  137.       }
  138.       else if( strcmp(fitf, "gg")==0 ) {
  139.                         //fgg->SetParameters(2600.,0*TDC_BIN,2*TDC_BIN, 100.,0*TDC_BIN,10*TDC_BIN);
  140.                         fgg->SetParameters(hp1d->GetMaximum(), hp1d->GetBinCenter(hp1d->GetMaximumBin()), 0.03,
  141.                                           hp1d->GetMaximum(), hp1d->GetBinCenter(hp1d->GetMaximumBin())/10., 0.1);
  142.                         hp1d->Fit(fgg, "0ql", "", -50*TDC_BIN, 50*TDC_BIN);
  143.                         hp1d->Fit(fgg, "ql", "", fgg->GetParameter(1)-2.5*fitw*fgg->GetParameter(2),
  144.                                                  fgg->GetParameter(1)+3.5*fitw*fgg->GetParameter(2));
  145.                                                  
  146.                         printf("Ch[%d] Sigma =%6.1lfps\n",i,fgg->GetParameter(2)*1000.);
  147. /*
  148.                   fg->SetRange(-50., 200.);
  149.                         fg->SetParameters(fgg->GetParameter(0), fgg->GetParameter(1), fgg->GetParameter(2));
  150.                         fg->SetLineColor(2);
  151.                         fg->DrawCopy("LSAME");
  152.                         fg->SetParameters(fgg->GetParameter(3), fgg->GetParameter(4), fgg->GetParameter(5));
  153.                         fg->SetLineColor(3);
  154.                         fg->DrawCopy("LSAME");
  155. */
  156.       } else {printf("Wrong fit function (parameter 4)!!!\n"); return;}
  157.                 }      
  158.                 sprintf(fullname, "ps/%s_cTDC.eps", fname);
  159.                 c[cc]->SaveAs(fullname);
  160.   }  
  161.  
  162. // TDCs -----------------------------------------------------------------------------------------------
  163.  
  164.   if( strchr(plopt, 't') != NULL ) {   
  165.         c[++cc] = new TCanvas("Raw TDC", fname, 0, 500, 600, 430);
  166.                        
  167.                 int i=0;
  168.                         gPad->SetLogy();
  169.                                 sprintf(hname, "htdc%d",i); hp1d = (TH1F *) dir->Get(hname);
  170.                                 hp1d->SetTitle("TDC; TDC [ns]; N");
  171.                                 (hp1d->GetXaxis())->SetRangeUser(-10,20);
  172.                                 hp1d->DrawCopy();
  173.                                 //if( hp1d->GetMaximum() < 10 ) continue;
  174.                                
  175.                           fg->SetParameters(hp1d->GetMaximum(), hp1d->GetBinCenter(hp1d->GetMaximumBin()), 0.25);        
  176.                           fg->SetLineColor(kRed);
  177.                           hp1d->Fit(fg,"0QL", "", fg->GetParameter(1)-0.5, fg->GetParameter(1)+0.5);
  178.                           hp1d->Fit(fg,"QL", "", fg->GetParameter(1)-0.8, fg->GetParameter(1)+0.8);
  179.  
  180.                     printf("Ch[%d] RAW Sigma =%6.1lfps\n",i,fg->GetParameter(2)*1000.);
  181.                
  182.                 sprintf(fullname, "ps/%s_TDC.eps", fname);
  183.                 c[cc]->SaveAs(fullname);
  184.   }
  185.  
  186.  
  187.         //rootfile->Close();
  188. }
  189.