Index: /trunk/MagicSoft/Mars/Changelog
===================================================================
--- /trunk/MagicSoft/Mars/Changelog	(revision 6870)
+++ /trunk/MagicSoft/Mars/Changelog	(revision 6871)
@@ -68,4 +68,10 @@
      - 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
 
 
Index: /trunk/MagicSoft/Mars/mimage/MHillas.cc
===================================================================
--- /trunk/MagicSoft/Mars/mimage/MHillas.cc	(revision 6870)
+++ /trunk/MagicSoft/Mars/mimage/MHillas.cc	(revision 6871)
@@ -118,4 +118,31 @@
 // --------------------------------------------------------------------------
 //
+//  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
 //
Index: /trunk/MagicSoft/Mars/mimage/MHillas.h
===================================================================
--- /trunk/MagicSoft/Mars/mimage/MHillas.h	(revision 6870)
+++ /trunk/MagicSoft/Mars/mimage/MHillas.h	(revision 6871)
@@ -47,4 +47,7 @@
     TVector2 GetMean() const;
 
+    Double_t GetBorderLine() const;
+    Double_t GetArea() const;
+
     Float_t GetCosDelta() const { return fCosDelta; }
     Float_t GetSinDelta() const { return fSinDelta; }