Changeset 2920 for trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc

Timestamp:

01/26/04 19:32:58 (21 years ago)

Author:

gaug

Message:

*** empty log message ***

File:

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

trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc

@@ -260 +260 @@
 {
 
-    gStyle->SetOptFit(0);
+    gStyle->SetOptFit(1);
     gStyle->SetOptStat(111111);
 
@@ -366 +366 @@
       break;
     case kEPolya:
-      break;
-
+      fSinglePheFit = new TF1("SinglePheFit",&fPolya,min,max,8);
+      break;
     case kEMichele:
       break;
@@ -381 +381 @@
   const Double_t mu_1_guess = mu_0_guess + 50.;
   const Double_t si_1_guess = si_0_guess + si_0_guess;
+  // Michele
+  const Double_t lambda_1cat_guess = 0.5;
+  const Double_t lambda_1dyn_guess = 0.5;
+  const Double_t mu_1cat_guess = mu_0_guess + 50.;
+  const Double_t mu_1dyn_guess = mu_0_guess + 20.;
+  const Double_t si_1cat_guess = si_0_guess + si_0_guess;
+  const Double_t si_1dyn_guess = si_0_guess;
+  // Polya
+  const Double_t excessPoisson_guess = 0.5;
+  const Double_t delta1_guess     = 8.;
+  const Double_t delta2_guess     = 5.;
+  const Double_t electronicAmp_guess  = 0.0025;
+
   //
   // Initialize boundaries and start parameters
@@ -389 +402 @@
     case kEPoisson4:
       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","Area");
+      fSinglePheFit->SetParNames("#lambda","#mu_{0}","#mu_{1}","#sigma_{0}","#sigma_{1}","Area");
       fSinglePheFit->SetParLimits(0,0.,1.);
       fSinglePheFit->SetParLimits(1,-2.,2.);
@@ -400 +413 @@
     case kEPoisson6:
       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","Area");
+      fSinglePheFit->SetParNames("#lambda","#mu_{0}","#mu_{1}","#sigma_{0}","#sigma_{1}","Area");
       fSinglePheFit->SetParLimits(0,0.,1.);
       fSinglePheFit->SetParLimits(1,min,(max-min)/1.5);
@@ -410 +423 @@
 
     case kEPolya:
+      fSinglePheFit->SetParameters(lambda_guess, excessPoisson_guess,
+                                   delta1_guess,delta2_guess,
+                                   electronicAmp_guess,
+                                   si_0_guess,
+                                   norm, mu_0_guess);
+      fSinglePheFit->SetParNames("#lambda","b_{tot}",
+                                 "#delta_{1}","#delta_{2}",
+                                 "amp_{e}","#sigma_{0}",
+                                 "Area", "#mu_{0}");
+      fSinglePheFit->SetParLimits(0,0.,1.);
+      fSinglePheFit->SetParLimits(1,0.,1.); 
+      fSinglePheFit->SetParLimits(2,6.,12.);    
+      fSinglePheFit->SetParLimits(3,3.,8.);    
+      fSinglePheFit->SetParLimits(4,0.,0.005);    
+      fSinglePheFit->SetParLimits(5,min,(max-min)/1.5);
+      fSinglePheFit->SetParLimits(6,norm-0.1,norm+0.1);
+      fSinglePheFit->SetParLimits(7,-35.,15.);
       break;
 
     case kEMichele:
+
+
       break;
 
@@ -450 +482 @@
       fSigma1Err = fSinglePheFit->GetParError(4);
       break;
+    case kEPolya:
+      fLambda =  fSinglePheFit->GetParameter(0);
+      fMu0    =  fSinglePheFit->GetParameter(7);
+      fMu1    = 0.;
+      fSigma0 =  fSinglePheFit->GetParameter(5);
+      fSigma1 = 0.;
+
+      fLambdaErr = fSinglePheFit->GetParError(0);
+      fMu0Err    = fSinglePheFit->GetParError(7);
+      fMu1Err    = 0.;
+      fSigma0Err = fSinglePheFit->GetParError(5);
+      fSigma1Err = 0.;
     default:
       break;
@@ -479 +523 @@
 
   // Perform the cross-check fitting only the pedestal:
-  fSinglePhePedFit = new TF1("GausPed","gaus",rmin,fMu0);
+  fSinglePhePedFit = new TF1("GausPed","gaus",rmin,0.);
   fHBlindPixelCharge->Fit(fSinglePhePedFit,opt);