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

Added linearity limits.

Browse Source
This commit is contained in:
redford_s
2019-01-25 14:13:43 +01:00
parent eab04d984b
commit 9d4949ad14
1 changed files with 48 additions and 0 deletions
Download Patch File Download Diff File Expand all files Collapse all files

48
DB_analysis.cpp
View File

@ -554,6 +554,46 @@ int main(int argc, char* argv[]) {
TF1* flat_g0 = new TF1("flat_g0","0",rangemin0, rangemax0); TF1* flat_g0 = new TF1("flat_g0","0",rangemin0, rangemax0);
TF1* flat_g1 = new TF1("flat_g1","0",rangemin1, rangemax1); 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->SetMarkerColor(kBlue);
grap0_norm->SetLineColor(kBlue); grap0_norm->SetLineColor(kBlue);
flat_g0->SetLineColor(kBlue); flat_g0->SetLineColor(kBlue);
@ -562,6 +602,10 @@ int main(int argc, char* argv[]) {
grap0_norm->GetYaxis()->SetTitleOffset(0.9); grap0_norm->GetYaxis()->SetTitleOffset(0.9);
grap0_norm->Draw("AP"); grap0_norm->Draw("AP");
flat_g0->Draw("same"); 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()); 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); grap0_norm->SetMinimum(-100);
@ -577,6 +621,10 @@ int main(int argc, char* argv[]) {
grap1_norm->GetYaxis()->SetTitleOffset(0.9); grap1_norm->GetYaxis()->SetTitleOffset(0.9);
grap1_norm->Draw("AP"); grap1_norm->Draw("AP");
flat_g1->Draw("same"); 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()); 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); grap1_norm->SetMinimum(-50);