Changeset 1574 for trunk/MagicSoft/Mars/manalysis

Timestamp:

11/04/02 10:06:08 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MBlindPixelCalc.cc (modified) (5 diffs)
• MBlindPixelCalc.h (modified) (2 diffs)
• MCameraSmooth.cc (modified) (1 diff)
• MCerPhotEvt.cc (modified) (1 diff)
• MCerPhotEvt.h (modified) (1 diff)
• MCompProbCalc.cc (modified) (2 diffs)
• MCompProbCalc.h (modified) (1 diff)
• MHillasExt.cc (modified) (1 diff)

trunk/MagicSoft/Mars/manalysis/MBlindPixelCalc.cc

@@ -93 +93 @@
 Bool_t MBlindPixelCalc::PreProcess (MParList *pList)
 {
-    fPixels = (MBlindPixels*)pList->FindCreateObj("MBlindPixels");
+    if (TESTBIT(fFlags, kUseBlindPixels))
+        fPixels = (MBlindPixels*)pList->FindObject("MBlindPixels");
+    else
+        fPixels = (MBlindPixels*)pList->FindCreateObj("MBlindPixels");
     if (!fPixels)
         return kFALSE; 
@@ -106 +109 @@
     fGeomCam = (MGeomCam*)pList->FindObject("MGeomCam");
     if (!fGeomCam)
-        *fLog << warn << dbginf << "No camera geometry available... can't ude interpolation." << endl;
+        *fLog << warn << dbginf << "No camera geometry available... can't use interpolation." << endl;
+
+    if (TESTBIT(fFlags, kUseBlindPixels))
+        return kTRUE;
 
     const UShort_t size = fPixelsID.GetSize();
@@ -162 +168 @@
         Int_t num = TESTBIT(fFlags, kUseCentralPixel) ? 1 : 0;
 
-        nphot[i] = TESTBIT(fFlags, kUseCentralPixel) ? (*fEvt)[id].GetNumPhotons() : 0;
-        perr[i]  = TESTBIT(fFlags, kUseCentralPixel) ? (*fEvt)[id].GetErrorPhot()  : 0;
+        nphot[i] = TESTBIT(fFlags, kUseCentralPixel) ? pix.GetNumPhotons() : 0;
+        perr[i]  = TESTBIT(fFlags, kUseCentralPixel) ? pix.GetErrorPhot()  : 0;
         for (int j=0; j<n; j++)
         {
@@ -171 +177 @@
                 continue;
 
-            nphot[i] += (*fEvt)[nid].GetNumPhotons();
-            perr[i]  += (*fEvt)[nid].GetErrorPhot();
-
+            const MCerPhotPix *evtpix = fEvt->GetPixById(nid);
+            if (evtpix)
+            {
+                nphot[i] += evtpix->GetNumPhotons();
+                perr[i]  += evtpix->GetErrorPhot();
+            }
             num++;
         }
@@ -253 +262 @@
 Bool_t MBlindPixelCalc::ReInit(MParList *pList)
 {
+    if (TESTBIT(fFlags, kUseBlindPixels))
+        return kTRUE;
+
     //
     // If pixels are given by the user, we are already done

trunk/MagicSoft/Mars/manalysis/MBlindPixelCalc.h

@@ -28 +28 @@
     {
         kUseInterpolation = 1,
-        kUseCentralPixel  = 2
+        kUseCentralPixel  = 2,
+        kUseBlindPixels   = 3
     };
 
@@ -46 +47 @@
         b ? SETBIT(fFlags, kUseCentralPixel) : CLRBIT(fFlags, kUseCentralPixel);
     }
+    void SetUseBlindPixels(Bool_t b=kTRUE)
+    {
+        b ? SETBIT(fFlags, kUseBlindPixels) : CLRBIT(fFlags, kUseBlindPixels);
+    }
 
     Bool_t PreProcess(MParList *pList);

trunk/MagicSoft/Mars/manalysis/MCameraSmooth.cc

@@ -111 +111 @@
             const Int_t n = gpix.GetNumNeighbors();
 
-            photons[i] = fUseCentralPixel ? (*fEvt)[id].GetNumPhotons() : 0;
-            errors[i]  = fUseCentralPixel ? (*fEvt)[id].GetErrorPhot()  : 0;
+            Int_t num  = fUseCentralPixel ? 1 : 0;
+            photons[i] = fUseCentralPixel ? pix.GetNumPhotons() : 0;
+            errors[i]  = fUseCentralPixel ? pix.GetErrorPhot()  : 0;
+
             for (int j=0; j<n; j++)
             {
                 const UShort_t nid = gpix.GetNeighbor(j);
 
-                photons[i] += (*fEvt)[nid].GetNumPhotons();
-                errors[i]  += (*fEvt)[nid].GetErrorPhot();
+                const MCerPhotPix *evtpix = fEvt->GetPixById(nid);
+                if (evtpix)
+                {
+                    photons[i] += evtpix->GetNumPhotons();
+                    errors[i]  += evtpix->GetErrorPhot();
+                }
+                num++;
             }
 
-            photons[i] /= fUseCentralPixel ? n+1 : n;
-            errors[i]  /= fUseCentralPixel ? n+1 : n;
+            photons[i] /= num;
+            errors[i]  /= num;
         }
 

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.cc

@@ -306 +306 @@
 }
 
+// --------------------------------------------------------------------------
+//
+// Return a pointer to the pixel with the requested id. NULL if it doesn't
+// exist.
+//
+MCerPhotPix *MCerPhotEvt::GetPixById(int id) const
+{
+    TIter Next(fPixels);
+    MCerPhotPix *pix = NULL;
+
+    while ((pix=(MCerPhotPix*)Next()))
+        if (pix->GetPixId()==id)
+            return pix;
+
+    return NULL;
+}

trunk/MagicSoft/Mars/manalysis/MCerPhotEvt.h

@@ -47 +47 @@
     MCerPhotPix &operator[](int i) const { return *(MCerPhotPix*)(fPixels->UncheckedAt(i)); }
 
+    MCerPhotPix *GetPixById(int id) const;
+
     void Reset();
 

trunk/MagicSoft/Mars/manalysis/MCompProbCalc.cc

@@ -77 +77 @@
 //  - MHCompProb
 //  - all data values which were used to build the MHCompProb
-//  - MHadroness
+//  - MHadronness
 //
 Bool_t MCompProbCalc::PreProcess(MParList *plist)
@@ -136 +136 @@
 //  - For all data members multiply the probabilities.
 //  - For normalization take the n-th root of the result.
-//  - This is the hadroness stored in the MHadroness container
+//  - This is the hadroness stored in the MHadronness container
 //
 Bool_t MCompProbCalc::Process()

trunk/MagicSoft/Mars/manalysis/MCompProbCalc.h

@@ -13 +13 @@
 {
 private:
-    MHadronness *fHadronness; //! Output container (Hadroness)
+    MHadronness *fHadronness; //! Output container (Hadronness)
 
     TList *fData;           //! List of MDataChains to be used

trunk/MagicSoft/Mars/manalysis/MHillasExt.cc

@@ -39 +39 @@
 //
 // WARNING: Before you can use fAsym, fM3Long and fM3Trans you must
-//          multiply by the sign of MHillasSrc::fCosAlphaDelta
+//          multiply by the sign of MHillasSrc::fCosDeltaAlpha
 //
 ////////////////////////////////////////////////////////////////////////////