WebSVN – f9daq – Diff – /sipmScan/examples/sipm.C

 #include "TH1F.h"
 #include "TH2F.h"
 #include "TCanvas.h"
 #include "TStyle.h"
 #include "TPad.h"
+#include "TPaveText.h"
+//#include "TLabel.h"
 #include "TF1.h"
 #include "TGraph.h"
 #include "TSpectrum.h"
 #include "stdio.h"
 #include "include/RTUtil.h"
 double getNoise(TH2F*, int, int);
-int sipm(char filename[256] = "test", char plopt[256]="all", int chXstart=0, int chXend=7, int chYstart=0, int chYend=7, bool debug = false)
+int sipm(char filename[256] = "test", char plopt[256]="all", int parameter1=0, int parameter2=7, int chYstart=0, int chYend=7, bool debug = false)
+{
   const int c_nChannels = 64;
-  const double c_xOffset = 0; // mm
-  const double c_yOffset = 0;
   int map[8][8]={{32,34,53,55,40,42,61,63},
                  {48,50,37,39,56,58,45,47},
                  {33,35,52,54,41,43,60,62},
                  {49,51,36,38,57,59,44,46},
           return(0);
+        }
         // set draw style
         RTSetStyle(gStyle);
+        // Print results to .pdf (option "p")
+        if(strstr(plopt, "p") != NULL) {
+          TCanvas* canvasFirst = new TCanvas("canvasFirst","",500,500);
+          TPaveText* text = new TPaveText(.05,.1,.95,.8);
+          text->AddText("SiPM Scan results");
+          canvasFirst->cd();
+          text->Draw();
+    canvasFirst->Print("result.pdf(");
+  }
   if( strstr(plopt, "all") != NULL ) {
     TCanvas *canvas2 = new TCanvas("canvas2","Hits x;;",2000,2000);
     TCanvas *canvas3 = new TCanvas("canvas3","Hits y;;",2000,2000);
     canvas2->Divide(8,8);
   if( strstr(plopt, "x") != NULL ) {
     TCanvas *canvas10 = new TCanvas("canvas10","Ch x;;",500,500);
     TH1F* h_hitsx;
     canvas10->cd();
+    int chPosition = map[parameter1][parameter2];
     char hname[128];
-    sprintf(hname, "hnhitsx%d", chXstart);
+    sprintf(hname, "hnhitsx%d", chPosition);
     h_hitsx = (TH1F*)rootfile->Get(hname);
     h_hitsx->Draw();
+    if(strstr(plopt, "p") != NULL) {
+      canvas10->Print("result.pdf");
+    }
+  }
   if( strstr(plopt, "y") != NULL ) {
     TCanvas *canvas11 = new TCanvas("canvas11","Ch x;;",500,500);
     TH1F* h_hitsy;
     canvas11->cd();
+    int chPosition = map[parameter1][parameter2];
     char hname[128];
-    sprintf(hname, "hnhitsy%d", chXstart);
+    sprintf(hname, "hnhitsy%d", chPosition);
     h_hitsy = (TH1F*)rootfile->Get(hname);
     h_hitsy->Draw();
+    if(strstr(plopt, "p") != NULL) {
+      canvas11->Print("result.pdf");
+    }
+  }
   if( strstr(plopt, "share") != NULL ) {
    /*TCanvas *canvas4 = new TCanvas("canvas1","canvas1",1000,1000);
    int nChannels = chYend-chYstart+1;
      h_hitsy->Draw();
    }*/
    TCanvas *canvas4 = new TCanvas("canvas4","canvas4",500,500);
    canvas4->cd();
-   for(int i=chXstart; i<=chXend; i++) {
+   for(int i=parameter1; i<=parameter2; i++) {
      TH1F* h_hitsx;
      char hname[128];
      int chPosition = map[i][chYstart];
      sprintf(hname, "hnhitsx%d", chPosition);
      h_hitsx = (TH1F*)rootfile->Get(hname);
      h_hitsx->SetTitle("Scan X;x [mm]; Entries");
      h_hitsx->GetYaxis()->SetTitleOffset(1.3);
      h_hitsx->SetStats(0);
-     if (i == chXstart)
+     if (i == parameter1)
       h_hitsx->Draw();
      else {
       h_hitsx->SetLineColor(i+1);
       h_hitsx->Draw("same");
+     }
    TCanvas *canvas5 = new TCanvas("canvas5","canvas5",500,500);
    canvas5->cd();
    for(int i=chYstart; i<=chYend; i++) {
      TH1F* h_hitsy;
      char hname[128];
-     int chPosition = map[chXstart][i];
+     int chPosition = map[parameter1][i];
      sprintf(hname, "hnhitsy%d", chPosition);
      h_hitsy = (TH1F*)rootfile->Get(hname);
      h_hitsy->SetTitle("Scan Y;y [mm]; Entries");
      h_hitsy->GetYaxis()->SetTitleOffset(1.3);
      h_hitsy->SetStats(0);
     canvas51->cd();
     // Get the filled histogram from scan
     TH1F* h_x;
     char name[32];
-    sprintf(name, "hnhitsx%d", chXstart);
+    sprintf(name, "hnhitsx%d", parameter1);
     h_x = (TH1F*)rootfile->Get(name);
     h_x->DrawCopy();
     TH1F* h_y;
-    sprintf(name, "hnhitsy%d", chXstart);
+    sprintf(name, "hnhitsy%d", parameter1);
     h_y = (TH1F*)rootfile->Get(name);
     // Create and fill corrected histogram
     Int_t binsX = h_x->GetXaxis()->GetNbins();
     if(debug) printf("nBins: %d\n", binsX);
     * Draws also a 2d scan of these channels.
     */
   if (strstr(plopt, "line") != NULL) {
     TCanvas* canvas6 = new TCanvas("canvas6","canvas6",500,500);
     canvas6->cd(0);
     gStyle->SetOptStat(0);
     TH2F* h0 = (TH2F*) rootfile->Get("h2d0");
           Int_t binsX = h0->GetXaxis()->GetNbins();
     Int_t minX  = h0->GetXaxis()->GetFirst();
     Int_t maxX  = h0->GetXaxis()->GetLast()+1;
     Int_t binsY = h0->GetYaxis()->GetNbins();
     Int_t minY  = h0->GetYaxis()->GetFirst();
     TH1F* h_line[8];
     for(int j=0; j<8; j++) {
       h_line[j] = new TH1F("h_line", "h_line", binsX, xLowUser, xUpUser);
+    }
-    for(int j=chXstart; j<=chXend; j++) {
+    for(int j=parameter1; j<=parameter2; j++) {
       int chPosition = map[j][chYstart];
       char hname[128];
       sprintf(hname, "h2d%d", chPosition);
       int histogram = j;
       h[histogram] = (TH2F *) rootfile->Get(hname);
           //signal -= noise;
           //signal /= 5*10000.0;
           double eta = -log(1 - signal);
           double x = xLowUser + k*(xUpUser-xLowUser)/double(binsX);
           //double y = l*(yUpUser-yLowUser)/double(binsY);
-          h_line[j]->Fill(x-c_xOffset, signal);
+          h_line[j]->Fill(x, signal);
         //}
+      }
-      if (j == chXstart) {
+      if (j == parameter1) {
         h_line[j]->SetTitle("SiPM#2 w/o noise subtraction;x[mm];Hits");
         //h_line[j]->GetYaxis()->SetRangeUser(-0.05, 0.3);
         //h_line[j]->GetYaxis()->SetRangeUser(-50, 2500);
         h_line[j]->Draw("");
+      }
     //! 2d scan
     TCanvas* canvas61 = new TCanvas("canvas61","canvas61",8*200,300);
     canvas61->cd();
     TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX,xLowUser,xUpUser, binsY,yLowUser,yUpUser);
-    for(int i=chXstart; i<=chXend; i++) {
+    for(int i=parameter1; i<=parameter2; i++) {
-      //int canvasPosition = nX*(i-chYstart)+(j-chXstart)+1;
+      //int canvasPosition = nX*(i-chYstart)+(j-parameter1)+1;
-      //int canvasPosition = nX*(chYend-i)+chXstart+1;
+      //int canvasPosition = nX*(chYend-i)+parameter1+1;
       //if (debug) printf("canvas %d\n",canvasPosition);
         int chPosition = map[i][chYstart];
         char hname[128];
               sprintf(hname, "h2d%d", chPosition);
               int histogram = i;
             double eta = -log(1 - p);
             //double x = k*(xUpUser-xLowUser)/double(binsX);
             //double y = l*(yUpUser-yLowUser)/double(binsY);
             double x = h[histogram]->GetXaxis()->GetBinCenter(k);
             double y = h[histogram]->GetYaxis()->GetBinCenter(l);
-            h_corrected->Fill(x-c_xOffset, y-c_yOffset, eta);
+            h_corrected->Fill(x, y, eta);
+          }
+        }
+    }
     h_corrected->SetTitle("SiPM#2 n_pe = - ln(P0);x[mm];y[mm]");
     h_corrected->GetZaxis()->SetRangeUser(-0.05,0.30);
     int nX = 3;
     int nY = 3;
     TCanvas* canvas7 = new TCanvas("canvas7","canvas7", nX*400,nY*400);
     printf("nx %d ny %d\n",nX,nY);
     canvas7->Divide(nX,nY);
-    for(int i=chXend-1; i<=chXend+1; i++) {
+    for(int i=parameter2-1; i<=parameter2+1; i++) {
-      for(int j=chXstart-1; j<=chXstart+1; j++) {
+      for(int j=parameter1-1; j<=parameter1+1; j++) {
-      //int canvasPosition = nX*(i-chYstart)+(j-chXstart)+1;
+      //int canvasPosition = nX*(i-chYstart)+(j-parameter1)+1;
-      int canvasPosition = nX*(chXend+1-i) + (j-chXstart+1) +1;
+      int canvasPosition = nX*(parameter2+1-i) + (j-parameter1+1) +1;
       if (debug) printf("canvas %d\n",canvasPosition);
         canvas7->cd(canvasPosition);
         char hname[128];
         int chPosition = map[j][i];
         sprintf(hname, "h2d%d", chPosition);
     // Number of photoelectrons - Poissonian correction
     TCanvas* canvas8 = new TCanvas("canvas8","canvas8", 1000,1000);
     canvas8->cd();
     gStyle->SetOptStat(0);
       char hname[128];
-      int chPosition = map[chXstart][chXend];
+      int chPosition = map[parameter1][parameter2];
       sprintf(hname, "h2d%d", chPosition);
       TH2F* h_2d = (TH2F*)rootfile->Get(hname);
       Int_t binsX = h_2d->GetXaxis()->GetNbins();
       Int_t minX  = h_2d->GetXaxis()->GetFirst();
       Double_t xUpUser   = h_2d->GetXaxis()->GetXmax();
       Double_t yLowUser  = h_2d->GetYaxis()->GetXmin();
       Double_t yUpUser =   h_2d->GetYaxis()->GetXmax();
       if (debug) printf("xLow %f xUp %f\n",xLowUser,xUpUser);
       TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, xLowUser, xUpUser, binsY, yLowUser, yUpUser);
+      double diffX = xUpUser - xLowUser;
+      double diffY = yUpUser - yLowUser;
+      //TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, 3.00, 3.00+diffX, binsY, 8.00, 8.00+diffY);
       //TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX, minX, maxX, binsY, minY, maxY);
       double noise = getNoise(h_2d, 1, 70);
       if(debug) printf("Noise = %f\n", noise);
       for(int k=minX; k<=maxX; k++) {
         for(int j=minY; j<=maxY; j++) {
           double signal = h_2d->GetBinContent(k,j); // detected
           //double p = ((signal - noise) > 1) ? (signal-noise) : 1;
           double p = signal - noise;
+          /*
-          p /= 10000.;
+          p /= 1000.;
           double p0 = 1.0 - p; // events with zero photons
           //double eta = (-log(p0) * p0 ) / (1-p0-0.00001);
           double eta = -log(p0); // constant of the poissonian statistics
           if (debug) printf("p=%f p0=%f log(p0)=%f eta=%f\n",p,p0,log(p0),eta);
+          */
           //double x = xLowUser + k*(xUpUser - xLowUser) / double(binsX);
           double x = h_2d->GetXaxis()->GetBinCenter(k);
           //double y = yLowUser + j*(yUpUser-yLowUser)/double(binsY);
           double y = h_2d->GetYaxis()->GetBinCenter(j);
           if (debug) printf("x=%f y=%f\n", x, y);
-          h_corrected->Fill(x + c_xOffset, y + c_yOffset, eta);
+          h_corrected->Fill(x, y, p);
+        }
+       }
-      h_corrected->SetTitle("n_pe = - ln P(0);x[mm];y[mm]");
+      //h_corrected->SetTitle("n_pe = - ln P(0);x[mm];y[mm]");
-      //h_corrected->SetTitle(";x[mm];y[mm]");
+      h_corrected->SetTitle(";x[mm];y[mm]");
-      gStyle->SetPalette(52,0);
+      //gStyle->SetPalette(52,0); // black and white for print
-      h_corrected->GetZaxis()->SetRangeUser(-0.05,0.7);
+      //h_corrected->GetZaxis()->SetRangeUser(-0.05,1);
-      h_corrected->SetContour(50);
+      //h_corrected->SetContour(100);
+      //gStyle->SetPalette(1);
+      //SetGS();
       h_corrected->Draw("colz");
       // collection efficiency
       int nPoints =0;
       double efficiency=0;
           efficiency += signal;
           nPoints++;
+        }
+      }
     printf("Signal sum = %f\n   # of points = %d\n",efficiency,nPoints);
+    if(strstr(plopt, "p") != NULL) {
+      canvas8->Print("result.pdf");
+    }
+  }
   /** Draws the sum of channel signals
     * Each channel is a 2d ('h2d') histogram
     * Suitable for 8x8 chs scan
     */
   if( strstr(plopt, "sum") != NULL ) {
-    int nX = chXend - chXstart + 1;
+    int nX = parameter2 - parameter1 + 1;
     int nY = chYend - chYstart + 1;
           TCanvas* canvas12 = new TCanvas("canvas12","c2",8*200, 8*200);
           canvas12->cd();
           gStyle->SetOptStat(0);
     TH2F* h_corrected = new TH2F("h_corrected","h_corrected",binsX,xLowUser,xUpUser, binsY,yLowUser,yUpUser);
     TH2F* h[nX*nY];
     // 2d histogram noise subtraction and poisson scaling
     for(int i=chYstart; i<=chYend; i++) {
-      for(int j=chXstart; j<=chXend; j++) {
+      for(int j=parameter1; j<=parameter2; j++) {
         int chPosition = map[j][i];
         char hname[128];
               sprintf(hname, "h2d%d", chPosition);
-              int histogram = nX*(i-chYstart)+(j-chXstart);
+              int histogram = nX*(i-chYstart)+(j-parameter1);
               h[histogram] = (TH2F *) rootfile->Get(hname);
               int noise = getNoise(h[histogram], 1, 170);
               if (debug) printf("noise: %d\n",noise);
               for(int k=minX; k<=maxX; k++) {
           for(int l=minY; l<=maxY; l++) {
             //double eta = (-log(p0) * p0 ) / (1-p0-0.00001);
             double eta = -log(p0);
             //printf("p=%f p0=%f log(p0)=%f eta=%f\n",p,p0,log(p0),eta);
             double x = k*(xUpUser-xLowUser)/double(binsX);
             double y = l*(yUpUser-yLowUser)/double(binsY);
-            h_corrected->Fill(x-c_xOffset, y-c_yOffset, eta);
+            h_corrected->Fill(x, y, eta);
+          }
+        }
+      }
+    }
   if (strstr(plopt, "beam") != NULL) {
   TCanvas* canvas9 = new TCanvas("canvas9","canvas9", 500,500);
   canvas9->cd();
   char hname[128];
-  sprintf(hname, "hnhitsx%d", 36);
+  sprintf(hname, "hnhitsx%d", 49);
   TH1F* h_laser = (TH1F*)rootfile->Get(hname);
   h_laser->Draw();
   h_laser->SetStats(1);
-  TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",17.02,17.30);
+  TF1* err = new TF1("err","[0]+[1]*TMath::Erf((x-[2])/[3])",parameter1,parameter2);
   err->SetParameter(0,2500);
   err->SetParameter(1, h_laser->GetMaximum());
   err->SetParameter(2, h_laser->GetBinCenter(h_laser->GetMaximumBin()));
   err->SetParameter(3, 0.001);
   h_laser->Fit(err,"qr");
       printf("\n");
+    }
     gStyle->SetOptStat(0);
     h_map->Draw("text");
+  }
+  if(strstr(plopt, "p") != NULL) {
+    TCanvas* canvasLast = new TCanvas("canvasLast","",500,500);
+    canvasLast->Print("result.pdf)");
+  }
   return(0);
+}
 /** Function calculates the noise from one channel

Subversion Repositories f9daq

(root)/sipmScan/examples/sipm.C – Rev 40 → 47