Changeset 8775 for trunk/MagicSoft/Mars/mmuon

Timestamp:

12/03/07 17:44:59 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc (modified) (5 diffs)

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $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 +70 @@
 // 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 +199 @@
 Bool_t MHSingleMuon::Fill(const MParContainer *par, const Stat_t w)
 {
+    fRelTimeMean  =  0;
+    fRelTimeSigma = -1;
+
     fHistPhi.Reset();
     fHistWidth.Reset();
@@ -247 +255 @@
     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 +275 @@
         // 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());