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

Timestamp:

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

Author:

tbretz

Message:

*** empty log message ***

File:

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

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;