detectors/JFCalibration
0
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

Consistent variable names.

Browse Source
This commit is contained in:
redford_s
2019-02-08 15:07:09 +01:00
parent c4ebbf4666
commit de0ef22d58
2 changed files with 252 additions and 252 deletions
Download Patch File Download Diff File Expand all files Collapse all files

292
BP_analysis.cpp
View File

@ -299,8 +299,6 @@ int main(int argc, char* argv[]) {
for (int i = 0; i < NCH; i++) { for (int i = 0; i < NCH; i++) {
if (i%100000==0) {cout << "another 100k" << endl;}
if (pixel_mask[i] == true) { if (pixel_mask[i] == true) {
vector<double> r0_adc; vector<double> r0_adc;
@ -336,30 +334,30 @@ int main(int argc, char* argv[]) {
} }
} }
TGraphErrors *grap0 = 0; TGraphErrors *grap_g0 = 0;
TGraphErrors *grap1 = 0; TGraphErrors *grap_g1 = 0;
TF1 *g0_fit = 0; TF1 *fit_g0 = 0;
TF1 *g1_fit = 0; TF1 *fit_g1 = 0;
double g0min = 0; double rangemin_g0 = 0;
double g0max = 0; double rangemax_g0 = 0;
double g1min = 0; double rangemin_g1 = 0;
double g1max = 0; double rangemax_g1 = 0;
// define graphs // define graphs
if (r0_adc.size() > 1) { if (r0_adc.size() > 1) {
grap0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
grap0->SetMarkerStyle(20); grap_g0->SetMarkerStyle(20);
grap0->SetMarkerColor(kBlue); grap_g0->SetMarkerColor(kBlue);
grap0->SetLineColor(kBlue); grap_g0->SetLineColor(kBlue);
} }
if (r1_adc.size() > 1) { if (r1_adc.size() > 1) {
grap1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
grap1->SetMarkerStyle(20); grap_g1->SetMarkerStyle(20);
grap1->SetMarkerColor(kGreen+2); grap_g1->SetMarkerColor(kGreen+2);
grap1->SetLineColor(kGreen+2); grap_g1->SetLineColor(kGreen+2);
} }
// plot the datapoints // plot the datapoints
@ -391,15 +389,15 @@ int main(int argc, char* argv[]) {
mapcanvas->SetLeftMargin(0.13); mapcanvas->SetLeftMargin(0.13);
mapcanvas->SetRightMargin(0.05); mapcanvas->SetRightMargin(0.05);
grap1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); grap_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]");
grap1->GetYaxis()->SetTitle("ADC [ADU]"); grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
grap1->GetYaxis()->SetTitleOffset(0.9); grap_g1->GetYaxis()->SetTitleOffset(0.9);
grap1->SetMinimum(1000); grap_g1->SetMinimum(1000);
grap1->SetMaximum(15000); grap_g1->SetMaximum(15000);
grap1->GetXaxis()->SetLimits(0,7200); grap_g1->GetXaxis()->SetLimits(0,7200);
grap1->Draw("AP"); grap_g1->Draw("AP");
if (r0_adc.size() > 1) { if (r0_adc.size() > 1) {
grap0->Draw("P"); grap_g0->Draw("P");
} }
mapcanvas->Update(); 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()); 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());
@ -411,25 +409,25 @@ int main(int argc, char* argv[]) {
// define fit ranges and fit // define fit ranges and fit
if (r0_adc.size() > 1) { if (r0_adc.size() > 1) {
g0min = *min_element(r0_filter.begin(),r0_filter.end()); rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end());
g0max = highestPointBeforeSwitching(r0_filter,r1_filter); rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter);
if (g0max > g0min) { if (rangemax_g0 > rangemin_g0) {
grap0->Fit("pol1","QRC","",g0min,g0max); grap_g0->Fit("pol1","QRC","",rangemin_g0,rangemax_g0);
g0_fit = (TF1*) grap0->GetFunction("pol1"); fit_g0 = (TF1*) grap_g0->GetFunction("pol1");
g0_fit->SetLineColor(kBlue); fit_g0->SetLineColor(kBlue);
g0_fit->SetParName(0,"G0 const"); fit_g0->SetParName(0,"G0 const");
g0_fit->SetParName(1,"G0 grad"); fit_g0->SetParName(1,"G0 grad");
g0_grad_hist->Fill(g0_fit->GetParameter(1)); g0_grad_hist->Fill(fit_g0->GetParameter(1));
g0_grader_hist->Fill(g0_fit->GetParError(1)); g0_grader_hist->Fill(fit_g0->GetParError(1));
g0_grad_map->Fill(i%NC,i/NC,g0_fit->GetParameter(1)); g0_grad_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1));
g0_grader_map->Fill(i%NC,i/NC,g0_fit->GetParError(1)); g0_grader_map->Fill(i%NC,i/NC,fit_g0->GetParError(1));
g0_const_hist->Fill(g0_fit->GetParameter(0)); g0_const_hist->Fill(fit_g0->GetParameter(0));
g0_conster_hist->Fill(g0_fit->GetParError(0)); g0_conster_hist->Fill(fit_g0->GetParError(0));
g0_const_map->Fill(i%NC,i/NC,g0_fit->GetParameter(0)); g0_const_map->Fill(i%NC,i/NC,fit_g0->GetParameter(0));
g0_conster_map->Fill(i%NC,i/NC,g0_fit->GetParError(0)); g0_conster_map->Fill(i%NC,i/NC,fit_g0->GetParError(0));
if ((i >= 58000 && i < 58000+10) || // bulk if ((i >= 58000 && i < 58000+10) || // bulk
(i >= 10 && i < 10+10) || // edge (i >= 10 && i < 10+10) || // edge
@ -457,16 +455,16 @@ int main(int argc, char* argv[]) {
mapcanvas->SetLeftMargin(0.13); mapcanvas->SetLeftMargin(0.13);
mapcanvas->SetRightMargin(0.05); mapcanvas->SetRightMargin(0.05);
grap0->GetXaxis()->SetTitle("Signal generator voltage [mV]"); grap_g0->GetXaxis()->SetTitle("Signal generator voltage [mV]");
grap0->GetYaxis()->SetTitle("ADC [ADU]"); grap_g0->GetYaxis()->SetTitle("ADC [ADU]");
grap0->GetYaxis()->SetTitleOffset(0.9); grap_g0->GetYaxis()->SetTitleOffset(0.9);
grap0->SetMinimum(1000); grap_g0->SetMinimum(1000);
grap0->SetMaximum(15000); grap_g0->SetMaximum(15000);
grap0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
grap0->Draw("AP"); grap_g0->Draw("AP");
g0_fit->Draw("same"); fit_g0->Draw("same");
mapcanvas->Update(); mapcanvas->Update();
TPaveStats *st0 = (TPaveStats*)grap0->FindObject("stats"); TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats");
st0->SetX1NDC(0.2); st0->SetX1NDC(0.2);
st0->SetX2NDC(0.54); st0->SetX2NDC(0.54);
st0->SetY1NDC(0.18); st0->SetY1NDC(0.18);
@ -478,53 +476,53 @@ int main(int argc, char* argv[]) {
vector<double> r0_adc_norm; vector<double> r0_adc_norm;
for (size_t j = 0; j < r0_adc.size(); j++) { for (size_t j = 0; j < r0_adc.size(); j++) {
r0_adc_norm.push_back(r0_adc[j] - g0_fit->Eval(r0_filter[j])); r0_adc_norm.push_back(r0_adc[j] - fit_g0->Eval(r0_filter[j]));
} }
TGraphErrors *grap0_norm = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcerr[0])); TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcerr[0]));
grap0_norm->SetMarkerColor(kBlue); norm_g0->SetMarkerColor(kBlue);
grap0_norm->SetLineColor(kBlue); norm_g0->SetLineColor(kBlue);
TF1* flat_g0 = new TF1("flat_g0","0",g0min,g0max); TF1* flat_g0 = new TF1("flat_g0","0",rangemin_g0,rangemax_g0);
flat_g0->SetLineColor(kBlue); flat_g0->SetLineColor(kBlue);
TF1* lin_g0_p1pc = new TF1("lin_g0_p1pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_p1pc = new TF1("lin_g0_p1pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_p1pc->SetParameter(0,g0_fit->GetParameter(0)/100.); lin_g0_p1pc->SetParameter(0,fit_g0->GetParameter(0)/100.);
lin_g0_p1pc->SetParameter(1,g0_fit->GetParameter(1)/100.); lin_g0_p1pc->SetParameter(1,fit_g0->GetParameter(1)/100.);
lin_g0_p1pc->SetLineColor(kRed); lin_g0_p1pc->SetLineColor(kRed);
TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_p05pc->SetParameter(0,g0_fit->GetParameter(0)/200.); lin_g0_p05pc->SetParameter(0,fit_g0->GetParameter(0)/200.);
lin_g0_p05pc->SetParameter(1,g0_fit->GetParameter(1)/200.); lin_g0_p05pc->SetParameter(1,fit_g0->GetParameter(1)/200.);
lin_g0_p05pc->SetLineColor(kOrange+1); lin_g0_p05pc->SetLineColor(kOrange+1);
TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_p02pc->SetParameter(0,g0_fit->GetParameter(0)/500.); lin_g0_p02pc->SetParameter(0,fit_g0->GetParameter(0)/500.);
lin_g0_p02pc->SetParameter(1,g0_fit->GetParameter(1)/500.); lin_g0_p02pc->SetParameter(1,fit_g0->GetParameter(1)/500.);
lin_g0_p02pc->SetLineColor(kOrange); lin_g0_p02pc->SetLineColor(kOrange);
TF1* lin_g0_m1pc = new TF1("lin_g0_m1pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_m1pc = new TF1("lin_g0_m1pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_m1pc->SetParameter(0,g0_fit->GetParameter(0)/-100.); lin_g0_m1pc->SetParameter(0,fit_g0->GetParameter(0)/-100.);
lin_g0_m1pc->SetParameter(1,g0_fit->GetParameter(1)/-100.); lin_g0_m1pc->SetParameter(1,fit_g0->GetParameter(1)/-100.);
lin_g0_m1pc->SetLineColor(kRed); lin_g0_m1pc->SetLineColor(kRed);
TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_m05pc->SetParameter(0,g0_fit->GetParameter(0)/-200.); lin_g0_m05pc->SetParameter(0,fit_g0->GetParameter(0)/-200.);
lin_g0_m05pc->SetParameter(1,g0_fit->GetParameter(1)/-200.); lin_g0_m05pc->SetParameter(1,fit_g0->GetParameter(1)/-200.);
lin_g0_m05pc->SetLineColor(kOrange+1); lin_g0_m05pc->SetLineColor(kOrange+1);
TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",g0min,g0max); TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_m02pc->SetParameter(0,g0_fit->GetParameter(0)/-500.); lin_g0_m02pc->SetParameter(0,fit_g0->GetParameter(0)/-500.);
lin_g0_m02pc->SetParameter(1,g0_fit->GetParameter(1)/-500.); lin_g0_m02pc->SetParameter(1,fit_g0->GetParameter(1)/-500.);
lin_g0_m02pc->SetLineColor(kOrange); lin_g0_m02pc->SetLineColor(kOrange);
grap0_norm->GetXaxis()->SetRangeUser(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); norm_g0->GetXaxis()->SetRangeUser(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
grap0_norm->SetMinimum(1.5*lin_g0_m02pc->Eval(g0max)); norm_g0->SetMinimum(1.5*lin_g0_m02pc->Eval(rangemax_g0));
grap0_norm->SetMaximum(1.5*lin_g0_p02pc->Eval(g0max)); norm_g0->SetMaximum(1.5*lin_g0_p02pc->Eval(rangemax_g0));
grap0_norm->GetXaxis()->SetTitle("Signal generator voltage [mV]"); norm_g0->GetXaxis()->SetTitle("Signal generator voltage [mV]");
grap0_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]");
grap0_norm->GetYaxis()->SetTitleOffset(0.9); norm_g0->GetYaxis()->SetTitleOffset(0.9);
grap0_norm->Draw("AP"); norm_g0->Draw("AP");
flat_g0->Draw("same"); flat_g0->Draw("same");
lin_g0_p1pc->Draw("same"); lin_g0_p1pc->Draw("same");
lin_g0_p05pc->Draw("same"); lin_g0_p05pc->Draw("same");
@ -532,39 +530,39 @@ int main(int argc, char* argv[]) {
lin_g0_m1pc->Draw("same"); lin_g0_m1pc->Draw("same");
lin_g0_m05pc->Draw("same"); lin_g0_m05pc->Draw("same");
lin_g0_m02pc->Draw("same"); lin_g0_m02pc->Draw("same");
grap0_norm->Draw("P"); 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()); 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)); mapcanvas->SaveAs((const char *)(savename));
mapcanvas->SetLeftMargin(0.1); mapcanvas->SetLeftMargin(0.1);
mapcanvas->SetRightMargin(0.13); mapcanvas->SetRightMargin(0.13);
delete grap0_norm; delete norm_g0;
} }
} }
} }
if (r1_adc.size() > 1) { if (r1_adc.size() > 1) {
g1min = lowestPointAfterSwitching(r1_filter,r0_filter); rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter);
g1max = *max_element(r1_filter.begin(),r1_filter.end()); rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end());
if (g1max > g1min) { if (rangemax_g1 > rangemin_g1) {
grap1->Fit("pol1","QRC","",g1min,g1max); grap_g1->Fit("pol1","QRC","",rangemin_g1,rangemax_g1);
g1_fit = (TF1*) grap1->GetFunction("pol1"); fit_g1 = (TF1*) grap_g1->GetFunction("pol1");
g1_fit->SetLineColor(kGreen+2); fit_g1->SetLineColor(kGreen+2);
g1_fit->SetParName(0,"G1 const"); fit_g1->SetParName(0,"G1 const");
g1_fit->SetParName(1,"G1 grad"); fit_g1->SetParName(1,"G1 grad");
g1_grad_hist->Fill(g1_fit->GetParameter(1)); g1_grad_hist->Fill(fit_g1->GetParameter(1));
g1_grader_hist->Fill(g1_fit->GetParError(1)); g1_grader_hist->Fill(fit_g1->GetParError(1));
g1_grad_map->Fill(i%NC,i/NC,g1_fit->GetParameter(1)); g1_grad_map->Fill(i%NC,i/NC,fit_g1->GetParameter(1));
g1_grader_map->Fill(i%NC,i/NC,g1_fit->GetParError(1)); g1_grader_map->Fill(i%NC,i/NC,fit_g1->GetParError(1));
g1_const_hist->Fill(g1_fit->GetParameter(0)); g1_const_hist->Fill(fit_g1->GetParameter(0));
g1_conster_hist->Fill(g1_fit->GetParError(0)); g1_conster_hist->Fill(fit_g1->GetParError(0));
g1_const_map->Fill(i%NC,i/NC,g1_fit->GetParameter(0)); g1_const_map->Fill(i%NC,i/NC,fit_g1->GetParameter(0));
g1_conster_map->Fill(i%NC,i/NC,g1_fit->GetParError(0)); g1_conster_map->Fill(i%NC,i/NC,fit_g1->GetParError(0));
if ((i >= 58000 && i < 58000+10) || // bulk if ((i >= 58000 && i < 58000+10) || // bulk
(i >= 10 && i < 10+10) || // edge (i >= 10 && i < 10+10) || // edge
@ -592,16 +590,16 @@ int main(int argc, char* argv[]) {
mapcanvas->SetLeftMargin(0.13); mapcanvas->SetLeftMargin(0.13);
mapcanvas->SetRightMargin(0.05); mapcanvas->SetRightMargin(0.05);
grap1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); grap_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]");
grap1->GetYaxis()->SetTitle("ADC [ADU]"); grap_g1->GetYaxis()->SetTitle("ADC [ADU]");
grap1->GetYaxis()->SetTitleOffset(0.9); grap_g1->GetYaxis()->SetTitleOffset(0.9);
grap1->SetMinimum(1000); grap_g1->SetMinimum(1000);
grap1->SetMaximum(15000); grap_g1->SetMaximum(15000);
grap1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
grap1->Draw("AP"); grap_g1->Draw("AP");
g1_fit->Draw("same"); fit_g1->Draw("same");
mapcanvas->Update(); mapcanvas->Update();
TPaveStats *st0 = (TPaveStats*)grap1->FindObject("stats"); TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats");
st0->SetX1NDC(0.6); st0->SetX1NDC(0.6);
st0->SetX2NDC(0.94); st0->SetX2NDC(0.94);
st0->SetY1NDC(0.18); st0->SetY1NDC(0.18);
@ -613,91 +611,91 @@ int main(int argc, char* argv[]) {
vector<double> r1_adc_norm; vector<double> r1_adc_norm;
for (size_t j = 0; j < r1_adc.size(); j++) { for (size_t j = 0; j < r1_adc.size(); j++) {
r1_adc_norm.push_back(r1_adc[j] - g1_fit->Eval(r1_filter[j])); r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j]));
} }
TGraphErrors *grap1_norm = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0])); TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcerr[0]));
grap1_norm->SetMarkerColor(kGreen+2); norm_g1->SetMarkerColor(kGreen+2);
grap1_norm->SetLineColor(kGreen+2); norm_g1->SetLineColor(kGreen+2);
TF1* flat_g1 = new TF1("flat_gi","0",g1min,g1max); TF1* flat_g1 = new TF1("flat_gi","0",rangemin_g1,rangemax_g1);
flat_g1->SetLineColor(kGreen+2); flat_g1->SetLineColor(kGreen+2);
TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",g1min,g1max); TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p02pc->SetParameter(0,g1_fit->GetParameter(0)/500.); lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.);
lin_g1_p02pc->SetParameter(1,g1_fit->GetParameter(1)/500.); lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.);
lin_g1_p02pc->SetLineColor(kOrange); lin_g1_p02pc->SetLineColor(kOrange);
TF1* lin_g1_p01pc = new TF1("lin_g1_p01pc","[0]+[1]*x",g1min,g1max); TF1* lin_g1_p01pc = new TF1("lin_g1_p01pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p01pc->SetParameter(0,g1_fit->GetParameter(0)/1000.); lin_g1_p01pc->SetParameter(0,fit_g1->GetParameter(0)/1000.);
lin_g1_p01pc->SetParameter(1,g1_fit->GetParameter(1)/1000.); lin_g1_p01pc->SetParameter(1,fit_g1->GetParameter(1)/1000.);
lin_g1_p01pc->SetLineColor(kYellow); lin_g1_p01pc->SetLineColor(kYellow);
TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",g1min,g1max); TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m02pc->SetParameter(0,g1_fit->GetParameter(0)/-500.); lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.);
lin_g1_m02pc->SetParameter(1,g1_fit->GetParameter(1)/-500.); lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.);
lin_g1_m02pc->SetLineColor(kOrange); lin_g1_m02pc->SetLineColor(kOrange);
TF1* lin_g1_m01pc = new TF1("lin_g1_m01pc","[0]+[1]*x",g1min,g1max); TF1* lin_g1_m01pc = new TF1("lin_g1_m01pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m01pc->SetParameter(0,g1_fit->GetParameter(0)/-1000.); lin_g1_m01pc->SetParameter(0,fit_g1->GetParameter(0)/-1000.);
lin_g1_m01pc->SetParameter(1,g1_fit->GetParameter(1)/-1000.); lin_g1_m01pc->SetParameter(1,fit_g1->GetParameter(1)/-1000.);
lin_g1_m01pc->SetLineColor(kYellow); lin_g1_m01pc->SetLineColor(kYellow);
grap1_norm->GetXaxis()->SetRangeUser(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); norm_g1->GetXaxis()->SetRangeUser(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
grap1_norm->SetMinimum(1.5*lin_g1_m02pc->Eval(g1min)); norm_g1->SetMinimum(1.5*lin_g1_m02pc->Eval(rangemin_g1));
grap1_norm->SetMaximum(1.5*lin_g1_p02pc->Eval(g1min)); norm_g1->SetMaximum(1.5*lin_g1_p02pc->Eval(rangemin_g1));
grap1_norm->GetXaxis()->SetTitle("Signal generator voltage [mV]"); norm_g1->GetXaxis()->SetTitle("Signal generator voltage [mV]");
grap1_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]");
grap1_norm->GetYaxis()->SetTitleOffset(0.9); norm_g1->GetYaxis()->SetTitleOffset(0.9);
grap1_norm->Draw("AP"); norm_g1->Draw("AP");
flat_g1->Draw("same"); flat_g1->Draw("same");
lin_g1_p02pc->Draw("same"); lin_g1_p02pc->Draw("same");
lin_g1_p01pc->Draw("same"); lin_g1_p01pc->Draw("same");
lin_g1_m02pc->Draw("same"); lin_g1_m02pc->Draw("same");
lin_g1_m01pc->Draw("same"); lin_g1_m01pc->Draw("same");
grap1_norm->Draw("P"); 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()); 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)); mapcanvas->SaveAs((const char *)(savename));
mapcanvas->SetLeftMargin(0.1); mapcanvas->SetLeftMargin(0.1);
mapcanvas->SetRightMargin(0.13); mapcanvas->SetRightMargin(0.13);
delete grap1_norm; delete norm_g1;
} }
} }
} }
// get ratio measurements // get ratio measurements
if (g0max > g0min && g1max > g1min) { if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
g0overg1_hist->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1er_hist->Fill(sqrt(pow(g0_fit->GetParError(1)/g0_fit->GetParameter(1),2) + pow(g1_fit->GetParError(1)/g1_fit->GetParameter(1),2))); g0overg1er_hist->Fill(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,g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
g0overg1er_map->Fill(i%NC,i/NC,sqrt(pow(g0_fit->GetParError(1)/g0_fit->GetParameter(1),2) + pow(g1_fit->GetParError(1)/g1_fit->GetParameter(1),2))); g0overg1er_map->Fill(i%NC,i/NC,sqrt(pow(fit_g0->GetParError(1)/fit_g0->GetParameter(1),2) + pow(fit_g1->GetParError(1)/fit_g1->GetParameter(1),2)));
if (isEdge(i)) { if (isEdge(i)) {
g0overg1_hist_isEdge->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
if (isInnerEdge(i)) { if (isInnerEdge(i)) {
g0overg1_hist_isInnerEdge->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isInnerEdge->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
if (isDouble(i)) { if (isDouble(i)) {
g0overg1_hist_isDouble->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
if (isNextToDouble(i)) { if (isNextToDouble(i)) {
g0overg1_hist_isNextToDouble->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isNextToDouble->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
if (isQuad(i)) { if (isQuad(i)) {
g0overg1_hist_isQuad->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isQuad->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
if (isBulk(i)) { if (isBulk(i)) {
g0overg1_hist_isBulk->Fill(g0_fit->GetParameter(1) / g1_fit->GetParameter(1)); g0overg1_hist_isBulk->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1));
} }
} }
delete grap0; delete grap_g0;
delete grap1; delete grap_g1;
} }
} }

212
DB_analysis.cpp
View File

@ -379,16 +379,18 @@ int main(int argc, char* argv[]) {
TH2F* range1minmap = new TH2F("range1minmap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5); TH2F* range1minmap = new TH2F("range1minmap","",NC,-0.5,NC-0.5,NR,-0.5,NR-0.5);
for (int i = 0; i < NCH; i++) { for (int i = 0; i < NCH; i++) {
if (pixel_mask[i] == true) { if (pixel_mask[i] == true) {
vector<double> r0_adc; vector<double> r0_adc;
vector<double> r0_adcer; vector<double> r0_adcer;
vector<double> r0_filter; vector<double> r0_filter;
vector<double> r0_filterer; vector<double> r0_ferr;
vector<double> r1_adc; vector<double> r1_adc;
vector<double> r1_adcer; vector<double> r1_adcer;
vector<double> r1_filter; vector<double> r1_filter;
vector<double> r1_filterer; vector<double> r1_ferr;
for (int j = 0; j < nsteps; j++) { for (int j = 0; j < nsteps; j++) {
@ -406,13 +408,13 @@ int main(int argc, char* argv[]) {
if (this_g0avg != 0 && this_g0avger > 0) { if (this_g0avg != 0 && this_g0avger > 0) {
r0_filter.push_back(current); r0_filter.push_back(current);
r0_filterer.push_back(0); r0_ferr.push_back(0);
r0_adc.push_back(this_g0avg); r0_adc.push_back(this_g0avg);
r0_adcer.push_back(this_g0avger); r0_adcer.push_back(this_g0avger);
} }
if (this_g1avg != 0 && this_g1avger > 0) { if (this_g1avg != 0 && this_g1avger > 0) {
r1_filter.push_back(current); r1_filter.push_back(current);
r1_filterer.push_back(0); r1_ferr.push_back(0);
r1_adc.push_back(this_g1avg); r1_adc.push_back(this_g1avg);
r1_adcer.push_back(this_g1avger); r1_adcer.push_back(this_g1avger);
} }
@ -421,24 +423,24 @@ int main(int argc, char* argv[]) {
TGraphErrors *grap_g0 = 0; TGraphErrors *grap_g0 = 0;
TGraphErrors *grap_g1 = 0; TGraphErrors *grap_g1 = 0;
TF1 *fit1 = 0; TF1 *fit_g0 = 0;
TF1 *fit2 = 0; TF1 *fit_g1 = 0;
double rangemin0 = 0; double rangemin_g0 = 0;
double rangemax0 = 0; double rangemax_g0 = 0;
double rangemin1 = 0; double rangemin_g1 = 0;
double rangemax1 = 0; double rangemax_g1 = 0;
// define graphs // define graphs
if (r0_adc.size() > 1) { if (r0_adc.size() > 1) {
grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_filterer[0]),&(r0_adcer[0])); grap_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcer[0]));
grap_g0->SetMarkerStyle(20); grap_g0->SetMarkerStyle(20);
grap_g0->SetMarkerColor(kBlue); grap_g0->SetMarkerColor(kBlue);
grap_g0->SetLineColor(kBlue); grap_g0->SetLineColor(kBlue);
} }
if (r1_adc.size() > 1) { if (r1_adc.size() > 1) {
grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_filterer[0]),&(r1_adcer[0])); grap_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcer[0]));
grap_g1->SetMarkerStyle(20); grap_g1->SetMarkerStyle(20);
grap_g1->SetMarkerColor(kGreen+2); grap_g1->SetMarkerColor(kGreen+2);
grap_g1->SetLineColor(kGreen+2); grap_g1->SetLineColor(kGreen+2);
@ -494,23 +496,23 @@ int main(int argc, char* argv[]) {
// define fit ranges and fit // define fit ranges and fit
if (r0_adc.size() > 1) { if (r0_adc.size() > 1) {
rangemin0 = *min_element(r0_filter.begin(),r0_filter.end()); rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end());
rangemax0 = highestPointBeforeSwitching(r0_filter,r1_filter); rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter);
range0maxhist->Fill(rangemax0); range0maxhist->Fill(rangemax_g0);
range0maxmap->Fill(i%NC,i/NC,rangemax0); range0maxmap->Fill(i%NC,i/NC,rangemax_g0);
if (rangemax0 > rangemin0) { if (rangemax_g0 > rangemin_g0) {
fit1 = new TF1("fit1","[0]+[1]*x",rangemin0,rangemax0); fit_g0 = new TF1("fit_g0","[0]+[1]*x",rangemin_g0,rangemax_g0);
fit1->SetParameter(0, 100.); fit_g0->SetParameter(0, 100.);
fit1->SetParameter(1, 1); fit_g0->SetParameter(1, 1);
fit1->SetLineColor(kBlue); fit_g0->SetLineColor(kBlue);
fit1->SetParName(0,"G0 const"); fit_g0->SetParName(0,"G0 const");
fit1->SetParName(1,"G0 grad"); fit_g0->SetParName(1,"G0 grad");
grap_g0->Fit(fit1,"QR",""); grap_g0->Fit(fit_g0,"QR","");
g0map->Fill(i%NC,i/NC,fit1->GetParameter(1)); g0map->Fill(i%NC,i/NC,fit_g0->GetParameter(1));
g0ermap->Fill(i%NC,i/NC,fit1->GetParError(1)); g0ermap->Fill(i%NC,i/NC,fit_g0->GetParError(1));
if ((i >= 58000 && i < 58000+10) || // bulk if ((i >= 58000 && i < 58000+10) || // bulk
(i >= 10 && i < 10+10) || // edge (i >= 10 && i < 10+10) || // edge
@ -545,7 +547,7 @@ int main(int argc, char* argv[]) {
grap_g0->GetYaxis()->SetTitleOffset(0.9); grap_g0->GetYaxis()->SetTitleOffset(0.9);
grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); grap_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
grap_g0->Draw("AP"); grap_g0->Draw("AP");
fit1->Draw("same"); fit_g0->Draw("same");
mapcanvas->Update(); mapcanvas->Update();
TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats");
st0->SetX1NDC(0.2); st0->SetX1NDC(0.2);
@ -559,41 +561,41 @@ int main(int argc, char* argv[]) {
vector<double> r0_adc_norm; vector<double> r0_adc_norm;
for (size_t j = 0; j < r0_adc.size(); j++) { for (size_t j = 0; j < r0_adc.size(); j++) {
r0_adc_norm.push_back(r0_adc[j] - fit1->Eval(r0_filter[j])); r0_adc_norm.push_back(r0_adc[j] - fit_g0->Eval(r0_filter[j]));
} }
TGraphErrors *grap0_norm = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_filterer[0]),&(r0_adcer[0])); TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcer[0]));
grap0_norm->SetMarkerColor(kBlue); norm_g0->SetMarkerColor(kBlue);
grap0_norm->SetLineColor(kBlue); norm_g0->SetLineColor(kBlue);
TF1* flat_g0 = new TF1("flat_g0","0",rangemin0, rangemax0); TF1* flat_g0 = new TF1("flat_g0","0",rangemin_g0, rangemax_g0);
flat_g0->SetLineColor(kBlue); flat_g0->SetLineColor(kBlue);
TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin0,rangemax0); TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_p02pc->SetParameter(0,fit1->GetParameter(0)/500.); lin_g0_p02pc->SetParameter(0,fit_g0->GetParameter(0)/500.);
lin_g0_p02pc->SetParameter(1,fit1->GetParameter(1)/500.); lin_g0_p02pc->SetParameter(1,fit_g0->GetParameter(1)/500.);
lin_g0_p02pc->SetLineColor(kOrange); lin_g0_p02pc->SetLineColor(kOrange);
TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin0,rangemax0); TF1* lin_g0_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_p05pc->SetParameter(0,fit1->GetParameter(0)/200.); lin_g0_p05pc->SetParameter(0,fit_g0->GetParameter(0)/200.);
lin_g0_p05pc->SetParameter(1,fit1->GetParameter(1)/200.); lin_g0_p05pc->SetParameter(1,fit_g0->GetParameter(1)/200.);
lin_g0_p05pc->SetLineColor(kOrange+1); lin_g0_p05pc->SetLineColor(kOrange+1);
TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin0,rangemax0); TF1* lin_g0_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_m02pc->SetParameter(0,fit1->GetParameter(0)/-500.); lin_g0_m02pc->SetParameter(0,fit_g0->GetParameter(0)/-500.);
lin_g0_m02pc->SetParameter(1,fit1->GetParameter(1)/-500.); lin_g0_m02pc->SetParameter(1,fit_g0->GetParameter(1)/-500.);
lin_g0_m02pc->SetLineColor(kOrange); lin_g0_m02pc->SetLineColor(kOrange);
TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin0,rangemax0); TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin_g0,rangemax_g0);
lin_g0_m05pc->SetParameter(0,fit1->GetParameter(0)/-200.); lin_g0_m05pc->SetParameter(0,fit_g0->GetParameter(0)/-200.);
lin_g0_m05pc->SetParameter(1,fit1->GetParameter(1)/-200.); lin_g0_m05pc->SetParameter(1,fit_g0->GetParameter(1)/-200.);
lin_g0_m05pc->SetLineColor(kOrange+1); lin_g0_m05pc->SetLineColor(kOrange+1);
grap0_norm->GetXaxis()->SetTitle("Tube current [mA]"); norm_g0->GetXaxis()->SetTitle("Tube current [mA]");
grap0_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]");
grap0_norm->GetYaxis()->SetTitleOffset(0.9); norm_g0->GetYaxis()->SetTitleOffset(0.9);
grap0_norm->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); norm_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end()));
grap0_norm->Draw("AP"); norm_g0->Draw("AP");
flat_g0->Draw("same"); flat_g0->Draw("same");
lin_g0_p02pc->Draw("same"); lin_g0_p02pc->Draw("same");
lin_g0_p05pc->Draw("same"); lin_g0_p05pc->Draw("same");
@ -601,12 +603,12 @@ int main(int argc, char* argv[]) {
lin_g0_m05pc->Draw("same"); lin_g0_m05pc->Draw("same");
sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g0norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g0norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SaveAs((const char *)(savename));
grap0_norm->SetMinimum(-100); norm_g0->SetMinimum(-100);
grap0_norm->SetMaximum(100); norm_g0->SetMaximum(100);
sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g0norm_zoom_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g0norm_zoom_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SaveAs((const char *)(savename));
delete grap0_norm; delete norm_g0;
delete flat_g0; delete flat_g0;
mapcanvas->SetLeftMargin(0.1); mapcanvas->SetLeftMargin(0.1);
@ -617,23 +619,23 @@ int main(int argc, char* argv[]) {
if (r1_adc.size() > 1) { if (r1_adc.size() > 1) {
rangemin1 = lowestPointAfterSwitching(r1_filter,r0_filter); rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter);
range1minhist->Fill(rangemin1); range1minhist->Fill(rangemin_g1);
range1minmap->Fill(i%NC,i/NC,rangemin1); range1minmap->Fill(i%NC,i/NC,rangemin_g1);
rangemax1 = *max_element(r1_filter.begin(),r1_filter.end()); rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end());
if (rangemax1 > rangemin1) { if (rangemax_g1 > rangemin_g1) {
fit2 = new TF1("fit2","[0]+[1]*x",rangemin1,rangemax1); fit_g1 = new TF1("fit_g1","[0]+[1]*x",rangemin_g1,rangemax_g1);
fit2->SetParameter(0, 1000.); fit_g1->SetParameter(0, 1000.);
fit2->SetParameter(1, -0.1); fit_g1->SetParameter(1, -0.1);
fit2->SetLineColor(kGreen+2); fit_g1->SetLineColor(kGreen+2);
fit2->SetParName(0,"G1 const"); fit_g1->SetParName(0,"G1 const");
fit2->SetParName(1,"G1 grad"); fit_g1->SetParName(1,"G1 grad");
grap_g1->Fit(fit2,"QR",""); grap_g1->Fit(fit_g1,"QR","");
g1map->Fill(i%NC,i/NC,fit2->GetParameter(1)); g1map->Fill(i%NC,i/NC,fit_g1->GetParameter(1));
g1ermap->Fill(i%NC,i/NC,fit2->GetParError(1)); g1ermap->Fill(i%NC,i/NC,fit_g1->GetParError(1));
if ((i >= 58000 && i < 58000+10) || // bulk if ((i >= 58000 && i < 58000+10) || // bulk
(i >= 10 && i < 10+10) || // edge (i >= 10 && i < 10+10) || // edge
@ -668,7 +670,7 @@ int main(int argc, char* argv[]) {
grap_g1->GetYaxis()->SetTitleOffset(0.9); grap_g1->GetYaxis()->SetTitleOffset(0.9);
grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); grap_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
grap_g1->Draw("AP"); grap_g1->Draw("AP");
fit2->Draw("same"); fit_g1->Draw("same");
mapcanvas->Update(); mapcanvas->Update();
TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats");
st0->SetX1NDC(0.6); st0->SetX1NDC(0.6);
@ -682,41 +684,41 @@ int main(int argc, char* argv[]) {
vector<double> r1_adc_norm; vector<double> r1_adc_norm;
for (size_t j = 0; j < r1_adc.size(); j++) { 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]));
} }
TGraphErrors *grap1_norm = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_filterer[0]),&(r1_adcer[0])); TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcer[0]));
grap1_norm->SetMarkerColor(kGreen+2); norm_g1->SetMarkerColor(kGreen+2);
grap1_norm->SetLineColor(kGreen+2); norm_g1->SetLineColor(kGreen+2);
TF1* flat_g1 = new TF1("flat_g1","0",rangemin1, rangemax1); TF1* flat_g1 = new TF1("flat_g1","0",rangemin_g1, rangemax_g1);
flat_g1->SetLineColor(kGreen+2); flat_g1->SetLineColor(kGreen+2);
TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin1,rangemax1); TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p02pc->SetParameter(0,fit2->GetParameter(0)/500.); lin_g1_p02pc->SetParameter(0,fit_g1->GetParameter(0)/500.);
lin_g1_p02pc->SetParameter(1,fit2->GetParameter(1)/500.); lin_g1_p02pc->SetParameter(1,fit_g1->GetParameter(1)/500.);
lin_g1_p02pc->SetLineColor(kOrange); lin_g1_p02pc->SetLineColor(kOrange);
TF1* lin_g1_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin1,rangemax1); TF1* lin_g1_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_p05pc->SetParameter(0,fit2->GetParameter(0)/200.); lin_g1_p05pc->SetParameter(0,fit_g1->GetParameter(0)/200.);
lin_g1_p05pc->SetParameter(1,fit2->GetParameter(1)/200.); lin_g1_p05pc->SetParameter(1,fit_g1->GetParameter(1)/200.);
lin_g1_p05pc->SetLineColor(kOrange+1); lin_g1_p05pc->SetLineColor(kOrange+1);
TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin1,rangemax1); TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m02pc->SetParameter(0,fit2->GetParameter(0)/-500.); lin_g1_m02pc->SetParameter(0,fit_g1->GetParameter(0)/-500.);
lin_g1_m02pc->SetParameter(1,fit2->GetParameter(1)/-500.); lin_g1_m02pc->SetParameter(1,fit_g1->GetParameter(1)/-500.);
lin_g1_m02pc->SetLineColor(kOrange); lin_g1_m02pc->SetLineColor(kOrange);
TF1* lin_g1_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin1,rangemax1); TF1* lin_g1_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin_g1,rangemax_g1);
lin_g1_m05pc->SetParameter(0,fit2->GetParameter(0)/-200.); lin_g1_m05pc->SetParameter(0,fit_g1->GetParameter(0)/-200.);
lin_g1_m05pc->SetParameter(1,fit2->GetParameter(1)/-200.); lin_g1_m05pc->SetParameter(1,fit_g1->GetParameter(1)/-200.);
lin_g1_m05pc->SetLineColor(kOrange+1); lin_g1_m05pc->SetLineColor(kOrange+1);
grap1_norm->GetXaxis()->SetTitle("Tube current [mA]"); norm_g1->GetXaxis()->SetTitle("Tube current [mA]");
grap1_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]");
grap1_norm->GetYaxis()->SetTitleOffset(0.9); norm_g1->GetYaxis()->SetTitleOffset(0.9);
grap1_norm->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); norm_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end()));
grap1_norm->Draw("AP"); norm_g1->Draw("AP");
flat_g1->Draw("same"); flat_g1->Draw("same");
lin_g1_p02pc->Draw("same"); lin_g1_p02pc->Draw("same");
lin_g1_p05pc->Draw("same"); lin_g1_p05pc->Draw("same");
@ -724,50 +726,50 @@ int main(int argc, char* argv[]) {
lin_g1_m05pc->Draw("same"); lin_g1_m05pc->Draw("same");
sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g1norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g1norm_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SaveAs((const char *)(savename));
grap1_norm->SetMinimum(-50); norm_g1->SetMinimum(-50);
grap1_norm->SetMaximum(50); norm_g1->SetMaximum(50);
sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g1norm_zoom_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str()); sprintf(savename,"plots/M%s/DirectBeam/graps_%s_%d_g1norm_zoom_M%s.png", module_str.c_str(), pixel_type.c_str(), i, module_str.c_str());
mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SaveAs((const char *)(savename));
mapcanvas->SetLeftMargin(0.1); mapcanvas->SetLeftMargin(0.1);
mapcanvas->SetRightMargin(0.13); mapcanvas->SetRightMargin(0.13);
delete grap1_norm; delete norm_g1;
delete flat_g1; delete flat_g1;
} }
} }
} }
// get ratio measurements // get ratio measurements
if (rangemax0 > rangemin0 && rangemax1 > rangemin1) { if (rangemax_g0 > rangemin_g0 && rangemax_g1 > rangemin_g1) {
g0overg1map->Fill(i%NC,i/NC,fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
g0overg1hist->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
g0overg1ermap->Fill(i%NC,i/NC,abs(fit1->GetParameter(1)/fit2->GetParameter(1))*sqrt(pow((fit1->GetParError(1)/fit1->GetParameter(1)),2) + pow((fit2->GetParError(1)/fit2->GetParameter(1)),2))); 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)));
g0overg1erhist->Fill(abs(fit1->GetParameter(1)/fit2->GetParameter(1))*sqrt(pow((fit1->GetParError(1)/fit1->GetParameter(1)),2) + pow((fit2->GetParError(1)/fit2->GetParameter(1)),2))); 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)));
if (isEdge(i)) { if (isEdge(i)) {
g0overg1hist_isEdge->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isEdge->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
if (isInnerEdge(i)) { if (isInnerEdge(i)) {
g0overg1hist_isInnerEdge->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isInnerEdge->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
if (isDouble(i)) { if (isDouble(i)) {
g0overg1hist_isDouble->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isDouble->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
if (isNextToDouble(i)) { if (isNextToDouble(i)) {
g0overg1hist_isNextToDouble->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isNextToDouble->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
if (isQuad(i)) { if (isQuad(i)) {
g0overg1hist_isQuad->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isQuad->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
if (isBulk(i)) { if (isBulk(i)) {
g0overg1hist_isBulk->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); g0overg1hist_isBulk->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1));
} }
} }
delete fit1; delete fit_g0;
delete fit2; delete fit_g1;
delete grap_g0; delete grap_g0;
delete grap_g1; delete grap_g1;