Changeset 6871

Timestamp:

03/21/05 16:37:00 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars

Files:

• Changelog (modified) (1 diff)
• mimage/MHillas.cc (modified) (1 diff)
• mimage/MHillas.h (modified) (1 diff)

trunk/MagicSoft/Mars/Changelog

@@ -68 +68 @@
      - for some output there was a replacement of pix.GetPixId()
        by hist.GetName() <undocumented>
+
+   * mimage/MHillas.[h,cc]:
+     - added new member function GetBorderLine returning the border
+       line correspondig to an ellipse with axis lnegth and width
+     - added new member function GetArea returning the area
+       correspondig to an ellipse with axis lnegth and width
 
 

trunk/MagicSoft/Mars/mimage/MHillas.cc

@@ -118 +118 @@
 // --------------------------------------------------------------------------
 //
+//  Analytical estimation of length of border line:
+//    U = pi*(a+b - sqrt(a*b))
+//
+//  GetBorderLine() [mm]
+//
+Double_t MHillas::GetBorderLine() const
+{
+    const Double_t a = fWidth+fLength;
+    const Double_t s = fWidth*fLength;
+
+    return TMath::Pi()*(1.5*a - TMath::Sqrt(s));
+}
+
+// --------------------------------------------------------------------------
+//
+//  Analytical estimation of length of border line:
+//    A = pi*a*b
+//
+//  GetArea() [mm^2]
+//
+Double_t MHillas::GetArea() const
+{
+    return TMath::Pi()*fWidth*fLength;
+}
+
+// --------------------------------------------------------------------------
+//
 // Print the hillas Parameters to *fLog
 //

trunk/MagicSoft/Mars/mimage/MHillas.h

@@ -47 +47 @@
     TVector2 GetMean() const;
 
+    Double_t GetBorderLine() const;
+    Double_t GetArea() const;
+
     Float_t GetCosDelta() const { return fCosDelta; }
     Float_t GetSinDelta() const { return fSinDelta; }