Changeset 12275

Timestamp:

10/26/11 10:53:20 (13 years ago)

Author:

kraehenb

Message:

Added more comments. Questions please to tpk@phys.ethz.ch.

File:

• fact/tools/rootmacros/gainanalysis.C (modified) (7 diffs)

fact/tools/rootmacros/gainanalysis.C

@@ -21 +21 @@
 {
 //******************************************************************************
-//Read a datafile and plot the DRS-calibrated data
-//ATTENTION: only works for ROI=1024
+//Read a file with G-APD singles, analyse them and plot the spectra
+//This source code is best read with a wide editor (>100 characters per line) or a horizontal scroll bar.
+//ATTENTION: this script only works for ROI=1024
 // (array indices of the calibration wrong otherwise)
 //******************************************************************************
@@ -92 +93 @@
         Double_t par_init[5+2*n_peaks] = {70,90,100,15,10,15,10,30,30};
         Double_t parread_tmp[5+2*n_peaks];
-        Double_t parread[data_px][5+2*n_peaks];
 //      TF1 *f_peaks = new TF1("total","[2]*TMath::Exp(-(x-[0])*(x-[0])/(2*[3]*[3])) + [4]*TMath::Exp(-(x-[0]-[1])*(x-[0]-[1])/(2*[5]*[5])) + [6]*TMath::Exp(-(x-[0]-2*[1])*(x-[0]-2*[1])/(2*[7]*[7])) + [8]*TMath::Exp(-(x-[0]-3*[1])*(x-[0]-3*[1])/(2*[9]*[9])) + [10]*TMath::Exp(-(x-[12])*(x-[12])/(2*[11]*[11]))",fit_start,fit_end);
 //      Double_t par_init[3+2*n_peaks] = {70,90,500000,15,50000,15,5000,15,500,15,800000,20,0};
@@ -120 +120 @@
         
 //-------------------------------------------
-//Oscilloscope
+//Oscilloscope (search for peaks and fill integral into the histograms)
 //-------------------------------------------
 //      float threshold = pedestal_mean+pedestal_rms;
@@ -174 +174 @@
 //      canv_spect->Update();
 //***********************
+        
+//-------------------------------------------
+//Fit the histograms and store the parameters into parread
+//parread_dist: peak difference (gain) per pixel, parread_x: x-Values for the TGraph
+//-------------------------------------------
+        Double_t parread[data_px][5+2*n_peaks];
         Double_t parread_dist[data_px];
         Double_t parread_x[data_px];
@@ -196 +202 @@
         }
         
+//-------------------------------------------
+//Plot the gain of all pixel
+//-------------------------------------------
         canv_gain->cd();
         TGraph* plot_gain = new TGraph(data_px,parread_x,parread_dist);
@@ -207 +216 @@
         
 //-------------------------------------------
-//Draw the data where the gain is out of limits
+//Draw the spectrum of pixel where the gain is out of limits
+//Drawn are the fitted function (f_peaks) and its components (f_peak[i], f_pedestal)
 //-------------------------------------------
         char temp = 'x';
@@ -238 +248 @@
                 }
         }
-//-------------------------------------------
-//Draw the data
+        
+//-------------------------------------------
+//Give the possibility to draw the spectrum of any pixel
 //-------------------------------------------
         canv_spect->cd();