Changeset 8625 for trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

Timestamp:

06/29/07 13:34:23 (17 years ago)

Author:

tbretz

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -1 +1 @@
 /* ======================================================================== *\
-! $Name: not supported by cvs2svn $:$Id: MHSingleMuon.cc,v 1.14 2006-10-18 10:35:52 tbretz Exp $
+! $Name: not supported by cvs2svn $:$Id: MHSingleMuon.cc,v 1.15 2007-06-29 12:29:47 tbretz Exp $
 ! --------------------------------------------------------------------------
 !
@@ -144 +144 @@
     bins.Apply(fHistWidth);
 
-    bins.SetEdges(101, -10, 10);
+    bins.SetEdges(101, -33, 33);   // +/- 33ns
     bins.Apply(fHistTime);
 }
@@ -221 +221 @@
 
         // use only the inner pixles. FIXME: This is geometry dependent
-        if(i>397)
+        if(gpix.GetAidx()>0)
             continue;
 
@@ -227 +227 @@
     }
     // Setup the function and perform the fit
-    TF1 g1("g1", "gaus", -10, 10);
+    TF1 g1("g1", "gaus");//, -fHistTime.GetXmin(), fHistTime.GetXmax());
 
     // Choose starting values as accurate as possible
     g1.SetParameter(0, fHistTime.GetMaximum());
     g1.SetParameter(1, 0);
-    g1.SetParameter(2, 0.2);
-
-    g1.SetParLimits(1, -0.5, 0.5);
-    g1.SetParLimits(2, 0, 1);
+    g1.SetParameter(2, 0.7); // FIXME! GetRMS instead???
+
+    // According to fMuonSearchPar->GetTimeRMS() identified muons
+    // do not have an arrival time rms>3
+    g1.SetParLimits(1, -1.7, 1.7);
+    g1.SetParLimits(2,  0,   3.4);
+
     // options : N  do not store the function, do not draw
     //           I  use integral of function in bin rather than value at bin center
     //           R  use the range specified in the function range
     //           Q  quiet mode
-    fHistTime.Fit(&g1, "QNRB");
+    fHistTime.Fit(&g1, "QNB");
 
     //    Double_t err;