Back to home page

sPhenix code displayed by LXR

 
 

    


Warning, file /analysis/LightFlavorRatios/yield_and_ratios/plot_fitsystematics.C was not indexed or was modified since last indexation (in which case cross-reference links may be missing, inaccurate or erroneous).

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_fitsystematics()
0014 {
0015   gStyle->SetOptStat(0);
0016 
0017   TFile* f = TFile::Open("fits.root");
0018 
0019   TH1F* Ks_sys = (TH1F*)f->Get("K_S0_fit_syserr_vspT");
0020   TH1F* Lambda_sys = (TH1F*)f->Get("Lambda0_fit_syserr_vspT");
0021   TH1F* ratio_sys = (TH1F*)f->Get("ratio_fit_syserr_vspT");
0022 
0023   Ks_sys->SetFillColor(kAzure);
0024   Lambda_sys->SetFillColor(kAzure);
0025   ratio_sys->SetFillColor(kAzure);
0026 
0027   Lambda_sys->GetYaxis()->SetRangeUser(0.,1.1*Lambda_sys->GetMaximum());
0028 
0029   TCanvas* c = new TCanvas("c","c",800,800);
0030   Ks_sys->Draw();
0031 
0032   TPaveText *pt;
0033   pt = new TPaveText(0.15,0.75,0.65,1., "NDC");
0034   pt->SetFillColor(0);
0035   pt->SetFillStyle(0);
0036   pt->SetTextFont(42);
0037   string label = "#it{#bf{sPHENIX}} Simulation";
0038   TText *pt_LaTex = pt->AddText(label.c_str());
0039   pt->SetBorderSize(0);
0040   pt->Draw();
0041   gPad->Modified();
0042 
0043   TPaveText *ptDate;
0044   ptDate = new TPaveText(0.67,0.79,0.85,0.95, "NDC");
0045   ptDate->SetFillColor(0);
0046   ptDate->SetFillStyle(0);
0047   ptDate->SetTextFont(42);
0048   std::string compilation_date = getDate();
0049   TText *pt_LaTexDate = ptDate->AddText(compilation_date.c_str());
0050   ptDate->SetBorderSize(0);
0051   ptDate->Draw();
0052   gPad->Modified();
0053 
0054   TCanvas* cl = new TCanvas("cl","cl",800,800);
0055   Lambda_sys->Draw();
0056   pt->Draw();
0057   ptDate->Draw();
0058   gPad->Modified();
0059 
0060   TCanvas* cr = new TCanvas("cr","cr",800,800);
0061   ratio_sys->Draw();
0062   pt->Draw();
0063   ptDate->Draw();
0064   gPad->Modified();
0065 
0066   c->SaveAs("Kshort_yield_sys.pdf");
0067   cl->SaveAs("Lambda_yield_sys.pdf");
0068   cr->SaveAs("ratio_yield_sys.pdf");
0069 }