data analysis scripts
You can not select more than 25 topics Topics must start with a letter or number, can include dashes ('-') and can be up to 35 characters long.
 
 
 
 

303 lines
14 KiB

{
gROOT->ProcessLine(".x lhcbStyle.C");
TCanvas c1_mg1_e1;
TGraphErrors * ge1_p1 = new TGraphErrors("inlistE255_p_bpm1_diff.txt","%lg %lg %lg");
TGraphErrors * ge1_p2 = new TGraphErrors("inlistE21_p_bpm1_diff.txt","%lg %lg %lg");
TGraphErrors * ge1_ic1 = new TGraphErrors("inlistE255_p_ic1_diff.txt","%lg %lg %lg");
TGraphErrors * ge1_ic2 = new TGraphErrors("inlistE21_p_ic1_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e1 = new TMultiGraph();
mg_e1->Add(ge1_p1,"pc"); ge1_p1->SetMarkerStyle(20);ge1_p1->SetLineStyle(3);
mg_e1->Add(ge1_p2,"pc"); ge1_p2->SetMarkerStyle(24);ge1_p2->SetLineStyle(3);
mg_e1->Add(ge1_ic1,"p"); ge1_ic1->SetMarkerStyle(29);
mg_e1->Add(ge1_ic2,"p"); ge1_ic2->SetMarkerStyle(30);
mg_e1->Draw("a");
mg_e1->SetTitle(" ");
mg_e1->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e1->GetYaxis()->SetTitle("#sqrt{2}#upoint #sigma(A_{1}-A_{2}) / (A_{1}+A_{2})");
// mg_e1->GetXaxis()->SetRangeUser(0,600);
//mg_e1->GetYaxis()->SetRangeUser(0,1);
mg_e1->SetMinimum(0.);
mg_e1->SetMaximum(0.18);
TLegend *mylegend = new TLegend(0.60,0.6,0.9,0.9);
mylegend->SetFillColor(0); // white background
mylegend->SetTextFont(22);
mylegend->SetBorderSize(0); // get rid of the box
mylegend->SetTextSize(0.045); // set text size
mylegend->AddEntry(ge1_p1,"BPM (221 MeV)","pl"); // options: p,l,f
mylegend->AddEntry(ge1_p2,"BPM (70 MeV)","pl"); // options: p,l,f
mylegend->AddEntry(ge1_ic1,"IC (221 MeV)","p"); // options: p,l,f
mylegend->AddEntry(ge1_ic2,"IC (70 MeV)","p"); // options: p,l,f
mylegend->Draw();
gPad->Modified();
c1_mg1_e1.SaveAs("figs/inlist_p_bpm_diff.pdf");
c1_mg1_e1.SaveAs("figs/inlist_p_bpm_diff.png");
c1_mg1_e1.SaveAs("figs/inlist_p_bpm_diff.C");
TCanvas c1_mg1_e2;
TGraphErrors * ge2_p1 = new TGraphErrors("inlistE255_c_bpm1_diff.txt","%lg %lg %lg");
TGraphErrors * ge2_p2 = new TGraphErrors("inlistE21_c_bpm1_diff.txt","%lg %lg %lg");
TGraphErrors * ge2_ic1 = new TGraphErrors("inlistE255_c_ic1_diff.txt","%lg %lg %lg");
TGraphErrors * ge2_ic2 = new TGraphErrors("inlistE21_c_ic1_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e2 = new TMultiGraph();
mg_e2->Add(ge2_p1,"pc"); ge2_p1->SetMarkerStyle(22);ge2_p1->SetLineStyle(3);
mg_e2->Add(ge2_p2,"pc"); ge2_p2->SetMarkerStyle(26);ge2_p2->SetLineStyle(3);
mg_e2->Add(ge2_ic1,"p"); ge2_ic1->SetMarkerStyle(29);
mg_e2->Add(ge2_ic2,"p"); ge2_ic2->SetMarkerStyle(30);
mg_e2->Draw("a");
mg_e2->SetTitle(" ");
mg_e2->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e2->GetYaxis()->SetTitle("#sqrt{2}#upoint #sigma(A_{1}-A_{2}) / (A_{1}+A_{2})");
// mg_e2->GetXaxis()->SetRangeUser(0,600);
//mg_e2->GetYaxis()->SetRangeUser(0,1);
mg_e2->SetMinimum(0.);
mg_e2->SetMaximum(0.1);
TLegend * mylegend2 = new TLegend(0.60,0.6,0.9,0.9);
mylegend2->SetFillColor(0); // white background
mylegend2->SetTextFont(22);
mylegend2->SetBorderSize(0); // get rid of the box
mylegend2->SetTextSize(0.045); // set text size
mylegend2->AddEntry(ge2_p1,"BPM (430 MeV/u)","pl"); // options: p,l,f
mylegend2->AddEntry(ge2_p2,"BPM (130 MeV/u)","pl"); // options: p,l,f
mylegend2->AddEntry(ge2_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
mylegend2->AddEntry(ge2_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegend2->Draw();
gPad->Modified();
c1_mg1_e2.SaveAs("figs/inlist_c_bpm_diff.pdf");
c1_mg1_e2.SaveAs("figs/inlist_c_bpm_diff.png");
c1_mg1_e2.SaveAs("figs/inlist_c_bpm_diff.C");
///////////////////////////////////////////////
TCanvas c1_mg1_e3;
TGraphErrors * ge3_p1 = new TGraphErrors("inlistE255_p_posfit_diff.txt","%lg %lg %lg");
TGraphErrors * ge3_p2 = new TGraphErrors("inlistE21_p_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge3_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge3_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e3 = new TMultiGraph();
mg_e3->Add(ge3_p1,"pc"); ge3_p1->SetMarkerStyle(20);ge3_p1->SetLineStyle(3);
mg_e3->Add(ge3_p2,"pc"); ge3_p2->SetMarkerStyle(24);ge3_p2->SetLineStyle(3);
// mg_e3->Add(ge3_ic1,"p"); ge3_ic1->SetMarkerStyle(29);
// mg_e3->Add(ge3_ic2,"p"); ge3_ic2->SetMarkerStyle(30);
mg_e3->Draw("a");
mg_e3->SetTitle(" ");
mg_e3->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e3->GetYaxis()->SetTitle("#sigma(X_{1}-X_{2})/#sqrt{2} / mm");
// mg_e3->GetXaxis()->SetRangeUser(0,600);
//mg_e3->GetYaxis()->SetRangeUser(0,1);
mg_e3->SetMinimum(0.);
// mg_e3->SetMaximum(0.1);
TLegend *mylegend3 = new TLegend(0.60,0.6,0.9,0.9);
mylegend3->SetFillColor(0); // white background
mylegend3->SetTextFont(22);
mylegend3->SetBorderSize(0); // get rid of the box
mylegend3->SetTextSize(0.045); // set text size
mylegend3->AddEntry(ge3_p1,"BPM (221 MeV)","pl"); // options: p,l,f
mylegend3->AddEntry(ge3_p2,"BPM (70 MeV)","pl"); // options: p,l,f
// mylegend3->AddEntry(ge3_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend3->AddEntry(ge3_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegend3->Draw();
gPad->Modified();
c1_mg1_e3.SaveAs("figs/inlist_p_posfit_diff.pdf");
c1_mg1_e3.SaveAs("figs/inlist_p_posfit_diff.png");
c1_mg1_e3.SaveAs("figs/inlist_p_possfit_diff.C");
TCanvas c1_mg1_e4;
TGraphErrors * ge4_p1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
TGraphErrors * ge4_p2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge4_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge4_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e4 = new TMultiGraph();
mg_e4->Add(ge4_p1,"pc"); ge4_p1->SetMarkerStyle(22);ge4_p1->SetLineStyle(3);
mg_e4->Add(ge4_p2,"pc"); ge4_p2->SetMarkerStyle(26);ge4_p2->SetLineStyle(3);
// mg_e4->Add(ge4_ic1,"p"); ge4_ic1->SetMarkerStyle(29);
// mg_e4->Add(ge4_ic2,"p"); ge4_ic2->SetMarkerStyle(30);
mg_e4->Draw("a");
mg_e4->SetTitle(" ");
mg_e4->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e4->GetYaxis()->SetTitle("#sigma(X_{1}-X_{2})/#sqrt{2} / mm");
// mg_e4->GetXaxis()->SetRangeUser(0,600);
//mg_e4->GetYaxis()->SetRangeUser(0,1);
mg_e4->SetMinimum(0.);
// mg_e4->SetMaximum(0.1);
TLegend *mylegend4 = new TLegend(0.60,0.6,0.9,0.9);
mylegend4->SetFillColor(0); // white background
mylegend4->SetTextFont(22);
mylegend4->SetBorderSize(0); // get rid of the box
mylegend4->SetTextSize(0.045); // set text size
mylegend4->AddEntry(ge4_p1,"BPM (430 MeV/u)","pl"); // options: p,l,f
mylegend4->AddEntry(ge4_p2,"BPM (130 MeV/u)","pl"); // options: p,l,f
// mylegend4->AddEntry(ge4_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend4->AddEntry(ge4_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegend4->Draw();
gPad->Modified();
c1_mg1_e4.SaveAs("figs/inlist_c_posfit_diff.pdf");
c1_mg1_e4.SaveAs("figs/inlist_c_posfit_diff.png");
c1_mg1_e4.SaveAs("figs/inlist_c_possfit_diff.C");
///////////////////////////////////////////////
TCanvas c1_mg1_e5;
TGraphErrors * ge5_p1 = new TGraphErrors("inlistE255_p_focusfit_diff.txt","%lg %lg %lg");
TGraphErrors * ge5_p2 = new TGraphErrors("inlistE21_p_focusfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge5_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge5_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e5 = new TMultiGraph();
mg_e5->Add(ge5_p1,"pc"); ge5_p1->SetMarkerStyle(20);ge5_p1->SetLineStyle(3);
mg_e5->Add(ge5_p2,"pc"); ge5_p2->SetMarkerStyle(24);ge5_p2->SetLineStyle(3);
// mg_e5->Add(ge5_ic1,"p"); ge5_ic1->SetMarkerStyle(29);
// mg_e5->Add(ge5_ic2,"p"); ge5_ic2->SetMarkerStyle(30);
mg_e5->Draw("a");
mg_e5->SetTitle(" ");
mg_e5->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
//mg_e5->GetYaxis()->SetTitle("#sqrt{2}#upoint#sigma(F_{1}-F_{2})/(F_{1}+F_{2})");
mg_e5->GetYaxis()->SetTitle("#sigma(F_{1}-F_{2})/#sqrt{2} / mm");
// mg_e5->GetXaxis()->SetRangeUser(0,600);
//mg_e5->GetYaxis()->SetRangeUser(0,1);
mg_e5->SetMinimum(0.);
// mg_e5->SetMaximum(0.1);
TLegend * mylegend5 = new TLegend(0.60,0.6,0.9,0.9);
mylegend5->SetFillColor(0); // white background
mylegend5->SetTextFont(22);
mylegend5->SetBorderSize(0); // get rid of the box
mylegend5->SetTextSize(0.045); // set text size
mylegend5->AddEntry(ge5_p1,"BPM (221 MeV)","pl"); // options: p,l,f
mylegend5->AddEntry(ge5_p2,"BPM (70 MeV)","pl"); // options: p,l,f
// mylegend5->AddEntry(ge5_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend5->AddEntry(ge5_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegend5->Draw();
gPad->Modified();
c1_mg1_e5.SaveAs("figs/inlist_p_focusfit_diff.pdf");
c1_mg1_e5.SaveAs("figs/inlist_p_focusfit_diff.png");
c1_mg1_e5.SaveAs("figs/inlist_p_focusfit_diff.C");
TCanvas c1_mg1_e6;
TGraphErrors * ge6_p1 = new TGraphErrors("inlistE255_c_focusfit_diff.txt","%lg %lg %lg");
TGraphErrors * ge6_p2 = new TGraphErrors("inlistE21_c_focusfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge6_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge6_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e6 = new TMultiGraph();
mg_e6->Add(ge6_p1,"pc"); ge6_p1->SetMarkerStyle(22);ge6_p1->SetLineStyle(3);
mg_e6->Add(ge6_p2,"pc"); ge6_p2->SetMarkerStyle(26);ge6_p2->SetLineStyle(3);
// mg_e6->Add(ge6_ic1,"p"); ge6_ic1->SetMarkerStyle(29);
// mg_e6->Add(ge6_ic2,"p"); ge6_ic2->SetMarkerStyle(30);
mg_e6->Draw("a");
mg_e6->SetTitle(" ");
mg_e6->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e6->GetYaxis()->SetTitle("#sqrt{2}#upoint#sigma(F_{1}-F_{2})/(F_{1}+F_{2})");
mg_e6->GetYaxis()->SetTitle("#sigma(F_{1}-F_{2})/#sqrt{2} / mm");
// mg_e6->GetXaxis()->SetRangeUser(0,600);
//mg_e6->GetYaxis()->SetRangeUser(0,1);
mg_e6->SetMinimum(0.);
// mg_e6->SetMaximum(0.1);
TLegend * mylegend6 = new TLegend(0.60,0.6,0.9,0.9);
mylegend6->SetFillColor(0); // white background
mylegend6->SetTextFont(22);
mylegend6->SetBorderSize(0); // get rid of the box
mylegend6->SetTextSize(0.045); // set text size
mylegend6->AddEntry(ge6_p1,"BPM (430 MeV/u)","pl"); // options: p,l,f
mylegend6->AddEntry(ge6_p2,"BPM (130 MeV/u)","pl"); // options: p,l,f
// mylegend6->AddEntry(ge6_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend6->AddEntry(ge6_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegend6->Draw();
gPad->Modified();
c1_mg1_e6.SaveAs("figs/inlist_c_focusfit_diff.pdf");
c1_mg1_e6.SaveAs("figs/inlist_c_focusfit_diff.png");
c1_mg1_e6.SaveAs("figs/inlist_c_focusfit_diff.C");
TCanvas c1_mg1_e7;
TGraphErrors * ge7_p1 = new TGraphErrors("inlistE255_c_peakfit.txt","%lg %lg %lg");
TGraphErrors * ge7_p2 = new TGraphErrors("inlistE21_c_peakfit.txt","%lg %lg %lg");
// TGraphErrors * ge7_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge7_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e7 = new TMultiGraph();
mg_e7->Add(ge7_p1,"pc"); ge7_p1->SetMarkerStyle(22);ge7_p1->SetLineStyle(1);
mg_e7->Add(ge7_p2,"pc"); ge7_p2->SetMarkerStyle(26);ge7_p2->SetLineStyle(1);
// mg_e7->Add(ge7_ic1,"p"); ge7_ic1->SetMarkerStyle(29);
// mg_e7->Add(ge7_ic2,"p"); ge7_ic2->SetMarkerStyle(30);
mg_e7->Draw("a");
mg_e7->SetTitle(" ");
mg_e7->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e7->GetYaxis()->SetTitle("SNR");
//mg_e7->GetYaxis()->SetTitle("#sigma(F_{1}-F_{2})/#sqrt{2}");
// mg_e7->GetXaxis()->SetRangeUser(0,600);
//mg_e7->GetYaxis()->SetRangeUser(0,1);
mg_e7->SetMinimum(1.);
gPad->SetLogy();
// mg_e7->SetMaximum(0.1);
TLegend *mylegende7 = new TLegend(0.60,0.2,0.9,0.5);
mylegende7->SetFillColor(0); // white background
mylegende7->SetTextFont(22);
mylegende7->SetBorderSize(0); // get rid of the box
mylegende7->SetTextSize(0.045); // set text size
mylegende7->AddEntry(ge7_p1,"BPM (430 MeV/u)","p"); // options: p,l,f
mylegende7->AddEntry(ge7_p2,"BPM (130 MeV/u)","p"); // options: p,l,f
// mylegend6->AddEntry(ge7_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend6->AddEntry(ge7_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegende7->Draw();
gPad->Modified();
c1_mg1_e7.SaveAs("figs/inlist_c_peakfit.pdf");
c1_mg1_e7.SaveAs("figs/inlist_c_peakfit.png");
c1_mg1_e7.SaveAs("figs/inlist_c_peakfit.C");
TCanvas c1_mg1_e8;
TGraphErrors * ge8_p1 = new TGraphErrors("inlistE255_p_peakfit.txt","%lg %lg %lg");
TGraphErrors * ge8_p2 = new TGraphErrors("inlistE21_p_peakfit.txt","%lg %lg %lg");
// TGraphErrors * ge8_c1 = new TGraphErrors("inlistE255_c_posfit_diff.txt","%lg %lg %lg");
// TGraphErrors * ge8_c2 = new TGraphErrors("inlistE21_c_posfit_diff.txt","%lg %lg %lg");
TMultiGraph * mg_e8 = new TMultiGraph();
mg_e8->Add(ge8_p1,"pc"); ge8_p1->SetMarkerStyle(20);ge8_p1->SetLineStyle(1);
mg_e8->Add(ge8_p2,"pc"); ge8_p2->SetMarkerStyle(24);ge8_p2->SetLineStyle(1);
// mg_e8->Add(ge8_ic1,"p"); ge8_ic1->SetMarkerStyle(29);
// mg_e8->Add(ge8_ic2,"p"); ge8_ic2->SetMarkerStyle(30);
mg_e8->Draw("a");
mg_e8->SetTitle(" ");
mg_e8->GetXaxis()->SetTitle("Intensity / 10^{6}s^{-1}");
mg_e8->GetYaxis()->SetTitle("SNR");
//mg_e8->GetYaxis()->SetTitle("#sigma(F_{1}-F_{2})/#sqrt{2}");
// mg_e8->GetXaxis()->SetRangeUser(0,600);
//mg_e8->GetYaxis()->SetRangeUser(0,1);
mg_e8->SetMinimum(1.);
gPad->SetLogy();
// mg_e8->SetMaximum(0.1);
TLegend * mylegende8 = new TLegend(0.60,0.2,0.9,0.5);
mylegende8->SetFillColor(0); // white background
mylegende8->SetTextFont(22);
mylegende8->SetBorderSize(0); // get rid of the box
mylegende8->SetTextSize(0.045); // set text size
mylegende8->AddEntry(ge8_p1,"BPM (221 MeV/u)","p"); // options: p,l,f
mylegende8->AddEntry(ge8_p2,"BPM (70 MeV/u)","p"); // options: p,l,f
// mylegend6->AddEntry(ge8_ic1,"IC (430 MeV/u)","p"); // options: p,l,f
// mylegend6->AddEntry(ge8_ic2,"IC (130 MeV/u)","p"); // options: p,l,f
mylegende8->Draw();
gPad->Modified();
c1_mg1_e8.SaveAs("figs/inlist_p_peakfit.pdf");
c1_mg1_e8.SaveAs("figs/inlist_p_peakfit.png");
c1_mg1_e8.SaveAs("figs/inlist_p_peakfit.C");
}