Changeset 2209 for trunk/MagicSoft/Mars/manalysis

Timestamp:

06/23/03 12:09:05 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MSigmabar.cc (modified) (2 diffs)
• MSigmabarCalc.cc (modified) (7 diffs)
• MSigmabarCalc.h (modified) (1 diff)

trunk/MagicSoft/Mars/manalysis/MSigmabar.cc

@@ -145 +145 @@
         }
 
-        const MGeomPix    &gpix = geom[id];
-        const MPedestalPix &pix =  ped[id];
+        const MGeomPix &gpix = geom[id];
 
         Int_t sector = (Int_t)(atan2(gpix.GetY(),gpix.GetX())*6 / (TMath::Pi()*2));
@@ -153 +152 @@
 
         // count only those pixels which have a sigma != 0.0
-        const Float_t sigma = pix.GetMeanRms();
+        const Float_t sigma = ped[id].GetMeanRms();
 
         if ( sigma <= 0 )

trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.cc

@@ -65 +65 @@
 }
 
-MSigmabarCalc::~MSigmabarCalc()
-{
-  // nothing special yet.
-}
-
 // --------------------------------------------------------------------------
 //
@@ -75 +70 @@
 //  if not create one and add them to the list
 //
-Bool_t MSigmabarCalc::PreProcess(MParList *pList)
+Int_t MSigmabarCalc::PreProcess(MParList *pList)
 {
     fCam = (MGeomCam*)pList->FindObject("MGeomCam");
@@ -91 +86 @@
     }
 
-    fSig = (MSigmabar*)pList->FindCreateObj("MSigmabar");
-    if (!fSig)
-    {
-        *fLog << err << "MSigmabar not found... aborting." << endl;
-        return kFALSE;
-    }
-
-    fSigParam = (MSigmabarParam*)pList->FindCreateObj("MSigmabarParam");
-    if (!fSigParam)
-    {
-        *fLog << err << "MSigmabarParam not found... aborting." << endl;
-        return kFALSE;
-    }
-
     fRun = (MRawRunHeader*)pList->FindObject("MRawRunHeader");
     if (!fRun)
@@ -111 +92 @@
         return kFALSE;
     }
-    
+
     // This is needed for determining min/max Theta
     fMcEvt = (MMcEvt*)pList->FindObject("MMcEvt");
     if (!fMcEvt)
-      {
-        *fLog << err << "MMcEvt not found... aborting." << endl;
-        return kFALSE;
-      }
+    {
+        *fLog << warn << "MMcEvt not found... aborting." << endl;
+        fThetaMin =  0;
+        fThetaMax = 90;
+    }
 
     fEvt = (MCerPhotEvt*)pList->FindObject("MCerPhotEvt");
@@ -126 +108 @@
         return kFALSE;
       }
+
+    fSig = (MSigmabar*)pList->FindCreateObj("MSigmabar");
+    if (!fSig)
+        return kFALSE;
+
+    fSigParam = (MSigmabarParam*)pList->FindCreateObj("MSigmabarParam");
+    if (!fSigParam)
+        return kFALSE;
+
 
     return kTRUE;
@@ -137 +128 @@
 // values. Preliminary.
 //
-Bool_t MSigmabarCalc::Process()
+Int_t MSigmabarCalc::Process()
 {
     const Double_t rc = fSig->Calc(*fCam, *fPed, *fEvt);
 
-    fSigmabarMax = TMath::Max(rc, fSigmabarMax);
-    fSigmabarMin = TMath::Min(rc, fSigmabarMin);
+    if (rc>fSigmabarMax) fSigmabarMax=rc;
+    if (rc<fSigmabarMin) fSigmabarMin=rc;
+
+    if (!fMcEvt)
+        return kTRUE;
 
     const Double_t theta = fMcEvt->GetTelescopeTheta()*kRad2Deg;
 
-    fThetaMax = TMath::Max(theta, fThetaMax);
-    fThetaMin = TMath::Min(theta, fThetaMin);
+    if (theta>fSigmabarMax) fSigmabarMax=theta;
+    if (theta<fSigmabarMax) fSigmabarMin=theta;
 
     return kTRUE;
@@ -169 +163 @@
 void MSigmabarCalc::Reset()
 {
-    fThetaMin = 200;    // there must be a function which gives me the hightest
-    fThetaMax = 0;      // value allowed for a certain type!
-    fSigmabarMin = 200;
+    if (!fMcEvt)
+        return;
+
+    fThetaMin    = FLT_MAX;
+    fThetaMax    = 0;
+    fSigmabarMin = FLT_MAX;
     fSigmabarMax = 0;
 }

trunk/MagicSoft/Mars/manalysis/MSigmabarCalc.h

@@ -48 +48 @@
     void Reset();
 
+    Int_t PreProcess(MParList *pList);
+    Int_t Process();
+
+    Bool_t ReInit(MParList *pList);
+
 public:
     MSigmabarCalc(const char *name=NULL, const char *title=NULL);
-    ~MSigmabarCalc();
-
-    Bool_t PreProcess(MParList *pList);
-    Bool_t ReInit(MParList *pList);
-    Bool_t Process();
 
     ClassDef(MSigmabarCalc, 0) // task for calculating sigmabar