Changeset 9850 for trunk/Mars/mimage

Timestamp:

08/12/10 09:34:13 (14 years ago)

Author:

tbretz

Message:

Adapted MHillasExt member functions. Fixed code to calculate weighted time spread.

Location:

trunk/Mars/mimage

Files:

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

trunk/Mars/mimage/MHillasExt.cc

@@ -121 +121 @@
     fSlopeLong  =  0;
     fSlopeTrans =  0;
+
+    fTimeSpread          = -1;
+    fTimeSpreadWeighted  = -1;
+    fSlopeSpread         = -1;
+    fSlopeSpreadWeighted = -1;
 }
 
@@ -202 +207 @@
 
         sumxw   += x*nphot;
-        sumyw   += x*nphot;
+        sumyw   += y*nphot;
         sumtw   += t*nphot;
 
@@ -250 +255 @@
     // Time spread
     //
-    fTimeSpread         = TMath::Sqrt(sumt2/cnt - sumt*sumt/cnt/cnt);
-    fTimeSpreadWeighted = TMath::Sqrt(sumt2w2 + (sumw2 - 2*sumtw2)*sumtw/sumw)/sumw;
+    fTimeSpread         = TMath::Sqrt(sumt2*cnt - sumt*sumt)/cnt;
+    fTimeSpreadWeighted = TMath::Sqrt(sumt2w2 + sumw2*sumtw*sumtw/sumw/sumw - 2*sumtw2*sumtw/sumw)/sumw;
 
     //
@@ -280 +285 @@
     const Double_t sumdt2w2  = sumt2w2 - 2*fSlopeLong*(c*sumtxw2 + s*sumtyw2) + fSlopeLong*fSlopeLong * (c*c*sumx2 + s*s*sumy2 + 2*c*s*sumxy) * sumw2;
 
-    fSlopeSpread         = TMath::Sqrt(sumdt2/cnt - sumdt*sumdt/cnt/cnt);
-    fSlopeSpreadWeighted = TMath::Sqrt(sumdt2w2 + (sumw2 - 2*sumdtw2)*sumdtw/sumw)/sumw;
-
-    //
-    // FIXME: Missing Time spread along slope
-    //
+    fSlopeSpread         = TMath::Sqrt(sumdt2*cnt - sumdt*sumdt)/cnt;
+    fSlopeSpreadWeighted = TMath::Sqrt(sumdt2w2 + sumw2*sumdtw*sumdtw/sumw/sumw - 2*sumdtw2*sumdtw/sumw)/sumw;
 
     //
@@ -314 +315 @@
 void MHillasExt::Set(const TArrayF &arr)
 {
-    if (arr.GetSize() != 5)
+    if (arr.GetSize() != 9)
         return;
 
@@ -322 +323 @@
     fSlopeLong  = arr.At(3);
     fSlopeTrans = arr.At(4);
+
+    fTimeSpread          = arr.At(5);
+    fTimeSpreadWeighted  = arr.At(6);
+    fSlopeSpread         = arr.At(7);
+    fSlopeSpreadWeighted = arr.At(8);
 }
 
@@ -341 +347 @@
     if (showtrans)
         *fLog << " - Slope Trans    [mm]  = " << fSlopeTrans << endl;
+    *fLog << " - Time Spread    [ns]  = " << fTimeSpread << endl;
+    *fLog << " - Slope Spread   [ns]  = " << fSlopeSpread << endl;
+    *fLog << " - Time Spread W  [ns]  = " << fTimeSpreadWeighted << endl;
+    *fLog << " - Slope Spread W [ns]  = " << fSlopeSpreadWeighted << endl;
 }
 
@@ -357 +367 @@
     *fLog << " - Slope Long     [mm]  = " << fSlopeLong*geom.GetConvMm2Deg()  << endl;
     *fLog << " - Slope Trans    [mm]  = " << fSlopeTrans*geom.GetConvMm2Deg() << endl;
+    *fLog << " - Slope Spread   [ns]  = " << fSlopeSpread << endl;
+    *fLog << " - Time Spread W  [ns]  = " << fTimeSpreadWeighted << endl;
+    *fLog << " - Slope Spread W [ns]  = " << fSlopeSpreadWeighted << endl;
 }
 

trunk/Mars/mimage/MHillasExt.h

@@ -35 +35 @@
     void Reset();
 
-    Float_t GetAsym() const               { return fAsym; }
-    Float_t GetM3Long() const             { return fM3Long; }
-    Float_t GetM3Trans() const            { return fM3Trans; }
-    Float_t GetSlopeLong() const          { return fSlopeLong; }
-    Float_t GetSlopeTrans() const         { return fSlopeTrans; }
-    Float_t GetTimeSpread() const         { return fTimeSpread; }
-    Float_t GetTimeSpreadWeighted() const { return fTimeSpreadWeighted; }
+    Float_t GetAsym() const                { return fAsym; }
+    Float_t GetM3Long() const              { return fM3Long; }
+    Float_t GetM3Trans() const             { return fM3Trans; }
+    Float_t GetSlopeLong() const           { return fSlopeLong; }
+    Float_t GetSlopeTrans() const          { return fSlopeTrans; }
+    Float_t GetTimeSpread() const          { return fTimeSpread; }
+    Float_t GetTimeSpreadWeighted() const  { return fTimeSpreadWeighted; }
+    Float_t GetSlopeSpread() const         { return fSlopeSpread; }
+    Float_t GetSlopeSpreadWeighted() const { return fSlopeSpreadWeighted; }
 
     Int_t Calc(const MGeomCam &geom, const MSignalCam &pix,