Changeset 2822

Timestamp:

01/15/04 17:31:50 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc

@@ -110 +110 @@
 
     fgSinglePheFitFunc = &gfKto4;
-    fgSinglePheFitNPar = 5;
+    fgSinglePheFitNPar = 6;
 }
 
@@ -352 +352 @@
   const Double_t mu_1_guess = mu_0_guess + 50.;
   const Double_t si_1_guess = si_0_guess + si_0_guess;
-
-  fSinglePheFit = new TF1("SinglePheFit",fgSinglePheFitFunc,rmin,rmax,fgSinglePheFitNPar+1);
+  const Double_t norm       = entries/gkSq2Pi;
+
+  fSinglePheFit = new TF1("SinglePheFit",fgSinglePheFitFunc,rmin,rmax,fgSinglePheFitNPar);
   //  fSinglePheFit = new TF1("SinglePheFit",fgSinglePheFitFunc,rmin,rmax,fgSinglePheFitNPar+1);
   //  fSinglePheFit->SetParameters(lambda_guess,mu_0_guess,mu_1_guess,si_0_guess,si_1_guess);
-  fSinglePheFit->SetParameters(lambda_guess,mu_0_guess,mu_1_guess,si_0_guess,si_1_guess,entries);
+  fSinglePheFit->SetParameters(lambda_guess,mu_0_guess,mu_1_guess,si_0_guess,si_1_guess,norm);
   //  fSinglePheFit->SetParNames("#lambda","#mu_0","#mu_1","#sigma_0","#sigma_1");
-  fSinglePheFit->SetParNames("#lambda","#mu_0","#mu_1","#sigma_0","#sigma_1","area");
+  fSinglePheFit->SetParNames("#lambda","#mu_0","#mu_1","#sigma_0","#sigma_1","Area");
   fSinglePheFit->SetParLimits(0,0.,1.);
   fSinglePheFit->SetParLimits(1,rmin,(rmax-rmin)/1.5);
@@ -364 +365 @@
   fSinglePheFit->SetParLimits(3,1.0,(rmax-rmin)/2.0);
   fSinglePheFit->SetParLimits(4,1.0,(rmax-rmin)/2.5);
+  fSinglePheFit->SetParLimits(5,norm-0.1,norm+0.1);
   //  fSinglePheFit->SetParLimits(5,entries/gkSq2Pi,entries/gkSq2Pi);
   //  fSinglePheFit->SetParLimits(5,0.,1.5*entries);
@@ -373 +375 @@
   //  fHBlindPixelCharge->Scale(gkSq2Pi*(float)bins/npx/entries);
   // fHBlindPixelCharge->Scale(gkSq2Pi/entries);
-  Float_t  norm   = entries/gkSq2Pi;
   //           norm  *= (Float_t)fSinglePheFit->GetNpx()/(Float_t)fBlindPixelChargenbins;
-  fSinglePheFit->SetParLimits(5,norm-0.1,norm+0.1);
 
   fHBlindPixelCharge->Fit(fSinglePheFit,opt);