systematics/trackMomScale/plot_sv_ks_result.C

0001 // Auto-generated by find_sv_smear_fracs.py
0002 // K0S: SV sim inv mass smeared to match data KFP mass.
0003 // Run with:  root -l -b -q plot_sv_ks_result.C
0004
0005 void plot_sv_ks_result() {
0006     gStyle->SetOptStat(0); gStyle->SetOptTitle(0);
0007
0008     const int N = 8;
0009     double avg_pt[]  = {0.2237, 0.2745, 0.3477, 0.4471, 0.5468, 0.6879, 0.8884, 1.3650};
0010     double sig_dk[]  = {3.6358, 3.9570, 4.5560, 5.1785, 5.8656, 6.3226, 7.4040, 9.0714};   // data KFP mass sigma  ← TARGET
0011     double sig_sv0[] = {3.2051, 2.8669, 2.5875, 2.6077, 2.5544, 2.6633, 2.9423, 3.1022};  // SV sim inv mass (no smear)
0012     double sig_m[]   = {3.6343, 3.9585, 4.5578, 5.1784, 5.8647, 6.3244, 7.3909, 9.0714};    // SV sim inv mass (matched smear)
0013     double fracs[]   = {0.7307, 1.1105, 1.5147, 1.7908, 2.1707, 2.2894, 2.8963, 3.4227};    // optimal pion smearing %
0014
0015     TCanvas *c1 = new TCanvas("c1","K0S SV smearing",950,620);
0016     c1->SetLeftMargin(0.12); c1->SetBottomMargin(0.13); c1->SetRightMargin(0.05);
0017
0018     auto sty = [](TGraph* g, Color_t col, Style_t mrk, Style_t ls=1){
0019         g->SetMarkerColor(col); g->SetLineColor(col);
0020         g->SetMarkerStyle(mrk); g->SetMarkerSize(1.6);
0021         g->SetLineWidth(2); g->SetLineStyle(ls);
0022     };
0023
0024     TGraph *gDK = new TGraph(N,avg_pt,sig_dk);
0025     TGraph *gS0 = new TGraph(N,avg_pt,sig_sv0);
0026     TGraph *gM  = new TGraph(N,avg_pt,sig_m);
0027
0028     sty(gDK, kBlack,  20, 1);
0029     sty(gS0, kBlue+1, 24, 2);
0030     sty(gM,  kRed+1,  22, 1);
0031
0032     TMultiGraph *mg = new TMultiGraph();
0033     mg->Add(gS0,"PL"); mg->Add(gM,"PL"); mg->Add(gDK,"PL");
0034     mg->Draw("A");
0035     mg->GetXaxis()->SetTitle("p_{T} (GeV/#it{c})");
0036     mg->GetYaxis()->SetTitle("#sigma_{K^{0}_{S}} (MeV/#it{c}^{2})");
0037     mg->GetYaxis()->SetRangeUser(0,12.2);
0038     mg->GetXaxis()->SetTitleSize(0.05);
0039     mg->GetYaxis()->SetTitleSize(0.05);
0040
0041     TLegend *leg = new TLegend(0.14,0.68,0.72,0.93);
0042     leg->SetBorderSize(0); leg->SetFillStyle(0); leg->SetTextSize(0.032);
0043     leg->AddEntry(gDK,"Data      (K_{S}^{0} mass, KFP) — target","lp");
0044     leg->AddEntry(gS0,"SV sim    (inv mass, no smear)","lp");
0045     leg->AddEntry(gM, "SV sim    (inv mass, matched smear)","lp");
0046     leg->Draw();
0047
0048     c1->SaveAs("plot_sv_ks_result.pdf");
0049     c1->SaveAs("plot_sv_ks_result.png");
0050     std::cout << "Saved plot_sv_ks_result.pdf/.png\n";
0051
0052     TCanvas *c2 = new TCanvas("c2","K0S pion smearing frac",900,550);
0053     c2->SetLeftMargin(0.13); c2->SetBottomMargin(0.13);
0054     TGraph *gF = new TGraph(N,avg_pt,fracs);
0055     sty(gF, kGreen+2, 23);
0056     gF->SetTitle(";p_{T} (GeV/#it{c});Optimal pion smearing (% of p_{T})");
0057     gF->GetYaxis()->SetRangeUser(0,15);
0058     gF->GetXaxis()->SetTitleSize(0.05); gF->GetYaxis()->SetTitleSize(0.05);
0059     gF->Draw("APL");
0060     c2->SaveAs("plot_sv_ks_frac.pdf");
0061     c2->SaveAs("plot_sv_ks_frac.png");
0062     std::cout << "Saved plot_sv_ks_frac.pdf/.png\n";
0063 }