From 0a83c6e01f127047bf85438de900aef5b2e306d5 Mon Sep 17 00:00:00 2001 From: redford_s Date: Mon, 13 Nov 2017 10:47:21 +0100 Subject: [PATCH] Adding the file to compare BP and DB data for a given module. --- .gitignore | 1 + BP_DB_Comp.cpp | 179 ++++++++++++++++++++++++++++++++++ run_JFModuleCalib_CuFBDCS.txt | 6 ++ 3 files changed, 186 insertions(+) create mode 100644 BP_DB_Comp.cpp diff --git a/.gitignore b/.gitignore index 33825e0..c5ad14d 100644 --- a/.gitignore +++ b/.gitignore @@ -5,6 +5,7 @@ JFMC_CuFluoPeakFit JFMC_CurrentSourceScan JFMC_DirectBeamScan dummyRootFileMaker +BP_DB_Comp data/ plots/ diff --git a/BP_DB_Comp.cpp b/BP_DB_Comp.cpp new file mode 100644 index 0000000..eef7a97 --- /dev/null +++ b/BP_DB_Comp.cpp @@ -0,0 +1,179 @@ +// to compare the results of backplane pulsing with direct beam + +#include "/afs/psi.ch/project/mythen/sophie/sls_detector_calibration/jungfrauCommonHeader.h" +#include "/afs/psi.ch/project/mythen/sophie/sls_detector_calibration/jungfrauCommonFunctions.h" + +#include "TFile.h" +#include "TH2F.h" +#include "TCanvas.h" + +int main(int argc, char* argv[]) { + jungfrauStyle(); + gROOT->SetBatch(1); + gStyle->SetOptFit(11); + + if (argc != 2) { + cout << "Correct usage:" << endl; + cout << "arg 1: specify module number" << endl; + exit(1); + } + + string module_str = argv[1]; + char savename[128]; + + TCanvas *mapcanvas = new TCanvas("mapcanvas","",150,10,800,400); + mapcanvas->SetLeftMargin(0.1); + mapcanvas->SetRightMargin(0.13); + mapcanvas->SetTopMargin(0.08); + mapcanvas->SetBottomMargin(0.15); + + TCanvas *c1 = new TCanvas("c1",""); + + // BP data + sprintf(savename,"../JFModuleCalib_CuFBPCS/data/M%s/BP_ratio_M%s.root", module_str.c_str(), module_str.c_str()); + TFile* saved_file2 = new TFile((const char *)(savename),"READ"); + TH2F* g0overg1_map_BP = (TH2F*)saved_file2->Get("g0overg1_map"); + + // DB data + sprintf(savename,"data/M%s/DB_ratio_M%s.root", module_str.c_str(), module_str.c_str()); + TFile* saved_file = new TFile((const char *)(savename),"READ"); + TH2F* g0overg1_map_DB = (TH2F*)saved_file->Get("g0overg1map"); + TH2F* g0overg1er_map_DB = (TH2F*)saved_file->Get("g0overg1ermap"); + + // difference and pull + TH1F* diff_g0overg1_hist = new TH1F("diff_g0overg1_hist","",100,-20,20); + TH2F* diff_g0overg1_map = new TH2F("diff_g0overg1_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH1F* pull_hist = new TH1F("pull_hist","",100,-20,20); + + // scatter plots + TH2F* scat = new TH2F("scat","",100,-50,0,100,-40,-10); + TH2F* scat_e = new TH2F("scat_e","",100,-50,0,100,-40,-10); + TH2F* scat_d = new TH2F("scat_d","",100,-50,0,100,-40,-10); + TH2F* scat_ntd = new TH2F("scat_ntd","",100,-50,0,100,-40,-10); + TH2F* scat_q = new TH2F("scat_q","",100,-50,0,100,-40,-10); + TH2F* scat_ie = new TH2F("scat_ie","",100,-50,0,100,-40,-10); + TH2F* scat_b = new TH2F("scat_b","",100,-50,0,100,-40,-10); + + double this_BP = 0; + double this_DB = 0; + + for (int i = 0; i < NCH; i ++) { + + this_BP = g0overg1_map_BP->GetBinContent((i%NC)+1,(i/NC)+1); + this_DB = g0overg1_map_DB->GetBinContent((i%NC)+1,(i/NC)+1); + + if (this_BP != 0 && this_DB != 0) { + diff_g0overg1_hist->Fill(this_DB - this_BP); + diff_g0overg1_map->SetBinContent((i%NC)+1,(i/NC)+1,this_DB - this_BP); + pull_hist->Fill((this_DB - this_BP) / g0overg1er_map_DB->GetBinContent((i%NC)+1,(i/NC)+1)); + + scat->Fill(this_DB, this_BP, 1); + if (isEdge(i)) { + scat_e->Fill(this_DB, this_BP, 1); + } + if (isDouble(i)) { + scat_d->Fill(this_DB, this_BP, 1); + } + if (isNextToDouble(i)) { + scat_ntd->Fill(this_DB, this_BP, 1); + } + if (isQuad(i)) { + scat_q->Fill(this_DB, this_BP, 1); + } + if (isInnerEdge(i)) { + scat_ie->Fill(this_DB, this_BP, 1); + } + if (isBulk(i)) { + scat_b->Fill(this_DB, this_BP, 1); + } + + } + } + + c1->cd(); + diff_g0overg1_hist->GetXaxis()->SetTitle("(G0 / G1)_{DB} - (G0 / G1)_{BP}"); + diff_g0overg1_hist->Draw(); + sprintf(savename,"plots/M%s/BP_DB_Comp/diff_g0overg1_hist.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + c1->SetLogy(); + diff_g0overg1_hist->Draw(); + sprintf(savename,"plots/M%s/BP_DB_Comp/diff_g0overg1_hist_log.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + c1->SetLogy(0); + + mapcanvas->cd(); + diff_g0overg1_map->GetXaxis()->SetTitle("Column"); + diff_g0overg1_map->GetYaxis()->SetTitle("Row"); + diff_g0overg1_map->GetYaxis()->SetTitleOffset(0.7); + diff_g0overg1_map->Draw("colz"); + diff_g0overg1_map->GetZaxis()->SetRangeUser(-20,20); + sprintf(savename,"plots/M%s/BP_DB_Comp/diff_g0overg1_map.png", module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + + c1->cd(); + pull_hist->GetXaxis()->SetTitle("(G0 / G1)_{DB} - (G0 / G1)_{BP} / #sigma_{(G0 / G1)_{DB}}"); + pull_hist->Fit("gaus"); + pull_hist->Draw(); + sprintf(savename,"plots/M%s/BP_DB_Comp/pull_hist.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + c1->SetRightMargin(0.2); + + // all + scat->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat->GetYaxis()->SetTitleOffset(1.3); + scat->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // edge pixels + scat_e->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_e->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_e->GetYaxis()->SetTitleOffset(1.3); + scat_e->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_e.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // double pixels + scat_d->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_d->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_d->GetYaxis()->SetTitleOffset(1.3); + scat_d->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_d.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // next to double pixels + scat_ntd->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_ntd->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_ntd->GetYaxis()->SetTitleOffset(1.3); + scat_ntd->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_ntd.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // quad pixels + scat_q->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_q->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_q->GetYaxis()->SetTitleOffset(1.3); + scat_q->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_q.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // inner edge pixels + scat_ie->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_ie->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_ie->GetYaxis()->SetTitleOffset(1.3); + scat_ie->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_ie.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + // bulk pixels + scat_b->GetXaxis()->SetTitle("(G0 / G1)_{DB}"); + scat_b->GetYaxis()->SetTitle("(G0 / G1)_{BP}"); + scat_b->GetYaxis()->SetTitleOffset(1.3); + scat_b->Draw("colz"); + sprintf(savename,"plots/M%s/BP_DB_Comp/scat_b.png", module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + +} diff --git a/run_JFModuleCalib_CuFBDCS.txt b/run_JFModuleCalib_CuFBDCS.txt index 201998b..cc7424b 100644 --- a/run_JFModuleCalib_CuFBDCS.txt +++ b/run_JFModuleCalib_CuFBDCS.txt @@ -5,6 +5,7 @@ Module 006 - ESRF ./JFMC_CuFluoPeak 006 /mnt/pcmoench_jungfrau_data/Module_006_161116 10us_500Hz_QS_allgain_pede 10us_500Hz_QS_G0_Cu ./JFMC_CuFluoPeakFit 006 ./JFMC_DirectBeamScan 006 /mnt/pcmoench_jungfrau_data/Module_006_161116 100us_500Hz_QS_allgain_pede_1055 100us_500Hz_QS_beam_scan_a +./BP_DB_Comp 006 ./JFMC_CurrentSourceScan 006 /mnt/pcmoench_jungfrau_data1_ib/Module_006_2017-08-21_CSBaselineTest CurrentSource_scan12 ./JFMC_CalibWriter 006 @@ -21,6 +22,7 @@ Module 021 - 1M top LCLS, PX ./JFMC_CuFluoPeak 021 /mnt/pcmoench_jungfrau_data/1M_021_022_150916 10us_200Hz_QS_allgain_pede__B 10us_200Hz_QS_G0_Cu__B ./JFMC_CuFluoPeakFit 021 ./JFMC_DirectBeamScan 021 /mnt/pcmoench_jungfrau_data/1M_021_022_150916 200us_200Hz_QS_allgain_pede_scan_d_B 200us_200Hz_QS_beam_scan_d_B +./BP_DB_Comp 021 ./JFMC_CurrentSourceScan 021 /mnt/pcmoench_jungfrau_data1_ib/1M_021_022_150916 Current_source_scan_a_B ./JFMC_CalibWriter 021 @@ -29,6 +31,7 @@ Module 022 - 1M bottom LCLS, PX ./JFMC_CuFluoPeak 022 /mnt/pcmoench_jungfrau_data/1M_021_022_150916 10us_200Hz_QS_allgain_pede__A 10us_200Hz_QS_G0_Cu__A ./JFMC_CuFluoPeakFit 022 ./JFMC_DirectBeamScan 022 /mnt/pcmoench_jungfrau_data/1M_021_022_150916 200us_200Hz_QS_allgain_pede_scan_d_A 200us_200Hz_QS_beam_scan_d_A +./BP_DB_Comp 022 ./JFMC_CurrentSourceScan 022 /mnt/pcmoench_jungfrau_data1_ib/Module_022_Calib CS_HS_NF_200Hz_scan1 ./JFMC_CalibWriter 022 @@ -49,6 +52,7 @@ Module 044 - added to the 1M to make the 1.5M, cSAXS ./JFMC_CuFluoPeak 044 /mnt/pcmoench_jungfrau_data/Module044 allpede_10us_500Hz_QS_G0_1500 CuFluo_10us_500Hz_QS_G0 ./JFMC_CuFluoPeakFit 044 ./JFMC_DirectBeamScan 044 /mnt/pcmoench_jungfrau_data/Module044 allpede_150us_500Hz_QS_1620 DB_170221_scan1 +./BP_DB_Comp 044 ./JFMC_CurrentSourceScan 044 /mnt/pcmoench_jungfrau_data1_ib/Module044 CS_170221_scan2 ./JFMC_CalibWriter 044 @@ -81,6 +85,7 @@ Module 066 - replacing 027 at ESRF ./JFMC_CuFluoPeak 066 /mnt/pcmoench_jungfrau_data1_ib/Module_066_2017-10-24_Calib pedeG0_1021 CuFluo_10us_500Hz ./JFMC_CuFluoPeakFit 066 ./JFMC_DirectBeamScan 066 /mnt/pcmoench_jungfrau_data1_ib/Module_066_2017-10-24_Calib DB_OriginalScan_200Hz_200us_pede DB_DoubleGranularity_500Hz_200us +./BP_DB_Comp 066 ./JFMC_CurrentSourceScan 066 /mnt/pcmoench_jungfrau_data1_ib/Module_066_2017-10-24_Calib CS_scan1 @@ -112,6 +117,7 @@ Module 094 ./JFMC_CuFluoPeak 094 /mnt/pcmoench_jungfrau_data1_ib/Module_094_2017-10-24_Calib pedeG0_1041 CuFluo_10us_500Hz ./JFMC_CuFluoPeakFit 094 ./JFMC_DirectBeamScan 094 /mnt/pcmoench_jungfrau_data1_ib/Module_094_2017-10-24_Calib DB_DoubleGranularity_500Hz_200us_pede DB_DoubleGranularity_500Hz_200us +./BP_DB_Comp 094 ./JFMC_CurrentSourceScan 094 /mnt/pcmoench_jungfrau_data1_ib/Module_094_2017-10-24_Calib CS_scan2