diff --git a/CuFluo_Ka_KaKb_Comp.cpp b/CuFluo_Ka_KaKb_Comp.cpp index 39f36d2..6a7b7ff 100644 --- a/CuFluo_Ka_KaKb_Comp.cpp +++ b/CuFluo_Ka_KaKb_Comp.cpp @@ -90,11 +90,10 @@ int main(int argc, char* argv[]) { int high_bin_peak = 0; if (gain_str == "HG0") { low_bin_peak = 701; - //high_bin_peak = 1101; // was this and is still this in the main fitter - high_bin_peak = 1151; + high_bin_peak = 1200; } else if (gain_str == "G0") { low_bin_peak = 301; - high_bin_peak = 601; + high_bin_peak = 651; } for (int j=1; j<9; j++) { @@ -233,11 +232,23 @@ int main(int argc, char* argv[]) { mypar_KaKb[0] = 0.0; mypar_KaKb[1] = 0.0; mypar_KaKb[2] = proj_peak->GetBinCenter(proj_peak->GetMaximumBin()); - mypar_KaKb[3] = 15.; + if (gain_str == "G0") { + mypar_KaKb[3] = 16.; + } else if (gain_str == "HG0") { + mypar_KaKb[3] = 29.; + } mypar_KaKb[4] = proj_peak->GetBinContent(proj_peak->GetMaximumBin()); - mypar_KaKb[5] = 0.18; - mypar_KaKb[6] = 1.121; - mypar_KaKb[7] = 0.11; + if (gain_str == "G0") { + mypar_KaKb[5] = 0.17; + } else if (gain_str == "HG0") { + mypar_KaKb[5] = 0.14; + } + mypar_KaKb[6] = 1.12; + if (gain_str == "G0") { + mypar_KaKb[7] = 0.12; + } else if (gain_str == "HG0") { + mypar_KaKb[7] = 0.14; + } Double_t emypar_KaKb[8]; energyCalibration *thiscalibration_KaKb = new energyCalibration();