Changeset 3125 for trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc
- Timestamp:
- 02/12/04 19:56:16 (21 years ago)
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/MagicSoft/Mars/mcalib/MHCalibrationBlindPixel.cc
r3120 r3125 687 687 { 688 688 689 *fLog << err << "min: " << min << endl; 690 *fLog << err << "max: " << max << endl; 691 689 692 // 690 693 // First guesses for the fit (should be as close to reality as possible, … … 692 695 // 693 696 const Stat_t entries = fHBlindPixelCharge->Integral("width"); 694 const Double_t lambda_guess = 0. 5;697 const Double_t lambda_guess = 0.1; 695 698 const Double_t maximum_bin = fHBlindPixelCharge->GetBinCenter(fHBlindPixelCharge->GetMaximumBin()); 696 699 const Double_t norm = entries/gkSq2Pi; … … 740 743 const Double_t electronicAmp_limit = fgBlindPixelElectronicAmpError; 741 744 745 *fLog << err << "pedestal: " << fMeanPedestal << endl; 746 *fLog << err << "sigma: " << fSigmaPedestal << endl; 747 742 748 // 743 749 // Initialize boundaries and start parameters … … 754 760 fSinglePheFit->SetParNames("#lambda","#mu_{0}","#mu_{1}","#sigma_{0}","#sigma_{1}","Area"); 755 761 756 fSinglePheFit->SetParLimits(0,0., 1.);762 fSinglePheFit->SetParLimits(0,0.,0.5); 757 763 758 764 if ((fMeanPedestal) && (fSigmaPedestal)) 759 765 fSinglePheFit->SetParLimits(1, 760 fMeanPedestal- 1.*fMeanPedestalErr,761 fMeanPedestal+ 1.*fMeanPedestalErr);766 fMeanPedestal-5.*fMeanPedestalErr, 767 fMeanPedestal+5.*fMeanPedestalErr); 762 768 else 763 769 fSinglePheFit->SetParLimits(1,-3.,0.); 764 770 765 fSinglePheFit->SetParLimits(2, (max-min)/2.,max);771 fSinglePheFit->SetParLimits(2,min,max); 766 772 767 773 if ((fMeanPedestal) && (fSigmaPedestal)) 768 774 fSinglePheFit->SetParLimits(3, 769 fSigmaPedestal- 3.*fSigmaPedestalErr,770 fSigmaPedestal+ 3.*fSigmaPedestalErr);775 fSigmaPedestal-5.*fSigmaPedestalErr, 776 fSigmaPedestal+5.*fSigmaPedestalErr); 771 777 else 772 778 fSinglePheFit->SetParLimits(3,1.0,(max-min)/2.0); 773 779 774 fSinglePheFit->SetParLimits(4, 1.0,(max-min));775 fSinglePheFit->SetParLimits(5,norm- 0.5,norm+0.5);780 fSinglePheFit->SetParLimits(4,0.,(max-min)); 781 fSinglePheFit->SetParLimits(5,norm-(0.5*norm),norm+(0.5*norm)); 776 782 break; 777 783 case kEPoisson5: … … 894 900 rmin = (rmin != 0.) ? rmin : fBlindPixelChargefirst; 895 901 rmax = (rmax != 0.) ? rmax : fBlindPixelChargelast; 896 897 902 if (!InitFit(rmin,rmax)) 898 903 return kFALSE; … … 956 961 { 957 962 958 Int_t nbins = 30;963 Int_t nbins = 25; 959 964 960 965 CutEdges(fHBlindPixelCharge,nbins); … … 1002 1007 1003 1008 } 1009
Note:
See TracChangeset
for help on using the changeset viewer.