detectors/JFCalibration
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Copied over best starting parameters from main fit file. Comparison results more stable now, but were not bad before.

Browse Source
This commit is contained in:
redford_s
2018-01-12 15:57:08 +01:00
parent 3f1e4497f4
commit 9b69367fc5
1 changed files with 18 additions and 7 deletions
Download Patch File Download Diff File Expand all files Collapse all files

25
CuFluo_Ka_KaKb_Comp.cpp
View File

@ -90,11 +90,10 @@ int main(int argc, char* argv[]) {
int high_bin_peak = 0; int high_bin_peak = 0;
if (gain_str == "HG0") { if (gain_str == "HG0") {
low_bin_peak = 701; low_bin_peak = 701;
//high_bin_peak = 1101; // was this and is still this in the main fitter high_bin_peak = 1200;
high_bin_peak = 1151;
} else if (gain_str == "G0") { } else if (gain_str == "G0") {
low_bin_peak = 301; low_bin_peak = 301;
high_bin_peak = 601; high_bin_peak = 651;
} }
for (int j=1; j<9; j++) { for (int j=1; j<9; j++) {
@ -233,11 +232,23 @@ int main(int argc, char* argv[]) {
mypar_KaKb[0] = 0.0; mypar_KaKb[0] = 0.0;
mypar_KaKb[1] = 0.0; mypar_KaKb[1] = 0.0;
mypar_KaKb[2] = proj_peak->GetBinCenter(proj_peak->GetMaximumBin()); 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[4] = proj_peak->GetBinContent(proj_peak->GetMaximumBin());
mypar_KaKb[5] = 0.18; if (gain_str == "G0") {
mypar_KaKb[6] = 1.121; mypar_KaKb[5] = 0.17;
mypar_KaKb[7] = 0.11; } 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]; Double_t emypar_KaKb[8];
energyCalibration *thiscalibration_KaKb = new energyCalibration(); energyCalibration *thiscalibration_KaKb = new energyCalibration();