diff --git a/CALIBRATION/shms_hgcer_calib/calibration.C b/CALIBRATION/shms_hgcer_calib/calibration.C index 4409e417a6d6a0f3bb00d32c3b961ae4b8ee11e3..3c7c23e0cdd9b66da7b5b6a6adaf067849938732 100644 --- a/CALIBRATION/shms_hgcer_calib/calibration.C +++ b/CALIBRATION/shms_hgcer_calib/calibration.C @@ -108,7 +108,7 @@ void calibration::SlaveBegin(TTree * /*tree*/) { ADC_min = -10; ADC_max = 200; - bins = 2*(abs(ADC_min) + abs(ADC_max)); + bins = 12*(abs(ADC_min) + abs(ADC_max)); } fPulseInt = new TH1F*[4]; @@ -655,7 +655,7 @@ void calibration::Terminate() nbins = (PulseInt[ipmt]->GetXaxis()->GetNbins()); //With the scale of ADC to NPE create a histogram that has the conversion applied - fscaled[ipmt] = new TH1F(Form("fscaled_PMT%d", ipmt+1), Form("Scaled ADC spectra for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 210 : 210, -1, fNGC ? 20 : 20); + fscaled[ipmt] = new TH1F(Form("fscaled_PMT%d", ipmt+1), Form("Scaled ADC spectra for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 1260 : 1260, -1, fNGC ? 20 : 20); //Fill this histogram bin by bin for (Int_t ibin=0; ibin<nbins; ibin++) @@ -679,7 +679,7 @@ void calibration::Terminate() fFullShow ? fscaled[ipmt]->Fit("Poisson","RQ") : fscaled[ipmt]->Fit("Poisson","RQN"); //Make and fill histogram with the background removed - fscaled_nobackground[ipmt] = new TH1F(Form("fscaled_nobackground_pmt%d", ipmt+1), Form("NPE spectra background removed for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 210 : 210, -1, fNGC ? 20 : 20); + fscaled_nobackground[ipmt] = new TH1F(Form("fscaled_nobackground_pmt%d", ipmt+1), Form("NPE spectra background removed for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 1260 : 1260, -1, fNGC ? 20 : 20); for (Int_t ibin=0; ibin<nbins; ibin++) { @@ -702,7 +702,7 @@ void calibration::Terminate() Gauss3->SetParLimits(8, 0.2, 0.5); fFullShow ? fscaled_nobackground[ipmt]->Fit("Gauss3","RQ") : fscaled_nobackground[ipmt]->Fit("Gauss3","RQN"); if (fFullShow) fscaled_nobackground[ipmt]->GetXaxis()->SetRangeUser(0,5); - if (fFullShow) fscaled_nobackground[ipmt]->GetYaxis()->SetRangeUser(0,0.2); + if (fFullShow) fscaled_nobackground[ipmt]->GetYaxis()->SetRangeUser(0,0.3); //Create a TGraphErrors to determine the spacing of the NPE y_npe[0] = Gauss3->GetParameter(1), y_npe[1] = Gauss3->GetParameter(4), y_npe[2] = Gauss3->GetParameter(7); @@ -726,7 +726,7 @@ void calibration::Terminate() Double_t xscale_mk2 = xscale * Gauss3->GetParameter(1); //Take this new xscale and repeat the exact same procedure as before - fscaled_mk2[ipmt] = new TH1F(Form("fhgc_scaled_mk2_PMT%d", ipmt+1), Form("Scaled ADC spectra for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 210 : 210, -1, fNGC ? 20 : 20); + fscaled_mk2[ipmt] = new TH1F(Form("fhgc_scaled_mk2_PMT%d", ipmt+1), Form("Scaled ADC spectra for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 1260 : 1260, -1, fNGC ? 20 : 20); //Fill this histogram bin by bin for (Int_t ibin=0; ibin<nbins; ibin++) @@ -750,7 +750,7 @@ void calibration::Terminate() fFullShow ? fscaled_mk2[ipmt]->Fit("Poisson","RQ"):fscaled_mk2[ipmt]->Fit("Poisson","RQN"); //Make and fill histogram with the background removed - fscaled_mk2_nobackground[ipmt] = new TH1F(Form("fscaled_mk2_nobackground_pmt%d", ipmt+1), Form("NPE spectra background removed for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 210 : 210, -1, fNGC ? 20 : 20); + fscaled_mk2_nobackground[ipmt] = new TH1F(Form("fscaled_mk2_nobackground_pmt%d", ipmt+1), Form("NPE spectra background removed for PMT%d; NPE; Normalized Counts",ipmt+1), fNGC ? 1260 : 1260, -1, fNGC ? 20 : 20); for (Int_t ibin=0; ibin<nbins; ibin++) {