Changeset 8625 for trunk/MagicSoft/Mars/mmuon

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mmuon

Files:

• MHSingleMuon.cc (modified) (4 diffs)
• MMuonSearchPar.cc (modified) (2 diffs)
• MMuonSearchPar.h (modified) (1 diff)

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;

trunk/MagicSoft/Mars/mmuon/MMuonSearchPar.cc

@@ -197 +197 @@
     // ----------------- Setup and call minuit -------------------
     const Float_t  delta = 30.;  // 3 mm (1/10 of an inner pixel size) Step to move.
-    const Double_t r     = geom.GetMaxRadius()*2;
+    //const Double_t r     = geom.GetMaxRadius()*2;
 
     // Save gMinuit
@@ -207 +207 @@
     minuit.Command("set nowarn"); // Switch off warning
     // Define Parameters
-    minuit.DefineParameter(0, "x",     x, delta, -r, r);
-    minuit.DefineParameter(1, "y",     y, delta, -r, r);
+    minuit.DefineParameter(0, "x",     x, delta,  0, 0);//-r, r);
+    minuit.DefineParameter(1, "y",     y, delta,  0, 0);//-r, r);
     minuit.DefineParameter(2, "r",     0, 1,      0, 0);
     minuit.DefineParameter(3, "sigma", 0, 1,      0, 0);

trunk/MagicSoft/Mars/mmuon/MMuonSearchPar.h

@@ -42 +42 @@
     Float_t GetCenterX()   const { return fCenterX; }
     Float_t GetCenterY()   const { return fCenterY; }
+    Float_t GetDist()      const { return TMath::Hypot(fCenterX, fCenterY); }
     Float_t GetTime()      const { return fTime; }
     Float_t GetTimeRms()   const { return fTimeRms; }