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

Timestamp:

04/06/04 13:41:56 (21 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mimage/MHillasExt.cc

@@ -47 +47 @@
 // fConc1    removed
 //
+// Version 3:
+// ----------
+// fMaxDist  added. Distance between center and most distant used pixel
+//
 //
 // WARNING: Before you can use fAsym, fM3Long and fM3Trans you must
@@ -96 +100 @@
     fM3Long  =  0;
     fM3Trans =  0;
+
+    fMaxDist = -1;
 }
 
@@ -109 +115 @@
     *fLog << " - 3.Moment Long  [mm]  = " << fM3Long  << endl;
     *fLog << " - 3.Moment Trans [mm]  = " << fM3Trans << endl;
+    *fLog << " - Max.Dist       [mm]  = " << fMaxDist << endl;
 }
 
@@ -123 +130 @@
     *fLog << " - 3.Moment Long  [deg] = " << fM3Long *geom.GetConvMm2Deg() << endl;
     *fLog << " - 3.Moment Trans [deg] = " << fM3Trans*geom.GetConvMm2Deg() << endl;
+    *fLog << " - Max.Dist       [deg] = " << fMaxDist*geom.GetConvMm2Deg() << endl;
 }
 
@@ -146 +154 @@
     const UInt_t npixevt = evt.GetNumPixels();
 
-    Int_t maxpixid = 0;
-    Float_t maxpix = 0;
+    Int_t maxpixid  = 0;
+    Float_t maxpix  = 0;
+    Float_t maxdist = 0;
 
     for (UInt_t i=0; i<npixevt; i++)
@@ -160 +169 @@
         const Double_t dx = gpix.GetX() - hil.GetMeanX();      // [mm]
         const Double_t dy = gpix.GetY() - hil.GetMeanY();      // [mm]
+
+        const Double_t dist = dx*dx+dy*dy;
+        if (dist>maxdist)
+            maxdist=dist;                                      // [mm^2]
 
         Double_t nphot = pix.GetNumPhotons();                  // [1]
@@ -197 +210 @@
     fM3Trans = m3y<0 ? -pow(-m3y, 1./3) : pow(m3y, 1./3);      // [mm]
 
+    fMaxDist = TMath::Sqrt(maxdist);                           // [mm]
+
     SetReadyToSave();