Changeset 7367

Timestamp:

09/27/05 14:40:03 (19 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mmuon

Files:

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

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -240 +240 @@
 
     //    Double_t err;
-    Double_t mean, meanerr, dummy;
+    Double_t sig, mean, dummy;
+    gMinuit->GetParameter(1, mean, dummy);  // get the mean  value
     gMinuit->GetParameter(2, sig,  dummy);  // get the sigma value
-    gMinuit->GetParameter(1, mean, dummy);  // get the sigma value
 
     for (Int_t i=0; i<entries; i++)

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.h

@@ -24 +24 @@
     Double_t fMargin;               //!
 
-    TProfile fHistPhi;   // Histogram of photon distribution along the arc.
-    TProfile fHistWidth; // Histogram of radial photon distribution of the arc.
+    TProfile fHistPhi;    // Histogram of photon distribution along the arc.
+    TProfile fHistWidth;  // Histogram of radial photon distribution of the arc.
     TH1F     fHistTime;   // Histogram of arrival time distribution along the arc.
 
@@ -50 +50 @@
 
 #endif
-
-

trunk/MagicSoft/Mars/mmuon/MMuonSearchPar.cc

@@ -38 +38 @@
 // infomation stored in
 //
-// MMuonCalibPar.
+//   MMuonCalibPar.
 //
 // The information will be available by using the task of
 //
-// MMuonCalibParCalc.
-//
+//   MMuonCalibParCalc.
+//
+// Version 2:
+// ----------
+//   + Float_t fTime;      // Mean arrival time of core pixels
+//   + Float_t fTimeRms;   // Rms of arrival time of core pixels
 //
 //  Input Containers:
@@ -90 +94 @@
     fCenterY   =  0;
     fTime      =  0;
+    fTimeRms   = -1;
 }
 
@@ -144 +149 @@
 {
     // ------- Make a temporaray copy of the signal ---------
+    // ------- and calculate arrival time parameters --------
     const Int_t n = geom.GetNumPixels();
 
@@ -153 +159 @@
     Int_t q=0;
 
-    Float_t fTime2=0;
+    Float_t mean=0;
+    Float_t sq  =0;
 
     for (int i=0; i<n; i++)
@@ -168 +175 @@
 
             //timing
-            if(pix.IsPixelCore())
-            {
-                fTime  += pix.GetArrivalTime();
-                fTime2 += pix.GetArrivalTime()*pix.GetArrivalTime();
-                q++;
-            }
+            if (!pix.IsPixelCore())
+                continue;
+
+            mean += pix.GetArrivalTime();
+            sq   += pix.GetArrivalTime()*pix.GetArrivalTime();
+            q++;
         }
     }
-    fTime    = fTime/q;
-    fTime2   = fTime2/q;
-    fTimeRms = TMath::Sqrt(fTime2-(fTime*fTime));
+
+    mean /= q;
+    sq   /= q;
+
+    fTime    = mean;
+    fTimeRms = TMath::Sqrt(sq-mean*mean);
 
     fSignal.Set(p);

trunk/MagicSoft/Mars/mmuon/MMuonSearchPar.h

@@ -21 +21 @@
     Float_t fCenterX;   // An estimated center position in X of the muon ring [mm]
     Float_t fCenterY;   // An estimated center position in Y of the muon ring [mm]
-    Float_t fTime;
-    Float_t fTimeRms;
+    Float_t fTime;      // Mean arrival time of core pixels
+    Float_t fTimeRms;   // Rms of arrival time of core pixels
 
     MArrayF fSignal;    //! Temporary storage for signal
@@ -57 +57 @@
     void   Print(const MGeomCam &geom, Option_t *opt=NULL) const;
 
-    ClassDef(MMuonSearchPar, 1) // Container to hold muon search parameters
+    ClassDef(MMuonSearchPar, 2) // Container to hold muon search parameters
 };