Changeset 4450 for trunk/MagicSoft/Mars/manalysis

Timestamp:

08/04/04 11:18:07 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MCameraData.cc (modified) (9 diffs)

trunk/MagicSoft/Mars/manalysis/MCameraData.cc

@@ -17 +17 @@
 !
 !   Author(s): Thomas Bretz, 10/2003 <mailto:tbretz@astro.uni-wuerzburg.de>
-!              Hendrik Bartko, 08/2004 <mailto:hbartko@mppmu.mpg.de>
-!
-!   Copyright: MAGIC Software Development, 2000-2003
+!   Author(s): Hendrik Bartko, 08/2004 <mailto:hbartko@mppmu.mpg.de>
+!
+!   Copyright: MAGIC Software Development, 2000-2004
 !
 !
@@ -27 +27 @@
 //
 // MCameraData
+//
+// This is a generalized class storing camera data. For example the cleaning
+// level for the image cleaning is one possibility.
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -56 +59 @@
 {
     fName  = name  ? name  : "MCameraData";
-    fTitle = title ? title : "Photons/PedRMS Information";
+    fTitle = title ? title : "Generalized storage container for camera contents";
 }
 
@@ -91 +94 @@
 // scales with the inverse square root of the pixel area.
 //
-
+// FIXME: Should the check noise<=0 be replaced by MBadPixels?
+//
 void MCameraData::CalcCleaningLevel(const MCerPhotEvt &evt, const MPedPhotCam &cam,
                                     const MGeomCam &geom)
@@ -138 +142 @@
 // inner and outer pixels into account.
 //
-
+// FIXME: Should the check noise<=0 be replaced by MBadPixels?
+//
 void MCameraData::CalcCleaningLevel2(const MCerPhotEvt &evt, const MPedPhotCam &cam,
-                                    const MGeomCam &geom)
-{
-    const Int_t n = geom.GetNumPixels();
-
-    fData.Set(n);
-    fData.Reset();
-
-    fValidity.Set(n);
-    fValidity.Reset();
-
-    const Int_t entries = evt.GetNumPixels();
-
-    const Float_t meannoise[2] = {cam.GetArea(0).GetRms(), cam.GetArea(1).GetRms()} ; // mean noise for the inner and the outer pixels
-
-    //
-    // check the number of all pixels against the noise level and
-    // set them to 'unused' state if necessary
-    // 
-    for (Int_t i=0; i<entries; i++)
-    {
-        const MCerPhotPix &pix = evt[i];
-
-        const Int_t idx = pix.GetPixId();
-        const Float_t noise = cam[idx].GetRms();
+                                     const MGeomCam &geom)
+{
+    const Int_t n = geom.GetNumPixels();
+
+    fData.Set(n);
+    fData.Reset();
+
+    fValidity.Set(n);
+    fValidity.Reset();
+
+    const Int_t entries   = evt.GetNumPixels();
+    const Float_t anoise0 = cam.GetArea(0).GetRms();
+    if (anoise0<=0)
+        return;
+
+    //
+    // check the number of all pixels against the noise level and
+    // set them to 'unused' state if necessary
+    // 
+    for (Int_t i=0; i<entries; i++)
+    {
+        const MCerPhotPix &pix = evt[i];
+
+        const Int_t   idx   = pix.GetPixId();
+        const Float_t noide = cam[idx].GetRms();
 
         if (noise<=0) // fData[idx]=0, fValidity[idx]=0
@@ -174 +180 @@
         // with the square root of the pixel area.
         // 
-
-        const UInt_t aidx = geom[idx].GetAidx();
-
-
-        if (meannoise[0] > 0 && meannoise[1] > 0)  
-        {
-          fData[idx] = pix.GetNumPhotons() * geom.GetPixRatio(idx) * meannoise[aidx] / meannoise[0] / noise;
-         
-        }
-        else fData[idx] = pix.GetNumPhotons() * geom.GetPixRatioSqrt(idx) / noise;
-
+        const UInt_t  aidx  = geom[idx].GetAidx();
+        const Float_t ratio = cam.GetArea(aidx).GetRms()/noise0;
+
+        fData[idx] = pix.GetNumPhotons() * geom.GetPixRatio(idx) * ratio / noise;
         fValidity[idx] = 1;
     }
@@ -231 +230 @@
 }
 
-
 // --------------------------------------------------------------------------
 //
@@ -238 +236 @@
 // inner pixel and the average pedestal RMS of the inner pixels (democratic
 // image cleaning, see TDAS 02-14).
-
-
+//
+// FIXME: Should the check noise<=0 be replaced by MBadPixels?
+//
 void MCameraData::CalcCleaningLevelDemocratic(const MCerPhotEvt &evt, const MPedPhotCam &cam,
                                               const MGeomCam &geom)
@@ -251 +250 @@
     fValidity.Reset();
 
-    const Int_t entries = evt.GetNumPixels();
-
-    const Float_t meannoise[2] = {cam.GetArea(0).GetRms(), cam.GetArea(1).GetRms()} ; // mean noise for the inner and the outer pixels
-
-        
-    //
-    // check the number of all pixels against the noise level and
-    // set them to 'unused' state if necessary
-    // 
-    for (Int_t i=0; i<entries; i++)
-    {
-      const MCerPhotPix &pix = evt[i];
-      
-      const Int_t idx = pix.GetPixId();
-      const Float_t noise = cam[idx].GetRms();
-      
-      if (noise<=0)
-        continue;
-        
-      //
-      // We calculate a correction factor which accounts for the     
-      // fact that pixels have different size (see TDAS 02-14).
-      // 
-      
-      
-      if (meannoise[0]>0) // fData[idx]=0, fValidity[idx]=0
-      {
-        fData[idx] = pix.GetNumPhotons() * geom.GetPixRatio(idx) / meannoise[0];
-        fValidity[idx] = 1;
-      }
-      else 
-        fValidity[idx] = 0;
-    }
-    
-}
-
-
-
-// --------------------------------------------------------------------------
-//
-// Returns, depending on the type flag:
-//
-//  0: Number of Photons*PixRatio
-//  1: Error*sqrt(PixRatio)
-//  2: Cleaning level = Num Photons*sqrt(PixRatio)/Error
-//  3: Number of Photons
-//  4: Error
+    const Int_t   entries = evt.GetNumPixels();
+    const Float_t noise0  = cam.GetArea(0).GetRms();
+    if (noise0<=0)
+        return;
+
+    //
+    // check the number of all pixels against the noise level and
+    // set them to 'unused' state if necessary
+    // 
+    for (Int_t i=0; i<entries; i++)
+    {
+        const MCerPhotPix &pix = evt[i];
+
+        const Int_t idx     = pix.GetPixId();
+        const Float_t noise = cam[idx].GetRms();
+
+        if (noise<=0)
+            continue;
+
+        //
+        // We calculate a correction factor which accounts for the
+        // fact that pixels have different size (see TDAS 02-14).
+        //
+        fData[idx] = pix.GetNumPhotons() * geom.GetPixRatio(idx) / noise0;
+        fValidity[idx] = 1;
+    }
+}
+
+// --------------------------------------------------------------------------
+//
+// Returns the contents of the pixel.
 //
 Bool_t MCameraData::GetPixelContent(Double_t &val, Int_t idx, const MGeomCam &cam, Int_t type) const