Changeset 6986 for trunk/MagicSoft/Mars/mmuon

Timestamp:

04/29/05 13:09:12 (20 years ago)

Author:

tbretz

Message:

*** empty log message ***

Location:

trunk/MagicSoft/Mars/mmuon

Files:

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

trunk/MagicSoft/Mars/mmuon/MHSingleMuon.cc

@@ -126 +126 @@
     }
 
-    MMuonSetup *setup = (MMuonSetup*)plist->FindObject("MMuonSetup");
+    MMuonSetup *setup = (MMuonSetup*)const_cast<MParList*>(plist)->FindCreateObj("MMuonSetup");
     if (!setup)
-    {
-        *fLog << warn << "MMuonSetup not found... abort." << endl;
-        return kFALSE;
-    }
+        return kFALSE;
+
     fMargin = setup->GetMargin()/fGeomCam->GetConvMm2Deg();
 
@@ -185 +183 @@
     //  Up to now, the error of pedestal is not taken into accout. This is not
     // of course correct. We will include this soon.
-    Double_t ADC2PhEl = 0.14;
-    Double_t Ffactor  = 1.4;
+    const Double_t Ffactor  = 1.4;
     for (Int_t i=0; i<fHistPhi.GetNbinsX()+1; i++)
     {
         const Float_t abs      = TMath::Abs(fHistPhi.GetBinContent(i));
-        const Float_t rougherr = TMath::Sqrt(abs/ADC2PhEl)*Ffactor;
+        const Float_t rougherr = TMath::Sqrt(abs)*Ffactor;
 
         fHistPhi.SetBinError(i, rougherr);
@@ -198 +195 @@
     {
         const Float_t abs      = TMath::Abs(fHistWidth.GetBinContent(i));
-        const Float_t rougherr = TMath::Sqrt(abs/ADC2PhEl)*Ffactor;
+        const Float_t rougherr = TMath::Sqrt(abs)*Ffactor;
 
         fHistWidth.SetBinError(i, rougherr);
@@ -298 +295 @@
             first = 0; // If maximum is on the left side of histogram
 
+    if (last-first<=3)
+        return kFALSE;
+
     // Now get the fit range
     const Float_t startfitval = fHistWidth.GetBinLowEdge(first+1);
@@ -316 +316 @@
     fHistWidth.Fit(&f1, "NQR");
 
-    if (last-first>3)
-        chi = f1.GetChisquare()/(last-first-3);
+    chi = f1.GetChisquare()/f1.GetNDF();
 
     Double_t err;

trunk/MagicSoft/Mars/mmuon/MMuonCalibParCalc.cc

@@ -187 +187 @@
 {
     // Calculation of Arc Phi etc...
-    const Float_t thresphi   = fMuonSetup->GetThresholdArcPhi()  /fGeomCam->GetConvMm2Deg();
-    const Float_t threswidth = fMuonSetup->GetThresholdArcWidth()/fGeomCam->GetConvMm2Deg();
+    const Float_t thresphi   = fMuonSetup->GetThresholdArcPhi();
+    const Float_t threswidth = fMuonSetup->GetThresholdArcWidth();
 
     Double_t peakphi, arcphi;

trunk/MagicSoft/Mars/mmuon/MMuonSetup.cc

@@ -45 +45 @@
 //
 MMuonSetup::MMuonSetup(const char *name, const char *title)
-    : fMargin(0.2), fThresholdArcPhi(0.1), fThresholdArcWidth(0.1)
+    : fMargin(0.2), fThresholdArcPhi(30), fThresholdArcWidth(30)
 {
     fName  = name  ? name  : "MMuonSetup";
@@ -55 +55 @@
 // Read the setup from a TEnv, eg:
 //   MMuonSetup.Margin:            0.2
-//   MMuonSetup.ThresholdArcPhi:   0.1
-//   MMuonSetup.ThresholdArcWidth: 0.1
+//   MMuonSetup.ThresholdArcPhi:   30
+//   MMuonSetup.ThresholdArcWidth: 30
 //
 Int_t MMuonSetup::ReadEnv(const TEnv &env, TString prefix, Bool_t print)

trunk/MagicSoft/Mars/mmuon/MMuonSetup.h

@@ -10 +10 @@
 private:
     Float_t fMargin;            // [deg] margin to evaluate muons. The defaut value is 60 mm, corresponding to 0.2 deg. This value can be changed by using the function of SetMargin
-    Float_t fThresholdArcPhi;   // [deg] The threshold value to define arc phi
-    Float_t fThresholdArcWidth; // [deg] The threshold value to define arc width
+    Float_t fThresholdArcPhi;   // [phe] The threshold value to define arc phi
+    Float_t fThresholdArcWidth; // [phe] The threshold value to define arc width
 
     Int_t ReadEnv(const TEnv &env, TString prefix, Bool_t print);