Changeset 7409 for trunk/MagicSoft/Mars/mfilter

Timestamp:

11/18/05 17:15:30 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mfilter

Files:

• MFEnergySlope.cc (modified) (4 diffs)
• MFEnergySlope.h (modified) (3 diffs)

trunk/MagicSoft/Mars/mfilter/MFEnergySlope.cc

@@ -18 +18 @@
 !   Author(s): Antonio Stamerra  02/2003 <mailto:antonio.stamerra@pi.infn.it>
 !
-!   Copyright: MAGIC Software Development, 2000-2003
+!   Copyright: MAGIC Software Development, 2000-2005
 !
 !
@@ -65 +65 @@
 // --------------------------------------------------------------------------
 //
+//  Default Constructor
+//
+MFEnergySlope::MFEnergySlope(const char *name, const char *title):
+    fNewSlope(-1), fMcMinEnergy(-1.), fMcMaxEnergy(-1.)
+{
+    fName  = name  ? name  : "MFEnergySlope";
+    fTitle = title ? title : "Filter to select energy with given slope";
+}
+
+// --------------------------------------------------------------------------
+//
 //     Constructor
 //
-MFEnergySlope::MFEnergySlope(const char *name, const char *title):
-  fNumSelectedEvts(0), fNewSlope(-1), fMcMinEnergy(-1.), fMcMaxEnergy(-1.)
-{
-  //    fContName = cname;
-  fName  = name  ? name  : "MFEnergySlope";
-  fTitle = title ? title : "Filter to select energy with given slope";
+MFEnergySlope::MFEnergySlope(Float_t slope, const char *name, const char *title):
+    fNewSlope(TMath::Abs(slope)), fMcMinEnergy(-1.), fMcMaxEnergy(-1.)
+{
+    fName  = name  ? name  : "MFEnergySlope";
+    fTitle = title ? title : "Filter to select energy with given slope";
 }
 
@@ -84 +94 @@
 Int_t MFEnergySlope::PreProcess(MParList *pList)
 {
- 
+    if (fNewSlope<0)
+    {
+        *fLog << err << "New slope still undefined... aborting." << endl;
+        return kFALSE;
+    }
+
+    fEvt = (MMcEvt*)pList->FindObject("MMcEvt");
+    if (!fEvt)
+    {
+        *fLog << err << "MMcEvt not found... aborting." << endl;
+        return kFALSE;
+    }
+
+    return kTRUE;
+}
+
+// --------------------------------------------------------------------------
+//
+//   Preprocess
+//  
+//  MC slope and min/max energy are read
+//  Normalization factor is computed
+//
+Bool_t MFEnergySlope::ReInit(MParList *pList)
+{
     MMcCorsikaRunHeader *runheader = (MMcCorsikaRunHeader*)pList->FindObject("MMcCorsikaRunHeader");
-
     if (!runheader)
-      {
-        *fLog << err << dbginf << fName << " [MMcCorsikaRunHeader] not found... aborting." << endl;
-        return kFALSE;
-      }
+    {
+        *fLog << err << "MMcCorsikaRunHeader not found... aborting." << endl;
+        return kFALSE;
+    }
+
     //
-    // Read info from the MC sample (it must be generated with 
+    // Read info from the MC sample (it must be generated with
     //   reflector ver.0.6 and camera ver. 0.6)
     //
-    fMcSlope = runheader->GetSlopeSpec();    
+    fMcSlope = runheader->GetSlopeSpec();
     if (fMcMinEnergy<0)
-      fMcMinEnergy = runheader->GetELowLim();
+        fMcMinEnergy = runheader->GetELowLim();
     if (fMcMinEnergy<0)
-      fMcMaxEnergy = runheader->GetEUppLim();
-
-    *fLog << inf;
-    *fLog << "MFEnergySlope::PreProcess: fetched MC info:" << endl;
-    *fLog << "  E Slope:     " << fMcSlope << endl;
-    *fLog << "  E Min:       " << fMcMinEnergy << endl;
-    *fLog << "  E Max:       " << fMcMaxEnergy << endl;
-    *fLog << "  New E Slope: " << fNewSlope << endl;
-    
+        fMcMaxEnergy = runheader->GetEUppLim();
+
     // Slope is used with positive values in the code
     if (fNewSlope < 0)
-      fNewSlope *= -1; 
+        fNewSlope *= -1;
     if (fMcSlope < 0)
-      fMcSlope *= -1;
-
-
-  // Set normalization on energy  
-    fN0 = pow(fNewSlope>fMcSlope?fMcMinEnergy:fMcMaxEnergy,fNewSlope-fMcSlope);
-
-  *fLog << "Normalization factor:" <<fN0 << endl;
-
-  //---
-    fEvt = (MMcEvt*)pList->FindObject("MMcEvt");
-    if (!fEvt)
-      {
-        *fLog << err << dbginf << fName << " [MMcEvt] not found... aborting." << endl;
-        return kFALSE;
-      }
+        fMcSlope *= -1;
+
+    // Calculate normalization factor
+    fN0 = TMath::Power(fNewSlope>fMcSlope ? fMcMinEnergy : fMcMaxEnergy, fNewSlope-fMcSlope);
+
+    // Print some infos
+    *fLog << inf;
+    *fLog << "Fetched MC info:" << endl;
+    *fLog << "  Change E Slope " << fMcSlope << " (" << fMcMinEnergy << " < E < ";
+    *fLog << fMcMaxEnergy << ") to " << fNewSlope << endl;
+    *fLog << "  Norm factor: " << fN0 << endl;
 
     return kTRUE;
@@ -144 +167 @@
 Int_t MFEnergySlope::Process()
 {
-  fResult = kTRUE;
-
-  // Energy slopes are the same: skip it
-  if (fNewSlope == fMcSlope)
+    fResult = kTRUE;
+
+    // Energy slopes are the same: skip it
+    if (fNewSlope == fMcSlope)
+        return kTRUE;
+
+    //  The value of the normalized spectrum is compared with a
+    //   random value in [0,1];
+    //   if the latter is higher the event is accepted
+    const Float_t energy = fEvt->GetEnergy();
+
+    /*
+     //
+     // If energy bounds different from MC ones are selected, then
+     // events outside these bounds are rejected, as their slope has
+     // not been changed.
+     //
+     if (energy > fMcMaxEnergy || energy < fMcMinEnergy)
+     {
+        fResult = kFALSE;
+        return kTRUE;
+     }
+     */
+
+    const Float_t Nexp = fN0 * pow(energy,fMcSlope-fNewSlope);
+    const Float_t Nrnd = gRandom->Uniform();
+
+    fResult = Nexp > Nrnd;
+
+    //if (fResult)
+    //    fNumSelectedEvts++;
+
     return kTRUE;
-  
-  //  The value of the normalized spectrum is compared with a 
-  //   random value in [0,1]; 
-  //   if the latter is higher the event is accepted
-  const Float_t energy = fEvt->GetEnergy();
-
-  /*
-  //
-  // If energy bounds different from MC ones are selected, then
-  // events outside these bounds are rejected, as their slope has
-  // not been changed.
-  //
-  if (energy > fMcMaxEnergy || energy < fMcMinEnergy)
-    {
-      fResult = kFALSE;
-      return kTRUE;
-    }
-  */
-
-  const Float_t Nexp = fN0 * pow(energy,fMcSlope-fNewSlope);
-  const Float_t Nrnd = gRandom->Uniform();
-
-  fResult = Nexp > Nrnd;
-
-  if (!fResult)
-      return kTRUE;
-
-  fNumSelectedEvts++;
-  return kTRUE;
-}
-
+}
+
+
+// --------------------------------------------------------------------------
+//
+//   MFEnergySlope.NewSlope: -2.8
+//
+Int_t MFEnergySlope::ReadEnv(const TEnv &env, TString prefix, Bool_t print)
+{
+    Bool_t rc = kFALSE;
+    if (IsEnvDefined(env, prefix, "NewSlope", print))
+    {
+        rc = kTRUE;
+        SetNewSlope(GetEnvValue(env, prefix, "NewSlope", fNewSlope));
+    }
+    return rc;
+}

trunk/MagicSoft/Mars/mfilter/MFEnergySlope.h

@@ -1 +1 @@
 #ifndef MARS_MFEnergySlope
 #define MARS_MFEnergySlope
-/////////////////////////////////////////////////////////////////////////////
-//                                                                         //
-// MFEnergySlope                                                           //
-//                                                                         //
-/////////////////////////////////////////////////////////////////////////////
 
 #ifndef MARS_MFilter
@@ -11 +6 @@
 #endif
 
+class MMcEvt;
 class MParList;
-class MMcEvt;
 class MMcCorsikaRunHeader;
 
@@ -18 +13 @@
 {
 private:
-    Int_t fNumSelectedEvts; // counter for number of selected events
+    //Int_t fNumSelectedEvts; // counter for number of selected events
 
-    MMcEvt *fEvt;           // Events used to determin energy slope
+    MMcEvt *fEvt;           //! Events used to determin energy slope
 
-    Bool_t fResult;         // Result returned by IsExpressionTrue
     Float_t fNewSlope;      // New slope set by user
+    Float_t fMcSlope;       //! Original energy slope from MC data
 
-    Float_t fMcSlope;       // Original energy slope from MC data
-    Float_t fMcMinEnergy;   // Starting energy of MC data
-    Float_t fMcMaxEnergy;   // Ending energy of MC data
+    Float_t fMcMinEnergy;   //! Starting energy of MC data
+    Float_t fMcMaxEnergy;   //! Ending energy of MC data
 
-    Float_t fN0;            // Normalization factor
+    Float_t fN0;            //! Normalization factor
 
-    Int_t PreProcess(MParList *pList);
-    Int_t Process();
+    Bool_t  fResult;        //! Result returned by IsExpressionTrue
+
+    // MTask
+    Int_t  PreProcess(MParList *pList);
+    Bool_t ReInit(MParList *pList);
+    Int_t  Process();
+
+    // MFilter
+    Bool_t IsExpressionTrue() const { return fResult; }
 
 public:
     MFEnergySlope(const char *name=NULL, const char *title=NULL);
+    MFEnergySlope(Float_t slope, const char *name=NULL, const char *title=NULL);
 
-    Bool_t IsExpressionTrue() const { return fResult; }
+    // Setter
+    void SetNewSlope(Float_t f)    { fNewSlope = TMath::Abs(f); }
+    void SetMcSlope(Float_t f)     { fMcSlope  = TMath::Abs(f); }
 
-    // Slope is used with positive values in the code 
-    void SetNewSlope(Float_t f) {fNewSlope = TMath::Abs(f);}
-    void SetMcSlope(Float_t f) {fMcSlope = TMath::Abs(f);}
+    void SetMcMinEnergy(Float_t f) { fMcMinEnergy = f; }
+    void SetMcMaxEnergy(Float_t f) { fMcMaxEnergy = f; }
 
-    void SetMcMinEnergy(Float_t f) {fMcMinEnergy = f;}
-    void SetMcMaxEnergy(Float_t f) {fMcMaxEnergy = f;}
+    // MParContainer
+    Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);
 
     ClassDef(MFEnergySlope, 0) // A Filter to select events with a given energy slope