void NuclearModification(
TString particle,
const int nbins,
Double_t xbins[],
Double_t xhbins[],
Double_t exl[],
Double_t exl2[],
Double_t exl0[],
Double_t yPercSigmapPbSystTotHigh[],
Double_t yPercSigmapPbSystTotLow[],
Double_t commonErrorP,
Double_t commonErrorN,
Double_t FFsysterror,
Double_t tagandprobcorrection[],
TCanvas *canvasSigma,
TCanvas *canvasRpA,
int PadNum
){
gROOT->SetStyle("Plain");
gStyle->SetOptTitle(0);
gStyle->SetOptStat(0);
TFile*filePPReference=new TFile(Form("../../../fonll/output%s.root",particle.Data()));
//###TFile*filePPReference=new TFile(Form("../../fonll/output%s.root",particle.Data()));
TGraphAsymmErrors*gaeBplusReference=(TGraphAsymmErrors*)filePPReference->Get(Form("gaeSigmaDecay%s",particle.Data()));
gaeBplusReference->SetName(Form("gae%sReference",particle.Data()));
TFile*filepPb=new TFile(Form("../Results%s/Sigma%s.root",particle.Data(),particle.Data()));
TH1F*hSigmapPbStat=(TH1F*)filepPb->Get("hPtSigma");
TH1F*hPt=(TH1F*)filepPb->Get("hPt");
TH1F*hEff=(TH1F*)filepPb->Get("hEff");
double scalingfactor=1e-6;
double yvalue,xvalue,yerrorhigh,yerrorlow;
for (int i=0;i<nbins;i++){
hSigmapPbStat->SetBinContent(i+1,scalingfactor*hSigmapPbStat->GetBinContent(i+1));
hSigmapPbStat->SetBinError(i+1,scalingfactor*hSigmapPbStat->GetBinError(i+1));
yvalue=-1.;
xvalue=-1.;
yerrorhigh=-1.;
yerrorlow=-1.;
gaeBplusReference->GetPoint(i,xvalue,yvalue);
yerrorhigh=gaeBplusReference->GetEYhigh()[i];
yerrorlow=gaeBplusReference->GetEYlow()[i];
gaeBplusReference->SetPoint(i,xvalue,yvalue*scalingfactor);
gaeBplusReference->SetPointEYhigh(i,yerrorhigh*scalingfactor);
gaeBplusReference->SetPointEYlow(i,yerrorlow*scalingfactor);
}
for (int i=0;i<nbins;i++){
hSigmapPbStat->SetBinContent(i+1,(1./tagandprobcorrection[i])*hSigmapPbStat->GetBinContent(i+1));
hSigmapPbStat->SetBinError(i+1,(1./tagandprobcorrection[i])*hSigmapPbStat->GetBinError(i+1));
}
Double_t yRefPP[nbins]; //value y reference
Double_t xRefPP[nbins]; //value x reference
Double_t yPPsystFONLLhigh[nbins]; //y err syst FONLL high
Double_t yPPsystFONLLlow[nbins]; //y err syst FONLL low
Double_t yPercPPsystFONLLhigh[nbins]; //y percentuale err syst FONLL high
Double_t yPercPPsystFONLLlow[nbins]; //y percentuale err syst FONLL low
Double_t ySigmapPb[nbins]; //value y pPb
Double_t xSigmapPb[nbins]; //value x pPb
Double_t ySigmapPbStat[nbins]; //y err stat pPb
Double_t yPercSigmapPbStat[nbins]; //y err stat pPb
Double_t yFONLL[nbins]; //1
Double_t yRpA[nbins]; //value y RpA
Double_t yRpAStat[nbins]; //y err stat RpA
Double_t yRpAsystFONLLhigh[nbins]; //y err syst FONLL RpA high
Double_t yRpAsystFONLLlow[nbins]; //y err syst FONLL RpA lzow
Double_t yPercRpAsystFONLLhigh[nbins]; //y percentuale err syst FONLL RpA high
Double_t yPercRpAsystFONLLlow[nbins]; //y percentuale err syst FONLL RpA low
Double_t ySigmapPbSystTotHigh[nbins]; //y percentuale err syst pPb TOT
Double_t ySigmapPbSystTotLow[nbins]; //y percentuale err syst pPb TOT
Double_t yPercRpPbSystTotHigh[nbins]; //y percentuale err syst RpPb TOT
Double_t yPercRpPbSystTotLow[nbins]; //y percentuale err syst RpPb TOT
Double_t yRpPbSystTotHigh[nbins]; //y percentuale err syst RpPb TOT
Double_t yRpPbSystTotLow[nbins]; //y percentuale err syst RpPb TOT
double x,y;
for (Int_t i=0;i<nbins;i++) {
gaeBplusReference->GetPoint(i,xRefPP[i],yRefPP[i]);
yPPsystFONLLhigh[i]=gaeBplusReference->GetEYhigh()[i];
yPPsystFONLLlow[i]=gaeBplusReference->GetEYlow()[i];
yPercPPsystFONLLhigh[i]=yPPsystFONLLhigh[i]/yRefPP[i];
yPercPPsystFONLLlow[i]=yPPsystFONLLlow[i]/yRefPP[i];
yPercPPsystFONLLhigh[i]=TMath::Sqrt(yPercPPsystFONLLhigh[i]*yPercPPsystFONLLhigh[i]+FFsysterror*FFsysterror);
yPercPPsystFONLLlow[i]=TMath::Sqrt(yPercPPsystFONLLlow[i]*yPercPPsystFONLLlow[i]+FFsysterror*FFsysterror);
}
for(Int_t i=0;i<nbins;i++) {
ySigmapPb[i]=hSigmapPbStat->GetBinContent(i+1);
ySigmapPbStat[i]=hSigmapPbStat->GetBinError(i+1);
yPercSigmapPbStat[i]=ySigmapPbStat[i]/ySigmapPb[i];
ySigmapPbSystTotHigh[i]=yPercSigmapPbSystTotHigh[i]*ySigmapPb[i];
ySigmapPbSystTotLow[i]=yPercSigmapPbSystTotLow[i]*ySigmapPb[i];
}
for(Int_t i=0;i<nbins;i++) {
yRpA[i]=ySigmapPb[i]/yRefPP[i];
yRpAStat[i]=ySigmapPbStat[i]/yRefPP[i];
yFONLL[i]=yRpA[i];
yPercRpAsystFONLLhigh[i]=(yPercPPsystFONLLlow[i]/(1-yPercPPsystFONLLlow[i]));
yPercRpAsystFONLLlow[i]=(yPercPPsystFONLLhigh[i]/(1+yPercPPsystFONLLhigh[i]));
yRpAsystFONLLhigh[i]=yPercRpAsystFONLLhigh[i]*yRpA[i];
yRpAsystFONLLlow[i]=yPercRpAsystFONLLlow[i]*yRpA[i];
yRpPbSystTotHigh[i]=yPercSigmapPbSystTotHigh[i]*yRpA[i];
yRpPbSystTotLow[i]=yPercSigmapPbSystTotLow[i]*yRpA[i];
std::cout << i << " , " << xbins[i] << " , " << ySigmapPb[i] << " , sta: " << ySigmapPbStat[i] << " , syslow: " << ySigmapPbSystTotLow[i] << " ,syshigh: " << ySigmapPbSystTotHigh[i] << std::endl;
std::cout << "FONLL: " << yRefPP[i] << " - " << yPPsystFONLLlow[i] << " + " << yPPsystFONLLhigh[i] << std::endl;
std::cout << i << " ####### " << xbins[i] << " , " << yRpA[i] << " , sta: " << yRpAStat[i] << " , syslow: " << yRpPbSystTotLow[i] << " ,syshigh: " << yRpPbSystTotHigh[i] << " ,FONLLlow: " << yRpAsystFONLLlow[i] << " , FONLLhigh: " << yRpAsystFONLLhigh[i] << std::endl;
TGraphAsymmErrors *gRpAstat = new TGraphAsymmErrors(nbins,xbins,yRpA,exl0,exl0,yRpAStat,yRpAStat);
TGraphAsymmErrors *gRpAsyst = new TGraphAsymmErrors(nbins,xbins,yRpA,exl,exl,yRpPbSystTotLow,yRpPbSystTotHigh);
TGraphAsymmErrors *gRpAsystFONLL = new TGraphAsymmErrors(nbins,xbins,yFONLL,exl,exl,yRpAsystFONLLlow,yRpAsystFONLLhigh);
}
TGraphAsymmErrors *gSigmasyst = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl2,exl2,ySigmapPbSystTotLow,ySigmapPbSystTotHigh);
//###TGraphAsymmErrors *gSigmasyst = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl,exl,ySigmapPbSystTotLow,ySigmapPbSystTotHigh);
gSigmasyst->SetTitle("Sigma syst uncertainty from pPb");
gSigmasyst->SetMarkerColor(1);
gSigmasyst->SetLineColor(1);
gSigmasyst->SetLineWidth(1);
gSigmasyst->SetMarkerStyle(21);
gSigmasyst->SetMarkerColor(1);
gSigmasyst->SetFillColor(kYellow-7);//5
gSigmasyst->SetFillStyle(1001);
TGraphAsymmErrors*gSigmasyst2=(TGraphAsymmErrors*)gSigmasyst->Clone();
gSigmasyst2->SetMarkerColor(1);
gSigmasyst2->SetMarkerStyle(25);
gSigmasyst2->SetFillColor(0);
gSigmasyst2->SetFillStyle(0);
gSigmasyst2->SetLineColor(1);//5
gSigmasyst2->SetLineStyle(1);
gSigmasyst2->SetLineWidth(1);
//###TGraphAsymmErrors *gSigmastat = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl,exl,ySigmapPbStat,ySigmapPbStat);
TGraphAsymmErrors *gSigmastat = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl0,exl0,ySigmapPbStat,ySigmapPbStat);
gSigmastat->SetTitle("Sigma stat uncertainty from pPb");
gSigmastat->SetMarkerColor(1);
gSigmastat->SetLineColor(1);
gSigmastat->SetLineWidth(2);//###
gSigmastat->SetMarkerStyle(21);
gSigmastat->SetMarkerColor(1);
gSigmastat->SetFillColor(0);
gSigmastat->SetFillStyle(0);
gSigmastat->SetFillStyle(0);
Double_t padcorrection;
if (particle=="Bplus") padcorrection=Bp_padratio/Bp_padratio;
else if (particle=="Bzero") padcorrection=Bp_padratio/B0_padratio;
else if (particle=="Bs") padcorrection=Bp_padratio/Bs_padratio;
std::cout << "padcorrection: " << padcorrection << std::endl;
canvasSigma->cd(PadNum);
//canvasSigma->Range(-1.989924,-0.2917772,25.49622,2.212202);
canvasSigma->SetFillColor(0);
canvasSigma->SetBorderMode(0);
canvasSigma->SetBorderSize(2);
//canvasSigma->SetLeftMargin(0.1451613);
//canvasSigma->SetRightMargin(0.05443548);
//canvasSigma->SetTopMargin(0.005);//0.08474576
//canvasSigma->SetBottomMargin(0.1165254);
canvasSigma->SetFrameBorderMode(0);
canvasSigma->SetFrameBorderMode(0);
canvasSigma->SetLogy();
//###TH2F* hempty=new TH2F("hempty","", 10, 0.1, 65., 10., 0.1, 1e3);
TH2F* hempty=new TH2F("hempty","", 10, 0.1, 65., 10., 0.3, 1e3);
hempty->GetXaxis()->SetTitle("p_{T} (GeV/c)");
//if(particle=="Bplus") hempty->GetYaxis()->SetTitle("d#sigma / dp_{T} (B^{+}) (pb GeV^{-1}c)");
//if(particle=="Bzero") hempty->GetYaxis()->SetTitle("d#sigma / dp_{T} (B^{0}) (pb GeV^{-1}c)");
//if(particle=="Bs") hempty->GetYaxis()->SetTitle("d#sigma / dp_{T} (B_{s}) (pb GeV^{-1}c)");
hempty->GetXaxis()->CenterTitle();
hempty->GetYaxis()->CenterTitle();
//###hempty->GetYaxis()->SetTitle("d#sigma / dp_{T}( #mub GeV^{-1}c)");
hempty->GetYaxis()->SetTitle("d#sigma / dp_{T}(#mub/(GeV/c))");
//###hempty->GetXaxis()->SetTitleOffset(1.0);//###1.0
if (particle=="Bplus") hempty->GetXaxis()->SetTitleOffset(0.90);//###1.0
else if (particle=="Bzero") hempty->GetXaxis()->SetTitleOffset(0.90);//###0.80
else if (particle=="Bs") hempty->GetXaxis()->SetTitleOffset(0.90);//###0.85
hempty->GetYaxis()->SetTitleOffset(1.0);//###1.3
hempty->GetXaxis()->SetTitleSize(0.070*padcorrection);//###0.055
hempty->GetYaxis()->SetTitleSize(0.070*padcorrection);//###0.055
hempty->GetXaxis()->SetTitleFont(42);
hempty->GetYaxis()->SetTitleFont(42);
hempty->GetXaxis()->SetLabelFont(42);
hempty->GetYaxis()->SetLabelFont(42);
hempty->GetXaxis()->SetLabelSize(0.060*padcorrection);//###0.055
hempty->GetYaxis()->SetLabelSize(0.060*padcorrection);//###0.055
if (particle=="Bplus") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.005
else if (particle=="Bzero") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.0001
else if (particle=="Bs") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.0005
//###hempty->GetXaxis()->SetLabelOffset(0.005);//###0.005
hempty->SetMaximum(2);
hempty->SetMinimum(0.);
hempty->Draw();
gaeBplusReference->SetMarkerColor(1);
gaeBplusReference->SetMarkerStyle(25);
gaeBplusReference->SetFillColor(kYellow-7);//5
gaeBplusReference->SetFillStyle(1001);
gaeBplusReference->SetLineColor(kAzure-3);
gaeBplusReference->SetLineWidth(1);
gSigmastat->SetMarkerColor(1);
gSigmastat->SetLineColor(1);
gSigmastat->SetLineWidth(2);
gSigmastat->SetMarkerStyle(21);
gSigmastat->SetMarkerColor(1);
// Histogram style
TH1D* hBplusReference = new TH1D("hBplusReference","",nbins,xhbins);
TH1D* hBplusReferenceEYhigh = new TH1D("hBplusReferenceEYhigh","",nbins,xhbins);
TH1D* hBplusReferenceEYlow = new TH1D("hBplusReferenceEYlow","",nbins,xhbins);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gaeBplusReference->GetPoint(i,xgae,ygae);
hBplusReference->SetBinContent(i+1,ygae);
hBplusReferenceEYhigh->SetBinContent(i+1,ygae+gaeBplusReference->GetEYhigh()[i]);
hBplusReferenceEYlow->SetBinContent(i+1,ygae-gaeBplusReference->GetEYlow()[i]);
std::cout << "### FONLL (" << i << ") : " << ygae << " + " << gaeBplusReference->GetEYhigh()[i] << " - " << gaeBplusReference->GetEYlow()[i] << std::endl;
}
//hBplusReference->SetMarkerColor(1);
//hBplusReference->SetMarkerStyle(25);
//hBplusReference->SetFillColor(kYellow-7);//5
//hBplusReference->SetFillStyle(1001);
hBplusReference->SetLineColor(kAzure-3);
hBplusReference->SetLineWidth(2);
hBplusReferenceEYhigh->SetLineColor(kAzure-3);
hBplusReferenceEYhigh->SetLineWidth(2);
hBplusReferenceEYhigh->SetLineStyle(2);
hBplusReferenceEYlow->SetLineColor(kAzure-3);
hBplusReferenceEYlow->SetLineWidth(2);
hBplusReferenceEYlow->SetLineStyle(2);
/*
hBplusReference->Draw("][,same");
hBplusReferenceEYhigh->Draw("][,same");
hBplusReferenceEYlow->Draw("][,same");
*/
/*
gaeBplusReference->Draw("2psame");//2same
TGraphAsymmErrors*gaeBplusReference2=(TGraphAsymmErrors*)gaeBplusReference->Clone();
gaeBplusReference2->SetMarkerColor(1);
gaeBplusReference2->SetMarkerStyle(25);
gaeBplusReference2->SetFillColor(0);
gaeBplusReference2->SetFillStyle(0);
gaeBplusReference2->SetLineColor(kAzure-3);
gaeBplusReference2->SetLineWidth(1);
gaeBplusReference2->Draw("2psame");//2same
*/
//coord. for legend for sigma in the B+ pannel
//###TLegend *legendSigma=new TLegend(0.468298,0.7045614,0.7678185,0.8757895,"");
//###TLegend *legendSigma=new TLegend(0.55,0.63,0.85,0.80,"");
TLegend *legendSigma=new TLegend(0.38,0.67,0.68,0.84,"");
legendSigma->SetBorderSize(0);
legendSigma->SetLineColor(0);
legendSigma->SetFillColor(0);
legendSigma->SetFillStyle(1001);
legendSigma->SetTextFont(42);
legendSigma->SetTextSize(0.060);//###0.045
/*
TBox *c = new TBox(3,1-commonErrorN,7,1+commonErrorP);
c->SetLineColor(5);
c->SetFillColor(5);
c->Draw();
TBox *c2 = c->Clone();
c2->SetLineColor(1);
c2->SetFillStyle(0);
c2->Draw();
*/
//TLegendEntry *ent_SigmapPb=legendSigma->AddEntry(gSigmastat,"pPb","pf");
TLegendEntry *ent_SigmapPb=legendSigma->AddEntry(gSigmasyst,"pPb","pf");
/*
ent_SigmapPb->SetTextFont(42);
ent_SigmapPb->SetLineColor(1);
ent_SigmapPb->SetMarkerColor(1);
ent_SigmapPb->SetFillColor(0);
ent_SigmapPb->SetFillStyle(0);
*/
ent_SigmapPb->SetTextFont(42);
ent_SigmapPb->SetMarkerColor(1);
ent_SigmapPb->SetMarkerStyle(25);
//ent_SigmapPb->SetFillColor(kYellow-7);//5
//ent_SigmapPb->SetFillStyle(1001);
ent_SigmapPb->SetLineColor(1);
ent_SigmapPb->SetLineWidth(0);
ent_SigmapPb->SetLineStyle(0);
/*
TLegendEntry *ent_Sigmapp=legendSigma->AddEntry(gaeBplusReference,"FONLL pp ref.","pf");
ent_Sigmapp->SetTextFont(42);
ent_Sigmapp->SetLineColor(kAzure-3);//5
ent_Sigmapp->SetLineStyle(1);
ent_Sigmapp->SetMarkerColor(1);
ent_Sigmapp->SetMarkerStyle(21);
*/
TLegendEntry *ent_Sigmapp=legendSigma->AddEntry(hBplusReference,"Scaled FONLL pp ref.","l");
ent_Sigmapp->SetTextFont(42);
ent_Sigmapp->SetLineColor(kAzure-3);//5
ent_Sigmapp->SetLineStyle(1);
ent_Sigmapp->SetMarkerColor(1);
ent_Sigmapp->SetMarkerStyle(21);
TLegendEntry *ent_Sigmapperr=legendSigma->AddEntry(hBplusReferenceEYhigh,"FONLL pp ref. uncert.","l");
ent_Sigmapperr->SetTextFont(42);
ent_Sigmapperr->SetLineColor(kAzure-3);//5
ent_Sigmapperr->SetLineStyle(2);
ent_Sigmapperr->SetMarkerColor(1);
ent_Sigmapperr->SetMarkerStyle(21);
//###gSigmasyst->SetFillColor(0);
//###gSigmasyst->SetFillStyle(0);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gSigmastat->GetPoint(i,xgae,ygae);
std::cout << "### cross section (" << i << ") : " << ygae << " stat: " << gSigmastat->GetEYhigh()[i] << " syst: " << gSigmasyst->GetEYhigh()[i] << std::endl;
}
gSigmasyst->Draw("2esame");//###2psame
gSigmasyst2->Draw("2esame");
hBplusReference->Draw("][,same");
hBplusReferenceEYhigh->Draw("][,same");
hBplusReferenceEYlow->Draw("][,same");
gSigmastat->SetFillColor(0);
gSigmastat->Draw("epsame");
/*
TBox *d = new TBox(0.1,1-commonErrorN,4,1+commonErrorP);
d->SetLineColor(1);
d->SetFillColor(0);
d->Draw();
*/
//if(PadNum==0 || PadNum==1){
if(PadNum==3 || PadNum==1 || PadNum==2){
legendSigma->Draw("same");
//TLatex * tlatex1=new TLatex(0.21,0.88801268,"CMS");
TLatex * tlatex1=new TLatex(0.19,0.86,"CMS");
tlatex1->SetNDC();
tlatex1->SetTextColor(1);
tlatex1->SetTextFont(62);//42
tlatex1->SetTextSize(0.07);//0.045
tlatex1->Draw();
/*
TLatex * tlatex12=new TLatex(0.35,0.86,"Preliminary");
tlatex12->SetNDC();
tlatex12->SetTextColor(1);
tlatex12->SetTextFont(52);//42
tlatex12->SetTextSize(0.05);//0.045
tlatex12->Draw();
*/
//TLatex * tlatexlumi=new TLatex(0.471371,0.88801268,"L = 34.8 nb^{-1} (pPb 5.02 TeV)");
//TLatex * tlatexlumi=new TLatex(0.471371,0.88801268,"34.6 nb^{-1} (pPb 5.02 TeV)");
TLatex * tlatexlumi;
//###TLatex * tlatexlumi=new TLatex(0.390,0.950,"34.6 nb^{-1} (pPb 5.02 TeV)");
if (PadNum==1) tlatexlumi = new TLatex(0.41,0.94,"34.6 nb^{-1} (pPb 5.02 TeV)");
else tlatexlumi = new TLatex(0.39,0.94,"34.6 nb^{-1} (pPb 5.02 TeV)");
tlatexlumi->SetNDC();
tlatexlumi->SetTextColor(1);
tlatexlumi->SetTextFont(42);
tlatexlumi->SetTextSize(0.06);//0.045
tlatexlumi->Draw();
}
//if(PadNum==2){
if(PadNum==3 || PadNum==1 || PadNum==2) {
//TLatex * tlatex4=new TLatex(0.60,0.73,"|y_{lab}| < 2.4");
TLatex * tlatex4=new TLatex(0.65,0.50,"|y_{lab}| < 2.4");
tlatex4->SetNDC();
tlatex4->SetTextColor(1);
tlatex4->SetTextFont(42);//42
tlatex4->SetTextSize(0.07);//0.045
tlatex4->Draw();
}
double GloUnc;
double xtl5=0.54;
if(PadNum==1) {GloUnc=Bp_commonErrorP*100;xtl5=0.54;}//###0.57
if(PadNum==2) {GloUnc=B0_commonErrorP*100;xtl5=0.47;}
if(PadNum==3) {GloUnc=Bs_commonErrorP*100;xtl5=0.42;}
xtl5=0.54;
//TLegend *legendSigma=new TLegend(0.55,0.63,0.85,0.80,"");
TLatex * tlatex5=new TLatex(xtl5,0.60,Form("Global uncert. %2.1f%%",GloUnc));
tlatex5->SetNDC();
tlatex5->SetTextColor(1);
tlatex5->SetTextFont(42);//42
tlatex5->SetTextSize(0.05*padcorrection);//0.045
tlatex5->Draw();
//###double xpos=0.8528226;
//###double ypos=0.6849894;
double xpos=0.90;//0.86
double ypos=0.85;//0.86
TString mypar;
if(particle=="Bplus") mypar="B^{+}";
if(particle=="Bzero") mypar="B^{0}";
if(particle=="Bs") mypar="B_{s}^{0}";
double xlgap=0.04;
/*
if(particle=="Bplus") xpos=(0.365816-xlgap)/0.365816;
if(particle=="Bzero") xpos=1-xlgap/(0.673101-0.365816);
if(particle=="Bs") xpos=0.94-xlgap/(1-0.673101);;
*/
TLatex * tlatex3=new TLatex(xpos,ypos,mypar.Data());
tlatex3->SetNDC();
tlatex3->SetTextColor(1);
tlatex3->SetTextFont(42);
tlatex3->SetTextSize(0.07*padcorrection);
tlatex3->Draw();
TGraphAsymmErrors *gRpAstat = new TGraphAsymmErrors(nbins,xbins,yRpA,exl0,exl0,yRpAStat,yRpAStat);
gRpAstat->SetTitle("RpA stat uncertainty from pPb");
gRpAstat->SetMarkerStyle(21);
gRpAstat->SetMarkerColor(1);
gRpAstat->SetLineColor(1);
gRpAstat->SetLineWidth(2);
gRpAstat->SetFillColor(0);
TGraphAsymmErrors *gRpAsyst = new TGraphAsymmErrors(nbins,xbins,yRpA,exl,exl,yRpPbSystTotLow,yRpPbSystTotHigh);
gRpAsyst->SetTitle("RpA syst uncertainty from pPb");
gRpAsyst->SetName("gRpAsyst");
/*
gRpAsyst->SetFillColor(0);
gRpAsyst->SetMarkerSize(0);
gRpAsyst->SetLineColor(1);
gRpAsyst->SetLineWidth(2);
gRpAsyst->SetFillStyle(0);
*/
gRpAsyst->SetFillColor(kYellow-7);//###0
//gRpAsyst->SetMarkerSize(0);
gRpAsyst->SetLineColor(1);
gRpAsyst->SetLineWidth(2);
gRpAsyst->SetFillStyle(1001);
gRpAsyst->SetMarkerStyle(21);
gRpAsyst->SetMarkerColor(1);
TGraphAsymmErrors*gRpAsyst2=(TGraphAsymmErrors*)gRpAsyst->Clone();
gRpAsyst2->SetMarkerColor(1);
gRpAsyst2->SetMarkerStyle(25);
gRpAsyst2->SetFillColor(0);
gRpAsyst2->SetFillStyle(0);
gRpAsyst2->SetLineColor(1);//5
gRpAsyst2->SetLineStyle(1);
gRpAsyst2->SetLineWidth(2);
TGraphAsymmErrors *gRpAsystFONLL = new TGraphAsymmErrors(nbins,xbins,yFONLL,exl,exl,yRpAsystFONLLlow,yRpAsystFONLLhigh);
gRpAsystFONLL->SetTitle("RpA syst uncertainty from FONLL reference");
/*
gRpAsystFONLL->SetFillColor(kYellow-7);//5
gRpAsystFONLL->SetLineColor(kAzure-3);//5
gRpAsystFONLL->SetMarkerColor(4);//kAzure-3);
TGraphAsymmErrors *gRpAsystFONLL2 = (TGraphAsymmErrors*)gRpAsystFONLL->Clone();
gRpAsystFONLL2->SetFillStyle(0);
gRpAsystFONLL2->SetLineColor(kAzure-3);//5
gRpAsystFONLL2->SetMarkerColor(4);//kAzure-3);
*/
gRpAsystFONLL->SetFillColor(0);//5,kYellow-7
gRpAsystFONLL->SetLineColor(kAzure-3);//5
gRpAsystFONLL->SetMarkerColor(4);//kAzure-3);
gRpAsystFONLL->SetLineStyle(1);
gRpAsystFONLL->SetLineWidth(2);
TGraphAsymmErrors*gRpAsystFONLL2=(TGraphAsymmErrors*)gRpAsystFONLL->Clone();
gRpAsystFONLL2->SetMarkerColor(1);
gRpAsystFONLL2->SetMarkerStyle(25);
gRpAsystFONLL2->SetFillColor(0);
gRpAsystFONLL2->SetFillStyle(0);
gRpAsystFONLL2->SetLineColor(kAzure-3);//5
gRpAsystFONLL2->SetLineStyle(1);
gRpAsystFONLL2->SetLineWidth(2);
canvasRpA->cd(PadNum);
//canvasRpA->Range(-1.989924,-0.2917772,25.49622,2.212202);
canvasRpA->SetFillColor(0);
canvasRpA->SetBorderMode(0);
canvasRpA->SetBorderSize(2);
/*
canvasRpA->SetLeftMargin(0.1451613);
canvasRpA->SetRightMargin(0.05443548);
canvasRpA->SetTopMargin(0.08474576);
canvasRpA->SetBottomMargin(0.1165254);
*/
canvasRpA->SetFrameBorderMode(0);
canvasRpA->SetFrameBorderMode(0);
hempty=new TH2F("hempty","",10,0.1, 62. ,30.,0.,3.0);
hempty->GetXaxis()->CenterTitle();
hempty->GetYaxis()->CenterTitle();
hempty->GetXaxis()->SetTitle("p_{T} (GeV/c)");
hempty->GetYaxis()->SetTitle("R^{FONLL}_{pA}");
//###hempty->GetXaxis()->SetTitleOffset(1.0);//###1.3
hempty->GetYaxis()->SetTitleOffset(1.0);//###1.1
hempty->GetXaxis()->SetTitleSize(0.070*padcorrection);//###0.055
hempty->GetYaxis()->SetTitleSize(0.070*padcorrection);//###0.055
hempty->GetXaxis()->SetTitleFont(42);
hempty->GetYaxis()->SetTitleFont(42);
hempty->GetXaxis()->SetLabelFont(42);
hempty->GetYaxis()->SetLabelFont(42);
hempty->GetXaxis()->SetLabelSize(0.060*padcorrection);//###0.055
hempty->GetYaxis()->SetLabelSize(0.060*padcorrection);//###0.055
if (particle=="Bplus") hempty->GetXaxis()->SetTitleOffset(0.90);//###1.0
else if (particle=="Bzero") hempty->GetXaxis()->SetTitleOffset(0.90);//###0.80
else if (particle=="Bs") hempty->GetXaxis()->SetTitleOffset(0.90);//###0.85
if (particle=="Bplus") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.005
else if (particle=="Bzero") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.0001
else if (particle=="Bs") hempty->GetXaxis()->SetLabelOffset(0.005);//###0.0005
//###hempty->SetMaximum(2);
hempty->SetMaximum(3.0);
hempty->SetMinimum(0.);
hempty->Draw();
TLine *l = new TLine(0,1, 65.,1);
l->SetLineStyle(2);
//TLine *line = new TLine(8.740882,1.017445,61,1.008586);
//###TLine *line = new TLine(4.740882,1.017445,61,1.008586);
TLine *line = new TLine(4.740882,1.00,61,1.00);
line->SetLineColor(kRed);//###
line->SetLineStyle(2);
line->SetLineWidth(2);
/*
gRpAsystFONLL->Draw("2same");
gRpAsystFONLL2->Draw("2same");
line->Draw();
gRpAsyst->Draw("2esame");
gRpAstat->Draw("psame");
*/
gRpAstat->SetMarkerStyle(21);
gRpAstat->SetLineColor(1);
gRpAstat->SetMarkerColor(1);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gRpAstat->GetPoint(i,xgae,ygae);
std::cout << "####### RpA (" << i << ") : " << ygae << " stat: " << gRpAstat->GetEYhigh()[i] << " syst: " << gRpAsyst->GetEYhigh()[i] << std::endl;
std::cout << "########### FONLL (" << i << ") : " << " + " << gRpAsystFONLL->GetEYhigh()[i] << " - " << gRpAsystFONLL->GetEYlow()[i] << std:: endl;
}
gRpAsystFONLL->Draw("2same");
gRpAsystFONLL2->Draw("2same");
//###line->Draw();
gRpAsyst->Draw("2esame");
gRpAsyst2->Draw("2esame");
gRpAstat->Draw("psame");
line->Draw();
TBox *a = new TBox(0.1,1-commonErrorN,4,1+commonErrorP);
a->SetLineColor(1);
a->SetFillColor(0);
a->Draw();
std::cout << "%%%%%%% RpA commonError " << " - " << commonErrorN << " + " << commonErrorP << std::endl;
TBox *b = new TBox(0.1,1-commonErrorN,4,1+commonErrorP);
b->SetLineColor(1);
b->SetFillColor(kGray);
b->Draw();
TBox *b2 = (TBox*)b->Clone();
b2->SetLineColor(1);
b2->SetFillStyle(0);
b2->Draw();
TLegend *legendRpA=new TLegend(0.22,0.60,0.51,0.83,"");
legendRpA->SetBorderSize(0);
legendRpA->SetLineColor(0);
legendRpA->SetFillColor(0);
legendRpA->SetFillStyle(1001);
legendRpA->SetTextFont(42);
legendRpA->SetTextSize(0.060);//###0.045
TLegendEntry *ent_RpAstat=legendRpA->AddEntry(gRpAstat,"R^{FONLL}_{pA}","p");
//###TLegendEntry *ent_RpAstat=legendRpA->AddEntry(gRpAstat,"pPb data","pf");
/*
ent_RpAstat->SetTextFont(42);
ent_RpAstat->SetLineColor(2);
ent_RpAstat->SetMarkerColor(2);
*/
ent_RpAstat->SetTextFont(42);
ent_RpAstat->SetMarkerColor(1);
ent_RpAstat->SetMarkerStyle(25);
//ent_RpAstat->SetFillColor(kYellow-7);//5
//ent_RpAstat->SetFillStyle(1001);
ent_RpAstat->SetLineColor(1);
ent_RpAstat->SetLineWidth(2);
ent_RpAstat->SetLineStyle(0);
TLegendEntry *ent_RpAsyst=legendRpA->AddEntry(gRpAsyst,"Syst. pPb data","f");
ent_RpAsyst->SetTextFont(42);
ent_RpAsyst->SetMarkerColor(1);
ent_RpAsyst->SetMarkerStyle(25);
//ent_RpAsyst->SetFillColor(kYellow-7);//5
//ent_RpAsyst->SetFillStyle(1001);
ent_RpAsyst->SetLineColor(1);
ent_RpAsyst->SetLineWidth(2);
ent_RpAsyst->SetLineStyle(0);
TLegendEntry *ent_RpAsystFONLL=legendRpA->AddEntry(gRpAsystFONLL,"Syst. FONLL pp ref.","f");
/*
ent_RpAsystFONLL->SetTextFont(42);
ent_RpAsystFONLL->SetLineColor(kAzure-3);//5
ent_RpAsystFONLL->SetLineStyle(1);
ent_RpAsystFONLL->SetLineWidth(1);
ent_RpAsystFONLL->SetMarkerColor(kYellow-7);//5
*/
ent_RpAsystFONLL->SetTextFont(42);
ent_RpAsystFONLL->SetLineColor(2);//###5
ent_RpAsystFONLL->SetLineStyle(1);
ent_RpAsystFONLL->SetMarkerColor(kYellow-7);//5
TLegendEntry *ent_RpAsystData=legendRpA->AddEntry(b,"Syst. int. lumi + B","f");
/*
ent_RpAsystData->SetTextFont(42);
ent_RpAsystData->SetLineColor(1);
ent_RpAsystData->SetMarkerColor(2);
*/
ent_RpAsystData->SetTextFont(42);
ent_RpAsystData->SetLineColor(2);
ent_RpAsystData->SetMarkerColor(2);
//if(PadNum==1||PadNum==3){
if(PadNum==3 || PadNum==1 || PadNum==2) {
legendRpA->Draw();
//TLatex * tlatex4=new TLatex(0.21,0.88801268,"CMS");
TLatex * tlatex4=new TLatex(0.19,0.85,"CMS");
tlatex4->SetNDC();
tlatex4->SetTextColor(1);
tlatex4->SetTextFont(62);//42
tlatex4->SetTextSize(0.07);//0.045
tlatex4->Draw();
/*
TLatex * tlatex42=new TLatex(0.35,0.86,"Preliminary");
tlatex42->SetNDC();
tlatex42->SetTextColor(1);
tlatex42->SetTextFont(52);//42
tlatex42->SetTextSize(0.05);//0.045
tlatex42->Draw();
*/
//TLatex * tlatex2=new TLatex(0.471371,0.88801268,"L = 34.8 nb^{-1} (pPb 5.02 TeV)");
//TLatex * tlatex2=new TLatex(0.471371,0.88801268,"34.6 nb^{-1} (pPb 5.02 TeV)");
//###TLatex * tlatex2=new TLatex(0.390,0.950,"34.6 nb^{-1} (pPb 5.02 TeV)");
TLatex * tlatex2;
if (PadNum==1) tlatex2 = new TLatex(0.41,0.94,"34.6 nb^{-1} (pPb 5.02 TeV)");
else tlatex2 = new TLatex(0.39,0.94,"34.6 nb^{-1} (pPb 5.02 TeV)");
tlatex2->SetNDC();
tlatex2->SetTextColor(1);
tlatex2->SetTextFont(42);
tlatex2->SetTextSize(0.06);//0.045
tlatex2->Draw();
}
//if(PadNum==2){
if(PadNum==3 || PadNum==1 || PadNum==2) {
TLatex * tlatex4=new TLatex(0.65,0.23,"|y_{lab}| < 2.4");
tlatex4->SetNDC();
tlatex4->SetTextColor(1);
tlatex4->SetTextFont(42);//42
tlatex4->SetTextSize(0.07);//0.045
tlatex4->Draw();
}
tlatex3->Draw();
TFile *fout=new TFile(Form("../Results%s/fileRpA%s.root",particle.Data(),particle.Data()),"recreate");
fout->cd();
gSigmasyst->SetName("gSigmasyst");
gSigmasyst->Write();
gaeBplusReference->Write();
hSigmapPbStat->Write();
}
void raaExpOpen_pt(const char* inputDir = "../macro_raa/outRoot", // the place where the input root files, with the histograms are
bool bSavePlots = true,
bool bDrawCh = true,
double xMax = 30 // x-axis range limit (ALICE D has x-range maximum value at 36, CMS has it at 30)
)
{
// set the style
gSystem->mkdir(Form("./figs/png"), kTRUE);
gSystem->mkdir(Form("./figs/pdf"), kTRUE);
setTDRStyle();
// read CMS graphs
TFile *pgRaaCms_pt = new TFile(Form("%s/makeRaa_pt.root",inputDir));
// ##################### HIGH PT ############################
TGraphErrors *pgCms = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsi");
TGraphErrors *pgCmsP = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsiP");
TGraphErrors *pgCmsSyst = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsiSyst");
pgCmsSyst->SetFillColorAlpha(kOrange-9,0.5);
pgCms->SetName("gNonPrJpsi");
TBox *lumi = (TBox*)pgRaaCms_pt->Get("lumi");
lumi->SetFillColor(kOrange-9);
lumi->SetFillStyle(1001);
lumi->SetX1(xMax-0.75); lumi->SetX2(xMax);
TBox *lumi_lowpt = (TBox*)lumi->Clone("lumi_lowpt");
lumi_lowpt->SetFillColor(kViolet-9);
lumi_lowpt->SetFillStyle(1001);
lumi_lowpt->SetX1(xMax-1.5); lumi_lowpt->SetX2(xMax-0.75);
// ##################### LOW PT ############################
TGraphErrors *pgCms_lowpt = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsi_pt365y1624");
TGraphErrors *pgCmsP_lowpt = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsiP_pt365y1624");
TGraphErrors *pgCmsSyst_lowpt = (TGraphErrors *)pgRaaCms_pt->Get("gNonPrJpsiSyst_pt365y1624");
pgCmsSyst_lowpt->SetFillColorAlpha(kViolet-9,0.5);
pgCms_lowpt->SetName("gNonPrJpsi_pt365y1624");
//--------------------------------------------------------------------
// *********** CMS Charged particle RAA
TGraphAsymmErrors *p8165_d5x1y1 = new TGraphAsymmErrors(p8165_d5x1y1_numpoints, p8165_d5x1y1_xval, p8165_d5x1y1_yval, p8165_d5x1y1_xerrminus, p8165_d5x1y1_xerrplus, p8165_d5x1y1_ystatminus, p8165_d5x1y1_ystatplus);
TGraphAsymmErrors *p8165_d5x1y1_syst = new TGraphAsymmErrors(p8165_d5x1y1_numpoints, p8165_d5x1y1_xval, p8165_d5x1y1_yval, p8165_d5x1y1_xerrminus, p8165_d5x1y1_xerrplus, p8165_d5x1y1_yerrminus, p8165_d5x1y1_yerrplus);
p8165_d5x1y1->SetName("cms_ch_raa");
// X-axis errors will be discarded
for (int i=0; i<p8165_d5x1y1_numpoints; i++) {
p8165_d5x1y1->SetPointEXlow(i,0);
p8165_d5x1y1->SetPointEXhigh(i,0);
p8165_d5x1y1_syst->SetPointEXlow(i,0.5);
p8165_d5x1y1_syst->SetPointEXhigh(i,0.5);
}
// *********** ALICE D RAA vs. pT
TGraphAsymmErrors *p9059_d15x1y1 = new TGraphAsymmErrors(p9059_d15x1y1_numpoints, p9059_d15x1y1_xval, p9059_d15x1y1_yval, p9059_d15x1y1_xerrminus, p9059_d15x1y1_xerrplus, p9059_d15x1y1_yerrminus, p9059_d15x1y1_yerrplus);
TGraphAsymmErrors *p9059_d15x1y1_syst = new TGraphAsymmErrors(p9059_d15x1y1_numpoints, p9059_d15x1y1_xval, p9059_d15x1y1_yval, p9059_d15x1y1_xerrminus, p9059_d15x1y1_xerrplus, p9059_d15x1y1_yerrminus, p9059_d15x1y1_yerrplus);
p9059_d15x1y1->SetName("alice_d_raa");
// X-axis errors will be discarded
for (int i=0; i<p9059_d15x1y1_numpoints; i++) {
p9059_d15x1y1->SetPointEXlow(i,0);
p9059_d15x1y1->SetPointEXhigh(i,0);
p9059_d15x1y1_syst->SetPointEXlow(i,0.5);
p9059_d15x1y1_syst->SetPointEXhigh(i,0.5);
if (xMax<=30 && (i+1)==p9059_d15x1y1_numpoints) {
double x,y;
p9059_d15x1y1->GetPoint(i,x,y);
p9059_d15x1y1->SetPoint(i,x+5,y);
p9059_d15x1y1_syst->GetPoint(i,x,y);
p9059_d15x1y1_syst->SetPoint(i,x+5,y);
}
}
// Style for graphs
p9059_d15x1y1->SetMarkerStyle(kOpenSquare);
p9059_d15x1y1->SetMarkerSize(1.3);
p9059_d15x1y1->SetMarkerColor(kBlack);
p9059_d15x1y1->SetLineColor(kBlack);
p9059_d15x1y1_syst->SetFillColorAlpha(kGray+1,0.5);
p8165_d5x1y1->SetMarkerStyle(20);
p8165_d5x1y1->SetMarkerSize(1.3);
p8165_d5x1y1->SetMarkerColor(kTeal+3);
p8165_d5x1y1->SetLineColor(kTeal+4);
p8165_d5x1y1_syst->SetFillColorAlpha(kTeal-1,0.5);
//---------------------------------------------------------
TLine *line = new TLine(0.,1,xMax,1);
line->SetLineStyle(1);
line->SetLineWidth(1);
TCanvas *pc = new TCanvas("pc","pc");
TF1 *f4 = new TF1("f4","1",0,xMax);
f4->SetLineWidth(1);
f4->GetXaxis()->SetTitle("p_{T} (GeV/c)");
f4->GetYaxis()->SetTitle("R_{AA}");
f4->GetXaxis()->SetRangeUser(0.0,xMax);
f4->GetYaxis()->SetRangeUser(0.0,1.5);
f4->GetXaxis()->CenterTitle(kTRUE);
f4->Draw();
lumi->Draw();
lumi_lowpt->Draw();
pc->Update();
pgCmsSyst->Draw("2"); // for drawing x-axis
p9059_d15x1y1_syst->Draw("2");
if (bDrawCh) p8165_d5x1y1_syst->Draw("2");
p9059_d15x1y1->Draw("pz");
if (bDrawCh) p8165_d5x1y1->Draw("pz");
pgCmsSyst->Draw("2");
pgCmsP->Draw("P");
pgCms->Draw("P");
pgCmsSyst_lowpt->Draw("2");
pgCmsP_lowpt->Draw("P");
pgCms_lowpt->Draw("P");
// additional info
CMS_lumi(pc,12003000,0);
TLegend *leg_cent = new TLegend(0.28,0.77,0.86,0.89,NULL,"brNDC"); // at top center
leg_cent->SetMargin(0.12);
leg_cent->SetBorderSize(0);
leg_cent->SetTextFont(132);
leg_cent->SetTextSize(0.03);
leg_cent->SetLineColor(1);
leg_cent->SetLineStyle(1);
leg_cent->SetLineWidth(1);
leg_cent->SetFillColor(19);
leg_cent->SetFillStyle(0);
TLegendEntry *entry_cent;
entry_cent=leg_cent->AddEntry("raab","Open beauty: Nonprompt J/#psi","");
entry_cent->SetTextFont(132);
entry_cent->SetTextSize(ltxSetTextSize3);
entry_cent=leg_cent->AddEntry("gNonPrJpsi_pt365y1624", "1.6 < |y| < 2.4","p");
entry_cent->SetTextFont(42);
entry_cent->SetTextSize(entrySize);
entry_cent=leg_cent->AddEntry("gNonPrJpsi", "|y| < 2.4","p");
entry_cent->SetTextFont(42);
entry_cent->SetTextSize(entrySize);
TLegend *leg_ch = new TLegend(0.28,0.685,0.86,0.765,NULL,"brNDC");
leg_ch->SetMargin(0.12);
leg_ch->SetBorderSize(0);
leg_ch->SetTextFont(132);
leg_ch->SetTextSize(0.03);
leg_ch->SetLineColor(1);
leg_ch->SetLineStyle(1);
leg_ch->SetLineWidth(1);
leg_ch->SetFillColor(19);
leg_ch->SetFillStyle(0);
if (bDrawCh) {
TLegendEntry *entry_ch;
entry_ch=leg_ch->AddEntry("cms_ch","Charged hadron","");
entry_ch->SetTextSize(ltxSetTextSize3);
entry_ch->SetTextFont(132);
entry_ch=leg_ch->AddEntry("cms_ch_raa","|#eta| < 1, Cent. 0-5%","p");
entry_ch->SetTextSize(entrySize);
entry_ch->SetTextFont(42);
}
TLegend *leg_alice;
if (bDrawCh) {
leg_alice = new TLegend(0.28,0.60,0.86,0.68,NULL,"brNDC");
} else {
leg_alice = new TLegend(0.28,0.685,0.86,0.765,NULL,"brNDC");
}
leg_alice->SetMargin(0.12);
leg_alice->SetBorderSize(0);
leg_alice->SetTextFont(132);
leg_alice->SetTextSize(0.03);
leg_alice->SetLineColor(1);
leg_alice->SetLineStyle(1);
leg_alice->SetLineWidth(1);
leg_alice->SetFillColor(kWhite);
leg_alice->SetFillStyle(1000);
TLegendEntry *entry_alice;
entry_alice=leg_alice->AddEntry("alice_d","Open charm: Prompt D (ALICE)","");
entry_alice->SetTextSize(ltxSetTextSize3);
entry_alice->SetTextFont(132);
entry_alice=leg_alice->AddEntry("alice_d_raa","|y| < 0.5, Cent. 0-10%","p");
entry_alice->SetTextSize(entrySize);
entry_alice->SetTextFont(42);
TLatex *lat = new TLatex();
lat->SetNDC();
lat->SetTextFont(42);
lat->SetTextSize(ltxSetTextSize2);
if (bDrawCh) lat->DrawLatex(0.63,0.53,"Cent. 0-100%");
else lat->DrawLatex(0.63,0.58,"Cent. 0-100%");
line->Draw();
leg_cent->Draw();
if (bDrawCh) leg_ch->Draw();
leg_alice->Draw();
gPad->RedrawAxis();
if(bSavePlots)
{
pc->SaveAs(Form("figs/pdf/raaExpOpen_pt_RaaCh%d.pdf",bDrawCh));
pc->SaveAs(Form("figs/png/raaExpOpen_pt_RaaCh%d.png",bDrawCh));
}
}
void categories_sig() {
gStyle->SetOptTitle(0);
gStyle->SetOptStat(0);
gStyle->SetCanvasColor(0);
gStyle->SetFrameBorderMode(0);
gStyle->SetPalette(1);
TText *text = new TText();
text->SetNDC();
text->SetTextSize(0.05);
TFile *file = TFile::Open("histograms_CMS-HGG_categories_sig.root");
file->cd();
TCanvas *c_bdtout = new TCanvas("c_bdtout","BDT output");
bdtout_all_tot = (TH1*)bdtout_cat0_tot->Clone();
bdtout_all_tot->Add(bdtout_cat1_tot);
bdtout_all_tot->Add(bdtout_cat2_tot);
bdtout_all_tot->Add(bdtout_cat3_tot);
cout << bdtout_all_tot->Integral() << endl;
cout << bdtout_all_tot->Integral(49,100)+bdtout_all_tot->GetBinContent(48)/2. << endl;
cout << bdtout_all_tot->GetBinLowEdge(49) << endl;
cout << bdtout_all_tot->GetBinLowEdge(48) << endl;
cout << (bdtout_all_tot->Integral(49,100)+bdtout_all_tot->GetBinContent(48)/2.)/bdtout_all_tot->Integral() << endl;
bdtout_passCiC_all_tot = (TH1*)bdtout_passCiC_cat0_tot->Clone();
bdtout_passCiC_all_tot->Add(bdtout_passCiC_cat1_tot);
bdtout_passCiC_all_tot->Add(bdtout_passCiC_cat2_tot);
bdtout_passCiC_all_tot->Add(bdtout_passCiC_cat3_tot);
bdtout_failpresel_cat0_tot->Rebin(2);
bdtout_failpresel_cat1_tot->Rebin(2);
bdtout_failpresel_cat2_tot->Rebin(2);
bdtout_failpresel_cat3_tot->Rebin(2);
bdtout_failpresel_all_tot = (TH1*)bdtout_failpresel_cat0_tot->Clone();
bdtout_failpresel_all_tot->Add(bdtout_failpresel_cat1_tot);
bdtout_failpresel_all_tot->Add(bdtout_failpresel_cat2_tot);
bdtout_failpresel_all_tot->Add(bdtout_failpresel_cat3_tot);
float nFail_m100180_bdtout005 = bdtout_all_tot->Integral(1,52) + bdtout_all_tot->GetBinContent(53)/2.;
float nFail_m100180 = bdtout_all_tot->Integral(1,100);
float frac_bdtout005 = nFail_m100180_bdtout005/nFail_m100180;
float frac_bdtout005_err = frac_err(nFail_m100180_bdtout005,nFail_m100180);
cout << bdtout_all_tot->Integral(1,52) << " " << bdtout_all_tot->GetBinContent(53)/2. << " " << bdtout_all_tot->Integral(1,100) << " " << bdtout_all_tot->GetBinLowEdge(53) << " " << frac_bdtout005 << " ± " << frac_bdtout005_err << endl;
bdtout_all_tot->SetLineColor(1);
bdtout_cat0_tot->SetLineColor(2);
bdtout_cat1_tot->SetLineColor(3);
bdtout_cat2_tot->SetLineColor(4);
bdtout_cat3_tot->SetLineColor(6);
bdtout_all_tot->SetLineWidth(2);
bdtout_cat0_tot->SetLineWidth(2);
bdtout_cat1_tot->SetLineWidth(2);
bdtout_cat2_tot->SetLineWidth(2);
bdtout_cat3_tot->SetLineWidth(2);
bdtout_cat0_tot->SetMarkerColor(2);
bdtout_cat1_tot->SetMarkerColor(3);
bdtout_cat2_tot->SetMarkerColor(4);
bdtout_cat3_tot->SetMarkerColor(6);
bdtout_all_tot->SetMarkerStyle(20);
bdtout_cat0_tot->SetMarkerStyle(20);
bdtout_cat1_tot->SetMarkerStyle(20);
bdtout_cat2_tot->SetMarkerStyle(20);
bdtout_cat3_tot->SetMarkerStyle(20);
bdtout_all_tot->SetMarkerSize(0.7);
bdtout_cat0_tot->SetMarkerSize(0.7);
bdtout_cat1_tot->SetMarkerSize(0.7);
bdtout_cat2_tot->SetMarkerSize(0.7);
bdtout_cat3_tot->SetMarkerSize(0.7);
bdtout_passCiC_all_tot->SetLineColor(1);
bdtout_passCiC_cat0_tot->SetLineColor(2);
bdtout_passCiC_cat1_tot->SetLineColor(3);
bdtout_passCiC_cat2_tot->SetLineColor(4);
bdtout_passCiC_cat3_tot->SetLineColor(6);
bdtout_passCiC_lowPt_cat0_tot->SetLineColor(2);
bdtout_passCiC_lowPt_cat1_tot->SetLineColor(3);
bdtout_passCiC_lowPt_cat2_tot->SetLineColor(4);
bdtout_passCiC_lowPt_cat3_tot->SetLineColor(6);
bdtout_passCiC_highPt_cat0_tot->SetLineColor(2);
bdtout_passCiC_highPt_cat1_tot->SetLineColor(3);
bdtout_passCiC_highPt_cat2_tot->SetLineColor(4);
bdtout_passCiC_highPt_cat3_tot->SetLineColor(6);
bdtout_failpresel_all_tot->SetLineColor(1);
bdtout_failpresel_cat0_tot->SetLineColor(2);
bdtout_failpresel_cat1_tot->SetLineColor(3);
bdtout_failpresel_cat2_tot->SetLineColor(4);
bdtout_failpresel_cat3_tot->SetLineColor(6);
bdtout_all_tot->GetXaxis()->SetTitle("di-photon MVA output");
float boundaries[4] = {-0.05,0.49,0.79,.91};
float max = bdtout_all_tot->GetMaximum();
TBox* box = new TBox(-1.,0.,boundaries[0],max*1.05);
box->SetFillColor(38);
box->SetFillStyle(3002);
bdtout_all_tot->Draw("hist");
box->Draw("hist,same");
bdtout_all_tot->Draw("hist,same");
bdtout_cat0_tot->Draw("hist,same");
bdtout_cat1_tot->Draw("hist,same");
bdtout_cat2_tot->Draw("hist,same");
bdtout_cat3_tot->Draw("hist,same");
TLegend *leg;
leg = new TLegend(.14,.6,.46,.87);
leg->SetBorderSize(0);
leg->SetFillColor(10);
leg->SetTextSize(.035);
leg->AddEntry(bdtout_all_tot,"All");
leg->AddEntry(bdtout_cat0_tot,"both EB, both R9>0.94");
leg->AddEntry(bdtout_cat1_tot,"both EB, !both R9>0.94");
leg->AddEntry(bdtout_cat2_tot,"!both EB, both R9>0.94");
leg->AddEntry(bdtout_cat3_tot,"!both EB, !both R9>0.94");
leg->Draw("hist");
TLine* line[4];
for (int i=0; i<4; i++) {
line[i] = new TLine(boundaries[i],0.,boundaries[i],max*1.05);
line[i]->SetLineColor(4);
line[i]->SetLineWidth(2);
line[i]->SetLineStyle(9);
line[i]->Draw("hist");
}
gPad->RedrawAxis();
c_bdtout->SaveAs("categories_sig.png");
TCanvas *c_bdtout_passCiC = new TCanvas("c_bdtout_passCiC","BDT output, pass CiC");
bdtout_passCiC_cat0_tot_clone = (TH1*)bdtout_passCiC_cat0_tot->Clone();
bdtout_passCiC_cat1_tot_clone = (TH1*)bdtout_passCiC_cat1_tot->Clone();
bdtout_passCiC_cat2_tot_clone = (TH1*)bdtout_passCiC_cat2_tot->Clone();
bdtout_passCiC_cat3_tot_clone = (TH1*)bdtout_passCiC_cat3_tot->Clone();
bdtout_passCiC_all_tot->SetLineWidth(2);
bdtout_passCiC_cat0_tot_clone->SetLineWidth(2);
bdtout_passCiC_cat1_tot_clone->SetLineWidth(2);
bdtout_passCiC_cat2_tot_clone->SetLineWidth(2);
bdtout_passCiC_cat3_tot_clone->SetLineWidth(2);
bdtout_passCiC_all_tot->GetXaxis()->SetTitle("di-photon MVA output");
max = bdtout_passCiC_all_tot->GetMaximum();
TBox* box_passCiC = new TBox(-1.,0.,boundaries[0],max*1.05);
box_passCiC->SetFillColor(38);
box_passCiC->SetFillStyle(3002);
bdtout_passCiC_all_tot->Draw("hist");
box_passCiC->Draw("hist,same");
bdtout_passCiC_all_tot->Draw("hist,same");
bdtout_passCiC_cat0_tot_clone->Draw("hist,same");
bdtout_passCiC_cat1_tot_clone->Draw("hist,same");
bdtout_passCiC_cat2_tot_clone->Draw("hist,same");
bdtout_passCiC_cat3_tot_clone->Draw("hist,same");
leg->Draw("hist");
TLine* line_passCiC[4];
for (int i=0; i<4; i++) {
line_passCiC[i] = new TLine(boundaries[i],0.,boundaries[i],max*1.05);
line_passCiC[i]->SetLineColor(4);
line_passCiC[i]->SetLineWidth(2);
line_passCiC[i]->SetLineStyle(9);
line_passCiC[i]->Draw("hist");
}
gPad->RedrawAxis();
c_bdtout_passCiC->SaveAs("categories_passCiC_sig.png");
/*
TCanvas *c_bdtout_failpresel = new TCanvas("c_bdtout_failpresel","BDT output, pass CiC, fail presel");
bdtout_failpresel_cat0_tot_clone = (TH1*)bdtout_failpresel_cat0_tot->Clone();
bdtout_failpresel_cat1_tot_clone = (TH1*)bdtout_failpresel_cat1_tot->Clone();
bdtout_failpresel_cat2_tot_clone = (TH1*)bdtout_failpresel_cat2_tot->Clone();
bdtout_failpresel_cat3_tot_clone = (TH1*)bdtout_failpresel_cat3_tot->Clone();
bdtout_failpresel_all_tot->SetLineWidth(2);
bdtout_failpresel_cat0_tot_clone->SetLineWidth(2);
bdtout_failpresel_cat1_tot_clone->SetLineWidth(2);
bdtout_failpresel_cat2_tot_clone->SetLineWidth(2);
bdtout_failpresel_cat3_tot_clone->SetLineWidth(2);
bdtout_failpresel_all_tot->GetXaxis()->SetTitle("di-photon MVA output");
max = bdtout_failpresel_all_tot->GetMaximum();
TBox* box_failpresel = new TBox(-1.,0.,boundaries[0],max*1.05);
box_failpresel->SetFillColor(38);
box_failpresel->SetFillStyle(3002);
bdtout_failpresel_all_tot->Draw("hist");
box_failpresel->Draw("hist,same");
bdtout_failpresel_all_tot->Draw("hist,same");
bdtout_failpresel_cat0_tot_clone->Draw("hist,same");
bdtout_failpresel_cat1_tot_clone->Draw("hist,same");
bdtout_failpresel_cat2_tot_clone->Draw("hist,same");
bdtout_failpresel_cat3_tot_clone->Draw("hist,same");
leg->Draw("hist");
TLine* line_failpresel[4];
for (int i=0; i<4; i++) {
line_failpresel[i] = new TLine(boundaries[i],0.,boundaries[i],max*1.05);
line_failpresel[i]->SetLineColor(4);
line_failpresel[i]->SetLineWidth(2);
line_failpresel[i]->SetLineStyle(9);
line_failpresel[i]->Draw("hist");
}
gPad->RedrawAxis();
c_bdtout_failpresel->SaveAs("categories_failpresel_sig.png");
*/
TCanvas *c_bdtout_compareCiC = new TCanvas("c_bdtout_compareCiC","BDT output: pass CiC supertight",1000,650);
c_bdtout_compareCiC->Divide(2,2);
bdtout_passCiC_cat0_tot->SetFillColor(2);
bdtout_passCiC_cat1_tot->SetFillColor(3);
bdtout_passCiC_cat2_tot->SetFillColor(4);
bdtout_passCiC_cat3_tot->SetFillColor(6);
bdtout_passCiC_cat0_tot->SetFillStyle(3002);
bdtout_passCiC_cat1_tot->SetFillStyle(3002);
bdtout_passCiC_cat2_tot->SetFillStyle(3002);
bdtout_passCiC_cat3_tot->SetFillStyle(3002);
bdtout_cat0_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_cat1_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_cat2_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_cat3_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_cat0_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_cat1_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_cat2_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_cat3_tot->GetXaxis()->SetTitleSize(0.05);
c_bdtout_compareCiC->cd(1);
box_passCiC_cat0 = (TBox*)box->Clone();
box_passCiC_cat0->SetY2(bdtout_cat0_tot->GetMaximum()*1.05);
bdtout_cat0_tot->Draw("hist");
box_passCiC_cat0->Draw("hist,same");
bdtout_cat0_tot->Draw("hist,same");
bdtout_passCiC_cat0_tot->Draw("hist,same");
TLine* line_passCiC_cat0[4];
for (int i=0; i<4; i++) {
line_passCiC_cat0[i] = (TLine*)line[i]->Clone();
line_passCiC_cat0[i]->SetY2(bdtout_cat0_tot->GetMaximum()*1.05);
line_passCiC_cat0[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC->cd(2);
box_passCiC_cat1 = (TBox*)box->Clone();
box_passCiC_cat1->SetY2(bdtout_cat1_tot->GetMaximum()*1.05);
bdtout_cat1_tot->Draw("hist");
box_passCiC_cat1->Draw("hist,same");
bdtout_cat1_tot->Draw("hist,same");
bdtout_passCiC_cat1_tot->Draw("hist,same");
TLine* line_passCiC_cat1[4];
for (int i=0; i<4; i++) {
line_passCiC_cat1[i] = (TLine*)line[i]->Clone();
line_passCiC_cat1[i]->SetY2(bdtout_cat1_tot->GetMaximum()*1.05);
line_passCiC_cat1[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC->cd(3);
box_passCiC_cat2 = (TBox*)box->Clone();
box_passCiC_cat2->SetY2(bdtout_cat2_tot->GetMaximum()*1.05);
bdtout_cat2_tot->Draw("hist");
box_passCiC_cat2->Draw("hist,same");
bdtout_cat2_tot->Draw("hist,same");
bdtout_passCiC_cat2_tot->Draw("hist,same");
TLine* line_passCiC_cat2[4];
for (int i=0; i<4; i++) {
line_passCiC_cat2[i] = (TLine*)line[i]->Clone();
line_passCiC_cat2[i]->SetY2(bdtout_cat2_tot->GetMaximum()*1.05);
line_passCiC_cat2[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC->cd(4);
box_passCiC_cat3 = (TBox*)box->Clone();
box_passCiC_cat3->SetY2(bdtout_cat3_tot->GetMaximum()*1.05);
bdtout_cat3_tot->Draw("hist");
box_passCiC_cat3->Draw("hist,same");
bdtout_cat3_tot->Draw("hist,same");
bdtout_passCiC_cat3_tot->Draw("hist,same");
TLine* line_passCiC_cat3[4];
for (int i=0; i<4; i++) {
line_passCiC_cat3[i] = (TLine*)line[i]->Clone();
line_passCiC_cat3[i]->SetY2(bdtout_cat3_tot->GetMaximum()*1.05);
line_passCiC_cat3[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC->SaveAs("categories_compareCiC_sig.png");
TCanvas *c_bdtout_lowPt = new TCanvas("c_bdtout_lowPt","BDT output: di-photon pT<40 GeV");
bdtout_lowPt_all_tot = (TH1*)bdtout_lowPt_cat0_tot->Clone();
bdtout_lowPt_all_tot->Add(bdtout_lowPt_cat1_tot);
bdtout_lowPt_all_tot->Add(bdtout_lowPt_cat2_tot);
bdtout_lowPt_all_tot->Add(bdtout_lowPt_cat3_tot);
bdtout_lowPt_all_tot->SetLineColor(1);
bdtout_lowPt_cat0_tot->SetLineColor(2);
bdtout_lowPt_cat1_tot->SetLineColor(3);
bdtout_lowPt_cat2_tot->SetLineColor(4);
bdtout_lowPt_cat3_tot->SetLineColor(6);
bdtout_lowPt_all_tot->SetLineWidth(2);
bdtout_lowPt_cat0_tot->SetLineWidth(2);
bdtout_lowPt_cat1_tot->SetLineWidth(2);
bdtout_lowPt_cat2_tot->SetLineWidth(2);
bdtout_lowPt_cat3_tot->SetLineWidth(2);
bdtout_lowPt_all_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_lowPt_all_tot->GetYaxis()->SetTitle("Events/0.02");
bdtout_lowPt_all_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_lowPt_all_tot->GetXaxis()->SetLabelSize(0.05);
bdtout_lowPt_all_tot->GetYaxis()->SetTitleSize(0.05);
bdtout_lowPt_all_tot->GetYaxis()->SetLabelSize(0.05);
max = bdtout_lowPt_all_tot->GetMaximum();
TBox* box_lowPt = new TBox(-1.,0.,boundaries[0],max*1.05);
box_lowPt->SetFillColor(38);
box_lowPt->SetFillStyle(3002);
bdtout_lowPt_all_tot->Draw("hist");
box_lowPt->Draw("hist,same");
bdtout_lowPt_all_tot->Draw("hist,same");
bdtout_lowPt_cat0_tot->Draw("hist,same");
bdtout_lowPt_cat1_tot->Draw("hist,same");
bdtout_lowPt_cat2_tot->Draw("hist,same");
bdtout_lowPt_cat3_tot->Draw("hist,same");
TLegend *leg_lowPt = (TLegend*)leg->Clone();
leg_lowPt->Clear();
leg_lowPt->AddEntry(bdtout_all_tot,"All (p_{T}(#gamma#gamma) < 40 GeV)");
leg_lowPt->AddEntry(bdtout_cat0_tot,"both EB, both R9>0.94");
leg_lowPt->AddEntry(bdtout_cat1_tot,"both EB, !both R9>0.94");
leg_lowPt->AddEntry(bdtout_cat2_tot,"!both EB, both R9>0.94");
leg_lowPt->AddEntry(bdtout_cat3_tot,"!both EB, !both R9>0.94");
leg_lowPt->Draw("hist");
TLine* line_lowPt[4];
for (int i=0; i<4; i++) {
line_lowPt[i] = new TLine(boundaries[i],0.,boundaries[i],max*1.05);
line_lowPt[i]->SetLineColor(4);
line_lowPt[i]->SetLineWidth(2);
line_lowPt[i]->SetLineStyle(9);
line_lowPt[i]->Draw("hist");
}
gPad->RedrawAxis();
c_bdtout_lowPt->SaveAs("categories_lowPt_sig.png");
TCanvas *c_bdtout_compareCiC_lowPt = new TCanvas("c_bdtout_compareCiC_lowPt","BDT output: di-photon pT<40 GeV, pass CiC supertight",1000,650);
c_bdtout_compareCiC_lowPt->Divide(2,2);
bdtout_passCiC_lowPt_cat0_tot->SetFillColor(2);
bdtout_passCiC_lowPt_cat1_tot->SetFillColor(3);
bdtout_passCiC_lowPt_cat2_tot->SetFillColor(4);
bdtout_passCiC_lowPt_cat3_tot->SetFillColor(6);
bdtout_passCiC_lowPt_cat0_tot->SetFillStyle(3002);
bdtout_passCiC_lowPt_cat1_tot->SetFillStyle(3002);
bdtout_passCiC_lowPt_cat2_tot->SetFillStyle(3002);
bdtout_passCiC_lowPt_cat3_tot->SetFillStyle(3002);
bdtout_lowPt_cat0_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_lowPt_cat1_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_lowPt_cat2_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_lowPt_cat3_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_lowPt_cat0_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_lowPt_cat1_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_lowPt_cat2_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_lowPt_cat3_tot->GetXaxis()->SetTitleSize(0.05);
c_bdtout_compareCiC_lowPt->cd(1);
box_passCiC_lowPt_cat0 = (TBox*)box->Clone();
box_passCiC_lowPt_cat0->SetY2(bdtout_lowPt_cat0_tot->GetMaximum()*1.05);
bdtout_lowPt_cat0_tot->Draw("hist");
box_passCiC_lowPt_cat0->Draw("hist,same");
bdtout_lowPt_cat0_tot->Draw("hist,same");
bdtout_passCiC_lowPt_cat0_tot->Draw("hist,same");
TLine* line_passCiC_lowPt_cat0[4];
for (int i=0; i<4; i++) {
line_passCiC_lowPt_cat0[i] = (TLine*)line[i]->Clone();
line_passCiC_lowPt_cat0[i]->SetY2(bdtout_lowPt_cat0_tot->GetMaximum()*1.05);
line_passCiC_lowPt_cat0[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_lowPt->cd(2);
box_passCiC_lowPt_cat1 = (TBox*)box->Clone();
box_passCiC_lowPt_cat1->SetY2(bdtout_lowPt_cat1_tot->GetMaximum()*1.05);
bdtout_lowPt_cat1_tot->Draw("hist");
box_passCiC_lowPt_cat1->Draw("hist,same");
bdtout_lowPt_cat1_tot->Draw("hist,same");
bdtout_passCiC_lowPt_cat1_tot->Draw("hist,same");
TLine* line_passCiC_lowPt_cat1[4];
for (int i=0; i<4; i++) {
line_passCiC_lowPt_cat1[i] = (TLine*)line[i]->Clone();
line_passCiC_lowPt_cat1[i]->SetY2(bdtout_lowPt_cat1_tot->GetMaximum()*1.05);
line_passCiC_lowPt_cat1[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_lowPt->cd(3);
box_passCiC_lowPt_cat2 = (TBox*)box->Clone();
box_passCiC_lowPt_cat2->SetY2(bdtout_lowPt_cat2_tot->GetMaximum()*1.05);
bdtout_lowPt_cat2_tot->Draw("hist");
box_passCiC_lowPt_cat2->Draw("hist,same");
bdtout_lowPt_cat2_tot->Draw("hist,same");
bdtout_passCiC_lowPt_cat2_tot->Draw("hist,same");
TLine* line_passCiC_lowPt_cat2[4];
for (int i=0; i<4; i++) {
line_passCiC_lowPt_cat2[i] = (TLine*)line[i]->Clone();
line_passCiC_lowPt_cat2[i]->SetY2(bdtout_lowPt_cat2_tot->GetMaximum()*1.05);
line_passCiC_lowPt_cat2[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_lowPt->cd(4);
box_passCiC_lowPt_cat3 = (TBox*)box->Clone();
box_passCiC_lowPt_cat3->SetY2(bdtout_lowPt_cat3_tot->GetMaximum()*1.05);
bdtout_lowPt_cat3_tot->Draw("hist");
box_passCiC_lowPt_cat3->Draw("hist,same");
bdtout_lowPt_cat3_tot->Draw("hist,same");
bdtout_passCiC_lowPt_cat3_tot->Draw("hist,same");
TLine* line_passCiC_lowPt_cat3[4];
for (int i=0; i<4; i++) {
line_passCiC_lowPt_cat3[i] = (TLine*)line[i]->Clone();
line_passCiC_lowPt_cat3[i]->SetY2(bdtout_lowPt_cat3_tot->GetMaximum()*1.05);
line_passCiC_lowPt_cat3[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_lowPt->SaveAs("categories_compareCiC_lowPt_sig.png");
TCanvas *c_bdtout_highPt = new TCanvas("c_bdtout_highPt","BDT output: di-photon pT>40 GeV");
bdtout_highPt_all_tot = (TH1*)bdtout_highPt_cat0_tot->Clone();
bdtout_highPt_all_tot->Add(bdtout_highPt_cat1_tot);
bdtout_highPt_all_tot->Add(bdtout_highPt_cat2_tot);
bdtout_highPt_all_tot->Add(bdtout_highPt_cat3_tot);
bdtout_highPt_all_tot->SetLineColor(1);
bdtout_highPt_cat0_tot->SetLineColor(2);
bdtout_highPt_cat1_tot->SetLineColor(3);
bdtout_highPt_cat2_tot->SetLineColor(4);
bdtout_highPt_cat3_tot->SetLineColor(6);
bdtout_highPt_all_tot->SetLineWidth(2);
bdtout_highPt_cat0_tot->SetLineWidth(2);
bdtout_highPt_cat1_tot->SetLineWidth(2);
bdtout_highPt_cat2_tot->SetLineWidth(2);
bdtout_highPt_cat3_tot->SetLineWidth(2);
bdtout_highPt_all_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_highPt_all_tot->GetYaxis()->SetTitle("Events/0.02");
bdtout_highPt_all_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_highPt_all_tot->GetXaxis()->SetLabelSize(0.05);
bdtout_highPt_all_tot->GetYaxis()->SetTitleSize(0.05);
bdtout_highPt_all_tot->GetYaxis()->SetLabelSize(0.05);
max = 1.2*bdtout_highPt_all_tot->GetMaximum();
bdtout_highPt_all_tot->SetMaximum(max);
TBox* box_highPt = new TBox(-1.,0.,boundaries[0],max);
box_highPt->SetFillColor(38);
box_highPt->SetFillStyle(3002);
bdtout_highPt_all_tot->Draw("hist");
box_highPt->Draw("hist,same");
bdtout_highPt_all_tot->Draw("hist,same");
bdtout_highPt_cat0_tot->Draw("hist,same");
bdtout_highPt_cat1_tot->Draw("hist,same");
bdtout_highPt_cat2_tot->Draw("hist,same");
bdtout_highPt_cat3_tot->Draw("hist,same");
TLegend *leg_highPt = (TLegend*)leg->Clone();
leg_highPt->Clear();
leg_highPt->AddEntry(bdtout_all_tot,"All (p_{T}(#gamma#gamma) > 40 GeV)");
leg_highPt->AddEntry(bdtout_cat0_tot,"both EB, both R9>0.94");
leg_highPt->AddEntry(bdtout_cat1_tot,"both EB, !both R9>0.94");
leg_highPt->AddEntry(bdtout_cat2_tot,"!both EB, both R9>0.94");
leg_highPt->AddEntry(bdtout_cat3_tot,"!both EB, !both R9>0.94");
leg_highPt->Draw("hist");
TLine* line_highPt[4];
for (int i=0; i<4; i++) {
line_highPt[i] = new TLine(boundaries[i],0.,boundaries[i],max);
line_highPt[i]->SetLineColor(4);
line_highPt[i]->SetLineWidth(2);
line_highPt[i]->SetLineStyle(9);
line_highPt[i]->Draw("hist");
}
gPad->RedrawAxis();
c_bdtout_highPt->SaveAs("categories_highPt_sig.png");
TCanvas *c_bdtout_compareCiC_highPt = new TCanvas("c_bdtout_compareCiC_highPt","BDT output: di-photon pT>40 GeV, pass CiC supertight",1000,650);
c_bdtout_compareCiC_highPt->Divide(2,2);
bdtout_passCiC_highPt_cat0_tot->SetFillColor(2);
bdtout_passCiC_highPt_cat1_tot->SetFillColor(3);
bdtout_passCiC_highPt_cat2_tot->SetFillColor(4);
bdtout_passCiC_highPt_cat3_tot->SetFillColor(6);
bdtout_passCiC_highPt_cat0_tot->SetFillStyle(3002);
bdtout_passCiC_highPt_cat1_tot->SetFillStyle(3002);
bdtout_passCiC_highPt_cat2_tot->SetFillStyle(3002);
bdtout_passCiC_highPt_cat3_tot->SetFillStyle(3002);
bdtout_highPt_cat0_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_highPt_cat1_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_highPt_cat2_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_highPt_cat3_tot->GetXaxis()->SetTitle("di-photon MVA output");
bdtout_highPt_cat0_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_highPt_cat1_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_highPt_cat2_tot->GetXaxis()->SetTitleSize(0.05);
bdtout_highPt_cat3_tot->GetXaxis()->SetTitleSize(0.05);
c_bdtout_compareCiC_highPt->cd(1);
box_passCiC_highPt_cat0 = (TBox*)box->Clone();
box_passCiC_highPt_cat0->SetY2(bdtout_highPt_cat0_tot->GetMaximum()*1.05);
bdtout_highPt_cat0_tot->Draw("hist");
box_passCiC_highPt_cat0->Draw("hist,same");
bdtout_highPt_cat0_tot->Draw("hist,same");
bdtout_passCiC_highPt_cat0_tot->Draw("hist,same");
TLine* line_passCiC_highPt_cat0[4];
for (int i=0; i<4; i++) {
line_passCiC_highPt_cat0[i] = (TLine*)line[i]->Clone();
line_passCiC_highPt_cat0[i]->SetY2(bdtout_highPt_cat0_tot->GetMaximum()*1.05);
line_passCiC_highPt_cat0[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_highPt->cd(2);
box_passCiC_highPt_cat1 = (TBox*)box->Clone();
box_passCiC_highPt_cat1->SetY2(bdtout_highPt_cat1_tot->GetMaximum()*1.05);
bdtout_highPt_cat1_tot->Draw("hist");
box_passCiC_highPt_cat1->Draw("hist,same");
bdtout_highPt_cat1_tot->Draw("hist,same");
bdtout_passCiC_highPt_cat1_tot->Draw("hist,same");
TLine* line_passCiC_highPt_cat1[4];
for (int i=0; i<4; i++) {
line_passCiC_highPt_cat1[i] = (TLine*)line[i]->Clone();
line_passCiC_highPt_cat1[i]->SetY2(bdtout_highPt_cat1_tot->GetMaximum()*1.05);
line_passCiC_highPt_cat1[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_highPt->cd(3);
box_passCiC_highPt_cat2 = (TBox*)box->Clone();
box_passCiC_highPt_cat2->SetY2(bdtout_highPt_cat2_tot->GetMaximum()*1.05);
bdtout_highPt_cat2_tot->Draw("hist");
box_passCiC_highPt_cat2->Draw("hist,same");
bdtout_highPt_cat2_tot->Draw("hist,same");
bdtout_passCiC_highPt_cat2_tot->Draw("hist,same");
TLine* line_passCiC_highPt_cat2[4];
for (int i=0; i<4; i++) {
line_passCiC_highPt_cat2[i] = (TLine*)line[i]->Clone();
line_passCiC_highPt_cat2[i]->SetY2(bdtout_highPt_cat2_tot->GetMaximum()*1.05);
line_passCiC_highPt_cat2[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_highPt->cd(4);
box_passCiC_highPt_cat3 = (TBox*)box->Clone();
box_passCiC_highPt_cat3->SetY2(bdtout_highPt_cat3_tot->GetMaximum()*1.05);
bdtout_highPt_cat3_tot->Draw("hist");
box_passCiC_highPt_cat3->Draw("hist,same");
bdtout_highPt_cat3_tot->Draw("hist,same");
bdtout_passCiC_highPt_cat3_tot->Draw("hist,same");
TLine* line_passCiC_highPt_cat3[4];
for (int i=0; i<4; i++) {
line_passCiC_highPt_cat3[i] = (TLine*)line[i]->Clone();
line_passCiC_highPt_cat3[i]->SetY2(bdtout_highPt_cat3_tot->GetMaximum()*1.05);
line_passCiC_highPt_cat3[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"!both EB, !both R9>0.94");
gPad->RedrawAxis();
c_bdtout_compareCiC_highPt->SaveAs("categories_compareCiC_highPt_sig.png");
TCanvas *c_ptVsBdtout = new TCanvas("c_ptVsBdtout","pT(#gamma#gamma) vs MVA output",1000,650);
c_ptVsBdtout->Divide(2,2);
pt_vs_bdtout_cat0_all = (TH2*)pt_vs_bdtout_cat0_ggh_m125_8TeV->Clone();
pt_vs_bdtout_cat1_all = (TH2*)pt_vs_bdtout_cat1_ggh_m125_8TeV->Clone();
pt_vs_bdtout_cat2_all = (TH2*)pt_vs_bdtout_cat2_ggh_m125_8TeV->Clone();
pt_vs_bdtout_cat3_all = (TH2*)pt_vs_bdtout_cat3_ggh_m125_8TeV->Clone();
pt_vs_bdtout_cat4_all = (TH2*)pt_vs_bdtout_cat4_ggh_m125_8TeV->Clone();
pt_vs_bdtout_cat0_all->Add(pt_vs_bdtout_cat0_vbf_m125_8TeV);
pt_vs_bdtout_cat1_all->Add(pt_vs_bdtout_cat1_vbf_m125_8TeV);
pt_vs_bdtout_cat2_all->Add(pt_vs_bdtout_cat2_vbf_m125_8TeV);
pt_vs_bdtout_cat3_all->Add(pt_vs_bdtout_cat3_vbf_m125_8TeV);
pt_vs_bdtout_cat4_all->Add(pt_vs_bdtout_cat4_vbf_m125_8TeV);
pt_vs_bdtout_cat0_all->Add(pt_vs_bdtout_cat0_wzh_m125_8TeV);
pt_vs_bdtout_cat1_all->Add(pt_vs_bdtout_cat1_wzh_m125_8TeV);
pt_vs_bdtout_cat2_all->Add(pt_vs_bdtout_cat2_wzh_m125_8TeV);
pt_vs_bdtout_cat3_all->Add(pt_vs_bdtout_cat3_wzh_m125_8TeV);
pt_vs_bdtout_cat4_all->Add(pt_vs_bdtout_cat4_wzh_m125_8TeV);
pt_vs_bdtout_cat0_all->Add(pt_vs_bdtout_cat0_tth_m125_8TeV);
pt_vs_bdtout_cat1_all->Add(pt_vs_bdtout_cat1_tth_m125_8TeV);
pt_vs_bdtout_cat2_all->Add(pt_vs_bdtout_cat2_tth_m125_8TeV);
pt_vs_bdtout_cat3_all->Add(pt_vs_bdtout_cat3_tth_m125_8TeV);
pt_vs_bdtout_cat4_all->Add(pt_vs_bdtout_cat4_tth_m125_8TeV);
pt_vs_bdtout_cat0_all->GetXaxis()->SetTitle("di-photon MVA output");
pt_vs_bdtout_cat1_all->GetXaxis()->SetTitle("di-photon MVA output");
pt_vs_bdtout_cat2_all->GetXaxis()->SetTitle("di-photon MVA output");
pt_vs_bdtout_cat3_all->GetXaxis()->SetTitle("di-photon MVA output");
pt_vs_bdtout_cat0_all->GetXaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat1_all->GetXaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat2_all->GetXaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat3_all->GetXaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat0_all->GetYaxis()->SetTitle("di-photon p_{T} (GeV)");
pt_vs_bdtout_cat1_all->GetYaxis()->SetTitle("di-photon p_{T} (GeV)");
pt_vs_bdtout_cat2_all->GetYaxis()->SetTitle("di-photon p_{T} (GeV)");
pt_vs_bdtout_cat3_all->GetYaxis()->SetTitle("di-photon p_{T} (GeV)");
pt_vs_bdtout_cat0_all->GetYaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat1_all->GetYaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat2_all->GetYaxis()->SetTitleSize(0.05);
pt_vs_bdtout_cat3_all->GetYaxis()->SetTitleSize(0.05);
c_ptVsBdtout->cd(1);
pt_vs_bdtout_cat0_all->Draw("colz");
box_ptVsBdtout = (TBox*)box->Clone();
box_ptVsBdtout->SetY2(200.);
box_ptVsBdtout->Draw("hist");
pt_vs_bdtout_cat0_all->Draw("colz,same");
TLine* line_ptVsBdtout[4];
for (int i=0; i<4; i++) {
line_ptVsBdtout[i] = (TLine*)line[i]->Clone();
line_ptVsBdtout[i]->SetY2(200.);
line_ptVsBdtout[i]->Draw("hist");
}
text->DrawText(0.15,0.75,"both EB, both R9>0.94");
gPad->RedrawAxis();
c_ptVsBdtout->cd(2);
pt_vs_bdtout_cat1_all->Draw("colz");
box_ptVsBdtout->Draw("hist");
pt_vs_bdtout_cat1_all->Draw("colz,same");
for (int i=0; i<4; i++) line_ptVsBdtout[i]->Draw("hist");
text->DrawText(0.15,0.75,"both EB, !both R9>0.94");
gPad->RedrawAxis();
c_ptVsBdtout->cd(3);
pt_vs_bdtout_cat2_all->Draw("colz");
box_ptVsBdtout->Draw("hist");
pt_vs_bdtout_cat2_all->Draw("colz,same");
for (int i=0; i<4; i++) line_ptVsBdtout[i]->Draw("hist");
text->DrawText(0.15,0.75,"!both EB, both R9>0.94");
gPad->RedrawAxis();
c_ptVsBdtout->cd(4);
pt_vs_bdtout_cat3_all->Draw("colz");
box_ptVsBdtout->Draw("hist");
pt_vs_bdtout_cat3_all->Draw("colz,same");
for (int i=0; i<4; i++) line_ptVsBdtout[i]->Draw("hist");
text->DrawText(0.15,0.75,"!both EB, !both R9>0.94");
gPad->RedrawAxis();
c_ptVsBdtout->SaveAs("ptVsBdtout_sig.png");
TCanvas *c_2D = new TCanvas("c_2D","min(R9) vs max(eta), sublead eta vs lead eta",1250,500);
c_2D->Divide(4,2);
minR9_vs_maxEta_cat2_all = (TH2*)minR9_vs_maxEta_cat2_ggh_m125_8TeV->Clone();
minR9_vs_maxEta_cat3_all = (TH2*)minR9_vs_maxEta_cat3_ggh_m125_8TeV->Clone();
minR9_vs_maxEta_cat4_all = (TH2*)minR9_vs_maxEta_cat4_ggh_m125_8TeV->Clone();
minR9_vs_maxEta_cat5_all = (TH2*)minR9_vs_maxEta_cat5_ggh_m125_8TeV->Clone();
minR9_vs_maxEta_cat6_all = (TH2*)minR9_vs_maxEta_cat6_ggh_m125_8TeV->Clone();
minR9_vs_maxEta_cat2_all->Add(minR9_vs_maxEta_cat2_vbf_m125_8TeV);
minR9_vs_maxEta_cat3_all->Add(minR9_vs_maxEta_cat3_vbf_m125_8TeV);
minR9_vs_maxEta_cat4_all->Add(minR9_vs_maxEta_cat4_vbf_m125_8TeV);
minR9_vs_maxEta_cat5_all->Add(minR9_vs_maxEta_cat5_vbf_m125_8TeV);
minR9_vs_maxEta_cat6_all->Add(minR9_vs_maxEta_cat6_vbf_m125_8TeV);
minR9_vs_maxEta_cat2_all->Add(minR9_vs_maxEta_cat2_wzh_m125_8TeV);
minR9_vs_maxEta_cat3_all->Add(minR9_vs_maxEta_cat3_wzh_m125_8TeV);
minR9_vs_maxEta_cat4_all->Add(minR9_vs_maxEta_cat4_wzh_m125_8TeV);
minR9_vs_maxEta_cat5_all->Add(minR9_vs_maxEta_cat5_wzh_m125_8TeV);
minR9_vs_maxEta_cat6_all->Add(minR9_vs_maxEta_cat6_wzh_m125_8TeV);
minR9_vs_maxEta_cat2_all->Add(minR9_vs_maxEta_cat2_tth_m125_8TeV);
minR9_vs_maxEta_cat3_all->Add(minR9_vs_maxEta_cat3_tth_m125_8TeV);
minR9_vs_maxEta_cat4_all->Add(minR9_vs_maxEta_cat4_tth_m125_8TeV);
minR9_vs_maxEta_cat5_all->Add(minR9_vs_maxEta_cat5_tth_m125_8TeV);
minR9_vs_maxEta_cat6_all->Add(minR9_vs_maxEta_cat6_tth_m125_8TeV);
eta2_vs_eta1_cat2_all = (TH2*)eta2_vs_eta1_cat2_ggh_m125_8TeV->Clone();
eta2_vs_eta1_cat3_all = (TH2*)eta2_vs_eta1_cat3_ggh_m125_8TeV->Clone();
eta2_vs_eta1_cat4_all = (TH2*)eta2_vs_eta1_cat4_ggh_m125_8TeV->Clone();
eta2_vs_eta1_cat5_all = (TH2*)eta2_vs_eta1_cat5_ggh_m125_8TeV->Clone();
eta2_vs_eta1_cat6_all = (TH2*)eta2_vs_eta1_cat6_ggh_m125_8TeV->Clone();
eta2_vs_eta1_cat2_all->Add(eta2_vs_eta1_cat2_vbf_m125_8TeV);
eta2_vs_eta1_cat3_all->Add(eta2_vs_eta1_cat3_vbf_m125_8TeV);
eta2_vs_eta1_cat4_all->Add(eta2_vs_eta1_cat4_vbf_m125_8TeV);
eta2_vs_eta1_cat5_all->Add(eta2_vs_eta1_cat5_vbf_m125_8TeV);
eta2_vs_eta1_cat6_all->Add(eta2_vs_eta1_cat6_vbf_m125_8TeV);
eta2_vs_eta1_cat2_all->Add(eta2_vs_eta1_cat2_wzh_m125_8TeV);
eta2_vs_eta1_cat3_all->Add(eta2_vs_eta1_cat3_wzh_m125_8TeV);
eta2_vs_eta1_cat4_all->Add(eta2_vs_eta1_cat4_wzh_m125_8TeV);
eta2_vs_eta1_cat5_all->Add(eta2_vs_eta1_cat5_wzh_m125_8TeV);
eta2_vs_eta1_cat6_all->Add(eta2_vs_eta1_cat6_wzh_m125_8TeV);
eta2_vs_eta1_cat2_all->Add(eta2_vs_eta1_cat2_tth_m125_8TeV);
eta2_vs_eta1_cat3_all->Add(eta2_vs_eta1_cat3_tth_m125_8TeV);
eta2_vs_eta1_cat4_all->Add(eta2_vs_eta1_cat4_tth_m125_8TeV);
eta2_vs_eta1_cat5_all->Add(eta2_vs_eta1_cat5_tth_m125_8TeV);
eta2_vs_eta1_cat6_all->Add(eta2_vs_eta1_cat6_tth_m125_8TeV);
minR9_vs_maxEta_cat2_all->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat3_all->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat4_all->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat5_all->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat6_all->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat2_all->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat3_all->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat4_all->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat5_all->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat6_all->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat2_all->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat3_all->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat4_all->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat5_all->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat6_all->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat2_all->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat3_all->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat4_all->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat5_all->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat6_all->GetYaxis()->SetTitleSize(0.05);
line_minr9 = new TLine(0.,0.94,2.5,0.94);
line_maxeta = new TLine(1.479,0.5,1.479,1.);
//line_minr9->SetLineColor(4);
line_minr9->SetLineWidth(2);
line_minr9->SetLineStyle(9);
//line_maxeta->SetLineColor(4);
line_maxeta->SetLineWidth(2);
line_maxeta->SetLineStyle(9);
c_2D->cd(1);
minR9_vs_maxEta_cat2_all->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D->cd(2);
minR9_vs_maxEta_cat3_all->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D->cd(3);
minR9_vs_maxEta_cat4_all->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D->cd(4);
minR9_vs_maxEta_cat5_all->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
//c_2D->cd(5);
//minR9_vs_maxEta_cat6_all->Draw("colz");
//line_minr9->Draw("hist");
//line_maxeta->Draw("hist");
eta2_vs_eta1_cat2_all->GetYaxis()->SetTitle("sublead #eta");
eta2_vs_eta1_cat3_all->GetYaxis()->SetTitle("sublead #eta");
eta2_vs_eta1_cat4_all->GetYaxis()->SetTitle("sublead #eta");
eta2_vs_eta1_cat5_all->GetYaxis()->SetTitle("sublead #eta");
eta2_vs_eta1_cat6_all->GetYaxis()->SetTitle("sublead #eta");
eta2_vs_eta1_cat2_all->GetXaxis()->SetTitle("lead #eta");
eta2_vs_eta1_cat3_all->GetXaxis()->SetTitle("lead #eta");
eta2_vs_eta1_cat4_all->GetXaxis()->SetTitle("lead #eta");
eta2_vs_eta1_cat5_all->GetXaxis()->SetTitle("lead #eta");
eta2_vs_eta1_cat6_all->GetXaxis()->SetTitle("lead #eta");
eta2_vs_eta1_cat2_all->GetXaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat3_all->GetXaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat4_all->GetXaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat5_all->GetXaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat6_all->GetXaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat2_all->GetYaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat3_all->GetYaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat4_all->GetYaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat5_all->GetYaxis()->SetTitleSize(0.05);
eta2_vs_eta1_cat6_all->GetYaxis()->SetTitleSize(0.05);
c_2D->cd(5);
eta2_vs_eta1_cat2_all->Draw("colz");
c_2D->cd(6);
eta2_vs_eta1_cat3_all->Draw("colz");
c_2D->cd(7);
eta2_vs_eta1_cat4_all->Draw("colz");
c_2D->cd(8);
eta2_vs_eta1_cat5_all->Draw("colz");
//c_2D->cd(10);
//eta2_vs_eta1_cat6_all->Draw("colz");
c_2D->SaveAs("categories_2D_sig.png");
/*
TCanvas *c_2D_v2 = new TCanvas("c_2D_v2","min(R9) vs max(eta), |lead eta| vs |deltaEeta|",1500,500);
c_2D_v2->Divide(5,2);
minR9_vs_maxEta_cat2_tot->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat3_tot->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat4_tot->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat5_tot->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat6_tot->GetXaxis()->SetTitle("max(#eta)");
minR9_vs_maxEta_cat2_tot->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat3_tot->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat4_tot->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat5_tot->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat6_tot->GetYaxis()->SetTitle("min(R9)");
minR9_vs_maxEta_cat2_tot->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat3_tot->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat4_tot->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat5_tot->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat6_tot->GetXaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat2_tot->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat3_tot->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat4_tot->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat5_tot->GetYaxis()->SetTitleSize(0.05);
minR9_vs_maxEta_cat6_tot->GetYaxis()->SetTitleSize(0.05);
c_2D_v2->cd(1);
minR9_vs_maxEta_cat2_tot->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D_v2->cd(2);
minR9_vs_maxEta_cat3_tot->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D_v2->cd(3);
minR9_vs_maxEta_cat4_tot->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D_v2->cd(4);
minR9_vs_maxEta_cat5_tot->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
c_2D_v2->cd(5);
minR9_vs_maxEta_cat6_tot->Draw("colz");
line_minr9->Draw("hist");
line_maxeta->Draw("hist");
eta1_vs_deltaEta_cat2_tot->GetYaxis()->SetTitle("|lead #eta|");
eta1_vs_deltaEta_cat3_tot->GetYaxis()->SetTitle("|lead #eta|");
eta1_vs_deltaEta_cat4_tot->GetYaxis()->SetTitle("|lead #eta|");
eta1_vs_deltaEta_cat5_tot->GetYaxis()->SetTitle("|lead #eta|");
eta1_vs_deltaEta_cat6_tot->GetYaxis()->SetTitle("|lead #eta|");
eta1_vs_deltaEta_cat2_tot->GetXaxis()->SetTitle("|#Delta#eta|");
eta1_vs_deltaEta_cat3_tot->GetXaxis()->SetTitle("|#Delta#eta|");
eta1_vs_deltaEta_cat4_tot->GetXaxis()->SetTitle("|#Delta#eta|");
eta1_vs_deltaEta_cat5_tot->GetXaxis()->SetTitle("|#Delta#eta|");
eta1_vs_deltaEta_cat6_tot->GetXaxis()->SetTitle("|#Delta#eta|");
eta1_vs_deltaEta_cat2_tot->GetXaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat3_tot->GetXaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat4_tot->GetXaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat5_tot->GetXaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat6_tot->GetXaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat2_tot->GetYaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat3_tot->GetYaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat4_tot->GetYaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat5_tot->GetYaxis()->SetTitleSize(0.05);
eta1_vs_deltaEta_cat6_tot->GetYaxis()->SetTitleSize(0.05);
c_2D_v2->cd(6);
eta1_vs_deltaEta_cat2_tot->Draw("colz");
c_2D_v2->cd(7);
eta1_vs_deltaEta_cat3_tot->Draw("colz");
c_2D_v2->cd(8);
eta1_vs_deltaEta_cat4_tot->Draw("colz");
c_2D_v2->cd(9);
eta1_vs_deltaEta_cat5_tot->Draw("colz");
c_2D_v2->cd(10);
eta1_vs_deltaEta_cat6_tot->Draw("colz");
c_2D_v2->SaveAs("categories_2D_v2_sig.png");
*/
}
void NuclearModificationY(){
gROOT->SetStyle("Plain");
gStyle->SetOptTitle(0);
gStyle->SetOptStat(0);
/*
gStyle->SetPadRightMargin(0.010);//###0.020
gStyle->SetPadLeftMargin(0.165);
gStyle->SetPadTopMargin(0.075);
gStyle->SetPadBottomMargin(0.145);
*/
TFile*filePPReference=new TFile(Form("../../../fonll/output%sY.root",particle.Data()));
//TFile*filePPReference=new TFile(Form("../../fonll/output%sY.root",particle.Data()));
TGraphAsymmErrors*gaeBplusReference=(TGraphAsymmErrors*)filePPReference->Get(Form("gaeSigmaDecay%s",particle.Data()));
gaeBplusReference->SetName(Form("gae%sReference",particle.Data()));
TFile*filepPb=new TFile(Form("../Results%s_y/Sigma%s.root",particle.Data(),particle.Data()));
TH1F*hSigmapPbStat=(TH1F*)filepPb->Get("hPtSigma");
double scalingfactor=1e-6;
double yvalue,xvalue,yerrorhigh,yerrorlow;
for (int i=0;i<nbins;i++){
hSigmapPbStat->SetBinContent(i+1,scalingfactor*hSigmapPbStat->GetBinContent(i+1));
hSigmapPbStat->SetBinError(i+1,scalingfactor*hSigmapPbStat->GetBinError(i+1));
yvalue=-1.;
xvalue=-1.;
yerrorhigh=-1.;
yerrorlow=-1.;
gaeBplusReference->GetPoint(i,xvalue,yvalue);
yerrorhigh=gaeBplusReference->GetEYhigh()[i];
yerrorlow=gaeBplusReference->GetEYlow()[i];
gaeBplusReference->SetPoint(i,xvalue,yvalue*scalingfactor);
gaeBplusReference->SetPointEYhigh(i,yerrorhigh*scalingfactor);
gaeBplusReference->SetPointEYlow(i,yerrorlow*scalingfactor);
}
for (int i=0;i<nbins;i++){
hSigmapPbStat->SetBinContent(i+1,(1./tagandprobcorrection[i])*(hSigmapPbStat->GetBinContent(i+1)));
hSigmapPbStat->SetBinError(i+1,(1./tagandprobcorrection[i])*(hSigmapPbStat->GetBinError(i+1)));
}
Double_t yRefPP[nbins]; //value y reference
Double_t xRefPP[nbins]; //value x reference
Double_t yPPsystFONLLhigh[nbins]; //y err syst FONLL high
Double_t yPPsystFONLLlow[nbins]; //y err syst FONLL low
Double_t yPercPPsystFONLLhigh[nbins]; //y percentuale err syst FONLL high
Double_t yPercPPsystFONLLlow[nbins]; //y percentuale err syst FONLL low
Double_t ySigmapPb[nbins]; //value y pPb
Double_t xSigmapPb[nbins]; //value x pPb
Double_t ySigmapPbStat[nbins]; //y err stat pPb
Double_t yPercSigmapPbStat[nbins]; //y err stat pPb
Double_t yFONLL[nbins]; //1
Double_t yRpA[nbins]; //value y RpA
Double_t yRpAStat[nbins]; //y err stat RpA
Double_t yRpAsystFONLLhigh[nbins]; //y err syst FONLL RpA high
Double_t yRpAsystFONLLlow[nbins]; //y err syst FONLL RpA lzow
Double_t yPercRpAsystFONLLhigh[nbins]; //y percentuale err syst FONLL RpA high
Double_t yPercRpAsystFONLLlow[nbins]; //y percentuale err syst FONLL RpA low
Double_t ySigmapPbSystTotHigh[nbins]; //y percentuale err syst pPb TOT
Double_t ySigmapPbSystTotLow[nbins]; //y percentuale err syst pPb TOT
Double_t yPercRpPbSystTotHigh[nbins]; //y percentuale err syst RpPb TOT
Double_t yPercRpPbSystTotLow[nbins]; //y percentuale err syst RpPb TOT
Double_t yRpPbSystTotHigh[nbins]; //y percentuale err syst RpPb TOT
Double_t yRpPbSystTotLow[nbins]; //y percentuale err syst RpPb TOT
Double_t yRFB[nbinsRFB];
Double_t yRFBStat[nbinsRFB];
Double_t yRFBSystTotHigh[nbinsRFB];
Double_t yRFBSystTotLow[nbinsRFB];
double x,y;
for (Int_t i=0;i<nbins;i++) {
gaeBplusReference->GetPoint(i,xRefPP[i],yRefPP[i]);
yPPsystFONLLhigh[i]=gaeBplusReference->GetEYhigh()[i];
yPPsystFONLLlow[i]=gaeBplusReference->GetEYlow()[i];
yPercPPsystFONLLhigh[i]=yPPsystFONLLhigh[i]/yRefPP[i];
yPercPPsystFONLLlow[i]=yPPsystFONLLlow[i]/yRefPP[i];
}
for(Int_t i=0;i<nbins;i++) {
ySigmapPb[i]=hSigmapPbStat->GetBinContent(i+1);
ySigmapPbStat[i]=hSigmapPbStat->GetBinError(i+1);
yPercSigmapPbStat[i]=ySigmapPbStat[i]/ySigmapPb[i];
ySigmapPbSystTotHigh[i]=yPercSigmapPbSystTotHigh[i]*ySigmapPb[i];
ySigmapPbSystTotLow[i]=yPercSigmapPbSystTotLow[i]*ySigmapPb[i];
}
for(Int_t i=0;i<nbins;i++) {
yRpA[i]=ySigmapPb[i]/yRefPP[i];
yRpAStat[i]=ySigmapPbStat[i]/yRefPP[i];
yFONLL[i]=yRpA[i];
yPercRpAsystFONLLhigh[i]=(yPercPPsystFONLLlow[i]/(1-yPercPPsystFONLLlow[i]));
yPercRpAsystFONLLlow[i]=(yPercPPsystFONLLhigh[i]/(1+yPercPPsystFONLLhigh[i]));
yRpAsystFONLLhigh[i]=yPercRpAsystFONLLhigh[i]*yRpA[i];
yRpAsystFONLLlow[i]=yPercRpAsystFONLLlow[i]*yRpA[i];
yRpPbSystTotHigh[i]=yPercSigmapPbSystTotHigh[i]*yRpA[i];
yRpPbSystTotLow[i]=yPercSigmapPbSystTotLow[i]*yRpA[i];
std::cout << i << " , " << xbins[i] << " , " << ySigmapPb[i] << " , sta: " << ySigmapPbStat[i] << " , syslow: " << ySigmapPbSystTotLow[i] << " ,syshigh: " << ySigmapPbSystTotHigh[i] << std::endl;
std::cout << "FONLL: " << yRefPP[i] << " - " << yPPsystFONLLlow[i] << " + " << yPPsystFONLLhigh[i] << std::endl;
std::cout << i << " ####### " << xbins[i] << " , " << yRpA[i] << " , sta: " << yRpAStat[i] << " , syslow: " << yRpPbSystTotLow[i] << " ,syshigh: " << yRpPbSystTotHigh[i] << " ,FONLLlow: " << yRpAsystFONLLlow[i] << " , FONLLhigh: " << yRpAsystFONLLhigh[i] << std::endl;
}
//RFB stuff
//central value
yRFB[0] = ySigmapPb[3]/ySigmapPb[2];
yRFB[1] = ySigmapPb[4]/ySigmapPb[1];
//Stat.
yRFBStat[0] = yRFB[0]*sqrt(pow(ySigmapPbStat[2]/ySigmapPb[2],2)+pow(ySigmapPbStat[3]/ySigmapPb[3],2));
yRFBStat[1] = yRFB[1]*sqrt(pow(ySigmapPbStat[4]/ySigmapPb[4],2)+pow(ySigmapPbStat[1]/ySigmapPb[1],2));
//Syst.
yRFBSystTotHigh[0] = yRFB[0]*yRFBSystTotHighRel[0];
yRFBSystTotLow[0] = yRFB[0]*yRFBSystTotLowRel[0];
yRFBSystTotHigh[1] = yRFB[1]*yRFBSystTotHighRel[1];
yRFBSystTotLow[1] = yRFB[1]*yRFBSystTotLowRel[1];
TGraphAsymmErrors *gSigmasyst = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl2,exl2,ySigmapPbSystTotLow,ySigmapPbSystTotHigh);
//###TGraphAsymmErrors *gSigmasyst = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl,exl,ySigmapPbSystTotLow,ySigmapPbSystTotHigh);
gSigmasyst->SetTitle("Sigma syst uncertainty from pPb");
gSigmasyst->SetMarkerColor(1);
gSigmasyst->SetLineColor(1);
gSigmasyst->SetLineWidth(3);
gSigmasyst->SetMarkerStyle(21);
gSigmasyst->SetMarkerColor(1);
gSigmasyst->SetFillColor(kYellow-7);//5
gSigmasyst->SetFillStyle(1001);
TGraphAsymmErrors*gSigmasyst2=(TGraphAsymmErrors*)gSigmasyst->Clone();
gSigmasyst2->SetMarkerColor(1);
gSigmasyst2->SetMarkerStyle(25);
gSigmasyst2->SetFillColor(0);
gSigmasyst2->SetFillStyle(0);
gSigmasyst2->SetLineColor(1);//5
gSigmasyst2->SetLineStyle(1);
gSigmasyst2->SetLineWidth(3);
//###TGraphAsymmErrors *gSigmastat = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl,exl,ySigmapPbStat,ySigmapPbStat);
TGraphAsymmErrors *gSigmastat = new TGraphAsymmErrors(nbins,xbins,ySigmapPb,exl0,exl0,ySigmapPbStat,ySigmapPbStat);
gSigmastat->SetTitle("Sigma stat uncertainty from pPb");
gSigmastat->SetMarkerColor(1);
gSigmastat->SetLineColor(1);
gSigmastat->SetLineWidth(3);
gSigmastat->SetMarkerStyle(21);
gSigmastat->SetMarkerColor(1);
gSigmastat->SetFillColor(0);
gSigmastat->SetFillStyle(0);
//###TCanvas *canvasSigma=new TCanvas("canvasSigma","canvasSigma",600,500);
TCanvas *canvasSigma=new TCanvas("canvasSigma","canvasSigma",500,500);
canvasSigma->cd();
canvasSigma->Range(-1.989924,-0.2917772,25.49622,2.212202);
canvasSigma->SetFillColor(0);
canvasSigma->SetBorderMode(0);
canvasSigma->SetBorderSize(2);
canvasSigma->SetLeftMargin(0.200);
canvasSigma->SetRightMargin(0.025);
canvasSigma->SetTopMargin(0.080);
canvasSigma->SetBottomMargin(0.150);
canvasSigma->SetFrameBorderMode(0);
TH2F* hempty=new TH2F("hempty","",10,-3.5,2.5,10,0,650.);
hempty->GetXaxis()->SetTitle("y_{CM}");
hempty->GetXaxis()->CenterTitle();
hempty->GetYaxis()->CenterTitle();
//###hempty->GetYaxis()->SetTitle("d#sigma / dy_{CM}(#mub GeV^{-1}c)");
hempty->GetYaxis()->SetTitle("d#sigma / dy_{CM}(#mub)");
hempty->GetXaxis()->SetTitleOffset(0.90);
hempty->GetYaxis()->SetTitleOffset(1.34);
hempty->GetXaxis()->SetTitleSize(0.070);//###0.045
hempty->GetYaxis()->SetTitleSize(0.070);//###0.045
hempty->GetXaxis()->SetTitleFont(42);
hempty->GetYaxis()->SetTitleFont(42);
hempty->GetXaxis()->SetLabelFont(42);
hempty->GetYaxis()->SetLabelFont(42);
hempty->GetXaxis()->SetLabelSize(0.060);//###0.037
hempty->GetYaxis()->SetLabelSize(0.060);//###0.037
hempty->SetMaximum(2);
hempty->SetMinimum(0.);
hempty->Draw();
// Histogram style
TH1D* hBplusReference = new TH1D("hBplusReference","",nbins,xhbins);
TH1D* hBplusReferenceEYhigh = new TH1D("hBplusReferenceEYhigh","",nbins,xhbins);
TH1D* hBplusReferenceEYlow = new TH1D("hBplusReferenceEYlow","",nbins,xhbins);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gaeBplusReference->GetPoint(i,xgae,ygae);
hBplusReference->SetBinContent(i+1,ygae);
hBplusReferenceEYhigh->SetBinContent(i+1,ygae+gaeBplusReference->GetEYhigh()[i]);
hBplusReferenceEYlow->SetBinContent(i+1,ygae-gaeBplusReference->GetEYlow()[i]);
}
//hBplusReference->SetMarkerColor(1);
//hBplusReference->SetMarkerStyle(25);
//hBplusReference->SetFillColor(kYellow-7);//5
//hBplusReference->SetFillStyle(1001);
hBplusReference->SetLineColor(kAzure-3);
hBplusReference->SetLineWidth(3);
hBplusReferenceEYhigh->SetLineColor(kAzure-3);
hBplusReferenceEYhigh->SetLineWidth(3);
hBplusReferenceEYhigh->SetLineStyle(2);
hBplusReferenceEYlow->SetLineColor(kAzure-3);
hBplusReferenceEYlow->SetLineWidth(3);
hBplusReferenceEYlow->SetLineStyle(2);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gaeBplusReference->GetPoint(i,xgae,ygae);
hBplusReference->SetBinContent(i+1,ygae);
hBplusReferenceEYhigh->SetBinContent(i+1,ygae+gaeBplusReference->GetEYhigh()[i]);
hBplusReferenceEYlow->SetBinContent(i+1,ygae-gaeBplusReference->GetEYlow()[i]);
std::cout << "### FONLL (" << i << ") : " << ygae << " + " << gaeBplusReference->GetEYhigh()[i] << " - " << gaeBplusReference->GetEYlow()[i] << std::endl;
}
for (int i=0;i<nbins;i++){
double xgae,ygae;
gSigmastat->GetPoint(i,xgae,ygae);
std::cout << "### cross section (" << i << ") : " << ygae << " stat: " << gSigmastat->GetEYhigh()[i] << " syst: " << gSigmasyst->GetEYhigh()[i] << std::endl;
}
/*
hBplusReference->Draw("][,same");
hBplusReferenceEYhigh->Draw("][,same");
hBplusReferenceEYlow->Draw("][,same");
*/
/*
gaeBplusReference->SetMarkerColor(1);
gaeBplusReference->SetMarkerStyle(25);
gaeBplusReference->SetFillColor(kYellow-7);//5
gaeBplusReference->SetFillStyle(1001);
gaeBplusReference->SetLineColor(kAzure-3);//5
gaeBplusReference->SetLineStyle(1);
gaeBplusReference->SetLineWidth(0);
gaeBplusReference->Draw("2psame");
TGraphAsymmErrors*gaeBplusReference2=(TGraphAsymmErrors*)gaeBplusReference->Clone();
gaeBplusReference2->SetMarkerColor(1);
gaeBplusReference2->SetMarkerStyle(25);
gaeBplusReference2->SetFillColor(0);
gaeBplusReference2->SetFillStyle(0);
gaeBplusReference2->SetLineColor(kAzure-3);//5
gaeBplusReference2->SetLineStyle(1);
gaeBplusReference2->SetLineWidth(1);
gaeBplusReference2->Draw("2psame");
*/
//gSigmasyst->SetFillColor(0);
//gSigmasyst->SetFillStyle(0);
gSigmasyst->Draw("2esame");//2same
gSigmasyst2->Draw("2esame");
hBplusReference->Draw("][,same");
hBplusReferenceEYhigh->Draw("][,same");
hBplusReferenceEYlow->Draw("][,same");
gSigmastat->SetMarkerColor(1);
gSigmastat->SetLineColor(1);
gSigmastat->SetLineWidth(3);
gSigmastat->SetMarkerStyle(21);
gSigmastat->SetMarkerColor(1);
gSigmastat->SetFillColor(0);
gSigmastat->Draw("epsame");
//###TLegend *legendSigma=new TLegend(0.233871,0.6701903,0.5322581,0.7526427,"");
//TLegend *legendSigma=new TLegend(0.40,0.60,0.70,0.78,"");
TLegend *legendSigma=new TLegend(0.23,0.56,0.57,0.77,"");
legendSigma->SetBorderSize(0);
legendSigma->SetLineColor(0);
legendSigma->SetFillColor(0);
legendSigma->SetFillStyle(1001);
legendSigma->SetTextFont(42);
legendSigma->SetTextSize(0.055);//###0.045
/*
TBox *c = new TBox(0.1,1-commonErrorN,4,1+commonErrorP);
c->SetLineColor(5);
c->SetFillColor(5);
c->Draw();
*/
gSigmastat->SetMarkerSize(1);
TLegendEntry *ent_SigmapPb=legendSigma->AddEntry(gSigmastat,"pPb data","p");
ent_SigmapPb->SetTextFont(42);
ent_SigmapPb->SetMarkerColor(1);
ent_SigmapPb->SetMarkerStyle(25);
ent_SigmapPb->SetLineColor(1);
ent_SigmapPb->SetLineWidth(2);
ent_SigmapPb->SetLineStyle(0);
TLegendEntry *ent_SigmapPbsyst=legendSigma->AddEntry(gSigmasyst,"Syst. pPb data","f");
ent_SigmapPbsyst->SetTextFont(42);
ent_SigmapPbsyst->SetMarkerColor(1);
ent_SigmapPbsyst->SetMarkerStyle(25);
ent_SigmapPbsyst->SetLineColor(1);
ent_SigmapPbsyst->SetLineWidth(2);
ent_SigmapPbsyst->SetLineStyle(0);
/*
///
TLegendEntry *ent_SigmapPb=legendSigma->AddEntry(gSigmasyst,"pPb","pf");
ent_SigmapPb->SetTextFont(42);
ent_SigmapPb->SetLineColor(1);
ent_SigmapPb->SetMarkerColor(1);
ent_SigmapPb->SetMarkerStyle(25);
//ent_SigmapPb->SetFillColor(kYellow-7);//5
//ent_SigmapPb->SetFillStyle(1001);
ent_SigmapPb->SetLineColor(1);
ent_SigmapPb->SetLineWidth(0);
ent_SigmapPb->SetLineStyle(0);
////
*/
/*
TLegendEntry *ent_Sigmapp=legendSigma->AddEntry(gaeBplusReference,"FONLL pp ref.","fp");
ent_Sigmapp->SetTextFont(42);
ent_Sigmapp->SetLineColor(kAzure-3);//5
ent_Sigmapp->SetLineStyle(1);
ent_Sigmapp->SetMarkerColor(1);
ent_Sigmapp->SetMarkerStyle(25);
*/
TLegendEntry *ent_Sigmapp=legendSigma->AddEntry(hBplusReference,"Scaled FONLL pp ref.","l");
ent_Sigmapp->SetTextFont(42);
ent_Sigmapp->SetLineColor(kAzure-3);//5
ent_Sigmapp->SetLineStyle(1);
ent_Sigmapp->SetMarkerColor(1);
ent_Sigmapp->SetMarkerStyle(21);
TLegendEntry *ent_Sigmapperr=legendSigma->AddEntry(hBplusReferenceEYhigh,"FONLL pp ref. uncert.","l");
ent_Sigmapperr->SetTextFont(42);
ent_Sigmapperr->SetLineColor(kAzure-3);//5
ent_Sigmapperr->SetLineStyle(2);
ent_Sigmapperr->SetMarkerColor(1);
ent_Sigmapperr->SetMarkerStyle(21);
legendSigma->Draw("psame");
/*
TBox *d = new TBox(0.1,1-commonErrorN,4,1+commonErrorP);
d->SetLineColor(1);
d->SetFillColor(0);
d->Draw();
*/
//###TLatex * tlatex1=new TLatex(0.23,0.84,"CMS");
TLatex * tlatex1=new TLatex(0.24,0.85,"CMS");
tlatex1->SetNDC();
tlatex1->SetTextColor(1);
tlatex1->SetTextFont(62);//###42
tlatex1->SetTextSize(0.07);//###0.045
tlatex1->Draw();
/*
TLatex * tlatex12=new TLatex(0.24,0.80,"Preliminary");
tlatex12->SetNDC();
tlatex12->SetTextColor(1);
tlatex12->SetTextFont(52);//###42
tlatex12->SetTextSize(0.05);//###0.045
tlatex12->Draw();
*/
TString mypar="B^{+}";
//###TLatex * tlatexlumi=new TLatex(0.671371,0.7801268,"L = 34.8 nb^{-1}");
TLatex * tlatexlumi=new TLatex(0.40,0.94,"34.6 nb^{-1} (pPb 5.02 TeV)");
tlatexlumi->SetNDC();
tlatexlumi->SetTextColor(1);
tlatexlumi->SetTextFont(42);
tlatexlumi->SetTextSize(0.06);//###0.045
tlatexlumi->Draw();
//###double xpos=0.8528226;
//###double ypos=0.6849894;
double xpos=0.85;
double ypos=0.75;
TLatex * tlatex3=new TLatex(xpos,ypos,mypar.Data());
tlatex3->SetNDC();
tlatex3->SetTextColor(1);
tlatex3->SetTextFont(42);
tlatex3->SetTextSize(0.07);//###0.06
tlatex3->Draw();
//###TLatex * tlatex4=new TLatex(0.51,0.84,"10 < p_{T}^{B} < 60 GeV/c");
TLatex * tlatex4=new TLatex(0.51,0.84,"10 < p_{T} < 60 GeV/c");
tlatex4->SetNDC();
tlatex4->SetTextColor(1);
tlatex4->SetTextFont(42);
tlatex4->SetTextSize(0.06);
tlatex4->Draw();
TLatex * tlatex5=new TLatex(0.62,0.20,Form("Global uncert. %2.1f%%",commonErrorP*100));
tlatex5->SetNDC();
tlatex5->SetTextColor(1);
tlatex5->SetTextFont(42);
tlatex5->SetTextSize(0.04);
tlatex5->Draw();
canvasSigma->SaveAs(Form("../Results%s_y/canvasSigma%s.pdf",particle.Data(),particle.Data()));
canvasSigma->SaveAs(Form("../Results%s_y/canvasSigma%s.png",particle.Data(),particle.Data()));
canvasSigma->SaveAs(Form("../Results%s_y/canvasSigma%s.eps",particle.Data(),particle.Data()));
TGraphAsymmErrors *gRpAstat = new TGraphAsymmErrors(nbins,xbins,yRpA,exl,exl,yRpAStat,yRpAStat);
gRpAstat->SetTitle("RpA stat uncertainty from pPb");
gRpAstat->SetMarkerStyle(21);
gRpAstat->SetMarkerColor(1);
gRpAstat->SetLineColor(1);
gRpAstat->SetLineWidth(3);
gRpAstat->SetFillColor(0);
TGraphAsymmErrors *gRpAsyst = new TGraphAsymmErrors(nbins,xbins,yRpA,exl,exl,yRpPbSystTotLow,yRpPbSystTotHigh);
gRpAsyst->SetTitle("RpA syst uncertainty from pPb");
gRpAsyst->SetFillColor(kYellow-7);//###0
//gRpAsyst->SetMarkerSize(0);
gRpAsyst->SetLineColor(1);
gRpAsyst->SetLineWidth(3);
gRpAsyst->SetFillStyle(1001);
gRpAsyst->SetMarkerStyle(21);
gRpAsyst->SetMarkerColor(1);
TGraphAsymmErrors*gRpAsyst2=(TGraphAsymmErrors*)gRpAsyst->Clone();
gRpAsyst2->SetMarkerColor(1);
gRpAsyst2->SetMarkerStyle(25);
gRpAsyst2->SetFillColor(0);
gRpAsyst2->SetFillStyle(0);
gRpAsyst2->SetLineColor(1);//5
gRpAsyst2->SetLineStyle(1);
gRpAsyst2->SetLineWidth(3);
TGraphAsymmErrors *gRpAsystFONLL = new TGraphAsymmErrors(nbins,xbins,yFONLL,exl,exl,yRpAsystFONLLlow,yRpAsystFONLLhigh);
gRpAsystFONLL->SetTitle("RpA syst uncertainty from FONLL reference");
gRpAsystFONLL->SetFillColor(0);//5,kYellow-7
gRpAsystFONLL->SetLineColor(kAzure-3);//5
gRpAsystFONLL->SetMarkerColor(4);//kAzure-3);
gRpAsystFONLL->SetLineStyle(2);//###1
gRpAsystFONLL->SetLineWidth(3);
TGraphAsymmErrors*gRpAsystFONLL2=(TGraphAsymmErrors*)gRpAsystFONLL->Clone();
gRpAsystFONLL2->SetMarkerColor(1);
gRpAsystFONLL2->SetMarkerStyle(25);
gRpAsystFONLL2->SetFillColor(0);
gRpAsystFONLL2->SetFillStyle(0);
gRpAsystFONLL2->SetLineColor(kAzure-3);//5
gRpAsystFONLL2->SetLineStyle(2);//###1
gRpAsystFONLL2->SetLineStyle(3);
gRpAsystFONLL2->SetLineWidth(3);
//###
TBox *gRpAsystFONLL_box0 = new TBox(xbins[0]-exl[0],yFONLL[0]-yRpAsystFONLLlow[0],xbins[0]+exl[0],yFONLL[0]+yRpAsystFONLLhigh[0]);
gRpAsystFONLL_box0->SetLineColor(kAzure-3);
gRpAsystFONLL_box0->SetLineWidth(3);
gRpAsystFONLL_box0->SetLineStyle(2);
gRpAsystFONLL_box0->SetFillColor(0);
gRpAsystFONLL_box0->SetFillStyle(0);
TBox *gRpAsystFONLL_box1 = new TBox(xbins[1]-exl[1],yFONLL[1]-yRpAsystFONLLlow[1],xbins[1]+exl[1],yFONLL[1]+yRpAsystFONLLhigh[1]);
gRpAsystFONLL_box1->SetLineColor(kAzure-3);
gRpAsystFONLL_box1->SetLineWidth(3);
gRpAsystFONLL_box1->SetLineStyle(2);
gRpAsystFONLL_box1->SetFillColor(0);
gRpAsystFONLL_box1->SetFillStyle(0);
TBox *gRpAsystFONLL_box2 = new TBox(xbins[2]-exl[2],yFONLL[2]-yRpAsystFONLLlow[2],xbins[2]+exl[2],yFONLL[2]+yRpAsystFONLLhigh[2]);
gRpAsystFONLL_box2->SetLineColor(kAzure-3);
gRpAsystFONLL_box2->SetLineWidth(3);
gRpAsystFONLL_box2->SetLineStyle(2);
gRpAsystFONLL_box2->SetFillColor(0);
gRpAsystFONLL_box2->SetFillStyle(0);
TBox *gRpAsystFONLL_box3 = new TBox(xbins[3]-exl[3],yFONLL[3]-yRpAsystFONLLlow[3],xbins[3]+exl[3],yFONLL[3]+yRpAsystFONLLhigh[3]);
gRpAsystFONLL_box3->SetLineColor(kAzure-3);
gRpAsystFONLL_box3->SetLineWidth(3);
gRpAsystFONLL_box3->SetLineStyle(2);
gRpAsystFONLL_box3->SetFillColor(0);
gRpAsystFONLL_box3->SetFillStyle(0);
TBox *gRpAsystFONLL_box4 = new TBox(xbins[4]-exl[4],yFONLL[4]-yRpAsystFONLLlow[4],xbins[4]+exl[4],yFONLL[4]+yRpAsystFONLLhigh[4]);
gRpAsystFONLL_box4->SetLineColor(kAzure-3);
gRpAsystFONLL_box4->SetLineWidth(3);
gRpAsystFONLL_box4->SetLineStyle(2);
gRpAsystFONLL_box4->SetFillColor(0);
gRpAsystFONLL_box4->SetFillStyle(0);
// Histogram style
TH1D* hRpAsystFONLL = new TH1D("hRpAsystFONLL","",nbins,xhbins);
TH1D* hRpAsystFONLLEYhigh = new TH1D("hRpAsystFONLLEYhigh","",nbins,xhbins);
TH1D* hRpAsystFONLLEYlow = new TH1D("hRpAsystFONLLEYlow","",nbins,xhbins);
for (int i=0;i<nbins;i++){
hRpAsystFONLL->SetBinContent(i+1,yFONLL[i]);
hRpAsystFONLLEYhigh->SetBinContent(i+1,yFONLL[i]+yRpAsystFONLLhigh[i]);
hRpAsystFONLLEYlow->SetBinContent(i+1,yFONLL[i]-yRpAsystFONLLlow[i]);
}
hRpAsystFONLLEYhigh->SetLineColor(kAzure-3);
hRpAsystFONLLEYhigh->SetLineWidth(3);
hRpAsystFONLLEYhigh->SetLineStyle(2);
hRpAsystFONLLEYhigh->SetFillColor(0);
hRpAsystFONLLEYhigh->SetFillStyle(0);
hRpAsystFONLLEYlow->SetLineColor(kAzure-3);
hRpAsystFONLLEYlow->SetLineWidth(3);
hRpAsystFONLLEYlow->SetLineStyle(2);
hRpAsystFONLLEYlow->SetFillColor(0);
hRpAsystFONLLEYlow->SetFillStyle(0);
hRpAsystFONLLEYlow->Draw("][,same");
//###TCanvas *canvasRpA=new TCanvas("canvasRpA","canvasRpA",600,500);
TCanvas *canvasRpA=new TCanvas("canvasRpA","canvasRpA",500,500);
canvasRpA->Range(-1.989924,-0.2917772,25.49622,2.212202);
canvasRpA->SetFillColor(0);
canvasRpA->SetBorderMode(0);
canvasRpA->SetBorderSize(2);
canvasRpA->SetLeftMargin(0.200);
canvasRpA->SetRightMargin(0.025);
canvasRpA->SetTopMargin(0.080);
canvasRpA->SetBottomMargin(0.150);
canvasRpA->SetFrameBorderMode(0);
canvasRpA->SetFrameBorderMode(0);
//###TLegend *legendRpA=new TLegend(0.23,0.58,0.57,0.77,"");
TLegend *legendRpA=new TLegend(0.23,0.56,0.57,0.77,"");
legendRpA->SetBorderSize(0);
legendRpA->SetLineColor(0);
legendRpA->SetFillColor(0);
legendRpA->SetFillStyle(1001);
legendRpA->SetTextFont(42);
legendRpA->SetTextSize(0.055);//###0.04
hempty=new TH2F("hempty","",4,-3.5,2.7,40.,0.,4.0);
hempty->GetXaxis()->SetTitle("y_{CM}");
/*
if(particle=="Bplus") hempty->GetYaxis()->SetTitle("R^{FONLL}_{pA} (B^{+})");
if(particle=="Bzero") hempty->GetYaxis()->SetTitle("R^{FONLL}_{pA} (B^{0})");
if(particle=="Bs") hempty->GetYaxis()->SetTitle("R^{FONLL}_{pA} (B_{s})");
*/
hempty->GetYaxis()->SetTitle("R^{FONLL}_{pA}");
hempty->GetXaxis()->CenterTitle();
hempty->GetYaxis()->CenterTitle();
hempty->GetXaxis()->SetTitleOffset(0.90);
hempty->GetYaxis()->SetTitleOffset(1.34);
hempty->GetXaxis()->SetTitleSize(0.070);//###0.045
hempty->GetYaxis()->SetTitleSize(0.070);//###0.045
hempty->GetXaxis()->SetTitleFont(42);
hempty->GetYaxis()->SetTitleFont(42);
hempty->GetXaxis()->SetLabelFont(42);
hempty->GetYaxis()->SetLabelFont(42);
hempty->GetXaxis()->SetLabelSize(0.060);//###0.040
hempty->GetYaxis()->SetLabelSize(0.060);//###0.040
hempty->SetMaximum(4.0);
hempty->SetMinimum(0.);
hempty->Draw();
TLine *l = new TLine(0,1,70,1);
l->SetLineStyle(2);
legendRpA->Draw();
gRpAstat->SetMarkerStyle(21);
gRpAstat->SetLineColor(1);
gRpAstat->SetMarkerColor(1);
for (int i=0;i<nbins;i++){
double xgae,ygae;
gRpAstat->GetPoint(i,xgae,ygae);
std::cout << "####### RpA (" << i << ") : " << ygae << " stat: " << gRpAstat->GetEYhigh()[i] << " syst: " << gRpAsyst->GetEYhigh()[i] << std::endl;
std::cout << "########### FONLL (" << i << ") : " << " + " << gRpAsystFONLL->GetEYhigh()[i] << " - " << gRpAsystFONLL->GetEYlow()[i] << std:: endl;
}
//###gRpAsystFONLL->Draw("e2same");//###2same
//###gRpAsystFONLL2->Draw("2esame");//###2same
/*
gRpAsystFONLL_box0->Draw("same");
gRpAsystFONLL_box1->Draw("same");
gRpAsystFONLL_box2->Draw("same");
gRpAsystFONLL_box3->Draw("same");
gRpAsystFONLL_box4->Draw("same");
*/
hRpAsystFONLLEYlow->Draw("][,same");
hRpAsystFONLLEYhigh->Draw("][,same");
gRpAsyst->Draw("2esame");
gRpAsyst2->Draw("2esame");
gRpAstat->Draw("psame");
std::cout << "%%%%%%% RpA commonError " << " - " << commonErrorN << " + " << commonErrorP << std::endl;
TBox *b = new TBox(-3.49,1-commonErrorN,-3.1,1+commonErrorP);
b->SetLineColor(1);
b->SetLineWidth(1);
b->SetLineStyle(1);
b->SetFillColor(kGray);
b->Draw();
TBox *b2=(TBox*)b->Clone();
b2->SetFillStyle(0);
b2->SetLineStyle(1);
b2->SetLineWidth(1);
b2->Draw();
TGraphAsymmErrors*gRpAstat2=(TGraphAsymmErrors*)gRpAstat->Clone();
gRpAstat2->SetTitle("RpA stat uncertainty from pPb");
gRpAstat2->SetMarkerStyle(21);
gRpAstat2->SetMarkerColor(1);
gRpAstat2->SetLineColor(1);
gRpAstat2->SetLineWidth(3);
gRpAstat2->SetFillColor(0);
gRpAstat2->SetMarkerSize(1.5);
TLegendEntry *ent_RpAstat=legendRpA->AddEntry(gRpAstat2,"R^{FONLL}_{pA}","p");
ent_RpAstat->SetTextFont(42);
ent_RpAstat->SetMarkerColor(1);
ent_RpAstat->SetMarkerStyle(25);
ent_RpAstat->SetMarkerSize(30);//###3
//ent_RpAstat->SetFillColor(kYellow-7);//5
//ent_RpAstat->SetFillStyle(1001);
ent_RpAstat->SetLineColor(1);
ent_RpAstat->SetLineWidth(3);
ent_RpAstat->SetLineStyle(0);
//TLegendEntry *ent_RpAsyst=legendRpA->AddEntry(gRpAsyst,"R^{FONLL}_{pA}","pf");
TLegendEntry *ent_RpAsyst=legendRpA->AddEntry(gRpAsyst,"Syst. pPb data","f");
ent_RpAsyst->SetTextFont(42);
ent_RpAsyst->SetMarkerColor(1);
ent_RpAsyst->SetMarkerStyle(25);
//ent_RpAsyst->SetFillColor(kYellow-7);//5
//ent_RpAsyst->SetFillStyle(1001);
ent_RpAsyst->SetLineColor(1);
ent_RpAsyst->SetLineWidth(3);
ent_RpAsyst->SetLineStyle(0);
TLegendEntry *ent_RpAsystFONLL=legendRpA->AddEntry(gRpAsystFONLL,"Syst. FONLL pp ref.","f");
ent_RpAsystFONLL->SetTextFont(42);
ent_RpAsystFONLL->SetLineColor(kAzure-3);//###5
ent_RpAsystFONLL->SetLineStyle(1);//###0
ent_RpAsystFONLL->SetLineWidth(3);
ent_RpAsystFONLL->SetMarkerColor(kYellow-7);//5
//###TLegendEntry *ent_RpAsystData=legendRpA->AddEntry(b,"Syst. int. lumi + BR","f");
TLegendEntry *ent_RpAsystData=legendRpA->AddEntry(b,"Syst. int. lumi + B","f");
ent_RpAsystData->SetTextFont(42);
ent_RpAsystData->SetLineColor(2);
ent_RpAsystData->SetMarkerColor(2);
tlatex1->Draw();
//###tlatex12->Draw();
tlatexlumi->Draw();
tlatex3->Draw();
tlatex4->Draw();
canvasRpA->SaveAs(Form("../Results%s_y/canvasRpA%s.pdf",particle.Data(),particle.Data()));
canvasRpA->SaveAs(Form("../Results%s_y/canvasRpA%s.png",particle.Data(),particle.Data()));
canvasRpA->SaveAs(Form("../Results%s_y/canvasRpA%s.eps",particle.Data(),particle.Data()));
canvasRpA->SaveAs(Form("../Results%s_y/canvasRpA%s.ps",particle.Data(),particle.Data()));
/*
TCanvas *canvasRFB=new TCanvas("canvasRFB","canvasRFB",500,500);
canvasRFB->cd();
canvasRFB->Range(-1.989924,-0.2917772,25.49622,2.212202);
canvasRFB->SetFillColor(0);
canvasRFB->SetBorderMode(0);
canvasRFB->SetBorderSize(2);
canvasRFB->SetLeftMargin(0.1451613);
canvasRFB->SetRightMargin(0.05443548);
canvasRFB->SetTopMargin(0.08474576);
canvasRFB->SetBottomMargin(0.1165254);
canvasRFB->SetFrameBorderMode(0);
canvasRFB->SetFrameBorderMode(0);
TLegend *legendRFB=new TLegend(0.2560484,0.6490486,0.5483871,0.7716702,"");
legendRFB->SetBorderSize(0);
legendRFB->SetLineColor(0);
legendRFB->SetFillColor(0);
legendRFB->SetFillStyle(1001);
legendRFB->SetTextFont(42);
legendRFB->SetTextSize(0.045);
TH2F* hempty=new TH2F("hempty","",4,-0.1,2.,10.,0.,2.5);
hempty->GetXaxis()->SetTitle("|y_{CM}|");
hempty->GetYaxis()->SetTitle("R_{FB}");
hempty->GetXaxis()->SetTitleOffset(1.1);
hempty->GetYaxis()->SetTitleOffset(1.3);
hempty->GetXaxis()->SetTitleSize(0.045);
hempty->GetYaxis()->SetTitleSize(0.045);
hempty->GetXaxis()->SetTitleFont(42);
hempty->GetYaxis()->SetTitleFont(42);
hempty->GetXaxis()->SetLabelFont(42);
hempty->GetYaxis()->SetLabelFont(42);
hempty->GetXaxis()->SetLabelSize(0.04);
hempty->GetYaxis()->SetLabelSize(0.04);
hempty->GetXaxis()->CenterTitle();
hempty->GetYaxis()->CenterTitle();
hempty->SetMaximum(3);
hempty->SetMinimum(0.);
hempty->Draw();
legendRFB->Draw();
TGraphAsymmErrors *gRFBsyst = new TGraphAsymmErrors(nbinsRFB,xbinsRFB,yRFB,exlRFB,exlRFB,yRFBSystTotLow,yRFBSystTotHigh);
gRFBsyst->SetTitle("gRFBsyst");
gRFBsyst->SetName("gRFBsyst");
gRFBsyst->SetMarkerColor(4);
gRFBsyst->SetLineColor(4);
gRFBsyst->SetLineWidth(2);
gRFBsyst->SetMarkerStyle(0);
gRFBsyst->SetMarkerColor(1);
TGraphAsymmErrors *gRFBstat = new TGraphAsymmErrors(nbinsRFB,xbinsRFB,yRFB,exlRFB,exlRFB,yRFBStat, yRFBStat);
gRFBstat->SetTitle("gRFBstat");
gRFBstat->SetName("gRFBstat");
gRFBstat->SetMarkerColor(1);
gRFBstat->SetLineColor(1);
gRFBstat->SetLineWidth(2);
gRFBstat->SetMarkerStyle(21);
gRFBstat->SetMarkerColor(1);
gRFBstat->SetFillColor(0);
gRFBstat->SetFillStyle(0);
gRFBsyst->SetFillColor(0);
gRFBsyst->SetFillStyle(0);
gRFBsyst->Draw("2same");
gRFBstat->Draw("psame");
TLegendEntry *ent_RFB=legendRFB->AddEntry(gRFBstat,"R_{FB} stat.","pl");
ent_RFB->SetTextFont(42);
ent_RFB->SetLineColor(2);
ent_RFB->SetMarkerColor(2);
TLegendEntry *ent_RFB=legendRFB->AddEntry(gRFBsyst,"R_{FB} syst.","l");
ent_RFB->SetTextFont(42);
ent_RFB->SetLineColor(4);
// ent_RFB->SetMarkerColor(2);
TLatex * tlatex1=new TLatex(0.1612903,0.8625793,"CMS pPb #sqrt{s_{NN}}= 5.02 TeV");
tlatex1->SetNDC();
tlatex1->SetTextColor(1);
tlatex1->SetTextFont(42);
tlatex1->SetTextSize(0.045);
tlatex1->Draw();
TLatex * tlatexlumi=new TLatex(0.671371,0.7801268,"L = 34.8 nb^{-1}");
tlatexlumi->SetNDC();
tlatexlumi->SetTextColor(1);
tlatexlumi->SetTextFont(42);
tlatexlumi->SetTextSize(0.045);
tlatexlumi->Draw();
//double xpos1=0.6528226;
// double ypos1=0.6849894;
TLatex * tlatex3=new TLatex(xpos,ypos,mypar.Data());
tlatex3->SetNDC();
tlatex3->SetTextColor(1);
tlatex3->SetTextFont(42);
tlatex3->SetTextSize(0.06);
tlatex3->Draw();
canvasRFB->SaveAs(Form("../Results%s_y/CanvasRDFB.pdf",particle.Data(),particle.Data()));
*/
// l->Draw();
// canvasRpA->SaveAs(Form("../Results%sY/canvasRpA%s.pdf",particle.Data(),particle.Data()));
// canvasRpA->SaveAs(Form("../Results%s_y/canvasRpA%s.pdf",particle.Data(),particle.Data()));
// TFile *fout=new TFile(Form("../Results%sY/fileRpA%s.root",particle.Data(),particle.Data()),"recreate");
TFile *fout=new TFile(Form("../Results%s_y/fileRpA%s.root",particle.Data(),particle.Data()),"recreate");
fout->cd();
gSigmasyst->SetName("gSigmasyst");
gSigmasyst->Write();
gaeBplusReference->Write();
hSigmapPbStat->Write();
}
