Changeset 13558

Timestamp:

05/04/12 16:27:24 (13 years ago)

Author:

Jens Buss

Message:

now histograms will be plotted again

File:

• fact/tools/rootmacros/PulseTemplates/FPulseOverlay.C (modified) (21 diffs)

fact/tools/rootmacros/PulseTemplates/FPulseOverlay.C

@@ -157 +157 @@
 
     //Root-File Objects
-//    TObjArray*  hList[sampleSetSize] = {NULL};
+//    TObjArray*  hList[pixelSetSize] = {NULL};
     TObjArray*  hRootList           = NULL;
 
@@ -175 +175 @@
 //void SaveHistograms( const char*, const char*, TObjArray*, int );
 
-void FillHistograms(Pixel*, vector<Region>*, int, int, TString, int, int);
+void FillHistograms(Pixel*, vector<Region>*, int, int, int , int );
 void DrawTestHistograms( int);
 void ProduceDebugHistos( vector<Region> *pZXings);
-
+bool  UseThisPulse( int, int, int, int, int, int);
 void UpdateCanvases( int, int);
 void DeletePixelCanvases( int, int );
@@ -190 +190 @@
         const char*     datafilename        = "/fact/raw/2011/11/09/20111109_006.fits.gz",
         const char*     drsfilename         = "/fact/raw/2011/11/09/20111109_003.drs.fits.gz",
-        const char*     OutRootFileName     = "test.root",
-        bool            ProduceGraphic      = false,
+        const char*     OutRootFileName     = "/home_nfs/isdc/jbbuss/analysis/FPulseTemplate/test.root",
+        bool            ProduceGraphic      = true,
         bool            spikeDebug          = false,
-        bool            debugPixel          = false,
+        bool            debugPixel          = true,
         bool            testmode            = false,
         bool            saveResults         = false,
         int             verbosityLevel      = 1,        // different verbosity levels can be implemented here
         int             firstevent          = 1,        // Hast to be between 1 and infinity NOT 0!!!!
-        int             nevents             = 100,
-        int             firstpixel          = 0,
-        int             npixel              = -1,
-        int             sampleSetSize       = 40,
-        int             maxPulseOrder       = 3,
+        int             nevents             = 1000,
+        int             firstpixel          = 15,
+        int             npixel              = 3,
+        int             pixelSetSize       = 40,
+        int             maxPulseOrder       = 4,
         int             AmplWindowWidth     = 14,       //Width of Window for selection of pulses to histograms
         TString         histoOptions        = "S",
-        float           GainMean            = 8,
+        float           GainMean            = 9,
         float           BSLMean             = -1,        //4 Histogramms will be drawn, decide how far pulsehights differ from eachother
         int             avg1                = 8,
@@ -275 +275 @@
             cName   ="cgpDistributions";
             cName   += pulse_order;
-            cTitle   ="Pulses of Order: ";
+            cTitle   ="Distributions of Pulses with Order of: ";
             cTitle   += pulse_order;
             cgpDistributions[pulse_order] = new TCanvas(cName,cTitle, 0,pulse_order*20,800,800);
@@ -281 +281 @@
             cName   ="cgpPixelPulses";
             cName   += pulse_order;
-            cTitle   ="Pulses of Order: ";
+            cTitle   ="Overlays of Pulses with Order of: ";
             cTitle   += pulse_order;
             cgpPixelPulses[pulse_order] = new TCanvas(cName,cTitle, 0,pulse_order*20,800,800);
@@ -405 +405 @@
 // Loop over Pixel Sets
 //-------------------------------------
-    for ( int firstPixelOfSample = firstpixel;
-          firstPixelOfSample < firstpixel + npixel;
-          firstPixelOfSample += sampleSetSize )
+    for ( int firstPixelOfSet = firstpixel;
+          firstPixelOfSet < firstpixel + npixel;
+          firstPixelOfSet += pixelSetSize )
     {
 
@@ -413 +413 @@
         {
             cout << "------------------------------------------------" << endl
-                 << "...processing Sample from Pixel "
-                 << firstPixelOfSample
+                 << "...processing Set from Pixel "
+                 << firstPixelOfSet
                  << " to Pixel "
-                 << firstPixelOfSample+sampleSetSize-1 << endl;
+                 << firstPixelOfSet+pixelSetSize-1 << endl;
         }
 
@@ -425 +425 @@
         {
             // Get an Event --> consists of 1440 Pixel ...erm....data
-            cout << endl << "processing event " << ev << endl;
             datafile->GetRow( ev - 1 );
 
@@ -431 +430 @@
             {
              cout << "-------------------------------------" << endl
-                  << "...processing Sample from Pixel "
-                  << firstPixelOfSample
+                  << "...processing Set from Pixel "
+                  << firstPixelOfSet
                   << " to Pixel "
-                  << firstPixelOfSample+sampleSetSize-1
+                  << firstPixelOfSet+pixelSetSize-1
                   << "... Event: " << CurrentEventID
                   << "/" << nevents << endl;
@@ -442 +441 @@
 // Loops over every Pixel of a Set of Pixels
 //--------------------------------------------------------------------
-            for ( int pixelID = firstPixelOfSample;
-                    pixelID < firstPixelOfSample + sampleSetSize
+            for ( int pixelID = firstPixelOfSet;
+                    pixelID < firstPixelOfSet + pixelSetSize
                     && pixelID < firstpixel + npixel;
                     pixelID++ )
@@ -450 +449 @@
                 {
                  cout << "-------------------------------------" << endl
-                      << "...processing Sample from Pixel "
-                      << firstPixelOfSample
+                      << "...processing Set from Pixel "
+                      << firstPixelOfSet
                       << " to Pixel "
-                      << firstPixelOfSample+sampleSetSize-1
+                      << firstPixelOfSet+pixelSetSize-1
                       << "... Event: " << CurrentEventID
                       << "/" << nevents << endl
@@ -550 +549 @@
                             AmplWindowWidth,
                             ev,
-                            histoOptions,
+//                            histoOptions,
                             maxPulseOrder,
                             verbosityLevel
@@ -573 +572 @@
                         if (debugPixel)
                         {
-                            for ( int order = 0;
-                                order < maxPulseOrder;
-                                order++
-                                )
-                            {
-                                pixel[pixelID]->DrawHistograms(
-                                            cgpPixelPulses[order],
+//                            for ( int order = 0;
+//                                order < maxPulseOrder;
+//                                order++
+//                                )
+//                            {
+                                pixel[firstPixelOfSet]->DrawHistograms(
+                                            cgpPixelPulses,
                                             pulsesCanvasFrameNrs
                                             );
-                            }
+//                            }
                             if (testmode)
                             {
@@ -605 +604 @@
                          << endl << "-------------------------------------"<< endl;
                 }
-            }//End of Loop over Sample
-//-------------------------------------
-// end of Loops over Sample
+            }//End of Loop over Set
+//-------------------------------------
+// end of Loops over Set
 //-------------------------------------
 
@@ -621 +620 @@
 //-------------------------------------
 
-        for ( int pixelID = firstPixelOfSample;
-              pixelID < firstPixelOfSample + sampleSetSize
+        for ( int pixelID = firstPixelOfSet;
+              pixelID < firstPixelOfSet + pixelSetSize
               && pixelID < firstpixel + npixel;
               pixelID++ )
         {
             //here is what happends at the end of each loop over all Events
-//            pixel[pixelID]->DrawHistograms(
-//                        cgpPixelPulses,
-//                        pulsesCanvasFrameNrs
-//                        );
-
-//            SaveHistograms(     //save histograms of generell results into output root file
-//                        OutRootFileName,
-//                        CreateSubDirName(pixel[pixelID]->mChid),
-//                        pixel[pixelID]->hList,
-//                        saveResults,
-//                        verbosityLevel
-//                        );
 
             if (ProduceGraphic)
             {
-            UpdateCanvases(
-                           verbosityLevel,
-                           maxPulseOrder
-                           );
+                pixel[pixelID]->DrawHistograms(
+                            cgpPixelPulses,
+                            pulsesCanvasFrameNrs
+                            );
+
+                UpdateCanvases(
+                               verbosityLevel,
+                               maxPulseOrder
+                               );
             }
 
@@ -717 +709 @@
                     verbosityLevel
                     );
-    }
-
-    delete datafile;
-    delete[] pixel;
+        delete      cgpTestHistos;
+
+        if (spikeDebug)
+            delete      cFiltered;
+    }
+
+    delete      datafile;
+    delete[]    pixel;
+
+
     return( 0 );
 }
@@ -813 +811 @@
         int             AmplWindowWidth,
         int             eventNumber,
-        TString         histoOptions,
         int             maxPulseOrder,
         int             verbosityLevel
@@ -819 +816 @@
 {
     if (verbosityLevel > 2) cout << endl << "...filling pulse histograms" ;
+    if (verbosityLevel > 3) cout << endl << "...EventNR " << eventNumber ;
     bool use_this_peak=false;
     int order_of_pulse=0;
@@ -831 +829 @@
             continue;
         }
+
         // Define axis range of Histogram
         int Left        = reg->maxPos - gPixelOverlayXaxisLeft;
@@ -843 +842 @@
 
 
-//        hTesthisto2->Fill( eventNumber, reg->slopeOfRisingEdge ) ;
-
         if (verbosityLevel > 2) cout << endl << "\t...choosing Histogram" ;
-
         //determine order of pulse and which histogram shall be filled
-        for(int order = 0; order < maxPulseOrder; order++)
+        for ( int order = 0; order < maxPulseOrder; order++)
         {
-            if (Ameas[reg->maxPos] >= ((gGainMean*(order+1)) - (AmplWindowWidth/2))
-                && Ameas[reg->maxPos] < ((gGainMean*(order+1)) + AmplWindowWidth/2))
-            {
-                if (verbosityLevel > 2) cout << "...#" << order ;
+            use_this_peak = UseThisPulse(
+                        reg->maxPos,
+                        order,
+                        AmplWindowWidth,
+                        gGainMean,
+                        maxPulseOrder,
+                        verbosityLevel
+                        );
+            if ( use_this_peak)
+            {
                 order_of_pulse = order;
-                use_this_peak = true;
                 break;
             }
-            else if (order >= (maxPulseOrder - 1))
-            {
-                use_this_peak = false;
-                if (verbosityLevel > 2) cout << "...NONE" << endl ;
+        }
+        //Fill histograms
+        if (use_this_peak)
+        {
+            //Histograms for Distributions
+            if (CurrentPixel->mOptions.Contains("S") )
+            {
+                if (verbosityLevel > 2)
+                        cout << endl << "\t...filling Edge Histogram" ;
+                CurrentPixel->hSlopeRisingEdge[order_of_pulse]->Fill( reg->slopeOfRisingEdge ) ;
             }
-            if (use_this_peak)
-            {
-                if ( histoOptions.Contains("S") )
-                {
-//                     if (verbosityLevel > 2)
-                         cout << endl << "\t...filling Edge Histogram" ;
-                    CurrentPixel->hSlopeRisingEdge[order]->Fill( reg->slopeOfRisingEdge ) ;
-                }
+
+            //Histograms for Maximum Overlay
+            for ( int pos = Left; pos < Right; pos++)
+            {
+                CurrentPixel->hMaxOverlay[order_of_pulse]->Fill( pos - (reg->maxPos), Ameas[pos]) ;
+                CurrentPixel->hMaxProfile[order_of_pulse]->Fill( pos - (reg->maxPos), Ameas[pos]) ;
             }
-        }
-        //Fill overlay und profile histograms
-        if (use_this_peak){
-            for ( int pos = Left; pos < Right; pos++){
-//                if ();
-               CurrentPixel->hMaxOverlay[order_of_pulse]->Fill( pos - (reg->maxPos), Ameas[pos]) ;
-               CurrentPixel->hMaxProfile[order_of_pulse]->Fill( pos - (reg->maxPos), Ameas[pos]) ;
-            }
-            for ( int pos = EdgeLeft; pos < EdgeRight; pos++){
+            //Histograms for Edge Overlay
+            for ( int pos = EdgeLeft; pos < EdgeRight; pos++)
+            {
                 CurrentPixel->hEdgeOverlay[order_of_pulse]->Fill( pos - (reg->halfRisingEdgePos), Ameas[pos]) ;
                 CurrentPixel->hEdgeProfile[order_of_pulse]->Fill( pos - (reg->halfRisingEdgePos), Ameas[pos]) ;
-
             }
+
             if (verbosityLevel > 2) cout << "...done" << endl;
         }
+
         //Breakout option
         if (breakout)
-        break;
+        {
+            break;
+        }
    }
     if (verbosityLevel > 2) cout << "...done" << endl;
@@ -895 +897 @@
 // end of FillHistograms
 //----------------------------------------------------------------------------
+
+bool UseThisPulse(
+        int             maxPos,
+        int             order,
+        int             AmplWindowWidth,
+        int             gGainMean,
+        int             maxPulseOrder,
+        int             verbosityLevel
+        )
+{
+    //determine if pulse is of given order
+bool use_this_peak = false;
+        if (Ameas[maxPos] >= ((gGainMean*(order+1)) - (AmplWindowWidth/2))
+            && Ameas[maxPos] < ((gGainMean*(order+1)) + AmplWindowWidth/2))
+        {
+            if (verbosityLevel > 3) cout << "...#" << order ;
+            use_this_peak = true;
+        }
+        else
+        {
+            if (verbosityLevel > 3)
+            {
+                if (order >= (maxPulseOrder - 1))
+                {
+                    cout << "...NONE" << endl ;
+                }
+            }
+            use_this_peak = false;
+        }
+
+return use_this_peak;
+}
+
 
 void