Changeset 3761 for trunk/MagicSoft/Mars/manalysis

Timestamp:

04/15/04 22:54:38 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/manalysis

Files:

• MPedPhotCalc.cc (modified) (8 diffs)
• MPedPhotCalc.h (modified) (1 diff)

trunk/MagicSoft/Mars/manalysis/MPedPhotCalc.cc

@@ -18 +18 @@
 !   Author(s): Josep Flix 1/2004 <mailto:jflix@ifae.es>
 !   Author(s): Javier Rico 1/2004 <mailto:jrico@ifae.es>
+!   Author(s): Markus Gaug 4/2004 <mailto:markus@ifae.es>
 !
 !   Copyright: MAGIC Software Development, 2000-2004
@@ -55 +56 @@
 #include "MPedPhotCam.h"
 
+#include "MBadPixelsPix.h"
+#include "MBadPixelsCam.h"
+
 ClassImp(MPedPhotCalc);
 
@@ -63 +67 @@
 //
 MPedPhotCalc::MPedPhotCalc(const char *name, const char *title)
+    : fPedestals(NULL), fCerPhot(NULL), fBadPixels(NULL)
 {
   fName  = name  ? name  : "MPedPhotCalc";
-  fTitle = title ? title : "Task to calculate pedestals from the charge computed from pedestal runs (in units of photons)";
+  fTitle = title ? title : "Task to calculate pedestals in units of photons";
 }
 
@@ -73 +78 @@
 //
 //  - MCerPhotEvt
+//  - MBadPixelsCam
 // 
 // The following output containers are also searched and created if
@@ -87 +93 @@
       *fLog << err << "MPedPhotCalc::PreProcess Error: MCerPhotEvt not found... aborting." << endl;
       return kFALSE;
+    }
+
+
+  fBadPixels = (MBadPixelsCam*)pList->FindObject("MBadPixelsCam");
+  if (!fBadPixels)
+    {
+      *fLog << warn << "MPedPhotCalc::PreProcess Warning: No MBadPixelsCam found." << endl;
     }
 
@@ -130 +143 @@
   }
 
+
   return kTRUE;
 }
@@ -142 +156 @@
     for(UInt_t i=0;i<fCerPhot->GetNumPixels();i++)
     {
-        const MCerPhotPix &pix = (*fCerPhot)[i];
+ 
+       const MCerPhotPix &pix = (*fCerPhot)[i];
+
+        if (fBadPixels)
+          {
+            const MBadPixelsPix &bad = (*fBadPixels)[i];
+            if (bad.IsBad())
+              continue;
+          }
 
         const Float_t nphot = pix.GetNumPhotons();
-        const Int_t idx     = pix.GetPixId();
-
-        fSumx[idx]  += nphot;
-        fSumx2[idx] += nphot*nphot;
+
+        fSumx[i]  += nphot;
+        fSumx2[i] += nphot*nphot;
     }
 
@@ -168 +189 @@
     for(Int_t i=0; i<npix; i++)
     {
-        const Float_t sum  = fSumx[i];
-        const Float_t sum2 = fSumx2[i];
-
-        const Float_t photped = sum/n;
-        const Float_t photrms = TMath::Sqrt((sum2-sum*sum/n)/(n-1.));
-
-        (*fPedestals)[i].Set(photped,photrms);
+
+      if (fBadPixels)
+        {
+          const MBadPixelsPix &bad = (*fBadPixels)[i];
+          if (bad.IsBad())
+            continue;
+        }
+
+      const Float_t sum  = fSumx[i];
+      const Float_t sum2 = fSumx2[i];
+      
+      const Float_t photped = sum/n;
+      const Float_t photrms = TMath::Sqrt((sum2-sum*sum/n)/(n-1.));
+      
+      (*fPedestals)[i].Set(photped,photrms);
     }
 

trunk/MagicSoft/Mars/manalysis/MPedPhotCalc.h

@@ -12 +12 @@
 class MPedPhotCam;
 class MCerPhotEvt;
-
+class MBadPixelsCam;
 class MPedPhotCalc : public MTask
 {
-    MPedPhotCam  *fPedestals;  // Pedestals of all pixels in the camera
-    MCerPhotEvt  *fCerPhot;
 
-    TArrayF fSumx;   // sum of values
-    TArrayF fSumx2;  // sum of squared values
+  MPedPhotCam   *fPedestals;  // Pedestals of all pixels in the camera
+  MCerPhotEvt   *fCerPhot;    // Calibrated Cherenkov events 
+  MBadPixelsCam *fBadPixels;  // Bad Pixels
+  
+  TArrayF fSumx;   // sum of values
+  TArrayF fSumx2;  // sum of squared values
 
-    Bool_t ReInit(MParList *pList);
+  Bool_t ReInit(MParList *pList);
 
-    Int_t PreProcess(MParList *pList);
-    Int_t Process();
-    Int_t PostProcess();
+  Int_t PreProcess(MParList *pList);
+  Int_t Process();
+  Int_t PostProcess();
+  
+public:
+  MPedPhotCalc(const char *name=NULL, const char *title=NULL);
 
-public:
-    MPedPhotCalc(const char *name=NULL, const char *title=NULL);
-
-    ClassDef(MPedPhotCalc, 0)//Task to calculate pedestals from the charge computed from pedestal runs (in units of photons)
+  ClassDef(MPedPhotCalc, 0) //Task to calculate pedestals in units of photons
 };