cascadeAnalysis/finalPlots/cutEffPlots.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 cutEffPlots()
0014 {
0015    TFile* cutEffFile = TFile::Open("/sphenix/user/rosstom/TrackingAnalysis/Physics_Val_TF/sim/pTComp/cutEfficiencyCorrections.root", "READ");
0016
0017    TH1F* h_pT_KS0 = (TH1F*)cutEffFile->Get("h_cutEff_pT_KS0");
0018    TH1F* h_pT_Xi = (TH1F*)cutEffFile->Get("h_cutEff_pT_Xi");
0019    TH1F* h_eta_KS0 = (TH1F*)cutEffFile->Get("h_cutEff_eta_KS0");
0020    TH1F* h_eta_Xi = (TH1F*)cutEffFile->Get("h_cutEff_eta_Xi");
0021    TH1F* h_phi_KS0 = (TH1F*)cutEffFile->Get("h_cutEff_phi_KS0");
0022    TH1F* h_phi_Xi = (TH1F*)cutEffFile->Get("h_cutEff_phi_Xi");
0023    TH1F* h_rap_KS0 = (TH1F*)cutEffFile->Get("h_cutEff_rap_KS0");
0024    TH1F* h_rap_Xi = (TH1F*)cutEffFile->Get("h_cutEff_rap_Xi");
0025
0026    TH1F* h_ratio_pT = (TH1F*)h_pT_Xi->Clone("h_ratio_pT");
0027    TH1F* h_ratio_eta = (TH1F*)h_eta_Xi->Clone("h_ratio_eta");
0028    TH1F* h_ratio_phi = (TH1F*)h_phi_Xi->Clone("h_ratio_phi");
0029    TH1F* h_ratio_rap = (TH1F*)h_rap_Xi->Clone("h_ratio_rap");
0030
0031    h_ratio_pT->Divide(h_pT_Xi,h_pT_KS0,1.,1.);
0032    h_ratio_eta->Divide(h_eta_Xi,h_eta_KS0,1.,1.);
0033    h_ratio_phi->Divide(h_phi_Xi,h_phi_KS0,1.,1.);
0034    h_ratio_rap->Divide(h_rap_Xi,h_rap_KS0,1.,1.);
0035
0036    //pT Figure
0037    h_ratio_pT->GetYaxis()->SetTitle("KFParticle Cut Eff. Ratio");
0038    h_ratio_pT->GetXaxis()->SetTitle("Mother p_{T} [GeV/c]");
0039    h_ratio_pT->SetMaximum(1.0);
0040    h_ratio_pT->SetMinimum(0.0);
0041    TCanvas *c1  = new TCanvas("myCanvas", "myCanvas",800,800);
0042    gPad->SetTopMargin(0.1);
0043    h_ratio_pT->Draw("EP");
0044
0045    TPaveText *pt;
0046    pt = new TPaveText(0.15,0.75,0.65,1., "NDC");
0047    pt->SetFillColor(0);
0048    pt->SetFillStyle(0);
0049    pt->SetTextFont(42);
0050    string label = "#it{#bf{sPHENIX}} Simulation";
0051    TText *pt_LaTex = pt->AddText(label.c_str());
0052    pt->SetBorderSize(0);
0053    pt->Draw();
0054    gPad->Modified();
0055
0056    TPaveText *ptDate;
0057    ptDate = new TPaveText(0.67,0.79,0.85,0.95, "NDC");
0058    ptDate->SetFillColor(0);
0059    ptDate->SetFillStyle(0);
0060    ptDate->SetTextFont(42);
0061    std::string compilation_date = getDate();
0062    TText *pt_LaTexDate = ptDate->AddText(compilation_date.c_str());
0063    ptDate->SetBorderSize(0);
0064    ptDate->Draw();
0065    gPad->Modified();
0066
0067    c1->SaveAs("cutEff_pT.png");
0068    c1->Close();
0069
0070    //eta Figure
0071    h_ratio_eta->GetYaxis()->SetTitle("KFParticle Cut Eff. Ratio");
0072    h_ratio_eta->GetXaxis()->SetTitle("Mother #eta");
0073    h_ratio_eta->SetMaximum(1.0);
0074    h_ratio_eta->SetMinimum(0.0);
0075    TCanvas *c2  = new TCanvas("myCanvas2", "myCanvas2",800,800);
0076    gPad->SetTopMargin(0.1);
0077    h_ratio_eta->Draw("EP");
0078
0079    pt->Draw();
0080    gPad->Modified();
0081
0082    ptDate->Draw();
0083    gPad->Modified();
0084
0085    c2->SaveAs("cutEff_eta.png");
0086    c2->Close();
0087
0088    //phi Figure
0089    h_ratio_phi->GetYaxis()->SetTitle("KFParticle Cut Eff. Ratio");
0090    h_ratio_phi->GetXaxis()->SetTitle("Mother #phi [rad]");
0091    h_ratio_phi->SetMaximum(1.0);
0092    h_ratio_phi->SetMinimum(0.0);
0093    TCanvas *c3  = new TCanvas("myCanvas3", "myCanvas3",800,800);
0094    gPad->SetTopMargin(0.1);
0095    h_ratio_phi->Draw("EP");
0096
0097    pt->Draw();
0098    gPad->Modified();
0099
0100    ptDate->Draw();
0101    gPad->Modified();
0102
0103    c3->SaveAs("cutEff_phi.png");
0104    c3->Close();
0105
0106    //rap Figure
0107    h_ratio_rap->GetYaxis()->SetTitle("KFParticle Cut Eff. Ratio");
0108    h_ratio_rap->GetXaxis()->SetTitle("Mother Rapidity");
0109    h_ratio_rap->SetMaximum(1.0);
0110    h_ratio_rap->SetMinimum(0.0);
0111    TCanvas *c4  = new TCanvas("myCanvas4", "myCanvas4",800,800);
0112    gPad->SetTopMargin(0.1);
0113    h_ratio_rap->Draw("EP");
0114
0115    pt->Draw();
0116    gPad->Modified();
0117
0118    ptDate->Draw();
0119    gPad->Modified();
0120
0121    c4->SaveAs("cutEff_rap.png");
0122    c4->Close();
0123 }