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

Timestamp:

02/27/13 12:25:30 (12 years ago)

Author:

Jens Buss

Message:

bootstapping>: bin values extracted from bincenter

File:

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

fact/tools/rootmacros/PulseTemplates/templateextractors.C

@@ -121 +121 @@
 //----------------------------------------------------------------------------
 
-Double_t MedianOfH1 (
-        TH1*            inputHisto
+double MedianOfH1 ( //THE PROBLEM MUST BE HERE!!!!!!!
+        TH1D*            inputHisto
         )
 {
    //compute the median for 1-d histogram h1
    Int_t nbins = inputHisto->GetXaxis()->GetNbins();
-   Double_t *x = new Double_t[nbins];
-   Double_t *y = new Double_t[nbins];
+   double *x = new double[nbins];
+   double *y = new double[nbins];
    for (Int_t i=0;i<nbins;i++) {
-      x[i] = inputHisto->GetXaxis()->GetBinCenter(i+1);
-      y[i] = inputHisto->GetBinContent(i+1);
+      x[i] = inputHisto->GetXaxis()->GetBinCenter(i);
+      y[i] = inputHisto->GetBinContent(i);
    }
-   Double_t median = TMath::Median(nbins,x,y);
+   double median = TMath::Median(nbins,x,y);
    delete [] x;
    delete [] y;
@@ -189 +189 @@
     double  Median[numIt];
     double  Max[numIt];
-    double  parameter[3];
-    double  parameterErr[3];
+    double  parameter[3]    = {0,0,0};
+    double  parameterErr[3] = {0,0,0};
+    for (int i = 0; i < numIt; i++)
+    {
+        Mean[i]     = 0;
+        Median[i]   = 0;
+        Max[i]      = 0;
+    }
 
 //    TCanvas test_canvas;
@@ -196 +202 @@
 //    test_canvas.cd(1);
 //    inputHisto->Draw();
+
     for (int i = 0; i < numIt; i++)
     {
 //        test_canvas.cd(2);
-        TH1* tempHisto = (TH1*)inputHisto->Clone("tempHisto");
+        TH1D* tempHisto = (TH1D*)inputHisto->Clone("tempHisto");
 //                    "",
 //                    "",
@@ -207 +214 @@
 //                    );
         tempHisto->Reset();
+        sample.SetSeed(sample.GetSeed()+1);
         sample.BootstrapTH1(inputHisto, tempHisto);
 
         Mean[i]     = tempHisto->GetMean();
         Median[i]   = MedianOfH1 ( tempHisto );
-
         Max[i]      = tempHisto->GetBinCenter( tempHisto->GetMaximumBin() );
 
-        //improved determination of maximum
+        //improved determination of modus
 //        TF1 gaus("fgaus", "gaus", Max[i]-10, Max[i]+10);
 //        tempHisto->Fit("fgaus", "WWQRN0");
@@ -540 +547 @@
     TH1F*   hOutputMeanHisto        = NULL;
     TH1F*   hOutputMedianHisto      = NULL;
-    TH1D*    hTempHisto              = NULL;
+    TH1D*    hTempHisto             = NULL;
     double  max_prop                = 0;
     double  median                  = 0;
@@ -600 +607 @@
     int slice_max  = hInputHisto->GetXaxis()->GetLast();
 
+    if (verbosityLevel > 2)
+    {
+        cout << "\t...looping over slice range " << slice_min
+             << " to " << slice_max << endl;
+    }
+
     for (Int_t slice=slice_min;slice<=slice_max;slice++)
     {
 
-        hTempHisto  = hInputHisto->ProjectionY("",slice,slice);
+        hTempHisto  = hInputHisto->ProjectionY("",slice,slice,"o");
 
         //containers for errors
-        double slMean[3];
-        double slError[3];
+        double slMean[3]    ={0,0,0};
+        double slError[3]   ={0,0,0};
 
         //calculate Errors with bootstrapping