LightFlavorRatios/swimming_correction/plot_cutsystematics.C

0001 std::string getDate()
0002 {
0003     std::time_t t = std::time(0);   // get time now
0004     std::tm* now = std::localtime(&t);
0005
0006     std::stringstream date;
0007     date << (now->tm_mon + 1) << '/'
0008          <<  now->tm_mday << '/'
0009          << (now->tm_year + 1900);
0010     return date.str();
0011 }
0012
0013 void plot_cutsystematics()
0014 {
0015   gStyle->SetOptStat(0);
0016
0017   TFile* f = TFile::Open("LamdbaKsCutEfficiency_200MeV_hists.root");
0018
0019   TH1F* hCentral = (TH1F*)f->Get("hEffRatio_pT");
0020   TH1F* hMax = (TH1F*)f->Get("hEffRatio_pT_max");
0021   TH1F* hMin = (TH1F*)f->Get("hEffRatio_pT_min");
0022
0023   TH1F* hSys = (TH1F*)f->Get("cuteff_syserr_pT");
0024
0025   hCentral->SetLineColor(kBlack);
0026   hCentral->SetMarkerColor(kBlack);
0027   hCentral->SetMarkerStyle(kFullCircle);
0028   hCentral->SetMarkerSize(0.7);
0029   hCentral->GetYaxis()->SetRangeUser(0.0, 1.05);
0030
0031   hMin->SetLineColor(kBlue);
0032   hMin->SetMarkerColor(kBlue);
0033   hMin->SetMarkerStyle(kFullCircle);
0034   hMin->SetMarkerSize(0.7);
0035
0036   hMax->SetLineColor(kBlue);
0037   hMax->SetMarkerColor(kBlue);
0038   hMax->SetMarkerStyle(kFullCircle);
0039   hMin->SetMarkerSize(0.7);
0040
0041   hSys->SetFillColor(kBlue);
0042
0043   TCanvas* c = new TCanvas("c","c",800,800);
0044   hCentral->Draw("E");
0045   hMin->Draw("E SAME");
0046   hMax->Draw("E SAME");
0047
0048   TPaveText *pt;
0049   pt = new TPaveText(0.15,0.75,0.65,1., "NDC");
0050   pt->SetFillColor(0);
0051   pt->SetFillStyle(0);
0052   pt->SetTextFont(42);
0053   string label = "#it{#bf{sPHENIX}} Simulation";
0054   TText *pt_LaTex = pt->AddText(label.c_str());
0055   pt->SetBorderSize(0);
0056   pt->Draw();
0057   gPad->Modified();
0058
0059   TPaveText *ptDate;
0060   ptDate = new TPaveText(0.67,0.79,0.85,0.95, "NDC");
0061   ptDate->SetFillColor(0);
0062   ptDate->SetFillStyle(0);
0063   ptDate->SetTextFont(42);
0064   std::string compilation_date = getDate();
0065   TText *pt_LaTexDate = ptDate->AddText(compilation_date.c_str());
0066   ptDate->SetBorderSize(0);
0067   ptDate->Draw();
0068   gPad->Modified();
0069
0070   TCanvas* cs = new TCanvas("cs","cs",800,800);
0071   hSys->Draw();
0072   pt->Draw();
0073   ptDate->Draw();
0074   gPad->Modified();
0075
0076   c->SaveAs("systematics_varying_cuteff.pdf");
0077   cs->SaveAs("cuteff_sys_pT.pdf");
0078 }