diff --git a/JFMC_CurrentSourceScan.cpp b/JFMC_CurrentSourceScan.cpp index 274f4df..0d398b3 100644 --- a/JFMC_CurrentSourceScan.cpp +++ b/JFMC_CurrentSourceScan.cpp @@ -29,7 +29,6 @@ int main(int argc, char* argv[]) { cout << "arg 1: specify module number" << endl; cout << "arg 2: specify data location" << endl; cout << "arg 3: specify data file" << endl; - cout << "eg: ./JFMC_CurrentSourceScan 006 ..." << endl; cout << " " << endl; exit(1); } @@ -575,18 +574,24 @@ int main(int argc, char* argv[]) { } } + TH2F *g0map = new TH2F("g0map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g0ermap = new TH2F("g0ermap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + 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 *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); + TH2F *g0overg1map = new TH2F("g0overg1map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); + TH2F *g0overg1ermap = new TH2F("g0overg1ermap","",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); - 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 *g0hist = new TH1F("g0hist","",100,10000,50000); + TH1F *g0erhist = new TH1F("g0erhist","",100,0,200); + TH1F *g1hist = 0; if (module_str == "032") { g1hist = new TH1F("g1hist","",100,-12000,-4000); @@ -603,11 +608,18 @@ int main(int argc, char* argv[]) { } TH1F *g2erhist = new TH1F("g2erhist","",100,0,10); + TH1F *g0overg1hist = new TH1F("g0overg1hist","",100,-50,0); + TH1F *g0overg1erhist = new TH1F("g0overg1erhist","",100,0,1); + 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 *g0overg1_isEdge = new TH1F("g0overg1_isEdge","",100,-50,0); + TH1F *g0overg1_isInnerEdge = new TH1F("g0overg1_isInnerEdge","",100,-50,0); + TH1F *g0overg1_isDouble = new TH1F("g0overg1_isDouble","",100,-50,0); + TH1F *g0overg1_isNextToDouble = new TH1F("g0overg1_isNextToDouble","",100,-50,0); + TH1F *g0overg1_isQuad = new TH1F("g0overg1_isQuad","",100,-50,0); + TH1F *g0overg1_isBulk = new TH1F("g0overg1_isBulk","",100,-50,0); TH1F *g1overg2_isEdge = new TH1F("g1overg2_isEdge","",100,11,15); TH1F *g1overg2_isInnerEdge = new TH1F("g1overg2_isInnerEdge","",100,11,15); @@ -689,7 +701,12 @@ int main(int argc, char* argv[]) { } } + mapcanvas->cd(); + mapcanvas->SetLeftMargin(0.13); + mapcanvas->SetRightMargin(0.05); + for (int i = 0; i < NCH; i++) { + if (pixel_mask[i] == 1) { vector r0_adc; @@ -702,10 +719,10 @@ int main(int argc, char* argv[]) { vector r1_adcerr; vector r1_ferr; - vector r3_adc; - vector r3_filter; - vector r3_adcerr; - vector r3_ferr; + vector r2_adc; + vector r2_filter; + vector r2_adcerr; + vector r2_ferr; // filter out points at zero and points with 0 error (mean taken from one entry in hist) int m_max = 0; @@ -734,104 +751,46 @@ int main(int argc, char* argv[]) { r1_ferr.push_back(0.); } if (adcmap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1) != 0 && adcermap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1) > 0) { - r3_adc.push_back(adcmap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1)); - r3_filter.push_back(filter[m]); - r3_adcerr.push_back(adcermap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1)); - r3_ferr.push_back(0.); + r2_adc.push_back(adcmap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1)); + r2_filter.push_back(filter[m]); + r2_adcerr.push_back(adcermap_avg_g2[m]->GetBinContent((i%NC)+1,(i/NC)+1)); + r2_ferr.push_back(0.); } } - if (r1_adc.size() > 1 && r3_adc.size() > 1) { + TGraphErrors *grap_g0 = 0; + TGraphErrors *grap_g1 = 0; + TGraphErrors *grap_g2 = 0; - TGraphErrors *grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); - TGraphErrors *grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); - TGraphErrors *grap_g2 = new TGraphErrors(r3_adc.size(),&(r3_filter[0]),&(r3_adc[0]),&(r3_ferr[0]),&(r3_adcerr[0])); + TF1 *fit_g0 = 0; + TF1 *fit_g1 = 0; + TF1 *fit_g2 = 0; - double rangemin2 = 0; - if (module_str == "032") { - rangemin2 = 0.75; - } else if (module_str == "008" || module_str == "040" || module_str == "044") { - rangemin2 = *min_element(r1_filter.begin(),r1_filter.end()); - } else if (module_str == "006" || module_str == "021" || module_str == "022" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "065" || module_str == "069" || module_str == "079" || module_str == "068" || module_str == "088" || module_str == "055" || module_str == "059" || module_str == "094" || module_str == "066" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "124" || module_str == "126") { - rangemin2 = 2.; - } - double rangemax2 = *max_element(r1_filter.begin(),r1_filter.end()); - TF1 *fit2 = new TF1("fit2","[0]+[1]*x",rangemin2,rangemax2); - fit2->SetParameter(0, 10000.); - fit2->SetParameter(1, -0.1); - fit2->SetLineColor(kGreen+2); - fit2->SetParName(0,"G1 const"); - fit2->SetParName(1,"G1 grad"); - grap_g1->Fit(fit2,"QR",""); + double rangemin_g0; + double rangemax_g0; + double rangemin_g1; + double rangemax_g1; + double rangemin_g2; + double rangemax_g2; - 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)); + if (r0_adc.size() > 1) { - double g1residsum = 0.; - double g1residn = 0.; - for (size_t it = 0; it < r1_adc.size(); it++) { - if (r1_filter[it] >= rangemin2) { - 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.; - } - } - diffFromLinG1map->Fill(i%NC,i/NC,g1residsum/g1residn); - - double rangemin3 = 0; - if (module_str == "032") { - rangemin3 = 5; - } else { - rangemin3 = *min_element(r3_filter.begin(),r3_filter.end()); - } - double rangemax3 = *max_element(r3_filter.begin(),r3_filter.end()); - TF1 *fit3 = new TF1("fit3","[0]+[1]*x",rangemin3, rangemax3); - fit3->SetParameter(0, 10000.); - fit3->SetParameter(1, -0.01); - fit3->SetLineColor(kRed); - fit3->SetParName(0,"G2 const"); - fit3->SetParName(1,"G2 grad"); - grap_g2->Fit(fit3,"QR+",""); + grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); - 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)); + rangemin_g0 = 0.15; + rangemax_g0 = *max_element(r0_filter.begin(),r0_filter.end()); + fit_g0 = new TF1("fit_g0","[0]+[1]*x",rangemin_g0, rangemax_g0); + fit_g0->SetParameter(0, 3000.); + fit_g0->SetParameter(1, 30000.); + fit_g0->SetLineColor(kBlue); + fit_g0->SetParName(0,"G0 const"); + fit_g0->SetParName(1,"G0 grad"); + grap_g0->Fit(fit_g0,"QR+",""); - double g2residsum = 0.; - double g2residn = 0.; - for (size_t it = 0; it < r3_adc.size(); 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.; - } - diffFromLinG2map->Fill(i%NC,i/NC,g2residsum/g2residn); - - 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)); - } - if (isInnerEdge(i)) { - g1overg2_isInnerEdge->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - } - if (isDouble(i)) { - g1overg2_isDouble->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - } - if (isNextToDouble(i)) { - g1overg2_isNextToDouble->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - } - if (isQuad(i)) { - g1overg2_isQuad->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - } - if (isBulk(i)) { - g1overg2_isBulk->Fill(fit2->GetParameter(1) / fit3->GetParameter(1)); - } + g0hist->Fill(fit_g0->GetParameter(1)); + g0map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)); + g0erhist->Fill(fit_g0->GetParError(1)); + g0ermap->Fill(i%NC,i/NC,fit_g0->GetParError(1)); if ((i >= 58000 && i < 58000+10) || // bulk (i >= 10 && i < 10+10) || // edge @@ -856,127 +815,406 @@ int main(int argc, char* argv[]) { pixel_type = "q"; } - mapcanvas->cd(); - mapcanvas->SetLeftMargin(0.13); - mapcanvas->SetRightMargin(0.05); - grap_g0->SetMarkerStyle(20); - grap_g1->SetMarkerStyle(20); - grap_g2->SetMarkerStyle(20); grap_g0->SetMarkerColor(kBlue); - grap_g1->SetMarkerColor(kGreen+2); - grap_g2->SetMarkerColor(kRed); grap_g0->SetLineColor(kBlue); - grap_g1->SetLineColor(kGreen+2); - grap_g2->SetLineColor(kRed); - - TF1 *fit2_e = new TF1("fit2_e","[0]+[1]*x",0.009,1.1E3); - TF1 *fit3_e = new TF1("fit3_e","[0]+[1]*x",0.009,1.1E3); - - fit2_e->SetParameters(fit2->GetParameter(0),fit2->GetParameter(1)); - fit2_e->SetLineColor(kGreen+2); - fit2_e->SetLineStyle(2); - fit3_e->SetParameters(fit3->GetParameter(0),fit3->GetParameter(1)); - fit3_e->SetLineColor(kRed); - fit3_e->SetLineStyle(2); - if (module_str == "021" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "124" || module_str == "126") { - grap_g1->GetXaxis()->SetLimits(0.009,5E2); - } else { - grap_g1->GetXaxis()->SetLimits(0.009,1.1E2); - } - grap_g1->SetMinimum(0); - grap_g1->SetMaximum(20000); - grap_g1->GetXaxis()->SetTitle("Integration time [#mus]"); - grap_g1->GetYaxis()->SetTitle("ADC [ADU]"); - grap_g1->GetYaxis()->SetTitleOffset(0.9); - grap_g1->Draw("AP"); - fit2->Draw("same"); + TF1 *extended_g0 = new TF1("extended_g0","[0]+[1]*x",0.009,1.1E3); + extended_g0->SetParameters(fit_g0->GetParameter(0),fit_g0->GetParameter(1)); + extended_g0->SetLineColor(kBlue); + extended_g0->SetLineStyle(2); + + grap_g0->SetMinimum(0); + grap_g0->SetMaximum(16000); + grap_g0->GetXaxis()->SetTitle("Integration time [#mus]"); + grap_g0->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g0->GetYaxis()->SetTitleOffset(0.9); + grap_g0->Draw("AP"); + fit_g0->Draw("same"); mapcanvas->Update(); - TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); + TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); st0->SetX1NDC(0.2); st0->SetX2NDC(0.54); st0->SetY1NDC(0.71); st0->SetY2NDC(0.90); st0->SetBorderSize(0); - st0->SetTextColor(kGreen+2); - grap_g2->Draw("P"); - fit3->Draw("same"); - mapcanvas->Update(); - TPaveStats *st1 = (TPaveStats*)grap_g2->FindObject("stats"); - st1->SetX1NDC(0.6); - st1->SetX2NDC(0.94); - st1->SetY1NDC(0.71); - st1->SetY2NDC(0.90); - st1->SetBorderSize(0); - st1->SetTextColor(kRed); - - grap_g0->Draw("P"); - fit2_e->Draw("same"); - fit3_e->Draw("same"); + st0->SetTextColor(kBlue); + extended_g0->Draw("same"); mapcanvas->SetLogx(); - sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g0_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + mapcanvas->SetLogx(0); + + vector r0_adc_norm; + for (size_t j = 0; j < r0_adc.size(); j++) { + r0_adc_norm.push_back(r0_adc[j] - fit_g0->Eval(r0_filter[j])); + } + + TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcerr[0])); + norm_g0->SetMarkerColor(kBlue); + norm_g0->SetLineColor(kBlue); + + TF1* flat_g0 = new TF1("flat_g0","0",rangemin_g0, rangemax_g0); + flat_g0->SetLineColor(kBlue); + + norm_g0->GetXaxis()->SetTitle("Integration time [#mus]"); + norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g0->GetYaxis()->SetTitleOffset(0.9); + norm_g0->SetMinimum(-100); + norm_g0->SetMaximum(100); + norm_g0->Draw("AP"); + flat_g0->Draw("same"); + sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g0norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + + delete norm_g0; + delete extended_g0; + delete flat_g0; + } + } + + if (r1_adc.size() > 1) { + + grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); + + rangemin_g1 = 0; + if (module_str == "032") { + rangemin_g1 = 0.75; + } else if (module_str == "008" || module_str == "040" || module_str == "044") { + rangemin_g1 = *min_element(r1_filter.begin(),r1_filter.end()); + } else if (module_str == "006" || module_str == "021" || module_str == "022" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "065" || module_str == "069" || module_str == "079" || module_str == "068" || module_str == "088" || module_str == "055" || module_str == "059" || module_str == "094" || module_str == "066" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "124" || module_str == "126") { + rangemin_g1 = 2.; + } + rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end()); + fit_g1 = new TF1("fit_g1","[0]+[1]*x",rangemin_g1,rangemax_g1); + fit_g1->SetParameter(0, 10000.); + fit_g1->SetParameter(1, -0.1); + fit_g1->SetLineColor(kGreen+2); + fit_g1->SetParName(0,"G1 const"); + fit_g1->SetParName(1,"G1 grad"); + grap_g1->Fit(fit_g1,"QR",""); + + g1hist->Fill(fit_g1->GetParameter(1)); + g1map->Fill(i%NC,i/NC,fit_g1->GetParameter(1)); + g1erhist->Fill(fit_g1->GetParError(1)); + g1ermap->Fill(i%NC,i/NC,fit_g1->GetParError(1)); + + if ((i >= 58000 && i < 58000+10) || // bulk + (i >= 10 && i < 10+10) || // edge + (i >= 1024+10 && i < 1024+10+10) || // inner edge + (i >= (256*1024)+10 && i < (256*1024)+10+10) || // double + (i >= (257*1024)+10 && i < (257*1024)+10+10) || // next to double + (i == (255*1024)+255) // quad + ) { + + string pixel_type = "x"; + if (i >= 58000 && i < 58000+10) { + pixel_type = "b"; + } else if (i >= 10 && i < 10+10) { + pixel_type = "e"; + } else if (i >= 1024+10 && i < 1024+10+10) { + pixel_type = "ie"; + } else if (i >= (256*1024)+10 && i < (256*1024)+10+10) { + pixel_type = "d"; + } else if (i >= (257*1024)+10 && i < (257*1024)+10+10) { + pixel_type = "ntd"; + } else if (i == (255*1024)+255) { + pixel_type = "q"; + } + + grap_g1->SetMarkerStyle(20); + grap_g1->SetMarkerColor(kGreen+2); + grap_g1->SetLineColor(kGreen+2); + + TF1 *extended_g1 = new TF1("extended_g1","[0]+[1]*x",0.009,1.1E3); + extended_g1->SetParameters(fit_g1->GetParameter(0),fit_g1->GetParameter(1)); + extended_g1->SetLineColor(kGreen+2); + extended_g1->SetLineStyle(2); + + grap_g1->SetMinimum(0); + grap_g1->SetMaximum(16000); + grap_g1->GetXaxis()->SetTitle("Integration time [#mus]"); + grap_g1->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g1->GetYaxis()->SetTitleOffset(0.9); + grap_g1->Draw("AP"); + fit_g1->Draw("same"); + mapcanvas->Update(); + TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); + st0->SetX1NDC(0.2); + st0->SetX2NDC(0.54); + st0->SetY1NDC(0.17); + st0->SetY2NDC(0.36); + st0->SetBorderSize(0); + st0->SetTextColor(kGreen+2); + extended_g1->Draw("same"); + mapcanvas->SetLogx(); + sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g1_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SetLogx(0); - // normalise vector r1_adc_norm; - vector r3_adc_norm; for (size_t j = 0; j < r1_adc.size(); j++) { - r1_adc_norm.push_back(r1_adc[j] - fit2->Eval(r1_filter[j])); + r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j])); } - for (size_t j = 0; j < r3_adc.size(); j++) { - r3_adc_norm.push_back(r3_adc[j] - fit3->Eval(r3_filter[j])); - } - TGraphErrors *grap1_norm = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0])); - TGraphErrors *grap2_norm = new TGraphErrors(r3_adc.size(),&(r3_filter[0]),&(r3_adc_norm[0]),&(r3_ferr[0]),&(r3_adcerr[0])); + + TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0])); + norm_g1->SetMarkerColor(kGreen+2); + norm_g1->SetLineColor(kGreen+2); - TF1* flat_g1 = new TF1("flat_g1","0",rangemin2, rangemax2); - TF1* flat_g2 = new TF1("flat_g2","0",rangemin3, rangemax3); - - grap1_norm->SetMarkerColor(kGreen+2); - grap1_norm->SetLineColor(kGreen+2); + TF1* flat_g1 = new TF1("flat_g1","0",rangemin_g1, rangemax_g1); flat_g1->SetLineColor(kGreen+2); - grap1_norm->GetXaxis()->SetTitle("Integration time [#mus]"); - grap1_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap1_norm->GetYaxis()->SetTitleOffset(0.9); - grap1_norm->SetMinimum(-100); - grap1_norm->SetMaximum(100); - grap1_norm->Draw("AP"); + + norm_g1->GetXaxis()->SetTitle("Integration time [#mus]"); + norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g1->GetYaxis()->SetTitleOffset(0.9); + norm_g1->SetMinimum(-100); + norm_g1->SetMaximum(100); + norm_g1->Draw("AP"); flat_g1->Draw("same"); sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g1norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - grap2_norm->SetMarkerColor(kRed); - grap2_norm->SetLineColor(kRed); + delete norm_g1; + delete extended_g1; + delete flat_g1; + } + } + + if (r2_adc.size() > 1) { + + grap_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc[0]),&(r2_ferr[0]),&(r2_adcerr[0])); + + rangemin_g2 = 0; + if (module_str == "032") { + rangemin_g2 = 5; + } else { + rangemin_g2 = *min_element(r2_filter.begin(),r2_filter.end()); + } + rangemax_g2 = *max_element(r2_filter.begin(),r2_filter.end()); + fit_g2 = new TF1("fit_g2","[0]+[1]*x",rangemin_g2, rangemax_g2); + fit_g2->SetParameter(0, 10000.); + fit_g2->SetParameter(1, -0.01); + fit_g2->SetLineColor(kRed); + fit_g2->SetParName(0,"G2 const"); + fit_g2->SetParName(1,"G2 grad"); + grap_g2->Fit(fit_g2,"QR+",""); + + g2hist->Fill(fit_g2->GetParameter(1)); + g2map->Fill(i%NC,i/NC,fit_g2->GetParameter(1)); + g2erhist->Fill(fit_g2->GetParError(1)); + g2ermap->Fill(i%NC,i/NC,fit_g2->GetParError(1)); + + if ((i >= 58000 && i < 58000+10) || // bulk + (i >= 10 && i < 10+10) || // edge + (i >= 1024+10 && i < 1024+10+10) || // inner edge + (i >= (256*1024)+10 && i < (256*1024)+10+10) || // double + (i >= (257*1024)+10 && i < (257*1024)+10+10) || // next to double + (i == (255*1024)+255) // quad + ) { + + string pixel_type = "x"; + if (i >= 58000 && i < 58000+10) { + pixel_type = "b"; + } else if (i >= 10 && i < 10+10) { + pixel_type = "e"; + } else if (i >= 1024+10 && i < 1024+10+10) { + pixel_type = "ie"; + } else if (i >= (256*1024)+10 && i < (256*1024)+10+10) { + pixel_type = "d"; + } else if (i >= (257*1024)+10 && i < (257*1024)+10+10) { + pixel_type = "ntd"; + } else if (i == (255*1024)+255) { + pixel_type = "q"; + } + + grap_g2->SetMarkerStyle(20); + grap_g2->SetMarkerColor(kRed); + grap_g2->SetLineColor(kRed); + + TF1 *extended_g2 = new TF1("extended_g2","[0]+[1]*x",0.009,1.1E3); + extended_g2->SetParameters(fit_g2->GetParameter(0),fit_g2->GetParameter(1)); + extended_g2->SetLineColor(kRed); + extended_g2->SetLineStyle(2); + + grap_g2->SetMinimum(0); + grap_g2->SetMaximum(16000); + grap_g2->GetXaxis()->SetTitle("Integration time [#mus]"); + grap_g2->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g2->GetYaxis()->SetTitleOffset(0.9); + grap_g2->Draw("AP"); + fit_g2->Draw("same"); + mapcanvas->Update(); + TPaveStats *st0 = (TPaveStats*)grap_g2->FindObject("stats"); + st0->SetX1NDC(0.2); + st0->SetX2NDC(0.54); + st0->SetY1NDC(0.17); + st0->SetY2NDC(0.36); + st0->SetBorderSize(0); + st0->SetTextColor(kRed); + extended_g2->Draw("same"); + mapcanvas->SetLogx(); + sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g2_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + mapcanvas->SetLogx(0); + + vector r2_adc_norm; + for (size_t j = 0; j < r2_adc.size(); j++) { + r2_adc_norm.push_back(r2_adc[j] - fit_g2->Eval(r2_filter[j])); + } + + TGraphErrors *norm_g2 = new TGraphErrors(r2_adc.size(),&(r2_filter[0]),&(r2_adc_norm[0]),&(r2_ferr[0]),&(r2_adcerr[0])); + norm_g2->SetMarkerColor(kRed); + norm_g2->SetLineColor(kRed); + + TF1* flat_g2 = new TF1("flat_g2","0",rangemin_g2, rangemax_g2); flat_g2->SetLineColor(kRed); - grap2_norm->GetXaxis()->SetTitle("Integration time [#mus]"); - grap2_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap2_norm->GetYaxis()->SetTitleOffset(0.9); - grap2_norm->SetMinimum(-100); - grap2_norm->SetMaximum(100); - grap2_norm->Draw("AP"); + + norm_g2->GetXaxis()->SetTitle("Integration time [#mus]"); + norm_g2->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g2->GetYaxis()->SetTitleOffset(0.9); + norm_g2->SetMinimum(-100); + norm_g2->SetMaximum(100); + norm_g2->Draw("AP"); flat_g2->Draw("same"); sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_g2norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - mapcanvas->SetLeftMargin(0.1); - mapcanvas->SetRightMargin(0.13); - delete fit2_e; - delete fit3_e; - delete grap1_norm; - delete grap2_norm; + delete norm_g2; + delete extended_g2; + delete flat_g2; } - - delete grap_g1; - delete grap_g2; - delete fit2; - delete fit3; - } + + if (r0_adc.size() > 1 && r1_adc.size() > 1) { + + g0overg1hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + g0overg1erhist->Fill(abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2))); + g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + g0overg1ermap->Fill(i%NC,i/NC,abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2))); + + if (isEdge(i)) { + g0overg1_isEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + if (isInnerEdge(i)) { + g0overg1_isInnerEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + if (isDouble(i)) { + g0overg1_isDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + if (isNextToDouble(i)) { + g0overg1_isNextToDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + if (isQuad(i)) { + g0overg1_isQuad->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + if (isBulk(i)) { + g0overg1_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + } + } + + + if (r1_adc.size() > 1 && r2_adc.size() > 1) { + + g1overg2hist->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + g1overg2erhist->Fill(abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2))); + g1overg2map->Fill(i%NC,i/NC,fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + g1overg2ermap->Fill(i%NC,i/NC,abs(fit_g1->GetParameter(1)/fit_g2->GetParameter(1))*sqrt(pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2) + pow((fit_g2->GetParError(1)/fit_g2->GetParameter(1)),2))); + + if (isEdge(i)) { + g1overg2_isEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + if (isInnerEdge(i)) { + g1overg2_isInnerEdge->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + if (isDouble(i)) { + g1overg2_isDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + if (isNextToDouble(i)) { + g1overg2_isNextToDouble->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + if (isQuad(i)) { + g1overg2_isQuad->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + if (isBulk(i)) { + g1overg2_isBulk->Fill(fit_g1->GetParameter(1) / fit_g2->GetParameter(1)); + } + } + + if (r0_adc.size() > 1 && r1_adc.size() > 1 && r2_adc.size() > 1) { + if ((i >= 58000 && i < 58000+10) || // bulk + (i >= 10 && i < 10+10) || // edge + (i >= 1024+10 && i < 1024+10+10) || // inner edge + (i >= (256*1024)+10 && i < (256*1024)+10+10) || // double + (i >= (257*1024)+10 && i < (257*1024)+10+10) || // next to double + (i == (255*1024)+255) // quad + ) { + + string pixel_type = "x"; + if (i >= 58000 && i < 58000+10) { + pixel_type = "b"; + } else if (i >= 10 && i < 10+10) { + pixel_type = "e"; + } else if (i >= 1024+10 && i < 1024+10+10) { + pixel_type = "ie"; + } else if (i >= (256*1024)+10 && i < (256*1024)+10+10) { + pixel_type = "d"; + } else if (i >= (257*1024)+10 && i < (257*1024)+10+10) { + pixel_type = "ntd"; + } else if (i == (255*1024)+255) { + pixel_type = "q"; + } + + grap_g0->SetMarkerStyle(20); + grap_g1->SetMarkerStyle(20); + grap_g2->SetMarkerStyle(20); + grap_g0->SetMarkerColor(kBlue); + grap_g1->SetMarkerColor(kGreen+2); + grap_g2->SetMarkerColor(kRed); + grap_g0->SetLineColor(kBlue); + grap_g1->SetLineColor(kGreen+2); + grap_g2->SetLineColor(kRed); + + grap_g0->GetFunction("fit_g0")->Delete(); + grap_g0->FindObject("stats")->Delete(); + grap_g1->GetFunction("fit_g1")->Delete(); + grap_g1->FindObject("stats")->Delete(); + grap_g2->GetFunction("fit_g2")->Delete(); + grap_g2->FindObject("stats")->Delete(); + + if (module_str == "021" || module_str == "035" || module_str == "039" || module_str == "045" || module_str == "048" || module_str == "050" || module_str == "057" || module_str == "062" || module_str == "064" || module_str == "067" || module_str == "080" || module_str == "086" || module_str == "100" || module_str == "102" || module_str == "103" || module_str == "104" || module_str == "108" || module_str == "109" || module_str == "110" || module_str == "113" || module_str == "119" || module_str == "124" || module_str == "126") { + grap_g1->GetXaxis()->SetLimits(0.009,5E2); + } else { + grap_g1->GetXaxis()->SetLimits(0.009,1.1E2); + } + grap_g1->SetMinimum(0); + grap_g1->SetMaximum(16000); + grap_g1->GetXaxis()->SetTitle("Integration time [#mus]"); + grap_g1->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g1->GetYaxis()->SetTitleOffset(0.9); + grap_g1->Draw("AP"); + grap_g2->Draw("P"); + grap_g0->Draw("P"); + mapcanvas->SetLogx(); + sprintf(savename,"plots/M%s/CurrentSource/pixel_%s_%d_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + mapcanvas->SetLogx(0); + + } + } + + delete fit_g0; + delete fit_g1; + delete fit_g2; + delete grap_g0; + delete grap_g1; + delete grap_g2; + } } + mapcanvas->SetLeftMargin(0.1); + mapcanvas->SetRightMargin(0.13); + c1->cd(); TPaveText *pave = new TPaveText(0.86,0.95,0.91,0.98,"blNDC"); @@ -985,6 +1223,13 @@ int main(int argc, char* argv[]) { pave->SetTextSize(0.06); pave->SetTextAlign(32); + g0hist->GetXaxis()->SetTitle("G0 gradient"); + g0hist->GetXaxis()->SetTitleOffset(1.1); + g0hist->GetYaxis()->SetTitleOffset(1.5); + g0hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0hist_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); + g1hist->GetXaxis()->SetTitle("G1 gradient"); g1hist->GetXaxis()->SetTitleOffset(1.1); g1hist->GetYaxis()->SetTitleOffset(1.5); @@ -999,6 +1244,38 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g2hist_M%s.png", module_str.c_str(), module_str.c_str()); c1->SaveAs((const char *)(savename)); + g0overg1hist->GetXaxis()->SetTitle("G0 / G1"); + g0overg1hist->GetXaxis()->SetTitleOffset(1.1); + g0overg1hist->GetYaxis()->SetTitleOffset(1.5); + g0overg1hist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1hist_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + g0overg1hist->Fit("gaus"); + g0overg1hist->Draw(); + c1->Update(); + TPaveText *pave2 = new TPaveText(0.6,0.8,0.94,0.94,"blNDC"); + pave2->SetBorderSize(0); + pave2->SetFillStyle(0); + pave2->SetTextSize(0.04); + pave2->SetTextAlign(32); + TF1* g0overg1hist_gaus = g0overg1hist->GetFunction("gaus"); + sprintf(savename,"Mean %0.4f #pm %0.4f", g0overg1hist_gaus->GetParameter(1), g0overg1hist_gaus->GetParError(1)); + pave2->AddText((const char *)(savename)); + sprintf(savename,"Sigma %0.4f #pm %0.4f", g0overg1hist_gaus->GetParameter(2), g0overg1hist_gaus->GetParError(2)); + pave2->AddText((const char *)(savename)); + pave2->Draw(); + g0overg1hist->SetStats(kFALSE); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1hist_fit_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); + + g0overg1erhist->GetXaxis()->SetTitle("G0 / G1 uncert"); + g0overg1erhist->GetXaxis()->SetTitleOffset(1.1); + g0overg1erhist->GetYaxis()->SetTitleOffset(1.5); + g0overg1erhist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1erhist_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); + g1overg2hist->GetXaxis()->SetTitle("G1 / G2"); g1overg2hist->GetXaxis()->SetTitleOffset(1.1); g1overg2hist->GetYaxis()->SetTitleOffset(1.5); @@ -1009,11 +1286,7 @@ int main(int argc, char* argv[]) { g1overg2hist->Fit("gaus"); g1overg2hist->Draw(); c1->Update(); - TPaveText *pave2 = new TPaveText(0.6,0.8,0.94,0.94,"blNDC"); - pave2->SetBorderSize(0); - pave2->SetFillStyle(0); - pave2->SetTextSize(0.04); - pave2->SetTextAlign(32); + pave2->Clear(); TF1* g1overg2hist_gaus = g1overg2hist->GetFunction("gaus"); sprintf(savename,"Mean %0.4f #pm %0.4f", g1overg2hist_gaus->GetParameter(1), g1overg2hist_gaus->GetParError(1)); pave2->AddText((const char *)(savename)); @@ -1031,6 +1304,13 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g1overg2erhist_M%s.png", module_str.c_str(), module_str.c_str()); c1->SaveAs((const char *)(savename)); + g0erhist->GetXaxis()->SetTitle("Fitted G0 uncert"); + g0erhist->GetXaxis()->SetTitleOffset(1.1); + g0erhist->GetYaxis()->SetTitleOffset(1.5); + g0erhist->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0erhist_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); + g1erhist->GetXaxis()->SetTitle("Fitted G1 uncert"); g1erhist->GetXaxis()->SetTitleOffset(1.1); g1erhist->GetYaxis()->SetTitleOffset(1.5); @@ -1045,26 +1325,39 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g2erhist_M%s.png", module_str.c_str(), module_str.c_str()); c1->SaveAs((const char *)(savename)); - 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_M%s.png", module_str.c_str(), module_str.c_str()); - c1->SaveAs((const char *)(savename)); + g0overg1_isEdge->SetLineColor(kBlue); + g0overg1_isInnerEdge->SetLineColor(kCyan); + g0overg1_isDouble->SetLineColor(kGreen+2); + g0overg1_isNextToDouble->SetLineColor(kRed); + g0overg1_isQuad->SetLineColor(kOrange); - 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_M%s.png", module_str.c_str(), module_str.c_str()); - c1->SaveAs((const char *)(savename)); + g0overg1_isEdge->Scale(1./g0overg1_isEdge->GetEntries()); + g0overg1_isInnerEdge->Scale(1./g0overg1_isInnerEdge->GetEntries()); + g0overg1_isDouble->Scale(1./g0overg1_isDouble->GetEntries()); + g0overg1_isNextToDouble->Scale(1./g0overg1_isNextToDouble->GetEntries()); + g0overg1_isQuad->Scale(1./g0overg1_isQuad->GetEntries()); + g0overg1_isBulk->Scale(1./g0overg1_isBulk->GetEntries()); - 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; + TLegend *leg1 = new TLegend(0.62,0.6,0.93,0.93); + leg1->AddEntry(g0overg1_isBulk, "Normal", "l"); + leg1->AddEntry(g0overg1_isDouble, "Double", "l"); + leg1->AddEntry(g0overg1_isNextToDouble, "Next to D", "l"); + leg1->AddEntry(g0overg1_isEdge, "Edge", "l"); + leg1->AddEntry(g0overg1_isInnerEdge, "Inner E", "l"); + + g0overg1_isDouble->GetXaxis()->SetTitle("G0 / G1"); + g0overg1_isDouble->GetYaxis()->SetTitle("Normalised"); + g0overg1_isDouble->GetYaxis()->SetTitleOffset(1.3); + g0overg1_isDouble->SetMinimum(0.0); + g0overg1_isDouble->SetMaximum(0.16); + g0overg1_isDouble->Draw(); + g0overg1_isEdge->Draw("same"); + g0overg1_isInnerEdge->Draw("same"); + g0overg1_isNextToDouble->Draw("same"); + g0overg1_isBulk->Draw("same"); + leg1->Draw("same"); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1_perType_M%s.png", module_str.c_str(), module_str.c_str()); + c1->SaveAs((const char *)(savename)); g1overg2_isEdge->SetLineColor(kBlue); g1overg2_isInnerEdge->SetLineColor(kCyan); @@ -1102,6 +1395,16 @@ int main(int argc, char* argv[]) { mapcanvas->cd(); + g0map->GetXaxis()->SetTitle("Column"); + g0map->GetYaxis()->SetTitle("Row"); + g0map->GetYaxis()->SetTitleOffset(0.7); + g0map->GetZaxis()->SetRangeUser(20000,40000); + g0map->Draw("colz"); + pave->AddText("G0 [ADU/#mus]"); + pave->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0map_M%s.png", module_str.c_str(), module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + g1map->GetXaxis()->SetTitle("Column"); g1map->GetYaxis()->SetTitle("Row"); g1map->GetYaxis()->SetTitleOffset(0.7); @@ -1111,6 +1414,7 @@ int main(int argc, char* argv[]) { g1map->GetZaxis()->SetRangeUser(-3000,0); } g1map->Draw("colz"); + pave->Clear(); pave->AddText("G1 [ADU/#mus]"); pave->Draw(); sprintf(savename,"plots/M%s/CurrentSource/g1map_M%s.png", module_str.c_str(), module_str.c_str()); @@ -1131,6 +1435,28 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g2map_M%s.png", module_str.c_str(), module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); + g0overg1map->GetXaxis()->SetTitle("Column"); + g0overg1map->GetYaxis()->SetTitle("Row"); + g0overg1map->GetYaxis()->SetTitleOffset(0.7); + g0overg1map->GetZaxis()->SetRangeUser(-50,0); + g0overg1map->Draw("colz"); + pave->Clear(); + pave->AddText("G0/G1"); + pave->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1map_M%s.png", module_str.c_str(), module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + + g0overg1ermap->GetXaxis()->SetTitle("Column"); + g0overg1ermap->GetYaxis()->SetTitle("Row"); + g0overg1ermap->GetYaxis()->SetTitleOffset(0.7); + g0overg1ermap->GetZaxis()->SetRangeUser(0,0.5); + g0overg1ermap->Draw("colz"); + pave->Clear(); + pave->AddText("#sigma_{G0/G1}"); + pave->Draw(); + sprintf(savename,"plots/M%s/CurrentSource/g0overg1ermap_M%s.png", module_str.c_str(), module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + g1overg2map->GetXaxis()->SetTitle("Column"); g1overg2map->GetYaxis()->SetTitle("Row"); g1overg2map->GetYaxis()->SetTitleOffset(0.7); @@ -1153,6 +1479,17 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g1overg2ermap_M%s.png", module_str.c_str(), module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); + g0ermap->GetXaxis()->SetTitle("Column"); + g0ermap->GetYaxis()->SetTitle("Row"); + g0ermap->GetYaxis()->SetTitleOffset(0.7); + g0ermap->Draw("colz"); + pave->Clear(); + pave->AddText("#sigma_{G0}"); + pave->Draw(); + g0ermap->GetZaxis()->SetRangeUser(0,200); + sprintf(savename,"plots/M%s/CurrentSource/g0ermap_M%s.png", module_str.c_str(), module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); + g1ermap->GetXaxis()->SetTitle("Column"); g1ermap->GetYaxis()->SetTitle("Row"); g1ermap->GetYaxis()->SetTitleOffset(0.7); @@ -1175,22 +1512,6 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/CurrentSource/g2ermap_M%s.png", module_str.c_str(), 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.04,0.04); - sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG1map_M%s.png", module_str.c_str(), module_str.c_str()); - mapcanvas->SaveAs((const char *)(savename)); - - diffFromLinG2map->GetXaxis()->SetTitle("Column"); - diffFromLinG2map->GetYaxis()->SetTitle("Row"); - diffFromLinG2map->GetYaxis()->SetTitleOffset(0.7); - diffFromLinG2map->Draw("colz"); - diffFromLinG2map->GetZaxis()->SetRangeUser(-0.1,0.1); - sprintf(savename,"plots/M%s/CurrentSource/diffFromLinG2map_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();