Changeset 8647 for trunk/MagicSoft/Mars/mimage/MNewImagePar.cc

Timestamp:

07/26/07 12:21:51 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mimage/MNewImagePar.cc (modified) (10 diffs)

trunk/MagicSoft/Mars/mimage/MNewImagePar.cc

@@ -32 +32 @@
 //    Float_t fLeakage1;             // (photons in most outer ring of pixels) over fSize
 //    Float_t fLeakage2;             // (photons in the 2 outer rings of pixels) over fSize
-//    Float_t fInnerLeakage1;        // (photons in most outer rings of inner pixels) over fInnerSize
-//    Float_t fInnerLeakage2;        // (photons in the 2 outer rings of inner pixels) over fInnerSize
-//    Float_t fInnerSize;            //
 //
 //    Float_t fConc;                 // [ratio] concentration ratio: sum of the two highest pixels / fSize
@@ -72 +69 @@
 //  - added fConcCore
 //
+// Version 6:
+// ----------
+//  - removed fInnerLeakage1
+//  - removed fInnerLeakage2
+//  - removed fInnerSize
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -109 +111 @@
     fLeakage1 = -1;
     fLeakage2 = -1;
-
-    fInnerLeakage1 = -1;
-    fInnerLeakage2 = -1;
-    fInnerSize     = -1;
 
     fConc     = -1;
@@ -139 +137 @@
     fCoreArea = 0;
 
-    fInnerSize = 0;
     fConcCore  = 0;
 
     Double_t edgepix1 = 0;
     Double_t edgepix2 = 0;
-
-    Double_t edgepixin1 = 0;
-    Double_t edgepixin2 = 0;
 
     Float_t maxpix1 = 0;                                 // [#phot]
@@ -157 +151 @@
     const Double_t rl = 1./(hillas.GetLength()*hillas.GetLength());
     const Double_t rw = 1./(hillas.GetWidth() *hillas.GetWidth());
-
-    const Bool_t ismagiclike =
-        geom.GetNumPixels() == 577 &&
-        geom.GetNumAreas()  == 2   &&
-        geom.GetPixRatio(396) > geom.GetPixRatio(397);
 
     UInt_t npix = evt.GetNumPixels();
@@ -245 +234 @@
            edgepix2 += nphot;
 
-        const Double_t ratio = geom.GetPixRatio(i);
-        if (TMath::Nint(ratio)==1) // Means this is a small (= inner pixel)
-        {
-            fInnerSize += nphot;
-
-            // Do calculation of "inner leakage" only for MAGIC-like geometry,
-            // i.e., 577 pixels, pixels of 2 different areas, inner part
-            // from pixel 0 to 396:
-            if (ismagiclike)
-            {
-                if(i > 270) // last two "rings" of inner pixels
-                {
-                    edgepixin2 += nphot;
-                    if(i > 330)  // last "ring" of inner pixels
-                        edgepixin1 += nphot;
-                }
-            }
-        }
-
         //
         // Now convert nphot from absolute number of photons or phe to signal
@@ -269 +239 @@
         // density:
         //
-        nphot *= ratio;
+        nphot *= geom.GetPixRatio(i);
 
         // Look for signal density in two highest pixels:
@@ -282 +252 @@
     }
 
-    fInnerLeakage1 = edgepixin1 / fInnerSize;
-    fInnerLeakage2 = edgepixin2 / fInnerSize;
-
     fLeakage1 = edgepix1 / hillas.GetSize();
     fLeakage2 = edgepix2 / hillas.GetSize();
@@ -318 +285 @@
     *fLog << " - Leakage1         [1] = " << fLeakage1      << endl;
     *fLog << " - Leakage2         [1] = " << fLeakage2      << endl;
-    *fLog << " - InnerLeakage1    [1] = " << fInnerLeakage1 << endl;
-    *fLog << " - InnerLeakage2    [1] = " << fInnerLeakage2 << endl;
-    *fLog << " - InnerSize      [phe] = " << fInnerSize     << endl;
     *fLog << " - Conc             [1] = " << fConc          << endl;
     *fLog << " - Conc1            [1] = " << fConc1         << endl;
@@ -342 +306 @@
     *fLog << " - Leakage1         [1] = " << fLeakage1      << endl;
     *fLog << " - Leakage2         [1] = " << fLeakage2      << endl;
-    *fLog << " - InnerLeakage1    [1] = " << fInnerLeakage1 << endl;
-    *fLog << " - InnerLeakage2    [1] = " << fInnerLeakage2 << endl;
-    *fLog << " - InnerSize      [phe] = " << fInnerSize     << endl;
     *fLog << " - Conc             [1] = " << fConc          << endl;
     *fLog << " - Conc1            [1] = " << fConc1         << endl;