diff --git a/BP_analysis.cpp b/BP_analysis.cpp index effc4b1..c4cc81f 100644 --- a/BP_analysis.cpp +++ b/BP_analysis.cpp @@ -299,8 +299,6 @@ int main(int argc, char* argv[]) { for (int i = 0; i < NCH; i++) { - if (i%100000==0) {cout << "another 100k" << endl;} - if (pixel_mask[i] == true) { vector r0_adc; @@ -336,30 +334,30 @@ int main(int argc, char* argv[]) { } } - TGraphErrors *grap0 = 0; - TGraphErrors *grap1 = 0; + TGraphErrors *grap_g0 = 0; + TGraphErrors *grap_g1 = 0; - TF1 *g0_fit = 0; - TF1 *g1_fit = 0; + TF1 *fit_g0 = 0; + TF1 *fit_g1 = 0; - double g0min = 0; - double g0max = 0; - double g1min = 0; - double g1max = 0; + double rangemin_g0 = 0; + double rangemax_g0 = 0; + double rangemin_g1 = 0; + double rangemax_g1 = 0; // define graphs if (r0_adc.size() > 1) { - grap0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc[0]),&(r0_ferr[0]),&(r0_adcerr[0])); - grap0->SetMarkerStyle(20); - grap0->SetMarkerColor(kBlue); - grap0->SetLineColor(kBlue); + 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) { - grap1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc[0]),&(r1_ferr[0]),&(r1_adcerr[0])); - grap1->SetMarkerStyle(20); - grap1->SetMarkerColor(kGreen+2); - grap1->SetLineColor(kGreen+2); + 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 @@ -391,15 +389,15 @@ int main(int argc, char* argv[]) { mapcanvas->SetLeftMargin(0.13); mapcanvas->SetRightMargin(0.05); - grap1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap1->GetYaxis()->SetTitle("ADC [ADU]"); - grap1->GetYaxis()->SetTitleOffset(0.9); - grap1->SetMinimum(1000); - grap1->SetMaximum(15000); - grap1->GetXaxis()->SetLimits(0,7200); - grap1->Draw("AP"); + 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) { - grap0->Draw("P"); + 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()); @@ -411,25 +409,25 @@ int main(int argc, char* argv[]) { // define fit ranges and fit if (r0_adc.size() > 1) { - g0min = *min_element(r0_filter.begin(),r0_filter.end()); - g0max = highestPointBeforeSwitching(r0_filter,r1_filter); + rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end()); + rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter); - if (g0max > g0min) { + if (rangemax_g0 > rangemin_g0) { - grap0->Fit("pol1","QRC","",g0min,g0max); - g0_fit = (TF1*) grap0->GetFunction("pol1"); - g0_fit->SetLineColor(kBlue); - g0_fit->SetParName(0,"G0 const"); - g0_fit->SetParName(1,"G0 grad"); + 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(g0_fit->GetParameter(1)); - g0_grader_hist->Fill(g0_fit->GetParError(1)); - g0_grad_map->Fill(i%NC,i/NC,g0_fit->GetParameter(1)); - g0_grader_map->Fill(i%NC,i/NC,g0_fit->GetParError(1)); - g0_const_hist->Fill(g0_fit->GetParameter(0)); - g0_conster_hist->Fill(g0_fit->GetParError(0)); - g0_const_map->Fill(i%NC,i/NC,g0_fit->GetParameter(0)); - g0_conster_map->Fill(i%NC,i/NC,g0_fit->GetParError(0)); + 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 @@ -457,16 +455,16 @@ int main(int argc, char* argv[]) { mapcanvas->SetLeftMargin(0.13); mapcanvas->SetRightMargin(0.05); - grap0->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap0->GetYaxis()->SetTitle("ADC [ADU]"); - grap0->GetYaxis()->SetTitleOffset(0.9); - grap0->SetMinimum(1000); - grap0->SetMaximum(15000); - grap0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); - grap0->Draw("AP"); - g0_fit->Draw("same"); + 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*)grap0->FindObject("stats"); + TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); st0->SetX1NDC(0.2); st0->SetX2NDC(0.54); st0->SetY1NDC(0.18); @@ -478,53 +476,53 @@ int main(int argc, char* argv[]) { vector r0_adc_norm; 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])); - grap0_norm->SetMarkerColor(kBlue); - grap0_norm->SetLineColor(kBlue); + 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",g0min,g0max); + 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",g0min,g0max); - lin_g0_p1pc->SetParameter(0,g0_fit->GetParameter(0)/100.); - lin_g0_p1pc->SetParameter(1,g0_fit->GetParameter(1)/100.); + 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",g0min,g0max); - lin_g0_p05pc->SetParameter(0,g0_fit->GetParameter(0)/200.); - lin_g0_p05pc->SetParameter(1,g0_fit->GetParameter(1)/200.); + 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",g0min,g0max); - lin_g0_p02pc->SetParameter(0,g0_fit->GetParameter(0)/500.); - lin_g0_p02pc->SetParameter(1,g0_fit->GetParameter(1)/500.); + 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",g0min,g0max); - lin_g0_m1pc->SetParameter(0,g0_fit->GetParameter(0)/-100.); - lin_g0_m1pc->SetParameter(1,g0_fit->GetParameter(1)/-100.); + 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",g0min,g0max); - lin_g0_m05pc->SetParameter(0,g0_fit->GetParameter(0)/-200.); - lin_g0_m05pc->SetParameter(1,g0_fit->GetParameter(1)/-200.); + 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",g0min,g0max); - lin_g0_m02pc->SetParameter(0,g0_fit->GetParameter(0)/-500.); - lin_g0_m02pc->SetParameter(1,g0_fit->GetParameter(1)/-500.); + 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); - grap0_norm->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)); - grap0_norm->SetMaximum(1.5*lin_g0_p02pc->Eval(g0max)); - grap0_norm->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap0_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap0_norm->GetYaxis()->SetTitleOffset(0.9); - grap0_norm->Draw("AP"); + 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"); @@ -532,39 +530,39 @@ int main(int argc, char* argv[]) { lin_g0_m1pc->Draw("same"); lin_g0_m05pc->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()); mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SetLeftMargin(0.1); mapcanvas->SetRightMargin(0.13); - delete grap0_norm; + delete norm_g0; } } } if (r1_adc.size() > 1) { - g1min = lowestPointAfterSwitching(r1_filter,r0_filter); - g1max = *max_element(r1_filter.begin(),r1_filter.end()); + rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter); + rangemax_g1 = *max_element(r1_filter.begin(),r1_filter.end()); - if (g1max > g1min) { + if (rangemax_g1 > rangemin_g1) { - grap1->Fit("pol1","QRC","",g1min,g1max); - g1_fit = (TF1*) grap1->GetFunction("pol1"); - g1_fit->SetLineColor(kGreen+2); - g1_fit->SetParName(0,"G1 const"); - g1_fit->SetParName(1,"G1 grad"); + 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(g1_fit->GetParameter(1)); - g1_grader_hist->Fill(g1_fit->GetParError(1)); - g1_grad_map->Fill(i%NC,i/NC,g1_fit->GetParameter(1)); - g1_grader_map->Fill(i%NC,i/NC,g1_fit->GetParError(1)); - g1_const_hist->Fill(g1_fit->GetParameter(0)); - g1_conster_hist->Fill(g1_fit->GetParError(0)); - g1_const_map->Fill(i%NC,i/NC,g1_fit->GetParameter(0)); - g1_conster_map->Fill(i%NC,i/NC,g1_fit->GetParError(0)); + 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 @@ -592,16 +590,16 @@ int main(int argc, char* argv[]) { mapcanvas->SetLeftMargin(0.13); mapcanvas->SetRightMargin(0.05); - grap1->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap1->GetYaxis()->SetTitle("ADC [ADU]"); - grap1->GetYaxis()->SetTitleOffset(0.9); - grap1->SetMinimum(1000); - grap1->SetMaximum(15000); - grap1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); - grap1->Draw("AP"); - g1_fit->Draw("same"); + 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*)grap1->FindObject("stats"); + TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); st0->SetX1NDC(0.6); st0->SetX2NDC(0.94); st0->SetY1NDC(0.18); @@ -613,91 +611,91 @@ int main(int argc, char* argv[]) { vector r1_adc_norm; 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])); - grap1_norm->SetMarkerColor(kGreen+2); - grap1_norm->SetLineColor(kGreen+2); + 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",g1min,g1max); + 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",g1min,g1max); - lin_g1_p02pc->SetParameter(0,g1_fit->GetParameter(0)/500.); - lin_g1_p02pc->SetParameter(1,g1_fit->GetParameter(1)/500.); + 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",g1min,g1max); - lin_g1_p01pc->SetParameter(0,g1_fit->GetParameter(0)/1000.); - lin_g1_p01pc->SetParameter(1,g1_fit->GetParameter(1)/1000.); + 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",g1min,g1max); - lin_g1_m02pc->SetParameter(0,g1_fit->GetParameter(0)/-500.); - lin_g1_m02pc->SetParameter(1,g1_fit->GetParameter(1)/-500.); + 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",g1min,g1max); - lin_g1_m01pc->SetParameter(0,g1_fit->GetParameter(0)/-1000.); - lin_g1_m01pc->SetParameter(1,g1_fit->GetParameter(1)/-1000.); + 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); - grap1_norm->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)); - grap1_norm->SetMaximum(1.5*lin_g1_p02pc->Eval(g1min)); - grap1_norm->GetXaxis()->SetTitle("Signal generator voltage [mV]"); - grap1_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap1_norm->GetYaxis()->SetTitleOffset(0.9); - grap1_norm->Draw("AP"); + 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"); - 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()); mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SetLeftMargin(0.1); mapcanvas->SetRightMargin(0.13); - delete grap1_norm; + delete norm_g1; } } } // 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)); - g0overg1er_hist->Fill(sqrt(pow(g0_fit->GetParError(1)/g0_fit->GetParameter(1),2) + pow(g1_fit->GetParError(1)/g1_fit->GetParameter(1),2))); + g0overg1_hist->Fill(fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + 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)); - 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))); + g0overg1_map->Fill(i%NC,i/NC,fit_g0->GetParameter(1) / fit_g1->GetParameter(1)); + 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)) { - 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)) { - 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)) { - 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)) { - 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)) { - 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)) { - 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 grap1; + delete grap_g0; + delete grap_g1; } } diff --git a/DB_analysis.cpp b/DB_analysis.cpp index 0799ecb..3d50ad5 100644 --- a/DB_analysis.cpp +++ b/DB_analysis.cpp @@ -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); for (int i = 0; i < NCH; i++) { + if (pixel_mask[i] == true) { vector r0_adc; vector r0_adcer; vector r0_filter; - vector r0_filterer; + vector r0_ferr; + vector r1_adc; vector r1_adcer; vector r1_filter; - vector r1_filterer; + vector r1_ferr; for (int j = 0; j < nsteps; j++) { @@ -406,13 +408,13 @@ int main(int argc, char* argv[]) { if (this_g0avg != 0 && this_g0avger > 0) { r0_filter.push_back(current); - r0_filterer.push_back(0); + r0_ferr.push_back(0); r0_adc.push_back(this_g0avg); r0_adcer.push_back(this_g0avger); } if (this_g1avg != 0 && this_g1avger > 0) { r1_filter.push_back(current); - r1_filterer.push_back(0); + r1_ferr.push_back(0); r1_adc.push_back(this_g1avg); r1_adcer.push_back(this_g1avger); } @@ -421,24 +423,24 @@ int main(int argc, char* argv[]) { TGraphErrors *grap_g0 = 0; TGraphErrors *grap_g1 = 0; - TF1 *fit1 = 0; - TF1 *fit2 = 0; + TF1 *fit_g0 = 0; + TF1 *fit_g1 = 0; - double rangemin0 = 0; - double rangemax0 = 0; - double rangemin1 = 0; - double rangemax1 = 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_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->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_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->SetMarkerColor(kGreen+2); grap_g1->SetLineColor(kGreen+2); @@ -494,23 +496,23 @@ int main(int argc, char* argv[]) { // define fit ranges and fit if (r0_adc.size() > 1) { - rangemin0 = *min_element(r0_filter.begin(),r0_filter.end()); - rangemax0 = highestPointBeforeSwitching(r0_filter,r1_filter); - range0maxhist->Fill(rangemax0); - range0maxmap->Fill(i%NC,i/NC,rangemax0); + rangemin_g0 = *min_element(r0_filter.begin(),r0_filter.end()); + rangemax_g0 = highestPointBeforeSwitching(r0_filter,r1_filter); + range0maxhist->Fill(rangemax_g0); + 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); - fit1->SetParameter(0, 100.); - fit1->SetParameter(1, 1); - fit1->SetLineColor(kBlue); - fit1->SetParName(0,"G0 const"); - fit1->SetParName(1,"G0 grad"); - grap_g0->Fit(fit1,"QR",""); + fit_g0 = new TF1("fit_g0","[0]+[1]*x",rangemin_g0,rangemax_g0); + fit_g0->SetParameter(0, 100.); + fit_g0->SetParameter(1, 1); + fit_g0->SetLineColor(kBlue); + fit_g0->SetParName(0,"G0 const"); + fit_g0->SetParName(1,"G0 grad"); + grap_g0->Fit(fit_g0,"QR",""); - g0map->Fill(i%NC,i/NC,fit1->GetParameter(1)); - g0ermap->Fill(i%NC,i/NC,fit1->GetParError(1)); + g0map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)); + g0ermap->Fill(i%NC,i/NC,fit_g0->GetParError(1)); if ((i >= 58000 && i < 58000+10) || // bulk (i >= 10 && i < 10+10) || // edge @@ -545,7 +547,7 @@ int main(int argc, char* argv[]) { 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->Draw("AP"); - fit1->Draw("same"); + fit_g0->Draw("same"); mapcanvas->Update(); TPaveStats *st0 = (TPaveStats*)grap_g0->FindObject("stats"); st0->SetX1NDC(0.2); @@ -559,41 +561,41 @@ int main(int argc, char* argv[]) { vector r0_adc_norm; 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])); - grap0_norm->SetMarkerColor(kBlue); - grap0_norm->SetLineColor(kBlue); + TGraphErrors *norm_g0 = new TGraphErrors(r0_adc.size(),&(r0_filter[0]),&(r0_adc_norm[0]),&(r0_ferr[0]),&(r0_adcer[0])); + norm_g0->SetMarkerColor(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); - TF1* lin_g0_p02pc = new TF1("lin_g0_p02pc","[0]+[1]*x",rangemin0,rangemax0); - lin_g0_p02pc->SetParameter(0,fit1->GetParameter(0)/500.); - lin_g0_p02pc->SetParameter(1,fit1->GetParameter(1)/500.); + 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_p05pc = new TF1("lin_g0_p05pc","[0]+[1]*x",rangemin0,rangemax0); - lin_g0_p05pc->SetParameter(0,fit1->GetParameter(0)/200.); - lin_g0_p05pc->SetParameter(1,fit1->GetParameter(1)/200.); + 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_m02pc = new TF1("lin_g0_m02pc","[0]+[1]*x",rangemin0,rangemax0); - lin_g0_m02pc->SetParameter(0,fit1->GetParameter(0)/-500.); - lin_g0_m02pc->SetParameter(1,fit1->GetParameter(1)/-500.); + 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); - TF1* lin_g0_m05pc = new TF1("lin_g0_m05pc","[0]+[1]*x",rangemin0,rangemax0); - lin_g0_m05pc->SetParameter(0,fit1->GetParameter(0)/-200.); - lin_g0_m05pc->SetParameter(1,fit1->GetParameter(1)/-200.); + 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); - grap0_norm->GetXaxis()->SetTitle("Tube current [mA]"); - grap0_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap0_norm->GetYaxis()->SetTitleOffset(0.9); - grap0_norm->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); - grap0_norm->Draw("AP"); + norm_g0->GetXaxis()->SetTitle("Tube current [mA]"); + norm_g0->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g0->GetYaxis()->SetTitleOffset(0.9); + norm_g0->GetXaxis()->SetLimits(*min_element(r0_filter.begin(),r0_filter.end()),*max_element(r0_filter.begin(),r0_filter.end())); + norm_g0->Draw("AP"); flat_g0->Draw("same"); lin_g0_p02pc->Draw("same"); lin_g0_p05pc->Draw("same"); @@ -601,12 +603,12 @@ int main(int argc, char* argv[]) { 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()); mapcanvas->SaveAs((const char *)(savename)); - grap0_norm->SetMinimum(-100); - grap0_norm->SetMaximum(100); + norm_g0->SetMinimum(-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()); mapcanvas->SaveAs((const char *)(savename)); - delete grap0_norm; + delete norm_g0; delete flat_g0; mapcanvas->SetLeftMargin(0.1); @@ -617,23 +619,23 @@ int main(int argc, char* argv[]) { if (r1_adc.size() > 1) { - rangemin1 = lowestPointAfterSwitching(r1_filter,r0_filter); - range1minhist->Fill(rangemin1); - range1minmap->Fill(i%NC,i/NC,rangemin1); - rangemax1 = *max_element(r1_filter.begin(),r1_filter.end()); + rangemin_g1 = lowestPointAfterSwitching(r1_filter,r0_filter); + range1minhist->Fill(rangemin_g1); + range1minmap->Fill(i%NC,i/NC,rangemin_g1); + 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); - fit2->SetParameter(0, 1000.); - fit2->SetParameter(1, -0.1); - fit2->SetLineColor(kGreen+2); - fit2->SetParName(0,"G1 const"); - fit2->SetParName(1,"G1 grad"); - grap_g1->Fit(fit2,"QR",""); + fit_g1 = new TF1("fit_g1","[0]+[1]*x",rangemin_g1,rangemax_g1); + fit_g1->SetParameter(0, 1000.); + fit_g1->SetParameter(1, -0.1); + fit_g1->SetLineColor(kGreen+2); + fit_g1->SetParName(0,"G1 const"); + fit_g1->SetParName(1,"G1 grad"); + grap_g1->Fit(fit_g1,"QR",""); - g1map->Fill(i%NC,i/NC,fit2->GetParameter(1)); - g1ermap->Fill(i%NC,i/NC,fit2->GetParError(1)); + g1map->Fill(i%NC,i/NC,fit_g1->GetParameter(1)); + g1ermap->Fill(i%NC,i/NC,fit_g1->GetParError(1)); if ((i >= 58000 && i < 58000+10) || // bulk (i >= 10 && i < 10+10) || // edge @@ -668,7 +670,7 @@ int main(int argc, char* argv[]) { 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->Draw("AP"); - fit2->Draw("same"); + fit_g1->Draw("same"); mapcanvas->Update(); TPaveStats *st0 = (TPaveStats*)grap_g1->FindObject("stats"); st0->SetX1NDC(0.6); @@ -682,41 +684,41 @@ int main(int argc, char* argv[]) { vector r1_adc_norm; for (size_t j = 0; j < r1_adc.size(); j++) { - r1_adc_norm.push_back(r1_adc[j] - fit2->Eval(r1_filter[j])); + r1_adc_norm.push_back(r1_adc[j] - fit_g1->Eval(r1_filter[j])); } - TGraphErrors *grap1_norm = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_filterer[0]),&(r1_adcer[0])); - grap1_norm->SetMarkerColor(kGreen+2); - grap1_norm->SetLineColor(kGreen+2); + TGraphErrors *norm_g1 = new TGraphErrors(r1_adc.size(),&(r1_filter[0]),&(r1_adc_norm[0]),&(r1_ferr[0]),&(r1_adcer[0])); + norm_g1->SetMarkerColor(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); - TF1* lin_g1_p02pc = new TF1("lin_g1_p02pc","[0]+[1]*x",rangemin1,rangemax1); - lin_g1_p02pc->SetParameter(0,fit2->GetParameter(0)/500.); - lin_g1_p02pc->SetParameter(1,fit2->GetParameter(1)/500.); + 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_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin1,rangemax1); - lin_g1_p05pc->SetParameter(0,fit2->GetParameter(0)/200.); - lin_g1_p05pc->SetParameter(1,fit2->GetParameter(1)/200.); + TF1* lin_g1_p05pc = new TF1("lin_g1_p05pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_p05pc->SetParameter(0,fit_g1->GetParameter(0)/200.); + lin_g1_p05pc->SetParameter(1,fit_g1->GetParameter(1)/200.); lin_g1_p05pc->SetLineColor(kOrange+1); - TF1* lin_g1_m02pc = new TF1("lin_g1_m02pc","[0]+[1]*x",rangemin1,rangemax1); - lin_g1_m02pc->SetParameter(0,fit2->GetParameter(0)/-500.); - lin_g1_m02pc->SetParameter(1,fit2->GetParameter(1)/-500.); + 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_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin1,rangemax1); - lin_g1_m05pc->SetParameter(0,fit2->GetParameter(0)/-200.); - lin_g1_m05pc->SetParameter(1,fit2->GetParameter(1)/-200.); + TF1* lin_g1_m05pc = new TF1("lin_g1_m05pc","[0]+[1]*x",rangemin_g1,rangemax_g1); + lin_g1_m05pc->SetParameter(0,fit_g1->GetParameter(0)/-200.); + lin_g1_m05pc->SetParameter(1,fit_g1->GetParameter(1)/-200.); lin_g1_m05pc->SetLineColor(kOrange+1); - grap1_norm->GetXaxis()->SetTitle("Tube current [mA]"); - grap1_norm->GetYaxis()->SetTitle("Normalised ADC [ADU]"); - grap1_norm->GetYaxis()->SetTitleOffset(0.9); - grap1_norm->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); - grap1_norm->Draw("AP"); + norm_g1->GetXaxis()->SetTitle("Tube current [mA]"); + norm_g1->GetYaxis()->SetTitle("Normalised ADC [ADU]"); + norm_g1->GetYaxis()->SetTitleOffset(0.9); + norm_g1->GetXaxis()->SetLimits(*min_element(r1_filter.begin(),r1_filter.end()),*max_element(r1_filter.begin(),r1_filter.end())); + norm_g1->Draw("AP"); flat_g1->Draw("same"); lin_g1_p02pc->Draw("same"); lin_g1_p05pc->Draw("same"); @@ -724,50 +726,50 @@ int main(int argc, char* argv[]) { 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()); mapcanvas->SaveAs((const char *)(savename)); - grap1_norm->SetMinimum(-50); - grap1_norm->SetMaximum(50); + norm_g1->SetMinimum(-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()); mapcanvas->SaveAs((const char *)(savename)); mapcanvas->SetLeftMargin(0.1); mapcanvas->SetRightMargin(0.13); - delete grap1_norm; + delete norm_g1; delete flat_g1; } } } // 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)); - g0overg1hist->Fill(fit1->GetParameter(1)/fit2->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))); - 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))); + g0overg1map->Fill(i%NC,i/NC,fit_g0->GetParameter(1)/fit_g1->GetParameter(1)); + g0overg1hist->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1)); + g0overg1ermap->Fill(i%NC,i/NC,abs(fit_g0->GetParameter(1)/fit_g1->GetParameter(1))*sqrt(pow((fit_g0->GetParError(1)/fit_g0->GetParameter(1)),2) + pow((fit_g1->GetParError(1)/fit_g1->GetParameter(1)),2))); + 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)) { - g0overg1hist_isEdge->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); + g0overg1hist_isEdge->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1)); } 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)) { - g0overg1hist_isDouble->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); + g0overg1hist_isDouble->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1)); } 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)) { - g0overg1hist_isQuad->Fill(fit1->GetParameter(1)/fit2->GetParameter(1)); + g0overg1hist_isQuad->Fill(fit_g0->GetParameter(1)/fit_g1->GetParameter(1)); } 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 fit2; + delete fit_g0; + delete fit_g1; delete grap_g0; delete grap_g1;