Removed the pixel mask. no_calib maps are now accurate, and things aren't hidden.
This commit is contained in:
redford_s
729
BP_analysis.cpp
729
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<double> r0_adc;
|
||||
vector<double> r0_filter;
|
||||
vector<double> r0_adcerr;
|
||||
vector<double> r0_ferr;
|
||||
|
||||
vector<double> r0_adc;
|
||||
vector<double> r0_filter;
|
||||
vector<double> r0_adcerr;
|
||||
vector<double> r0_ferr;
|
||||
vector<double> r1_adc;
|
||||
vector<double> r1_filter;
|
||||
vector<double> r1_adcerr;
|
||||
vector<double> r1_ferr;
|
||||
|
||||
vector<double> r1_adc;
|
||||
vector<double> r1_filter;
|
||||
vector<double> r1_adcerr;
|
||||
vector<double> 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<double> 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<double> 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<double> 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<double> 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);
|
||||
|
||||
Reference in New Issue
Block a user