From 630559cdfe9d8a2b1f8cf0a52df19f3b5844e2f4 Mon Sep 17 00:00:00 2001 From: redford_s Date: Fri, 22 Feb 2019 17:21:51 +0100 Subject: [PATCH] Removed the pixel mask. no_calib maps are now accurate, and things aren't hidden. --- BP_analysis.cpp | 729 ++++++++++++++++++++++-------------------------- 1 file changed, 336 insertions(+), 393 deletions(-) diff --git a/BP_analysis.cpp b/BP_analysis.cpp index c5c4620..d421fd8 100644 --- a/BP_analysis.cpp +++ b/BP_analysis.cpp @@ -5,7 +5,6 @@ #include "../sls_detector_calibration/jungfrauFile.C" #include "../sls_detector_calibration/jungfrauPedestal.C" -#include "../sls_detector_calibration/jungfrauPixelMask.C" #include "TGraphErrors.h" #include "TF1.h" @@ -35,44 +34,6 @@ int main(int argc, char* argv[]) { jungfrauFile *thisfile = new jungfrauFile(); - jungfrauPixelMask *pixelMaskObject = new jungfrauPixelMask(); - bool pixel_mask[NCH]; - pixelMaskObject->initialisePixelMask(pixel_mask); - if (module_str == "006") { - pixelMaskObject->maskChip(6, pixel_mask); - pixelMaskObject->maskSupercolumn(4,1, pixel_mask); - pixelMaskObject->maskSupercolumn(4,2, pixel_mask); - } else if (module_str == "029") { - pixelMaskObject->maskChip(2, pixel_mask); - } else if (module_str == "035") { - pixelMaskObject->maskSupercolumn(2, 2, pixel_mask); - } else if (module_str == "039") { - pixelMaskObject->maskChip(7, pixel_mask); - } else if (module_str == "064") { - pixelMaskObject->maskSupercolumn(3,4,pixel_mask); - } else if (module_str == "065") { - pixelMaskObject->maskChip(5,pixel_mask); - } else if (module_str == "067") { - pixelMaskObject->maskChip(1,pixel_mask); - pixelMaskObject->maskChip(2,pixel_mask); - pixelMaskObject->maskChip(3,pixel_mask); - pixelMaskObject->maskChip(4,pixel_mask); - } else if (module_str == "086") { - pixelMaskObject->maskChip(8,pixel_mask); - } else if (module_str == "094") { - pixelMaskObject->maskChip(6,pixel_mask); - } else if (module_str == "121") { - pixelMaskObject->maskChip(1,pixel_mask); - } else if (module_str == "124") { - pixelMaskObject->maskChip(1,pixel_mask); - } else if (module_str == "125") { - pixelMaskObject->maskChip(7,pixel_mask); - } else if (module_str == "223") { - pixelMaskObject->maskChip(2,pixel_mask); - } else if (module_str == "233") { - pixelMaskObject->maskChip(7,pixel_mask); - } - jungfrauPedestal *pedestalObjectG0 = new jungfrauPedestal(); pedestalObjectG0->pedestalSetNFrames(100); // using 100 frames, rolling window jungfrauPedestal *pedestalObjectG1 = new jungfrauPedestal(); @@ -115,7 +76,6 @@ int main(int argc, char* argv[]) { g1_grader_hist = new TH1F("g1_grader_hist","",100,0,0.0001); } - TH2F* pixelmask_map = new TH2F("pixelmask_map","",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* g0_grad_map = new TH2F("g0_grad_map","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); @@ -216,15 +176,13 @@ int main(int argc, char* argv[]) { } for (int i = 0; i < NCH; i++) { - if (pixel_mask[i] == true) { - avg_adcG0_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG0->pedestalOfChannel(i)); - avg_adcG0er_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG0->semOfChannel(i)); + avg_adcG0_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG0->pedestalOfChannel(i)); + avg_adcG0er_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG0->semOfChannel(i)); - avg_adcG1_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG1->pedestalOfChannel(i)); - avg_adcG1er_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG1->semOfChannel(i)); + avg_adcG1_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG1->pedestalOfChannel(i)); + avg_adcG1er_map[framecounter/100]->SetBinContent((i%NC)+1,(i/NC)+1,pedestalObjectG1->semOfChannel(i)); - } } pedestalObjectG0->pedestalClear(); @@ -307,67 +265,257 @@ int main(int argc, char* argv[]) { for (int i = 0; i < NCH; i++) { - if (pixel_mask[i] == true) { + vector r0_adc; + vector r0_filter; + vector r0_adcerr; + vector r0_ferr; - vector r0_adc; - vector r0_filter; - vector r0_adcerr; - vector r0_ferr; + vector r1_adc; + vector r1_filter; + vector r1_adcerr; + vector r1_ferr; - vector r1_adc; - vector r1_filter; - vector r1_adcerr; - vector r1_ferr; + for (int j = 0; j < 220; j++) { - for (int j = 0; j < 220; j++) { + double this_g0avg = avg_adcG0_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); + double this_g1avg = avg_adcG1_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); + double this_g0avger = avg_adcG0er_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); + double this_g1avger = avg_adcG1er_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); - double this_g0avg = avg_adcG0_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); - double this_g1avg = avg_adcG1_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); - double this_g0avger = avg_adcG0er_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); - double this_g1avger = avg_adcG1er_map[j]->GetBinContent((i%NC)+1,(i/NC)+1); + if (this_g0avg != 0) { + r0_filter.push_back(xs[j]); + r0_ferr.push_back(0.); + r0_adc.push_back(this_g0avg); + r0_adcerr.push_back(this_g0avger); + } + if (this_g1avg != 0) { + r1_filter.push_back(xs[j]); + r1_ferr.push_back(0.); + r1_adc.push_back(this_g1avg); + r1_adcerr.push_back(this_g1avger); + } + } + + TGraphErrors *grap_g0 = 0; + TGraphErrors *grap_g1 = 0; + + TF1 *fit_g0 = 0; + TF1 *fit_g1 = 0; + + double rangemin_g0 = 0; + double rangemax_g0 = 0; + double rangemin_g1 = 0; + double rangemax_g1 = 0; + + // define graphs + if (r0_adc.size() > 1) { + grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); + grap_g0->SetMarkerStyle(20); + grap_g0->SetMarkerColor(kBlue); + grap_g0->SetLineColor(kBlue); + } + + if (r1_adc.size() > 1) { + grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); + grap_g1->SetMarkerStyle(20); + grap_g1->SetMarkerColor(kGreen+2); + grap_g1->SetLineColor(kGreen+2); + } + + // plot the datapoints + if (r1_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_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); + grap_g1->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g1->GetYaxis()->SetTitleOffset(0.9); + grap_g1->SetMinimum(1000); + grap_g1->SetMaximum(15000); + grap_g1->GetXaxis()->SetLimits(0,7200); + grap_g1->Draw("AP"); + if (r0_adc.size() > 1) { + grap_g0->Draw("P"); + } + mapcanvas->Update(); + sprintf(savename,"plots/M%s/BackplanePulsing/pixel_%s_%d_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + mapcanvas->SaveAs((const char *)(savename)); - if (this_g0avg != 0) { - r0_filter.push_back(xs[j]); - r0_ferr.push_back(0.); - r0_adc.push_back(this_g0avg); - r0_adcerr.push_back(this_g0avger); - } - if (this_g1avg != 0) { - r1_filter.push_back(xs[j]); - r1_ferr.push_back(0.); - r1_adc.push_back(this_g1avg); - r1_adcerr.push_back(this_g1avger); - } } + } - TGraphErrors *grap_g0 = 0; - TGraphErrors *grap_g1 = 0; + // define fit ranges and fit + if (r0_adc.size() > 1) { - TF1 *fit_g0 = 0; - TF1 *fit_g1 = 0; + rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end()); + rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter); - double rangemin_g0 = 0; - double rangemax_g0 = 0; - double rangemin_g1 = 0; - double rangemax_g1 = 0; + if (rangemax_g0 > rangemin_g0) { - // define graphs - if (r0_adc.size() > 1) { - grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); - grap_g0->SetMarkerStyle(20); - grap_g0->SetMarkerColor(kBlue); - grap_g0->SetLineColor(kBlue); + grap_g0->Fit("pol1","QRC","",rangemin_g0,rangemax_g0); + fit_g0 = (TF1*) grap_g0->GetFunction("pol1"); + fit_g0->SetLineColor(kBlue); + fit_g0->SetParName(0,"G0 const"); + fit_g0->SetParName(1,"G0 grad"); + + g0_grad_hist->Fill(fit_g0->GetParameter(1)); + g0_grader_hist->Fill(fit_g0->GetParError(1)); + g0_grad_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)); + g0_grader_map->Fill(i%NC,i/NC,fit_g0->GetParError(1)); + g0_const_hist->Fill(fit_g0->GetParameter(0)); + g0_conster_hist->Fill(fit_g0->GetParError(0)); + g0_const_map->Fill(i%NC,i/NC,fit_g0->GetParameter(0)); + g0_conster_map->Fill(i%NC,i/NC,fit_g0->GetParError(0)); + + 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->GetXaxis()->SetTitle("Signal generator voltage [mV]"); + grap_g0->GetYaxis()->SetTitle("ADC [ADU]"); + grap_g0->GetYaxis()->SetTitleOffset(0.9); + grap_g0->SetMinimum(1000); + grap_g0->SetMaximum(15000); + grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); + grap_g0->Draw("AP"); + fit_g0->Draw("same"); + mapcanvas->Update(); + TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); + st0->SetX1NDC(0.2); + st0->SetX2NDC(0.54); + st0->SetY1NDC(0.18); + st0->SetY2NDC(0.37); + st0->SetBorderSize(0); + st0->SetTextColor(kBlue); + sprintf(savename,"plots/M%s/BackplanePulsing/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)); + + 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); + + TF1* lin_g0_p1pc = new TF1("lin_g0_p1pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_p1pc->SetParameter(0,fit_g0->GetParameter(0)/100.); + lin_g0_p1pc->SetParameter(1,fit_g0->GetParameter(1)/100.); + lin_g0_p1pc->SetLineColor(kRed); + + TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_p05pc->SetParameter(0,fit_g0->GetParameter(0)/200.); + lin_g0_p05pc->SetParameter(1,fit_g0->GetParameter(1)/200.); + lin_g0_p05pc->SetLineColor(kOrange+1); + + TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_p02pc->SetParameter(0,fit_g0->GetParameter(0)/500.); + lin_g0_p02pc->SetParameter(1,fit_g0->GetParameter(1)/500.); + lin_g0_p02pc->SetLineColor(kOrange); + + TF1* lin_g0_m1pc = new TF1("lin_g0_m1pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_m1pc->SetParameter(0,fit_g0->GetParameter(0)/-100.); + lin_g0_m1pc->SetParameter(1,fit_g0->GetParameter(1)/-100.); + lin_g0_m1pc->SetLineColor(kRed); + + TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_m05pc->SetParameter(0,fit_g0->GetParameter(0)/-200.); + lin_g0_m05pc->SetParameter(1,fit_g0->GetParameter(1)/-200.); + lin_g0_m05pc->SetLineColor(kOrange+1); + + TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin_g0,rangemax_g0); + lin_g0_m02pc->SetParameter(0,fit_g0->GetParameter(0)/-500.); + lin_g0_m02pc->SetParameter(1,fit_g0->GetParameter(1)/-500.); + lin_g0_m02pc->SetLineColor(kOrange); + + norm_g0->GetXaxis()->SetRangeUser(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); + norm_g0->SetMinimum(1.5*lin_g0_m02pc->Eval(rangemax_g0)); + norm_g0->SetMaximum(1.5*lin_g0_p02pc->Eval(rangemax_g0)); + norm_g0->GetXaxis()->SetTitle("Signal generator voltage [mV]"); + norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g0->GetYaxis()->SetTitleOffset(0.9); + norm_g0->Draw("AP"); + flat_g0->Draw("same"); + lin_g0_p1pc->Draw("same"); + lin_g0_p05pc->Draw("same"); + lin_g0_p02pc->Draw("same"); + lin_g0_m1pc->Draw("same"); + lin_g0_m05pc->Draw("same"); + lin_g0_m02pc->Draw("same"); + norm_g0->Draw("P"); + sprintf(savename,"plots/M%s/BackplanePulsing/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; + } } + } - if (r1_adc.size() > 1) { - grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); - grap_g1->SetMarkerStyle(20); - grap_g1->SetMarkerColor(kGreen+2); - grap_g1->SetLineColor(kGreen+2); - } + if (r1_adc.size() > 1) { - // plot the datapoints - if (r1_adc.size() > 1) { + rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter); + rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end()); + + if (rangemax_g1 > rangemin_g1) { + + grap_g1->Fit("pol1","QRC","",rangemin_g1,rangemax_g1); + fit_g1 = (TF1*) grap_g1->GetFunction("pol1"); + fit_g1->SetLineColor(kGreen+2); + fit_g1->SetParName(0,"G1 const"); + fit_g1->SetParName(1,"G1 grad"); + + g1_grad_hist->Fill(fit_g1->GetParameter(1)); + g1_grader_hist->Fill(fit_g1->GetParError(1)); + g1_grad_map->Fill(i%NC,i/NC,fit_g1->GetParameter(1)); + g1_grader_map->Fill(i%NC,i/NC,fit_g1->GetParError(1)); + g1_const_hist->Fill(fit_g1->GetParameter(0)); + g1_conster_hist->Fill(fit_g1->GetParError(0)); + g1_const_map->Fill(i%NC,i/NC,fit_g1->GetParameter(0)); + g1_conster_map->Fill(i%NC,i/NC,fit_g1->GetParError(0)); if ((i >= 58000 && i < 58000+10) || // bulk (i >= 10 && i < 10+10) || // edge @@ -397,301 +545,108 @@ int main(int argc, char* argv[]) { grap_g1->GetYaxis()->SetTitleOffset(0.9); grap_g1->SetMinimum(1000); grap_g1->SetMaximum(15000); - grap_g1->GetXaxis()->SetLimits(0,7200); + grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); grap_g1->Draw("AP"); - if (r0_adc.size() > 1) { - grap_g0->Draw("P"); - } + fit_g1->Draw("same"); mapcanvas->Update(); - sprintf(savename,"plots/M%s/BackplanePulsing/pixel_%s_%d_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); + TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); + st0->SetX1NDC(0.6); + st0->SetX2NDC(0.94); + st0->SetY1NDC(0.18); + st0->SetY2NDC(0.37); + st0->SetBorderSize(0); + st0->SetTextColor(kGreen+2); + sprintf(savename,"plots/M%s/BackplanePulsing/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)); - } - } - - // define fit ranges and fit - if (r0_adc.size() > 1) { - - rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end()); - rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter); - - if (rangemax_g0 > rangemin_g0) { - - grap_g0->Fit("pol1","QRC","",rangemin_g0,rangemax_g0); - fit_g0 = (TF1*) grap_g0->GetFunction("pol1"); - fit_g0->SetLineColor(kBlue); - fit_g0->SetParName(0,"G0 const"); - fit_g0->SetParName(1,"G0 grad"); - - g0_grad_hist->Fill(fit_g0->GetParameter(1)); - g0_grader_hist->Fill(fit_g0->GetParError(1)); - g0_grad_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)); - g0_grader_map->Fill(i%NC,i/NC,fit_g0->GetParError(1)); - g0_const_hist->Fill(fit_g0->GetParameter(0)); - g0_conster_hist->Fill(fit_g0->GetParError(0)); - g0_const_map->Fill(i%NC,i/NC,fit_g0->GetParameter(0)); - g0_conster_map->Fill(i%NC,i/NC,fit_g0->GetParError(0)); - - 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->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap_g0->GetYaxis()->SetTitle("ADC [ADU]"); - grap_g0->GetYaxis()->SetTitleOffset(0.9); - grap_g0->SetMinimum(1000); - grap_g0->SetMaximum(15000); - grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); - grap_g0->Draw("AP"); - fit_g0->Draw("same"); - mapcanvas->Update(); - TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); - st0->SetX1NDC(0.2); - st0->SetX2NDC(0.54); - st0->SetY1NDC(0.18); - st0->SetY2NDC(0.37); - st0->SetBorderSize(0); - st0->SetTextColor(kBlue); - sprintf(savename,"plots/M%s/BackplanePulsing/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)); - - 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); - - TF1* lin_g0_p1pc = new TF1("lin_g0_p1pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_p1pc->SetParameter(0,fit_g0->GetParameter(0)/100.); - lin_g0_p1pc->SetParameter(1,fit_g0->GetParameter(1)/100.); - lin_g0_p1pc->SetLineColor(kRed); - - TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_p05pc->SetParameter(0,fit_g0->GetParameter(0)/200.); - lin_g0_p05pc->SetParameter(1,fit_g0->GetParameter(1)/200.); - lin_g0_p05pc->SetLineColor(kOrange+1); - - TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_p02pc->SetParameter(0,fit_g0->GetParameter(0)/500.); - lin_g0_p02pc->SetParameter(1,fit_g0->GetParameter(1)/500.); - lin_g0_p02pc->SetLineColor(kOrange); - - TF1* lin_g0_m1pc = new TF1("lin_g0_m1pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_m1pc->SetParameter(0,fit_g0->GetParameter(0)/-100.); - lin_g0_m1pc->SetParameter(1,fit_g0->GetParameter(1)/-100.); - lin_g0_m1pc->SetLineColor(kRed); - - TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_m05pc->SetParameter(0,fit_g0->GetParameter(0)/-200.); - lin_g0_m05pc->SetParameter(1,fit_g0->GetParameter(1)/-200.); - lin_g0_m05pc->SetLineColor(kOrange+1); - - TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin_g0,rangemax_g0); - lin_g0_m02pc->SetParameter(0,fit_g0->GetParameter(0)/-500.); - lin_g0_m02pc->SetParameter(1,fit_g0->GetParameter(1)/-500.); - lin_g0_m02pc->SetLineColor(kOrange); - - norm_g0->GetXaxis()->SetRangeUser(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); - norm_g0->SetMinimum(1.5*lin_g0_m02pc->Eval(rangemax_g0)); - norm_g0->SetMaximum(1.5*lin_g0_p02pc->Eval(rangemax_g0)); - norm_g0->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - norm_g0->GetYaxis()->SetTitleOffset(0.9); - norm_g0->Draw("AP"); - flat_g0->Draw("same"); - lin_g0_p1pc->Draw("same"); - lin_g0_p05pc->Draw("same"); - lin_g0_p02pc->Draw("same"); - lin_g0_m1pc->Draw("same"); - lin_g0_m05pc->Draw("same"); - lin_g0_m02pc->Draw("same"); - norm_g0->Draw("P"); - sprintf(savename,"plots/M%s/BackplanePulsing/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; + vector r1_adc_norm; + for (size_t j = 0; j < r1_adc.size(); j++) { + r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j])); } + + 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_gi","0",rangemin_g1,rangemax_g1); + flat_g1->SetLineColor(kGreen+2); + + TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.); + lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.); + lin_g1_p02pc->SetLineColor(kOrange); + + TF1* lin_g1_p01pc = new TF1("lin_g1_p01pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_p01pc->SetParameter(0,fit_g1->GetParameter(0)/1000.); + lin_g1_p01pc->SetParameter(1,fit_g1->GetParameter(1)/1000.); + lin_g1_p01pc->SetLineColor(kYellow); + + TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.); + lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.); + lin_g1_m02pc->SetLineColor(kOrange); + + TF1* lin_g1_m01pc = new TF1("lin_g1_m01pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_m01pc->SetParameter(0,fit_g1->GetParameter(0)/-1000.); + lin_g1_m01pc->SetParameter(1,fit_g1->GetParameter(1)/-1000.); + lin_g1_m01pc->SetLineColor(kYellow); + + norm_g1->GetXaxis()->SetRangeUser(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); + norm_g1->SetMinimum(1.5*lin_g1_m02pc->Eval(rangemin_g1)); + norm_g1->SetMaximum(1.5*lin_g1_p02pc->Eval(rangemin_g1)); + norm_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); + norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g1->GetYaxis()->SetTitleOffset(0.9); + norm_g1->Draw("AP"); + flat_g1->Draw("same"); + lin_g1_p02pc->Draw("same"); + lin_g1_p01pc->Draw("same"); + lin_g1_m02pc->Draw("same"); + lin_g1_m01pc->Draw("same"); + norm_g1->Draw("P"); + sprintf(savename,"plots/M%s/BackplanePulsing/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)); + + delete norm_g1; } } - - if (r1_adc.size() > 1) { - - rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter); - rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end()); - - if (rangemax_g1 > rangemin_g1) { - - grap_g1->Fit("pol1","QRC","",rangemin_g1,rangemax_g1); - fit_g1 = (TF1*) grap_g1->GetFunction("pol1"); - fit_g1->SetLineColor(kGreen+2); - fit_g1->SetParName(0,"G1 const"); - fit_g1->SetParName(1,"G1 grad"); - - g1_grad_hist->Fill(fit_g1->GetParameter(1)); - g1_grader_hist->Fill(fit_g1->GetParError(1)); - g1_grad_map->Fill(i%NC,i/NC,fit_g1->GetParameter(1)); - g1_grader_map->Fill(i%NC,i/NC,fit_g1->GetParError(1)); - g1_const_hist->Fill(fit_g1->GetParameter(0)); - g1_conster_hist->Fill(fit_g1->GetParError(0)); - g1_const_map->Fill(i%NC,i/NC,fit_g1->GetParameter(0)); - g1_conster_map->Fill(i%NC,i/NC,fit_g1->GetParError(0)); - - 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->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap_g1->GetYaxis()->SetTitle("ADC [ADU]"); - grap_g1->GetYaxis()->SetTitleOffset(0.9); - grap_g1->SetMinimum(1000); - grap_g1->SetMaximum(15000); - grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); - grap_g1->Draw("AP"); - fit_g1->Draw("same"); - mapcanvas->Update(); - TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); - st0->SetX1NDC(0.6); - st0->SetX2NDC(0.94); - st0->SetY1NDC(0.18); - st0->SetY2NDC(0.37); - st0->SetBorderSize(0); - st0->SetTextColor(kGreen+2); - sprintf(savename,"plots/M%s/BackplanePulsing/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)); - - vector r1_adc_norm; - for (size_t j = 0; j < r1_adc.size(); j++) { - r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j])); - } - - 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_gi","0",rangemin_g1,rangemax_g1); - flat_g1->SetLineColor(kGreen+2); - - TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1); - lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.); - lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.); - lin_g1_p02pc->SetLineColor(kOrange); - - TF1* lin_g1_p01pc = new TF1("lin_g1_p01pc","[0]+[1]*x",rangemin_g1,rangemax_g1); - lin_g1_p01pc->SetParameter(0,fit_g1->GetParameter(0)/1000.); - lin_g1_p01pc->SetParameter(1,fit_g1->GetParameter(1)/1000.); - lin_g1_p01pc->SetLineColor(kYellow); - - TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1); - lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.); - lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.); - lin_g1_m02pc->SetLineColor(kOrange); - - TF1* lin_g1_m01pc = new TF1("lin_g1_m01pc","[0]+[1]*x",rangemin_g1,rangemax_g1); - lin_g1_m01pc->SetParameter(0,fit_g1->GetParameter(0)/-1000.); - lin_g1_m01pc->SetParameter(1,fit_g1->GetParameter(1)/-1000.); - lin_g1_m01pc->SetLineColor(kYellow); - - norm_g1->GetXaxis()->SetRangeUser(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); - norm_g1->SetMinimum(1.5*lin_g1_m02pc->Eval(rangemin_g1)); - norm_g1->SetMaximum(1.5*lin_g1_p02pc->Eval(rangemin_g1)); - norm_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - norm_g1->GetYaxis()->SetTitleOffset(0.9); - norm_g1->Draw("AP"); - flat_g1->Draw("same"); - lin_g1_p02pc->Draw("same"); - lin_g1_p01pc->Draw("same"); - lin_g1_m02pc->Draw("same"); - lin_g1_m01pc->Draw("same"); - norm_g1->Draw("P"); - sprintf(savename,"plots/M%s/BackplanePulsing/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)); - - delete norm_g1; - } - } - } - - // get ratio measurements - if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) { - - double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1); - double this_g0overg1er = sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2)); - - g0overg1_map->Fill(i%NC,i/NC,this_g0overg1); - g0overg1_hist->Fill(this_g0overg1); - g0overg1er_map->Fill(i%NC,i/NC,this_g0overg1er); - g0overg1er_hist->Fill(this_g0overg1er); - - if (isEdge(i)) { - g0overg1_hist_isEdge->Fill(this_g0overg1); - } - if (isInnerEdge(i)) { - g0overg1_hist_isInnerEdge->Fill(this_g0overg1); - } - if (isDouble(i)) { - g0overg1_hist_isDouble->Fill(this_g0overg1); - } - if (isNextToDouble(i)) { - g0overg1_hist_isNextToDouble->Fill(this_g0overg1); - } - if (isQuad(i)) { - g0overg1_hist_isQuad->Fill(this_g0overg1); - } - if (isBulk(i)) { - g0overg1_hist_isBulk->Fill(this_g0overg1); - } - } else { - no_g0g1_calib_map->Fill(i%NC,i/NC,1); - } - - delete grap_g0; - delete grap_g1; } + + // get ratio measurements + if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) { + + double this_g0overg1 = fit_g0->GetParameter(1)/fit_g1->GetParameter(1); + double this_g0overg1er = sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2)); + + g0overg1_map->Fill(i%NC,i/NC,this_g0overg1); + g0overg1_hist->Fill(this_g0overg1); + g0overg1er_map->Fill(i%NC,i/NC,this_g0overg1er); + g0overg1er_hist->Fill(this_g0overg1er); + + if (isEdge(i)) { + g0overg1_hist_isEdge->Fill(this_g0overg1); + } + if (isInnerEdge(i)) { + g0overg1_hist_isInnerEdge->Fill(this_g0overg1); + } + if (isDouble(i)) { + g0overg1_hist_isDouble->Fill(this_g0overg1); + } + if (isNextToDouble(i)) { + g0overg1_hist_isNextToDouble->Fill(this_g0overg1); + } + if (isQuad(i)) { + g0overg1_hist_isQuad->Fill(this_g0overg1); + } + if (isBulk(i)) { + g0overg1_hist_isBulk->Fill(this_g0overg1); + } + } else { + no_g0g1_calib_map->Fill(i%NC,i/NC,1); + } + + delete grap_g0; + delete grap_g1; } mapcanvas->SetLeftMargin(0.1); @@ -940,18 +895,6 @@ int main(int argc, char* argv[]) { sprintf(savename,"plots/M%s/BackplanePulsing/g0overg1er_map_M%s.png", module_str.c_str(), module_str.c_str()); mapcanvas->SaveAs((const char *)(savename)); - for (int i = 0; i < NCH; i++) { - if (pixel_mask[i] == false) { - pixelmask_map->Fill(i%NC,i/NC,1); - } - } - pixelmask_map->GetXaxis()->SetTitle("Column"); - pixelmask_map->GetYaxis()->SetTitle("Row"); - pixelmask_map->GetYaxis()->SetTitleOffset(0.7); - pixelmask_map->Draw("colz"); - sprintf(savename,"plots/M%s/BackplanePulsing/pixelmask_map_M%s.png", module_str.c_str(), module_str.c_str()); - 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);