Changeset 2092 for trunk/MagicSoft/Mars/mimage/MHillasExt.cc

Timestamp:

05/07/03 16:33:22 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mimage/MHillasExt.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mimage/MHillasExt.cc

@@ -32 +32 @@
 //
 //    extended image parameters
+//
+//
+// Version 1:
+// ----------
 // fConc     ratio of sum of two highest pixels over fSize
 // fConc1    ratio of highest pixel over fSize
@@ -38 +42 @@
 // fM3Trans  third moment along minor axis
 //
+// Version 2:
+// ----------
+// fConc     removed
+// fConc1    removed
+//
+//
 // WARNING: Before you can use fAsym, fM3Long and fM3Trans you must
 //          multiply by the sign of MHillasSrc::fCosDeltaAlpha
@@ -43 +53 @@
 ////////////////////////////////////////////////////////////////////////////
 /*
- // fLeakage1 ratio : (photons in most outer ring of pixels) over fSize
- // fLeakage2 ratio : (photons in the 2 outer rings of pixels) over fSize
- //
  // fAsymna   d/(d na) of ( sum(x*q^na)/sum(q^na), sum(y*q^na)/sum(q^na) )
  //           projected onto the major axis
@@ -155 +162 @@
             continue;                                          // [1]
         }
-
-        /*
-         //
-         //  power na for calculating fAsymna;
-         //  the value 1.5 was suggested by Thomas Schweizer
-         //
-         Double_t na = 1.5;
-
-         //
-         // get sums for calculating fAsymna
-         // the outer pixels are 4 times as big (in area)
-         // as the inner pixels !
-         //
-         const Double_t dummy = pow(nphot, na)/r;
-
-         sna +=     dummy;
-         xna += dzx*dummy;
-
-         sna1 += sna/nphot;
-         xna1 += xna/nphot;
-
-         //
-         // forward-backward asymmetry
-         //
-         fb += dzx<0 ? -nphot: nphot;
-         */
     }
 
@@ -186 +167 @@
 
     fAsym = (hil.GetMeanX()-maxp.GetX())*hil.GetCosDelta() +
-            (hil.GetMeanY()-maxp.GetY())*hil.GetSinDelta();    // [mm]
-
-    /*
-     fAsym0    =       fb / GetSize();
-     fAsymna   = na * (sna*xna1 - sna1*xna) / (sna*sna);
-     */
+            (hil.GetMeanY()-maxp.GetY())*hil.GetSinDelta();            // [mm]
 
     //