Changeset 1999 for trunk/MagicSoft/Mars

Timestamp:

04/24/03 11:09:17 (22 years ago)

Author:

wittek

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• manalysis/MPadSchweizer.cc (modified) (3 diffs)
• manalysis/MPadSchweizer.h (modified) (4 diffs)
• mfilter/MFCT1SelFinal.cc (modified) (1 diff)
• mhistmc/MHMcCT1CollectionArea.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/Changelog

@@ -1 +1 @@
                                                  -*-*- END OF LINE -*-*-
+
+ 2003/04/24: Wolfgang Wittek
+
+    * mfilter/MFCT1SelFinal.cc
+      - small correction in output
+
+    * mhistmc/MHMcCT1CollectionArea.cc
+      - plot all 3 histograms in DrawClone()
+
+    * manalysis/MPadSchweizer.[h,cc]
+      - cleanup as proposed by Thomas
+      - new member function SetHistograms() to set the pointers to the
+        histograms to be used in the padding
+
+
 
  2003/04/24: Antonio Stamerra

trunk/MagicSoft/Mars/manalysis/MPadSchweizer.cc

@@ -69 +69 @@
 #include "TH2.h"
 #include "TH3.h"
-#include "TProfile.h"
 #include "TRandom.h"
 #include "TCanvas.h"
@@ -93 +92 @@
 // Default constructor. 
 //
-MPadSchweizer::MPadSchweizer(const char *name, const char *title,
-                             TH2D *fHist2, TH3D *fHist3, TH3D *fHist3Diff) : 
-               fRunType(0), fGroup(0)
+MPadSchweizer::MPadSchweizer(const char *name, const char *title) 
 {
   fName  = name  ? name  : "MPadSchweizer";
   fTitle = title ? title : "Task for the padding (Schweizer)";
 
-  fHSigmaTheta    = fHist2;
-  fHSigmaPixTheta = fHist3;
-  fHDiffPixTheta  = fHist3Diff;
-
-  fHSigmaTheta->SetDirectory(NULL);
-  fHSigmaPixTheta->SetDirectory(NULL);
-  fHDiffPixTheta->SetDirectory(NULL);
-
-  Print();
+  //fHSigmaTheta    = fHist2;
+  //fHSigmaPixTheta = fHist3;
+  //fHDiffPixTheta  = fHist3Diff;
+
+  //fHSigmaTheta->SetDirectory(NULL);
+  //fHSigmaPixTheta->SetDirectory(NULL);
+  //fHDiffPixTheta->SetDirectory(NULL);
+
+  //Print();
 }
 
@@ -122 +119 @@
 // --------------------------------------------------------------------------
 //
+// Set the references to the histograms to be used in the padding
+//
+//
+void MPadSchweizer::SetHistograms(TH2D *fHist2, TH3D *fHist3, TH3D *fHist3Diff)
+{
+  fHSigmaTheta    = fHist2;
+  fHSigmaPixTheta = fHist3;
+  fHDiffPixTheta  = fHist3Diff;
+
+  fHSigmaTheta->SetDirectory(NULL);
+  fHSigmaPixTheta->SetDirectory(NULL);
+  fHDiffPixTheta->SetDirectory(NULL);
+
+  Print();
+}
+
+// --------------------------------------------------------------------------
+//
 // Set the option for the padding 
 //
+//  There are 2 options for the padding :                                  //
+//                                                                         //
+//  1) fPadFlag = 1 :                                                      //
+//     Generate first a Sigmabar using the 2D-histogram Sigmabar vs. Theta //
+//     (fHSigmaTheta). Then generate a pedestal sigma for each pixel using //
+//     the 3D-histogram Theta, pixel no., Sigma^2-Sigmabar^2               //
+//     (fHDiffPixTheta).                                                   //
+//                                                                         //
+//     This is the preferred option as it takes into account the           //
+//     correlations between the Sigma of a pixel and Sigmabar.             //
+//                                                                         //
+//  2) fPadFlag = 2 :                                                      //
+//     Generate a pedestal sigma for each pixel using the 3D-histogram     //
+//     Theta, pixel no., Sigma (fHSigmaPixTheta).                          //
+//                                                                         //
 void MPadSchweizer::SetPadFlag(Int_t padflag)
 {

trunk/MagicSoft/Mars/manalysis/MPadSchweizer.h

@@ -14 +14 @@
 #include "TH2.h"
 #include "TH3.h"
-#include "TProfile.h"
 
 
@@ -21 +20 @@
 class MPedestalCam;
 class MMcEvt;
-class MPedestalCam;
 class MSigmabar;
 class MParList;
@@ -63 +61 @@
 
 public:
-    MPadSchweizer(const char *name, const char *title,
-                  TH2D *fHist2, TH3D *fHist3, TH3D *fHist3Diff);
+    MPadSchweizer(const char *name=NULL, const char *title=NULL);
     ~MPadSchweizer();
+
+    void SetHistograms(TH2D     *fHist2=NULL, TH3D *fHist3=NULL, 
+                       TH3D *fHist3Diff=NULL);
 
     Bool_t PreProcess(MParList *pList);
@@ -72 +72 @@
     
     void SetPadFlag(Int_t padflag);
-    void SetRunType(Int_t runtype) { fRunType =  runtype; }
-    void SetGroup(Int_t group)     { fGroup   =  group; }
 
     ClassDef(MPadSchweizer, 0)    // task for the padding (Schweizer)

trunk/MagicSoft/Mars/mfilter/MFCT1SelFinal.cc

@@ -83 +83 @@
 void MFCT1SelFinal::SetCuts(Float_t hadmax, Float_t alphamax, Float_t distmax) 
 { 
-  fHadronnessMax = hadmax; 
+  fHadronnessMax =   hadmax; 
   fAlphaMax      = alphamax;
-  fDistMax       = distmax; 
-
-  *fLog << inf << "MFCT1SelFinal cut values : fHadronnessMax, fAlphaMax = "
-        << fHadronnessMax << ",  " << fAlphaMax << ",  fDistMax = " << fDistMax 
-<<  endl;
+  fDistMax       =  distmax; 
+
+  *fLog << inf << "MFCT1SelFinal cut values : fHadronnessMax, fAlphaMax, fDistMax = "
+        << fHadronnessMax << ",  " << fAlphaMax << ",  " << fDistMax <<  endl;
 }
 

trunk/MagicSoft/Mars/mhistmc/MHMcCT1CollectionArea.cc

@@ -183 +183 @@
 TObject *MHMcCT1CollectionArea::DrawClone(Option_t* option) const
 {
-  TCanvas *c = MH::MakeDefCanvas(fHistCol);
+  TCanvas &c = *MakeDefCanvas("CollArea", "Collection area plots", 600, 600);
+  c.Divide(2,2);
 
   //
@@ -190 +191 @@
   gROOT->SetSelectedPad(NULL);
 
+  c.cd(1);
+  fHistCol->SetDirectory(NULL);
   fHistCol->DrawCopy(option);
 
-  c->Modified();
-  c->Update();
-
-  return c;
+  c.cd(2);
+  fHistSel->SetDirectory(NULL);
+  fHistSel->DrawCopy(option);
+
+  c.cd(3);
+  fHistAll->SetDirectory(NULL);
+  fHistAll->DrawCopy(option);
+
+
+  c.Modified();
+  c.Update();
+
+  return &c;
 }