Index: fact/tools/rootmacros/PulseTemplates/templateextractors.C
===================================================================
--- fact/tools/rootmacros/PulseTemplates/templateextractors.C	(revision 14755)
+++ fact/tools/rootmacros/PulseTemplates/templateextractors.C	(revision 14760)
@@ -525,4 +525,5 @@
         TString         overlayType,
         int             pulse_order,
+        int             numBootstrapIt,// = 10, //number of iterations for bootstrapping
         int             verbosityLevel
         )
@@ -609,5 +610,5 @@
 
         //calculate Errors with bootstrapping
-        CalculateErrorsWithBootstrapping(hTempHisto, 10, slMean, slError);
+        CalculateErrorsWithBootstrapping(hTempHisto, numBootstrapIt, slMean, slError);
 
         if (verbosityLevel > 3)
@@ -617,11 +618,11 @@
 
         //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,6 +635,6 @@
             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,8 +642,8 @@
 
 
-        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;