Skip to content

Commit

Permalink
Add the compare code for PID validation
Browse files Browse the repository at this point in the history
  • Loading branch information
Biao Zhang committed Jun 4, 2021
1 parent eaee852 commit 93ce951
Showing 1 changed file with 202 additions and 0 deletions.
202 changes: 202 additions & 0 deletions analysis/SOverBCompvsPID.C
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,202 @@
void SOverBCompvsPID() {

TH1D *sig_combine;
TH1D *sig_nopid;
TH1D *sig_tof;
TH1D *sig_rich;

Int_t nBinsMerge = 9;
Double_t xAxis1[10] = {0,0.5,1,2,3,4,5,7,10,15};

TFile *filecombine = TFile::Open("combine/foutputJpsitoee.root","READ");
cout<<"ok"<<endl;
sig_combine = new TH1D(*((TH1D*)filecombine->Get("hBkgPerEvent")));

TFile *filenopid = TFile::Open("nopid/foutputJpsitoee.root","READ");
cout<<"ok"<<endl;
sig_nopid = new TH1D(*((TH1D*)filenopid->Get("hBkgPerEvent")));

TFile *filetof = TFile::Open("tof/foutputJpsitoee.root","READ");
sig_tof = new TH1D(*((TH1D*)filetof->Get("hBkgPerEvent")));

TFile *filerich = TFile::Open("rich/foutputJpsitoee.root","READ");
sig_rich = new TH1D(*((TH1D*)filerich->Get("hBkgPerEvent")));



int textfont = 42;
float textsize = 0.047;
float labelsize = 0.04;

TLatex* lat = new TLatex();
lat->SetNDC();
lat->SetTextFont(textfont);
lat->SetTextSize(textsize);

TLatex* lat2 = new TLatex();
lat2->SetNDC();
lat2->SetTextFont(textfont);
lat2->SetTextSize(textsize*1.0);

TLatex* lat3 = new TLatex();
lat3->SetNDC();
lat3->SetTextFont(textfont);
lat3->SetTextSize(textsize*1.0);

TCanvas *canvas3a = new TCanvas("canvas3a", "",443,100,700,700);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
gStyle->SetLegendBorderSize(0);
gStyle->SetLegendFont(textfont);
gStyle->SetLegendTextSize(textsize);

canvas3a->Range(0,0,1,1);
canvas3a->SetFillColor(0);
canvas3a->SetBorderMode(0);
canvas3a->SetBorderSize(2);
canvas3a->SetFrameBorderMode(0);


// ------------>Primitives in pad: canvas3a_1
TPad *canvas3a_1 = new TPad("canvas3a_1", "canvas3a_1",0.01,0.3,0.99,0.99);
canvas3a_1->Draw();
canvas3a_1->cd();
canvas3a_1->Range(-0.3750001,-2248.875,13.375,20249.88);
canvas3a_1->SetFillColor(0);
canvas3a_1->SetBorderMode(0);
canvas3a_1->SetBorderSize(2);
//canvas3a_1->SetGridx();
//canvas3a_1->SetGridy();
canvas3a_1->SetLeftMargin(0.15);
canvas3a_1->SetFrameBorderMode(0);
canvas3a_1->SetFrameBorderMode(0);

cout<<"ok"<<endl;

//sig_sc3->SetMarkerSize(1.25);
//sig_sc3->SetMinimum(2.27);
//sig_sc3->SetMaximum(2.30);
sig_combine->SetEntries(7);
sig_combine->SetStats(0);
sig_combine->SetLineWidth(2);
sig_combine->SetLineColor(1);
//sig_sc3->SetMarkerColor(1);
//sig_sc3->SetMarkerStyle(21);
sig_combine->GetXaxis()->SetTitle("#it{p}_{T}(J/#psi) [GeV/#it{c}]");
sig_combine->GetXaxis()->SetNdivisions(511);
sig_combine->GetXaxis()->SetRangeUser(0,15);
sig_combine->GetXaxis()->SetLabelFont(42);
sig_combine->GetXaxis()->SetLabelSize(0.04);
sig_combine->GetXaxis()->SetTitleSize(0.04);
sig_combine->GetXaxis()->SetTitleOffset(1.0);
sig_combine->GetYaxis()->SetRangeUser(0.000000003,10000);
sig_combine->GetXaxis()->SetTitleFont(42);

// sig_sc3->GetYaxis()->SetTitle("#frac{d^{2}#sigma}{dy dp_{T}} [#frac{#mu b}{GeV/c}]");
//sig_combine->GetYaxis()->SetTitle("Signal/Background (3#sigma)");
//sig_combine->GetYaxis()->SetTitle("Reconstruction Efficiency");
//sig_combine->GetYaxis()->SetTitle("Significance(3#sigma)");
sig_combine->GetYaxis()->SetTitle("Bkg/nEv (3#sigma)");
sig_combine->GetYaxis()->SetLabelFont(42);
sig_combine->GetYaxis()->SetTitleFont(42);
sig_combine->GetYaxis()->SetLabelSize(0.05);
sig_combine->GetYaxis()->SetTitleSize(0.07);
sig_combine->GetYaxis()->SetTitleOffset(0.9);
sig_combine->GetYaxis()->SetDecimals(kTRUE);

canvas3a_1->SetLogy();
sig_combine->Draw("");

sig_nopid->SetLineWidth(2);
sig_nopid->SetLineColor(2);
sig_tof->SetLineWidth(2);
sig_tof->SetLineColor(3);
sig_rich->SetLineWidth(2);
sig_rich->SetLineColor(4);
// sig_sc2->SetMarkerStyle(21);
// sig_sc2->SetMarkerSize(1.25);
//sig_sc2->SetMarkerColor(2);
sig_nopid->Draw("same");
sig_tof->Draw("same");
sig_rich->Draw("same");
lat->DrawLatex(0.2,0.83,"ALICE3 projection, pp #sqrt{s}= 14 TeV");
lat2->DrawLatex(0.2,0.76,"Pythia 8, mode = 2, N_{ev}=210 10^{12}");
lat3->DrawLatex(0.2,0.69,"J/#psi #rightarrow e^{+}e^{-}, |y| < 1.44, BR=5.94%");


/////////////////////////////////////////////////////////////////////// RATIO /////////////////////////////////////////////////////////////////////

TLegend *leg = new TLegend(0.6,0.15,0.9,0.40,0,"NDC");
leg->AddEntry(sig_combine,"TOF + RICH","lp");
leg->AddEntry(sig_nopid,"w/o PID","lp");
leg->AddEntry(sig_tof,"TOF","lp");
leg->AddEntry(sig_rich,"RICH","lp");
leg->SetFillStyle(0);
leg->SetBorderSize(0);
leg->SetFillColor(0);
//leg->SetTextSize(22);

leg->Draw();

canvas3a_1->Modified();
canvas3a->cd();


// ------------>Primitives in pad: canvas3a_2
TPad *canvas3a_2 = new TPad("canvas3a_2", "canvas3a_2",0.01,0.01,0.99,0.365);
canvas3a_2->Draw();
canvas3a_2->cd();
canvas3a_2->Range(-0.3750001,-0.7452094,13.375,0.6383566);
canvas3a_2->SetFillColor(0);
canvas3a_2->SetBorderMode(0);
canvas3a_2->SetBorderSize(2);
//canvas3a_2->SetGridx();
canvas3a_2->SetGridy();
canvas3a_2->SetLeftMargin(0.15);
canvas3a_2->SetBottomMargin(0.2);
canvas3a_2->SetFrameBorderMode(0);
canvas3a_2->SetFrameBorderMode(0);

TH1F *SigRatio = new TH1F("SigRatio","",9,xAxis1);
SigRatio->Divide(sig_combine,sig_nopid,1.,1.);

//SigRatio->SetMarkerStyle(1);
//SigRatio->SetMarkerSize(1);
SigRatio->SetMinimum(0.99);
SigRatio->SetMaximum(1.01);
SigRatio->SetEntries(7);
SigRatio->SetStats(0);
SigRatio->SetLineWidth(2);
SigRatio->SetLineColor(4);
SigRatio->GetXaxis()->SetTitle("p_{T}(J/#psi) [GeV/c]");
//SigRatio->GetXaxis()->SetNdivisions(511);
SigRatio->GetXaxis()->SetLabelFont(42);
SigRatio->GetXaxis()->SetLabelSize(0.085);
SigRatio->GetXaxis()->SetTitleSize(0.095);
SigRatio->GetXaxis()->SetTitleOffset(1.0);
// SigRatio->GetXaxis()->SetRangeUser(0,500);
SigRatio->GetYaxis()->SetRangeUser(-0.0000001,0.15);
SigRatio->GetXaxis()->SetTitleFont(42);
SigRatio->GetYaxis()->SetTitle("(TOF + RICH) / w/o PID");
SigRatio->GetYaxis()->SetNdivisions(106);
SigRatio->GetYaxis()->SetLabelFont(42);
SigRatio->GetYaxis()->SetLabelSize(0.085);
SigRatio->GetYaxis()->SetTitleSize(0.08);
SigRatio->GetYaxis()->SetTitleOffset(0.7);
SigRatio->GetYaxis()->SetTitleFont(42);
SigRatio->GetYaxis()->SetDecimals(kTRUE);
SigRatio->Draw("");

TLine *line = new TLine(0,1,15.,1);
line->SetLineColor(1);
line->SetLineWidth(4);
line->SetLineStyle(2);
line->Draw("same");

canvas3a_2->Modified();
canvas3a->cd();
canvas3a->Modified();
canvas3a->cd();
canvas3a->SetSelected(canvas3a);

}

0 comments on commit 93ce951

Please sign in to comment.