Index: trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc
===================================================================
--- trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc	(revision 8713)
+++ trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc	(revision 8775)
@@ -1,4 +1,4 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MHSingleMuon.cc,v 1.16 2007-08-06 16:50:45 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MHSingleMuon.cc,v 1.17 2007-12-03 17:44:59 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -70,8 +70,13 @@
 // your results!!
 //
-// Inputcontainer:
+// InputContainer:
 //   - MGeomCam
 //   - MMuonSearchPar
 //
+//
+// Class Version 2:
+// ----------------
+//   + Double_t fRelTimeMean;   // Result of the gaus fit to the arrival time
+//   + Double_t fRelTimeSigma;  // Result of the gaus fit to the arrival time
 //
 ////////////////////////////////////////////////////////////////////////////
@@ -194,4 +199,7 @@
 Bool_t MHSingleMuon::Fill(const MParContainer *par, const Stat_t w)
 {
+    fRelTimeMean  =  0;
+    fRelTimeSigma = -1;
+
     fHistPhi.Reset();
     fHistWidth.Reset();
@@ -247,12 +255,11 @@
     fHistTime.Fit(&g1, "QNB");
 
-    //    Double_t err;
-    Double_t sig, mean, dummy;
-    gMinuit->GetParameter(1, mean, dummy);  // get the mean  value
-    gMinuit->GetParameter(2, sig,  dummy);  // get the sigma value
+    Double_t dummy;
+    gMinuit->GetParameter(1, fRelTimeMean,   dummy);  // get the mean  value
+    gMinuit->GetParameter(2, fRelTimeSigma,  dummy);  // get the sigma value
 
     // The mean arrival time which was subtracted before will
     // be added again, now
-    const Double_t tm0 = fMuonSearchPar->GetTime()+mean;
+    const Double_t tm0 = fMuonSearchPar->GetTime()+fRelTimeMean;
 
     for (Int_t i=0; i<entries; i++)
@@ -268,5 +275,5 @@
         // if the signal is not near the estimated circle, it is ignored.
         if (TMath::Abs(dist-fMuonSearchPar->GetRadius())<fMargin &&
-            TMath::Abs(pix.GetArrivalTime()-tm0) < 2*sig)
+            TMath::Abs(pix.GetArrivalTime()-tm0) < 2*fRelTimeSigma)
         {
             fHistPhi.Fill(TMath::ATan2(dx, dy)*TMath::RadToDeg(), pix.GetNumPhotons());