From 9d4949ad14eb45efa9a36e800e701141f90bd6e7 Mon Sep 17 00:00:00 2001 From: redford_s Date: Fri, 25 Jan 2019 14:13:43 +0100 Subject: [PATCH] Added linearity limits. --- DB_analysis.cpp | 48 ++++++++++++++++++++++++++++++++++++++++++++++++ 1 file changed, 48 insertions(+) diff --git a/DB_analysis.cpp b/DB_analysis.cpp index 78dd5e8..5aede56 100644 --- a/DB_analysis.cpp +++ b/DB_analysis.cpp @@ -554,6 +554,46 @@ int main(int argc, char* argv[]) { TF1* flat_g0 = new TF1("flat_g0","0",rangemin0, rangemax0); TF1* flat_g1 = new TF1("flat_g1","0",rangemin1, rangemax1); + 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.); + 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.); + 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.); + 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.); + lin_g0_m05pc->SetLineColor(kOrange+1); + + 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.); + 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.); + 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.); + 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.); + lin_g1_m05pc->SetLineColor(kOrange+1); + grap0_norm->SetMarkerColor(kBlue); grap0_norm->SetLineColor(kBlue); flat_g0->SetLineColor(kBlue); @@ -562,6 +602,10 @@ int main(int argc, char* argv[]) { grap0_norm->GetYaxis()->SetTitleOffset(0.9); grap0_norm->Draw("AP"); flat_g0->Draw("same"); + lin_g0_p02pc->Draw("same"); + lin_g0_p05pc->Draw("same"); + lin_g0_m02pc->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()); mapcanvas->SaveAs((const char *)(savename)); grap0_norm->SetMinimum(-100); @@ -577,6 +621,10 @@ int main(int argc, char* argv[]) { grap1_norm->GetYaxis()->SetTitleOffset(0.9); grap1_norm->Draw("AP"); flat_g1->Draw("same"); + lin_g1_p02pc->Draw("same"); + lin_g1_p05pc->Draw("same"); + lin_g1_m02pc->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()); mapcanvas->SaveAs((const char *)(savename)); grap1_norm->SetMinimum(-50);