Changeset 4322 for trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

Timestamp:

06/22/04 14:24:03 (20 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc (modified) (6 diffs)

trunk/MagicSoft/Mars/mcalib/MCalibrationChargeCalc.cc

@@ -986 +986 @@
       area[0] = i;
 
-      TH1D *hist = camphes.ProjectionS(TArrayI(),area);
+      TH1D *hist = camphes.ProjectionS(TArrayI(),area,"_py",100);
       hist->Fit("gaus","Q");
       const Float_t mean  = hist->GetFunction("gaus")->GetParameter(1);
@@ -992 +992 @@
       const Int_t   ndf   = hist->GetFunction("gaus")->GetNDF();
 
-      if (ndf < 1)
+      if (ndf < 2)
         {
           *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " 
@@ -1004 +1004 @@
       const Double_t prob = hist->GetFunction("gaus")->GetProb();
 
-      if (prob < 0.0005)
+      if (prob < 0.001)
         {
           *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the number of photo-electrons " 
@@ -1394 +1394 @@
       area[0] = i;
 
-      TH1D *hist = camffactor.ProjectionS(TArrayI(),area);
+      TH1D *hist = camffactor.ProjectionS(TArrayI(),area,"_py",100);
       hist->Fit("gaus","Q");
       const Float_t mean  = hist->GetFunction("gaus")->GetParameter(1);
@@ -1400 +1400 @@
       const Int_t   ndf   = hist->GetFunction("gaus")->GetNDF();
 
-      if (ndf < 1)
+      if (ndf < 2)
         {
           *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor " 
@@ -1412 +1412 @@
       const Double_t prob = hist->GetFunction("gaus")->GetProb();
 
-      if (prob < 0.0005)
+      if (prob < 0.001)
         {
           *fLog << warn << GetDescriptor() << ": Cannot use a Gauss fit to the F-Factor "