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

Added maps where no calib is found.

Browse Source
This commit is contained in:
redford_s
2019-02-08 16:55:25 +01:00
parent d8515bcaa6
commit 8e9dadc17a
3 changed files with 47 additions and 1 deletions
Download Patch File Download Diff File Expand all files Collapse all files

22
CS_analysis.cpp
View File

@ -196,6 +196,9 @@ int main(int argc, char* argv[]) {
TH2F* pedeRMSG0 = new TH2F("pedeRMSG0","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5);
TH2F* no_g0g1_calib_map = new TH2F("no_g0g1_calib_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5);
TH2F* no_g1g2_calib_map = new TH2F("no_g1g2_calib_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5);
char savename[128];
int createHistoFile = 0;
@ -1414,9 +1417,10 @@ int main(int argc, char* argv[]) {
if (isBulk(i)) {
g0overg1_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
}
} else {
no_g0g1_calib_map->Fill(i%NC,i/NC,1);
}
if (rangemax_g1 > rangemin_g1 && rangemax_g2 > rangemin_g2) {
g1overg2hist->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
@ -1442,6 +1446,8 @@ int main(int argc, char* argv[]) {
if (isBulk(i)) {
g1overg2_isBulk->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1));
}
} else {
no_g1g2_calib_map->Fill(i%NC,i/NC,1);
}
delete fit_g0;
@ -1788,6 +1794,20 @@ int main(int argc, char* argv[]) {
mapcanvas->SaveAs((const char *)(savename));
}
no_g0g1_calib_map->GetXaxis()->SetTitle("Column");
no_g0g1_calib_map->GetYaxis()->SetTitle("Row");
no_g0g1_calib_map->GetYaxis()->SetTitleOffset(0.7);
no_g0g1_calib_map->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/no_g0g1_calib_map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
no_g1g2_calib_map->GetXaxis()->SetTitle("Column");
no_g1g2_calib_map->GetYaxis()->SetTitle("Row");
no_g1g2_calib_map->GetYaxis()->SetTitleOffset(0.7);
no_g1g2_calib_map->Draw("colz");
sprintf(savename,"plots/M%s/CurrentSource/no_g1g2_calib_map_M%s.png", module_str.c_str(), module_str.c_str());
mapcanvas->SaveAs((const char *)(savename));
sprintf(savename,"data/M%s/CS_ratio_M%s.root", module_str.c_str(), module_str.c_str());
TFile* saved_file2 = new TFile((const char *)(savename),"RECREATE");
g1overg2map->Write();