Changeset 1434 for trunk/MagicSoft/Mars/manalysis/MHillasSrc.cc

Timestamp:

07/25/02 09:14:36 (22 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/manalysis/MHillasSrc.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/manalysis/MHillasSrc.cc

@@ -19 +19 @@
 !   Author(s): Harald Kornmayer 1/2001
 !   Author(s): Rudolf Bock     10/2001 <mailto:Rudolf.Bock@cern.ch>
+!   Author(s): Wolfgang Wittek 06/2002 <mailto:wittek@mppmu.mpg.de>
 !
 !   Copyright: MAGIC Software Development, 2000-2002
@@ -32 +33 @@
 //
 //    source-dependent image parameters
-// fAlpha    angle between major axis and line source-to-center
-// fDist     distance from source to center of ellipse
+//
+// Version 1:
+// ----------
+//  fAlpha          angle between major axis and line source-to-center
+//  fDist           distance from source to center of ellipse
+//
+// Version 2:
+// ----------
+//  fHeadTail
+//
+// Version 3:
+// ----------
+//  fCosDeltaLenth  cosine of angle between d and a, where
+//                   - d is the vector from the source position to the
+//                     center of the ellipse
+//                   - a is a vector along the main axis of the ellipse
 //
 /////////////////////////////////////////////////////////////////////////////
@@ -57 +72 @@
 }
 
+void MHillasSrc::Reset()
+{
+    fDist          = -1;
+    fAlpha         =  0;
+    fHeadTail      =  0;
+    fCosDeltaAlpha =  0;
+}
+
 // --------------------------------------------------------------------------
 //
-//  calculation of source-dependent parameters
+//  Calculation of source-dependent parameters
+//  In case you don't call Calc from within an eventloop make sure, that
+//  you call the Reset member function before.
 //
-void MHillasSrc::Calc(const MHillas *hillas)
+Bool_t MHillasSrc::Calc(const MHillas *hillas)
 {
-    const Double_t mx = hillas->GetMeanX();        // [mm]
-    const Double_t my = hillas->GetMeanY();        // [mm]
+    fHillas = hillas;
 
-    const Double_t sx = mx - fSrcPos->GetX();      // [mm]
-    const Double_t sy = my - fSrcPos->GetY();      // [mm]
+    const Double_t mx   = GetMeanX();            // [mm]
+    const Double_t my   = GetMeanY();            // [mm]
 
-    const Double_t sd = sin(hillas->GetDelta());   // [1]
-    const Double_t cd = cos(hillas->GetDelta());   // [1]
+    const Double_t sx   = mx - fSrcPos->GetX();  // [mm]
+    const Double_t sy   = my - fSrcPos->GetY();  // [mm]
 
+    const Double_t sd   = sin(GetDelta());       // [1]
+    const Double_t cd   = cos(GetDelta());       // [1]
 
-    fHeadTail = cd*sx + sd*sy;                     // [mm]
-    fDist     = sqrt(sx*sx + sy*sy);               // [mm]
-    fAlpha    = atan((cd*sy - sd*sx)/fHeadTail);   // [rad]
-    fAlpha   *= kRad2Deg;                          // [deg]
+    const Double_t tand = tan(GetDelta());       // [1]
 
-    fHillas   = hillas;
+    fHeadTail = cd*sx + sd*sy;                   // [mm]
+
+    //
+    // Distance from source position to center of ellipse.
+    // If the distance is 0 distance, Alpha is not specified.
+    // The calculation has failed and returnes kFALSE.
+    //
+    Double_t dist = sqrt(sx*sx + sy*sy);         // [mm]
+
+    if (dist==0)
+    {
+        *fLog << warn << GetDescriptor() << ": Event has Dist==0... skipped." << endl;
+        return kFALSE;
+    }
+
+    fDist = dist;
+                                                 // [mm]
+    //
+    // Calculate Alpha and Cosda = cos(d,a)
+    // The sign of Cosda will be used for quantities containing 
+    // a head-tail information
+    //
+    const Double_t arg = (sy-tand*sx) / (fDist*sqrt(tand*tand+1));
+
+    fAlpha         = asin(arg)*kRad2Deg;        // [deg]
+    fCosDeltaAlpha = fHeadTail/fDist;           // [1]
 
     SetReadyToSave();
+
+    return kTRUE;
 } 
 
@@ -87 +137 @@
     *fLog << all;
     *fLog << "Source dependant Image Parameters (" << GetName() << ")" << endl;
-    *fLog << " - Dist     = " << fDist     << " mm"  << endl;
-    *fLog << " - Alpha    = " << fAlpha    << " deg" << endl;
-    *fLog << " - HeadTail = " << fHeadTail << " mm"  << endl;
+    *fLog << " - Dist          [mm]  = " << fDist << endl;
+    *fLog << " - Alpha         [deg] = " << fAlpha << endl;
+    *fLog << " - HeadTail      [mm]  = " << fHeadTail << endl;
+    *fLog << " - CosDeltaAlpha       = " << fCosDeltaAlpha << endl;
 }
 
@@ -96 +147 @@
 // overloaded MParContainer to read MHillasSrc from an ascii file
 //
+/*
 void MHillasSrc::AsciiRead(ifstream &fin)
 {
@@ -102 +154 @@
     fin >> fHeadTail;
 }
-
+*/
 // -----------------------------------------------------------------------
 //