Changeset 3125 for trunk/MagicSoft/Mars/mcalib

Timestamp:

02/12/04 19:56:16 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

• trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc (modified) (7 diffs)

trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc

@@ -687 +687 @@
 {
   
+  *fLog << err << "min: " << min << endl;
+  *fLog << err << "max: " << max << endl;  
+
   //
   // First guesses for the fit (should be as close to reality as possible, 
@@ -692 +695 @@
   //
   const Stat_t   entries      = fHBlindPixelCharge->Integral("width");
-  const Double_t lambda_guess = 0.5;
+  const Double_t lambda_guess = 0.1;
   const Double_t maximum_bin  = fHBlindPixelCharge->GetBinCenter(fHBlindPixelCharge->GetMaximumBin());
   const Double_t norm         = entries/gkSq2Pi;
@@ -740 +743 @@
   const Double_t electronicAmp_limit  = fgBlindPixelElectronicAmpError;
 
+  *fLog << err << "pedestal: " << fMeanPedestal << endl;
+  *fLog << err << "sigma: "    << fSigmaPedestal << endl;  
+
   //
   // Initialize boundaries and start parameters
@@ -754 +760 @@
       fSinglePheFit->SetParNames("#lambda","#mu_{0}","#mu_{1}","#sigma_{0}","#sigma_{1}","Area");
 
-      fSinglePheFit->SetParLimits(0,0.,1.);
+      fSinglePheFit->SetParLimits(0,0.,0.5);
 
       if ((fMeanPedestal) && (fSigmaPedestal))
         fSinglePheFit->SetParLimits(1,
-                                    fMeanPedestal-1.*fMeanPedestalErr,
-                                    fMeanPedestal+1.*fMeanPedestalErr);
+                                    fMeanPedestal-5.*fMeanPedestalErr,
+                                    fMeanPedestal+5.*fMeanPedestalErr);
       else
         fSinglePheFit->SetParLimits(1,-3.,0.);
 
-      fSinglePheFit->SetParLimits(2,(max-min)/2.,max);
+      fSinglePheFit->SetParLimits(2,min,max);
 
       if ((fMeanPedestal) && (fSigmaPedestal))
         fSinglePheFit->SetParLimits(3,
-                                    fSigmaPedestal-3.*fSigmaPedestalErr,
-                                    fSigmaPedestal+3.*fSigmaPedestalErr);
+                                    fSigmaPedestal-5.*fSigmaPedestalErr,
+                                    fSigmaPedestal+5.*fSigmaPedestalErr);
       else
         fSinglePheFit->SetParLimits(3,1.0,(max-min)/2.0);
 
-      fSinglePheFit->SetParLimits(4,1.0,(max-min));
-      fSinglePheFit->SetParLimits(5,norm-0.5,norm+0.5);
+      fSinglePheFit->SetParLimits(4,0.,(max-min));
+      fSinglePheFit->SetParLimits(5,norm-(0.5*norm),norm+(0.5*norm));
       break;
     case kEPoisson5:
@@ -894 +900 @@
   rmin = (rmin != 0.) ? rmin : fBlindPixelChargefirst;
   rmax = (rmax != 0.) ? rmax : fBlindPixelChargelast;
-
   if (!InitFit(rmin,rmax))
     return kFALSE;
@@ -956 +961 @@
 {
 
-  Int_t nbins = 30;
+  Int_t nbins = 25;
 
   CutEdges(fHBlindPixelCharge,nbins);
@@ -1002 +1007 @@
 
 }
+