Changeset 5365 for trunk/MagicSoft/Mars/mfilter

Timestamp:

11/09/04 20:49:06 (20 years ago)

Author:

moralejo

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mfilter/MFSize.cc (modified) (2 diffs)

trunk/MagicSoft/Mars/mfilter/MFSize.cc

@@ -59 +59 @@
 //     Constructor
 //
-MFSize::MFSize(const char *name, const char *title): fNumSelectedEvts(0)
+MFSize::MFSize(const char *name, const char *title): fNumSelectedEvts(0), fProb(NULL)
 {
   fName  = name  ? name  : "MFSize";
   fTitle = title ? title : "Filter to select events according to Size";
-
-  //
-  // Default fProb histogram: 20 bins in the log10(Size) range [2,4]
-  // Values are calculated to make the Size distribution of MC gammas
-  // in the standard MC sample to look like the distribution for MC 
-  // hadrons (p+He)
-  //
-  fProb = new TH1F("ProbVsLog10Size", "", 20, 2., 4.);
-
-  Float_t frac[20] = {0.144928, 0.171911, 0.18763, 0.209461, 0.253219, 
-                      0.305425, 0.384593, 0.485204, 0.612452, 0.708138, 
-                      0.754522, 0.728028, 0.774046, 0.791422, 0.808775, 
-                      0.896842, 1., 1., 1., 1.};
-  
-  for (Int_t i = 0; i < 20; i++)
-    fProb->SetBinContent(i+1, frac[i]);
 }
 
@@ -114 +98 @@
 Int_t MFSize::PreProcess(MParList *pList)
 {
+  if (!fProb)
+    {
+      *fLog << err << dbginf << fName <<
+        " Histogram fProb not set. You must call MFSize::SetProb(TH1F*) before the event loop." << endl;
+      return kFALSE;
+    }  
+
   fHillas = (MHillas*)pList->FindObject(AddSerialNumber("MHillas"));
     if (!fHillas)