Changeset 14760 for fact/tools/rootmacros/PulseTemplates/templateextractors.C

Timestamp:

12/21/12 15:32:48 (12 years ago)

Author:

Jens Buss

Message:

add number of bootstrapp iterations

File:

• fact/tools/rootmacros/PulseTemplates/templateextractors.C (modified) (5 diffs)

fact/tools/rootmacros/PulseTemplates/templateextractors.C

@@ -525 +525 @@
         TString         overlayType,
         int             pulse_order,
+        int             numBootstrapIt,// = 10, //number of iterations for bootstrapping
         int             verbosityLevel
         )
@@ -609 +610 @@
 
         //calculate Errors with bootstrapping
-        CalculateErrorsWithBootstrapping(hTempHisto, 10, slMean, slError);
+        CalculateErrorsWithBootstrapping(hTempHisto, numBootstrapIt, slMean, slError);
 
         if (verbosityLevel > 3)
@@ -617 +618 @@
 
         //Get maximum of slice's distribution
-        max_prop    = hTempHisto->GetBinCenter( hTempHisto->GetMaximumBin() );
-//                max_prop    = slMean[0];
+//        max_prop    = hTempHisto->GetBinCenter( hTempHisto->GetMaximumBin() );
+                max_prop    = slMean[0];
 
         //improve result by< fitting gaussian to slices distribution
-        TF1 gaus("fgaus", "gaus", max_prop-30, max_prop+30);
-        hTempHisto->Fit("fgaus", "QRN0");
-        max_prop        = gaus.GetParameter(1);
+//        TF1 gaus("fgaus", "gaus", max_prop-30, max_prop+30);
+//        hTempHisto->Fit("fgaus", "QRN0");
+//        max_prop        = gaus.GetParameter(1);
 
         //calculate error of max prop
@@ -634 +635 @@
             cout << "\t\t...calculating Median of slice " << slice << endl;
         }
-        median      = MedianOfH1(hTempHisto);
-//        median      = slMean[1];
+//        median      = MedianOfH1(hTempHisto);
+        median      = slMean[1];
         median_err  = slError[1]; //error from bootstraping
 
@@ -641 +642 @@
 
 
-        mean        = hTempHisto->GetMean();
-//        mean        = slMean[2];
-        mean_err    = hTempHisto->GetRMS(); //RMS of slice
-//        mean_err    = slError[2];         //error from bootstraping
+//        mean        = hTempHisto->GetMean();
+        mean        = slMean[2];
+//        mean_err    = hTempHisto->GetRMS()/hTempHisto->GetEntries(); //RMS of slice
+        mean_err    = slError[2];         //error from bootstraping
 
         if (verbosityLevel > 4) cout << "\t\t\t\t MaxProb of Slice " << slice << ": " << max_prop << endl;