Changeset 14063 for fact/tools/rootmacros/PulseTemplates/pixel.C

Timestamp:

06/04/12 18:52:53 (12 years ago)

Author:

Jens Buss

Message:

added function for shifting histograms in y-direction, histogramms for chi2test on pulseshape

File:

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

fact/tools/rootmacros/PulseTemplates/pixel.C

@@ -1 +1 @@
 #include <iostream>
+#include <TMath.h>
+#include <TF1.h>
 
 #include "rootfilehandler.h"
@@ -52 +54 @@
 
     hList           = new TList();
+
+    if (mOptions.Contains("C") )
+    {
+        hChi2EdgetoMax = new TH1F*[mMaxPulseOrder];
+    }
+
     if (options.Contains("L"))
     {
@@ -60 +68 @@
        BookPixelHistos();
     }
+
+    if (mOptions.Contains("C"))
+    {
+        BookDistributionHistos();
+    }
+
     BookTemplateHistos();
     BookEdgeTemplateHistos();
@@ -227 +241 @@
     {
         pixelCanvas[pulse_order]->cd( histoFrameNr[0] );
+        hMaxOverlay[pulse_order]->Draw("COLZ");
+
+        pixelCanvas[pulse_order]->cd( histoFrameNr[1] );
         hPixelMax[pulse_order]->Draw();
 
-        pixelCanvas[pulse_order]->cd( histoFrameNr[1] );
+        pixelCanvas[pulse_order]->cd( histoFrameNr[2] );
         hPixelMedian[pulse_order]->Draw();
 
-        pixelCanvas[pulse_order]->cd( histoFrameNr[2] );
+        pixelCanvas[pulse_order]->cd( histoFrameNr[3] );
         hPixelMean[pulse_order]->Draw();
+//        hPixelMax[pulse_order]->Draw("SAME");
+//        hPixelMedian[pulse_order]->Draw("SAME");
 
     }
@@ -249 +268 @@
     for (int pulse_order = 0; pulse_order < mMaxPulseOrder; pulse_order++)
     {
-        pixelCanvas[pulse_order]->cd( histoFrameNr[0] );
+        pixelCanvas[pulse_order]->cd( histoFrameNr[4] );
+        hEdgeOverlay[pulse_order]->Draw("COLZ");
+
+        pixelCanvas[pulse_order]->cd( histoFrameNr[5] );
         hPixelEdgeMax[pulse_order]->Draw();
 
-        pixelCanvas[pulse_order]->cd( histoFrameNr[1] );
+        pixelCanvas[pulse_order]->cd( histoFrameNr[6] );
         hPixelEdgeMedian[pulse_order]->Draw();
 
-        pixelCanvas[pulse_order]->cd( histoFrameNr[2] );
+        pixelCanvas[pulse_order]->cd( histoFrameNr[7] );
         hPixelEdgeMean[pulse_order]->Draw();
+//        hPixelEdgeMax[pulse_order]->Draw("SAME");
+//        hPixelEdgeMedian[pulse_order]->Draw("SAME");
 
     }
@@ -526 +550 @@
             }
 
+            if (mOptions.Contains("C"))
+            {
+                if (mVerbosityLevel > 3) cout << "\t...booking " << HistoName("hSlopeRisingEdge", order) << endl;
+                hChi2EdgetoMax[order] = new TH1F (
+                            HistoName("hChi2EdgetoMax", order),
+                            HistoTitle("Distribution of CHI2 comparison of Edge and Max", order),
+                            600,
+                            -299,
+                            300
+                            );
+                hChi2EdgetoMax[order]->SetAxisRange(
+                                    -1,
+                                    7,
+                                    "X");
+                hChi2EdgetoMax[order]->GetXaxis()->SetTitle( "p-Value" );
+                hChi2EdgetoMax[order]->GetYaxis()->SetTitle( "counts" );
+                hChi2EdgetoMax[order]->SetStats(mStats);
+//                hList->Add( hChi2EdgetoMax[order] );
+            }
+
             if (mVerbosityLevel > 2) cout << "...done" << endl;
         }
@@ -1061 +1105 @@
 
 
-
+void
+Pixel::ShiftHistoInY(
+        TH1*                histo,
+        float               shift
+        )
+{
+    int min_bin = histo->GetXaxis()->GetFirst();
+    int max_bin = histo->GetXaxis()->GetLast();
+
+    for (int bin = min_bin; bin <= max_bin; bin++)
+    {
+        histo->AddBinContent( bin, shift);
+    }
+
+}
+
+