From 8d65c4fe4e6d2fc1a760b1a434073648b0fc8532 Mon Sep 17 00:00:00 2001 From: redford_s Date: Tue, 4 Apr 2017 10:39:22 +0200 Subject: [PATCH] Better plot names and scales, to be consistent with other analyses. --- JFMC_CurrentSourceScan.cpp | 273 ++++++++++++++++++------------------- 1 file changed, 133 insertions(+), 140 deletions(-) diff --git a/JFMC_CurrentSourceScan.cpp b/JFMC_CurrentSourceScan.cpp index 4bf7473..97dc56c 100644 --- a/JFMC_CurrentSourceScan.cpp +++ b/JFMC_CurrentSourceScan.cpp @@ -355,40 +355,40 @@ int main(int argc, char* argv[]) { cout << "frame_counter " << frame_counter << endl; - TH1F *g1_avg; - if (module_str == "032") { - g1_avg = new TH1F("g1_avg","",100,-12000,-4000); - } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { - g1_avg = new TH1F("g1_avg","",100,-3000,0); - } - TH1F *g2_avg; - if (module_str == "032") { - g2_avg = new TH1F("g2_avg","",100,-1000,-200); - } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { - g2_avg = new TH1F("g2_avg","",100,-300,0); - } + TH2F *g1map = new TH2F("g1map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g1ermap = new TH2F("g1ermap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH2F *g1map_avg = new TH2F("g1map_avg","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH2F *g2map_avg = new TH2F("g2map_avg","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - - TH1F *g1overg2 = new TH1F("g1overg2","",100,11,15); - TH1F *g1overg2er = new TH1F("g1overg2er","",100,0,2); - - TH2F *g1overg2_map = new TH2F("g1overg2_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH2F *g1overg2er_map = new TH2F("g1overg2er_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g2map = new TH2F("g2map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g2ermap = new TH2F("g2ermap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH1F *g1_fit_er = new TH1F("g1_fit_er","",100,0,50); - TH1F *g2_fit_er = new TH1F("g2_fit_er","",100,0,10); - - TH2F *g1map_fit_er = new TH2F("g1map_fit_er","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH2F *g2map_fit_er = new TH2F("g2map_fit_er","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g1overg2map = new TH2F("g1overg2map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g1overg2ermap = new TH2F("g1overg2ermap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); - TH1F *diffFromLinG1 = new TH1F("diffFromLinG1","",100,-0.05,0.05); - TH1F *diffFromLinG2 = new TH1F("diffFromLinG2","",100,-0.5,0.5); - TH2F *diffFromLinG1map = new TH2F("diffFromLinG1map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); TH2F *diffFromLinG2map = new TH2F("diffFromLinG2map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH1F *g1hist; + if (module_str == "032") { + g1hist = new TH1F("g1hist","",100,-12000,-4000); + } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { + g1hist = new TH1F("g1hist","",100,-3000,0); + } + TH1F *g1erhist = new TH1F("g1erhist","",100,0,50); + + TH1F *g2hist; + if (module_str == "032") { + g2hist = new TH1F("g2hist","",100,-1000,-200); + } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { + g2hist = new TH1F("g2hist","",100,-300,0); + } + TH1F *g2erhist = new TH1F("g2erhist","",100,0,10); + + TH1F *g1overg2hist = new TH1F("g1overg2hist","",100,11,15); + TH1F *g1overg2erhist = new TH1F("g1overg2erhist","",100,0,1); + + TH1F *diffFromLinG1hist = new TH1F("diffFromLinG1hist","",100,-0.05,0.05); + TH1F *diffFromLinG2hist = new TH1F("diffFromLinG2hist","",100,-0.1,0.1); + TH1F *g1overg2_isEdge = new TH1F("g1overg2_isEdge","",100,11,15); TH1F *g1overg2_isInnerEdge = new TH1F("g1overg2_isInnerEdge","",100,11,15); TH1F *g1overg2_isDouble = new TH1F("g1overg2_isDouble","",100,11,15); @@ -466,16 +466,16 @@ int main(int argc, char* argv[]) { fit2->SetParName(1,"G1 grad"); grap_g1->Fit(fit2,"QR",""); - g1_avg->Fill(fit2->GetParameter(1)); - g1map_avg->Fill(i%NC,i/NC,fit2->GetParameter(1)); - g1_fit_er->Fill(fit2->GetParError(1)); - g1map_fit_er->Fill(i%NC,i/NC,fit2->GetParError(1)); + g1hist->Fill(fit2->GetParameter(1)); + g1map->Fill(i%NC,i/NC,fit2->GetParameter(1)); + g1erhist->Fill(fit2->GetParError(1)); + g1ermap->Fill(i%NC,i/NC,fit2->GetParError(1)); double g1residsum = 0.; double g1residn = 0.; for (size_t it = 0; it < r1_adc.size(); it++) { if (r1_filter[it] >= rangemin2) { - diffFromLinG1->Fill((r1_adc[it] - fit2->Eval(r1_filter[it])) / r1_adc[it]); + diffFromLinG1hist->Fill((r1_adc[it] - fit2->Eval(r1_filter[it])) / r1_adc[it]); g1residsum = g1residsum + ((r1_adc[it] - fit2->Eval(r1_filter[it])) / r1_adc[it]); g1residn = g1residn + 1.; } @@ -501,15 +501,15 @@ int main(int argc, char* argv[]) { fit3->SetParName(1,"G2 grad"); grap_g2->Fit(fit3,"QR+",""); - g2_avg->Fill(fit3->GetParameter(1)); - g2map_avg->Fill(i%NC,i/NC,fit3->GetParameter(1)); - g2_fit_er->Fill(fit3->GetParError(1)); - g2map_fit_er->Fill(i%NC,i/NC,fit3->GetParError(1)); + g2hist->Fill(fit3->GetParameter(1)); + g2map->Fill(i%NC,i/NC,fit3->GetParameter(1)); + g2erhist->Fill(fit3->GetParError(1)); + g2ermap->Fill(i%NC,i/NC,fit3->GetParError(1)); double g2residsum = 0.; double g2residn = 0.; for (size_t it = 0; it < r3_adc.size(); it++) { - diffFromLinG2->Fill((r3_adc[it] - fit3->Eval(r3_filter[it])) / r3_adc[it]); + diffFromLinG2hist->Fill((r3_adc[it] - fit3->Eval(r3_filter[it])) / r3_adc[it]); g2residsum = g2residsum + ((r3_adc[it] - fit3->Eval(r3_filter[it])) / r3_adc[it]); g2residn = g2residn + 1.; } @@ -517,10 +517,10 @@ int main(int argc, char* argv[]) { } if (r1_adc.size() > 1 && r3_adc.size() > 1) { - g1overg2->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - g1overg2er->Fill(abs(fit2->GetParameter(1)/fit3->GetParameter(1))*sqrt(pow((fit2->GetParError(1)/fit2->GetParameter(1)),2) + pow((fit3->GetParError(1)/fit3->GetParameter(1)),2))); - g1overg2_map->Fill(i%NC,i/NC,fit2->GetParameter(1) / fit3->GetParameter(1)); - g1overg2er_map->Fill(i%NC,i/NC,abs(fit2->GetParameter(1)/fit3->GetParameter(1))*sqrt(pow((fit2->GetParError(1)/fit2->GetParameter(1)),2) + pow((fit3->GetParError(1)/fit3->GetParameter(1)),2))); + g1overg2hist->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); + g1overg2erhist->Fill(abs(fit2->GetParameter(1)/fit3->GetParameter(1))*sqrt(pow((fit2->GetParError(1)/fit2->GetParameter(1)),2) + pow((fit3->GetParError(1)/fit3->GetParameter(1)),2))); + g1overg2map->Fill(i%NC,i/NC,fit2->GetParameter(1) / fit3->GetParameter(1)); + g1overg2ermap->Fill(i%NC,i/NC,abs(fit2->GetParameter(1)/fit3->GetParameter(1))*sqrt(pow((fit2->GetParError(1)/fit2->GetParameter(1)),2) + pow((fit3->GetParError(1)/fit3->GetParameter(1)),2))); if (isEdge(i)) { g1overg2_isEdge->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); @@ -645,68 +645,68 @@ int main(int argc, char* argv[]) { c1->cd(); - g1_avg->GetXaxis()->SetTitle("G1 gradient"); - g1_avg->GetXaxis()->SetTitleOffset(1.1); - g1_avg->GetYaxis()->SetTitleOffset(1.5); - g1_avg->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g1_avg.png", module_str.c_str()); + g1hist->GetXaxis()->SetTitle("G1 gradient"); + g1hist->GetXaxis()->SetTitleOffset(1.1); + g1hist->GetYaxis()->SetTitleOffset(1.5); + g1hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g1hist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - g2_avg->GetXaxis()->SetTitle("G2 gradient"); - g2_avg->GetXaxis()->SetTitleOffset(1.1); - g2_avg->GetYaxis()->SetTitleOffset(1.5); - g2_avg->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g2_avg.png", module_str.c_str()); + g2hist->GetXaxis()->SetTitle("G2 gradient"); + g2hist->GetXaxis()->SetTitleOffset(1.1); + g2hist->GetYaxis()->SetTitleOffset(1.5); + g2hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g2hist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - g1overg2->GetXaxis()->SetTitle("G1 / G2"); - g1overg2->GetXaxis()->SetTitleOffset(1.1); - g1overg2->GetYaxis()->SetTitleOffset(1.5); - g1overg2->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g1overg2.png", module_str.c_str()); + g1overg2hist->GetXaxis()->SetTitle("G1 / G2"); + g1overg2hist->GetXaxis()->SetTitleOffset(1.1); + g1overg2hist->GetYaxis()->SetTitleOffset(1.5); + g1overg2hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g1overg2hist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - g1overg2er->GetXaxis()->SetTitle("G1 / G2 uncert"); - g1overg2er->GetXaxis()->SetTitleOffset(1.1); - g1overg2er->GetYaxis()->SetTitleOffset(1.5); - g1overg2er->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g1overg2er.png", module_str.c_str()); + g1overg2erhist->GetXaxis()->SetTitle("G1 / G2 uncert"); + g1overg2erhist->GetXaxis()->SetTitleOffset(1.1); + g1overg2erhist->GetYaxis()->SetTitleOffset(1.5); + g1overg2erhist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g1overg2erhist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - g1_fit_er->GetXaxis()->SetTitle("Fitted G1 uncert"); - g1_fit_er->GetXaxis()->SetTitleOffset(1.1); - g1_fit_er->GetYaxis()->SetTitleOffset(1.5); - g1_fit_er->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g1er.png", module_str.c_str()); + g1erhist->GetXaxis()->SetTitle("Fitted G1 uncert"); + g1erhist->GetXaxis()->SetTitleOffset(1.1); + g1erhist->GetYaxis()->SetTitleOffset(1.5); + g1erhist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g1erhist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - g2_fit_er->GetXaxis()->SetTitle("Fitted G2 uncert"); - g2_fit_er->GetXaxis()->SetTitleOffset(1.1); - g2_fit_er->GetYaxis()->SetTitleOffset(1.5); - g2_fit_er->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_g2er.png", module_str.c_str()); + g2erhist->GetXaxis()->SetTitle("Fitted G2 uncert"); + g2erhist->GetXaxis()->SetTitleOffset(1.1); + g2erhist->GetYaxis()->SetTitleOffset(1.5); + g2erhist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g2erhist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - diffFromLinG1->GetXaxis()->SetTitle("Normalised residual G1"); - diffFromLinG1->GetYaxis()->SetTitle(""); - diffFromLinG1->GetXaxis()->SetTitleOffset(1.1); - diffFromLinG1->GetYaxis()->SetTitleOffset(1.5); - diffFromLinG1->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_diffFromLinG1.png", module_str.c_str()); + diffFromLinG1hist->GetXaxis()->SetTitle("Normalised residual G1"); + diffFromLinG1hist->GetYaxis()->SetTitle(""); + diffFromLinG1hist->GetXaxis()->SetTitleOffset(1.1); + diffFromLinG1hist->GetYaxis()->SetTitleOffset(1.5); + diffFromLinG1hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1hist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - diffFromLinG2->GetXaxis()->SetTitle("Normalised residual G2"); - diffFromLinG2->GetYaxis()->SetTitle(""); - diffFromLinG2->GetXaxis()->SetTitleOffset(1.1); - diffFromLinG2->GetYaxis()->SetTitleOffset(1.5); - diffFromLinG2->Draw(); - sprintf(savename,"plots/M%s/CurrentSource/module_diffFromLinG2.png", module_str.c_str()); + diffFromLinG2hist->GetXaxis()->SetTitle("Normalised residual G2"); + diffFromLinG2hist->GetYaxis()->SetTitle(""); + diffFromLinG2hist->GetXaxis()->SetTitleOffset(1.1); + diffFromLinG2hist->GetYaxis()->SetTitleOffset(1.5); + diffFromLinG2hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2hist.png", module_str.c_str()); c1->SaveAs((const char *)(savename)); - cout << "G1 underflow " << diffFromLinG1->GetBinContent(0) << endl; - cout << "G2 underflow " << diffFromLinG2->GetBinContent(0) << endl; - cout << "G1 overflow " << diffFromLinG1->GetBinContent(101) << endl; - cout << "G2 overflow " << diffFromLinG2->GetBinContent(101) << endl; + cout << "G1 underflow " << diffFromLinG1hist->GetBinContent(0) << endl; + cout << "G2 underflow " << diffFromLinG2hist->GetBinContent(0) << endl; + cout << "G1 overflow " << diffFromLinG1hist->GetBinContent(101) << endl; + cout << "G2 overflow " << diffFromLinG2hist->GetBinContent(101) << endl; g1overg2_isEdge->SetLineColor(kBlue); g1overg2_isInnerEdge->SetLineColor(kCyan); @@ -744,80 +744,73 @@ int main(int argc, char* argv[]) { mapcanvas->cd(); - g1map_avg->GetXaxis()->SetTitle("Column"); - g1map_avg->GetYaxis()->SetTitle("Row"); - g1map_avg->GetYaxis()->SetTitleOffset(0.7); + g1map->GetXaxis()->SetTitle("Column"); + g1map->GetYaxis()->SetTitle("Row"); + g1map->GetYaxis()->SetTitleOffset(0.7); if (module_str == "032") { - g1map_avg->GetZaxis()->SetRangeUser(-12000,-4000); + g1map->GetZaxis()->SetRangeUser(-12000,-4000); } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { - g1map_avg->GetZaxis()->SetRangeUser(-3000,0); + g1map->GetZaxis()->SetRangeUser(-3000,0); } - g1map_avg->Draw("colz"); - sprintf(savename,"plots/M%s/CurrentSource/module_g1map_avg.png", module_str.c_str()); + g1map->Draw("colz"); + sprintf(savename,"plots/M%s/CurrentSource/g1map.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g2map_avg->GetXaxis()->SetTitle("Column"); - g2map_avg->GetYaxis()->SetTitle("Row"); - g2map_avg->GetYaxis()->SetTitleOffset(0.7); + g2map->GetXaxis()->SetTitle("Column"); + g2map->GetYaxis()->SetTitle("Row"); + g2map->GetYaxis()->SetTitleOffset(0.7); if (module_str == "032") { - g2map_avg->GetZaxis()->SetRangeUser(-800,-400); + g2map->GetZaxis()->SetRangeUser(-800,-400); } else if (module_str == "022" || module_str == "021" || module_str == "006" || module_str == "008" || module_str == "044") { - g2map_avg->GetZaxis()->SetRangeUser(-300,0); + g2map->GetZaxis()->SetRangeUser(-200,-50); } - g2map_avg->Draw("colz"); - sprintf(savename,"plots/M%s/CurrentSource/module_g2map_avg.png", module_str.c_str()); + g2map->Draw("colz"); + sprintf(savename,"plots/M%s/CurrentSource/g2map.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g1overg2_map->GetXaxis()->SetTitle("Column"); - g1overg2_map->GetYaxis()->SetTitle("Row"); - g1overg2_map->GetYaxis()->SetTitleOffset(0.7); - g1overg2_map->GetZaxis()->SetRangeUser(10,15); - g1overg2_map->Draw("colz"); - sprintf(savename,"plots/M%s/CurrentSource/module_g1overg2_map.png", module_str.c_str()); + g1overg2map->GetXaxis()->SetTitle("Column"); + g1overg2map->GetYaxis()->SetTitle("Row"); + g1overg2map->GetYaxis()->SetTitleOffset(0.7); + g1overg2map->GetZaxis()->SetRangeUser(11,15); + g1overg2map->Draw("colz"); + sprintf(savename,"plots/M%s/CurrentSource/g1overg2map.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g1overg2_map->GetXaxis()->SetRangeUser(0,256); - g1overg2_map->GetYaxis()->SetRangeUser(240,280); - g1overg2_map->Draw("colz"); - sprintf(savename,"plots/M%s/CurrentSource/module_g1overg2_map_zoom.png", module_str.c_str()); + g1overg2map->GetXaxis()->SetRangeUser(0,256); + g1overg2map->GetYaxis()->SetRangeUser(240,280); + g1overg2map->Draw("colz"); + sprintf(savename,"plots/M%s/CurrentSource/g1overg2map_zoom.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g1overg2er_map->GetXaxis()->SetTitle("Column"); - g1overg2er_map->GetYaxis()->SetTitle("Row"); - g1overg2er_map->GetYaxis()->SetTitleOffset(0.7); - g1overg2er_map->GetZaxis()->SetRangeUser(0,2); - g1overg2er_map->Draw("colz"); - sprintf(savename,"plots/M%s/CurrentSource/module_g1overg2er_map.png", module_str.c_str()); + g1overg2ermap->GetXaxis()->SetTitle("Column"); + g1overg2ermap->GetYaxis()->SetTitle("Row"); + g1overg2ermap->GetYaxis()->SetTitleOffset(0.7); + g1overg2ermap->GetZaxis()->SetRangeUser(0,0.5); + g1overg2ermap->Draw("colz"); + sprintf(savename,"plots/M%s/CurrentSource/g1overg2ermap.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g1map_fit_er->GetXaxis()->SetTitle("Column"); - g1map_fit_er->GetYaxis()->SetTitle("Row"); - g1map_fit_er->GetYaxis()->SetTitleOffset(0.7); - g1map_fit_er->Draw("colz"); - g1map_fit_er->GetZaxis()->SetRangeUser(0,50); - /* - if (module_str == "032") { - g1map_fit_er->GetZaxis()->SetRangeUser(0,100); - } else if (module_str == "006" || module_str == "008" || module_str == "022" || module_str == "021") { - g1map_fit_er->GetZaxis()->SetRangeUser(0,2000); - } - */ - sprintf(savename,"plots/M%s/CurrentSource/module_g1map_er.png", module_str.c_str()); + g1ermap->GetXaxis()->SetTitle("Column"); + g1ermap->GetYaxis()->SetTitle("Row"); + g1ermap->GetYaxis()->SetTitleOffset(0.7); + g1ermap->Draw("colz"); + g1ermap->GetZaxis()->SetRangeUser(0,50); + sprintf(savename,"plots/M%s/CurrentSource/g1ermap.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - g2map_fit_er->GetXaxis()->SetTitle("Column"); - g2map_fit_er->GetYaxis()->SetTitle("Row"); - g2map_fit_er->GetYaxis()->SetTitleOffset(0.7); - g2map_fit_er->Draw("colz"); - g2map_fit_er->GetZaxis()->SetRangeUser(0,10); - sprintf(savename,"plots/M%s/CurrentSource/module_g2map_er.png", module_str.c_str()); + g2ermap->GetXaxis()->SetTitle("Column"); + g2ermap->GetYaxis()->SetTitle("Row"); + g2ermap->GetYaxis()->SetTitleOffset(0.7); + g2ermap->Draw("colz"); + g2ermap->GetZaxis()->SetRangeUser(0,5); + sprintf(savename,"plots/M%s/CurrentSource/g2ermap.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); diffFromLinG1map->GetXaxis()->SetTitle("Column"); diffFromLinG1map->GetYaxis()->SetTitle("Row"); diffFromLinG1map->GetYaxis()->SetTitleOffset(0.7); diffFromLinG1map->Draw("colz"); - diffFromLinG1map->GetZaxis()->SetRangeUser(-0.01,0.01); + diffFromLinG1map->GetZaxis()->SetRangeUser(-0.04,0.04); sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1map.png", module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); @@ -831,8 +824,8 @@ int main(int argc, char* argv[]) { 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"); - g1overg2_map->Write(); - g1overg2er_map->Write(); + g1overg2map->Write(); + g1overg2ermap->Write(); saved_file2->Close(); }