void OverlayAnalysis::PostDrawMultiGraph(GraphName graphName, TMultiGraph *pMultiGraph) const
{
if(NULL == pMultiGraph)
return;
// base customize
pMultiGraph->GetXaxis()->SetTitle("Distance between showers [cm]");
pMultiGraph->GetXaxis()->SetLabelFont(42);
pMultiGraph->GetXaxis()->SetTitleSize(0.05);
pMultiGraph->GetXaxis()->SetTitleOffset(1.);
pMultiGraph->GetXaxis()->SetTitleFont(42);
pMultiGraph->GetXaxis()->SetRangeUser(std::max(static_cast<int>(m_startDistance) - 5, 0), m_endDistance + 5);
pMultiGraph->GetYaxis()->SetTitle("");
pMultiGraph->GetYaxis()->SetLabelFont(42);
pMultiGraph->GetYaxis()->SetTitleSize(0.045);
pMultiGraph->GetYaxis()->SetTitleOffset(1.3);
pMultiGraph->GetYaxis()->SetTitleFont(42);
pMultiGraph->GetYaxis()->SetLabelSize(0.035);
TPaveText *pt = new TPaveText(0.3, 0.2, 0.93, 0.3, "tbNDC");
pt->SetTextSize(0.05);
pt->SetTextColor(kGray+2);
pt->SetFillColor(0);
pt->SetLineWidth(0);
pt->SetBorderSize(0);
pt->AddText("CALICE SDHCAL");
pt->SetTextFont(62);
pt->Draw();
// plot per plot customize
switch(graphName)
{
case N_PFOS:
pMultiGraph->GetYaxis()->SetTitle("<N_{pfos}>");
pMultiGraph->GetYaxis()->SetRangeUser(1, 3);
break;
case NEUTRAL_PURITY:
pMultiGraph->GetYaxis()->SetTitle("#rho_{neutral}");
pMultiGraph->GetYaxis()->SetRangeUser(0, 1);
break;
case NEUTRAL_EFFICIENCY:
pMultiGraph->GetYaxis()->SetTitle("#varepsilon_{neutral}");
pMultiGraph->GetYaxis()->SetRangeUser(0, 1);
break;
case NEUTRAL_RECOVER_PROBA:
pMultiGraph->GetYaxis()->SetTitle("P_{n>0}");
pMultiGraph->GetYaxis()->SetRangeUser(0, 1);
break;
case NEUTRAL_ENERGY_DIFFERENCE_EFFICIENT:
pMultiGraph->GetYaxis()->SetTitle("<E_{n,rec} - E_{n,meas}>, n>0");
pMultiGraph->GetYaxis()->SetRangeUser(-5, 5);
break;
}
}
void banner4Plot (){
TPaveText* pt ;
pt = new TPaveText(.14,0.91,.25,.94,"NDC");
pt->AddText("Preliminary");
pt->SetFillColor(0);
pt->SetTextSize(0.035);
pt->SetFillStyle(0);
pt->SetLineColor(0);
pt->SetLineWidth(0);
pt->SetMargin(0);
pt->SetShadowColor(0);
pt->Draw();
}
TPaveText* getHeader(double lumi, TString channelName)
{
TPaveText* pt = new TPaveText(0.18,0.75,0.18,0.88,"brNDC");
pt->SetBorderSize(1);
pt->SetTextFont(42);
pt->SetTextSize(0.04);
pt->SetLineColor(0);
pt->SetLineStyle(1);
pt->SetLineWidth(1);
pt->SetFillColor(0);
pt->SetFillStyle(1001);
pt->SetTextAlign(12);
pt->AddText("CMS Preliminary");
pt->AddText(Form("%.1f fb^{-1} at #sqrt{s} = 7 TeV", lumi));
if ( channelName != "" ) pt->AddText(channelName);
return pt;
}
void BackgroundPrediction(std::string pname,int rebin_factor,int model_number = 0,int imass=750, bool plotBands = false)
{
rebin = rebin_factor;
std::string fname = std::string("../fitFilesMETPT34/") + pname + std::string("/histos_bkg.root");
stringstream iimass ;
iimass << imass;
std::string dirName = "info_"+iimass.str()+"_"+pname;
gStyle->SetOptStat(000000000);
gStyle->SetPadGridX(0);
gStyle->SetPadGridY(0);
setTDRStyle();
gStyle->SetPadGridX(0);
gStyle->SetPadGridY(0);
gStyle->SetOptStat(0000);
writeExtraText = true; // if extra text
extraText = "Preliminary"; // default extra text is "Preliminary"
lumi_13TeV = "2.7 fb^{-1}"; // default is "19.7 fb^{-1}"
lumi_7TeV = "4.9 fb^{-1}"; // default is "5.1 fb^{-1}"
double ratio_tau=-1;
TFile *f=new TFile(fname.c_str());
TH1F *h_mX_CR_tau=(TH1F*)f->Get("distribs_18_10_1")->Clone("CR_tau");
TH1F *h_mX_SR=(TH1F*)f->Get("distribs_18_10_0")->Clone("The_SR");
double maxdata = h_mX_SR->GetMaximum();
double nEventsSR = h_mX_SR->Integral(600,4000);
ratio_tau=(h_mX_SR->GetSumOfWeights()/(h_mX_CR_tau->GetSumOfWeights()));
//double nEventsSR = h_mX_SR->Integral(600,4000);
std::cout<<"ratio tau "<<ratio_tau<<std::endl;
TH1F *h_SR_Prediction;
TH1F *h_SR_Prediction2;
if(blind) {
h_SR_Prediction2 = (TH1F*)h_mX_CR_tau->Clone("h_SR_Prediction2");
h_mX_CR_tau->Rebin(rebin);
h_mX_CR_tau->SetLineColor(kBlack);
h_SR_Prediction=(TH1F*)h_mX_CR_tau->Clone("h_SR_Prediction");
} else {
h_SR_Prediction2=(TH1F*)h_mX_SR->Clone("h_SR_Prediction2");
h_mX_SR->Rebin(rebin);
h_mX_SR->SetLineColor(kBlack);
h_SR_Prediction=(TH1F*)h_mX_SR->Clone("h_SR_Prediction");
}
h_SR_Prediction->SetMarkerSize(0.7);
h_SR_Prediction->GetYaxis()->SetTitleOffset(1.2);
h_SR_Prediction->Sumw2();
/*TFile *f_sig = new TFile((dirName+"/w_signal_"+iimass.str()+".root").c_str());
RooWorkspace* xf_sig = (RooWorkspace*)f_sig->Get("Vg");
RooAbsPdf *xf_sig_pdf = (RooAbsPdf *)xf_sig->pdf((std::string("signal_fixed_")+pname).c_str());
RooWorkspace w_sig("w");
w_sig.import(*xf_sig_pdf,RooFit::RenameVariable((std::string("signal_fixed_")+pname).c_str(),(std::string("signal_fixed_")+pname+std::string("low")).c_str()),RooFit::RenameAllVariablesExcept("low","x"));
xf_sig_pdf = w_sig.pdf((std::string("signal_fixed_")+pname+std::string("low")).c_str());
RooArgSet* biasVars = xf_sig_pdf->getVariables();
TIterator *it = biasVars->createIterator();
RooRealVar* var = (RooRealVar*)it->Next();
while (var) {
var->setConstant(kTRUE);
var = (RooRealVar*)it->Next();
}
*/
RooRealVar x("x", "m_{X} (GeV)", SR_lo, SR_hi);
RooRealVar nBackground((std::string("bg_")+pname+std::string("_norm")).c_str(),"nbkg",h_mX_SR->GetSumOfWeights());
RooRealVar nBackground2((std::string("alt_bg_")+pname+std::string("_norm")).c_str(),"nbkg",h_mX_SR->GetSumOfWeights());
std::string blah = pname;
//pname=""; //Antibtag=tag to constrain b-tag to the anti-btag shape
/* RooRealVar bg_p0((std::string("bg_p0_")+pname).c_str(), "bg_p0", 4.2, 0, 200.);
RooRealVar bg_p1((std::string("bg_p1_")+pname).c_str(), "bg_p1", 4.5, 0, 300.);
RooRealVar bg_p2((std::string("bg_p2_")+pname).c_str(), "bg_p2", 0.000047, 0, 10.1);
RooGenericPdf bg_pure = RooGenericPdf((std::string("bg_pure_")+blah).c_str(),"(pow(1-@0/13000,@1)/pow(@0/13000,@2+@3*log(@0/13000)))",RooArgList(x,bg_p0,bg_p1,bg_p2));
*/
RooRealVar bg_p0((std::string("bg_p0_")+pname).c_str(), "bg_p0", 0., -1000, 200.);
RooRealVar bg_p1((std::string("bg_p1_")+pname).c_str(), "bg_p1", -13, -1000, 1000.);
RooRealVar bg_p2((std::string("bg_p2_")+pname).c_str(), "bg_p2", -1.4, -1000, 1000.);
bg_p0.setConstant(kTRUE);
//RooGenericPdf bg_pure = RooGenericPdf((std::string("bg_pure_")+blah).c_str(),"(pow(@0/13000,@1+@2*log(@0/13000)))",RooArgList(x,bg_p1,bg_p2));
RooGenericPdf bg = RooGenericPdf((std::string("bg_")+blah).c_str(),"(pow(@0/13000,@1+@2*log(@0/13000)))",RooArgList(x,bg_p1,bg_p2));
/*TF1* biasFunc = new TF1("biasFunc","(0.63*x/1000-1.45)",1350,3600);
TF1* biasFunc2 = new TF1("biasFunc2","TMath::Min(2.,2.3*x/1000-3.8)",1350,3600);
double bias_term_s = 0;
if ((imass > 2450 && blah == "antibtag") || (imass > 1640 && blah == "btag")) {
if (blah == "antibtag") {
bias_term_s = 2.7*biasFunc->Eval(imass);
} else {
bias_term_s = 2.7*biasFunc2->Eval(imass);
}
bias_term_s/=nEventsSR;
}
RooRealVar bias_term((std::string("bias_term_")+blah).c_str(), "bias_term", 0., -bias_term_s, bias_term_s);
//bias_term.setConstant(kTRUE);
RooAddPdf bg((std::string("bg_")+blah).c_str(), "bg_all", RooArgList(*xf_sig_pdf, bg_pure), bias_term);
*/
string name_output = "CR_RooFit_Exp";
std::cout<<"Nevents "<<nEventsSR<<std::endl;
RooDataHist pred("pred", "Prediction from SB", RooArgList(x), h_SR_Prediction);
RooFitResult *r_bg=bg.fitTo(pred, RooFit::Minimizer("Minuit2"), RooFit::Range(SR_lo, SR_hi), RooFit::SumW2Error(kTRUE), RooFit::Save());
//RooFitResult *r_bg=bg.fitTo(pred, RooFit::Range(SR_lo, SR_hi), RooFit::Save());
//RooFitResult *r_bg=bg.fitTo(pred, RooFit::Range(SR_lo, SR_hi), RooFit::Save(),RooFit::SumW2Error(kTRUE));
std::cout<<" --------------------- Building Envelope --------------------- "<<std::endl;
//std::cout<< "bg_p0_"<< pname << " param "<<bg_p0.getVal() << " "<<bg_p0.getError()<<std::endl;
std::cout<< "bg_p1_"<< pname << " param "<<bg_p1.getVal() << " "<<100*bg_p1.getError()<<std::endl;
std::cout<< "bg_p2_"<< pname << " param "<<bg_p2.getVal() << " "<<100*bg_p2.getError()<<std::endl;
//std::cout<< "bias_term_"<< blah << " param 0 "<<bias_term_s<<std::endl;
RooPlot *aC_plot=x.frame();
pred.plotOn(aC_plot, RooFit::MarkerColor(kPink+2));
if (!plotBands) {
bg.plotOn(aC_plot, RooFit::VisualizeError(*r_bg, 2), RooFit::FillColor(kYellow));
bg.plotOn(aC_plot, RooFit::VisualizeError(*r_bg, 1), RooFit::FillColor(kGreen));
}
bg.plotOn(aC_plot, RooFit::LineColor(kBlue));
//pred.plotOn(aC_plot, RooFit::LineColor(kBlack), RooFit::MarkerColor(kBlack));
TGraph* error_curve[5]; //correct error bands
TGraphAsymmErrors* dataGr = new TGraphAsymmErrors(h_SR_Prediction->GetNbinsX()); //data w/o 0 entries
for (int i=2; i!=5; ++i) {
error_curve[i] = new TGraph();
}
error_curve[2] = (TGraph*)aC_plot->getObject(1)->Clone("errs");
int nPoints = error_curve[2]->GetN();
error_curve[0] = new TGraph(2*nPoints);
error_curve[1] = new TGraph(2*nPoints);
error_curve[0]->SetFillStyle(1001);
error_curve[1]->SetFillStyle(1001);
error_curve[0]->SetFillColor(kGreen);
error_curve[1]->SetFillColor(kYellow);
error_curve[0]->SetLineColor(kGreen);
error_curve[1]->SetLineColor(kYellow);
if (plotBands) {
RooDataHist pred2("pred2", "Prediction from SB", RooArgList(x), h_SR_Prediction2);
error_curve[3]->SetFillStyle(1001);
error_curve[4]->SetFillStyle(1001);
error_curve[3]->SetFillColor(kGreen);
error_curve[4]->SetFillColor(kYellow);
error_curve[3]->SetLineColor(kGreen);
error_curve[4]->SetLineColor(kYellow);
error_curve[2]->SetLineColor(kBlue);
error_curve[2]->SetLineWidth(3);
double binSize = rebin;
for (int i=0; i!=nPoints; ++i) {
double x0,y0, x1,y1;
error_curve[2]->GetPoint(i,x0,y0);
RooAbsReal* nlim = new RooRealVar("nlim","y0",y0,-100000,100000);
//double lowedge = x0 - (SR_hi - SR_lo)/double(2*nPoints);
//double upedge = x0 + (SR_hi - SR_lo)/double(2*nPoints);
double lowedge = x0 - binSize/2.;
double upedge = x0 + binSize/2.;
x.setRange("errRange",lowedge,upedge);
RooExtendPdf* epdf = new RooExtendPdf("epdf","extpdf",bg, *nlim,"errRange");
// Construct unbinned likelihood
RooAbsReal* nll = epdf->createNLL(pred2,NumCPU(2));
// Minimize likelihood w.r.t all parameters before making plots
RooMinimizer* minim = new RooMinimizer(*nll);
minim->setMinimizerType("Minuit2");
minim->setStrategy(2);
minim->setPrintLevel(-1);
minim->migrad();
minim->hesse();
RooFitResult* result = minim->lastMinuitFit();
double errm = nlim->getPropagatedError(*result);
//std::cout<<x0<<" "<<lowedge<<" "<<upedge<<" "<<y0<<" "<<nlim->getVal()<<" "<<errm<<std::endl;
error_curve[0]->SetPoint(i,x0,(y0-errm));
error_curve[0]->SetPoint(2*nPoints-i-1,x0,y0+errm);
error_curve[1]->SetPoint(i,x0,(y0-2*errm));
error_curve[1]->SetPoint(2*nPoints-i-1,x0,(y0+2*errm));
error_curve[3]->SetPoint(i,x0,-errm/sqrt(y0));
error_curve[3]->SetPoint(2*nPoints-i-1,x0,errm/sqrt(y0));
error_curve[4]->SetPoint(i,x0,-2*errm/sqrt(y0));
error_curve[4]->SetPoint(2*nPoints-i-1,x0,2*errm/sqrt(y0));
}
int npois = 0;
dataGr->SetMarkerSize(1.0);
dataGr->SetMarkerStyle (20);
const double alpha = 1 - 0.6827;
for (int i=0; i!=h_SR_Prediction->GetNbinsX(); ++i){
if (h_SR_Prediction->GetBinContent(i+1) > 0) {
int N = h_SR_Prediction->GetBinContent(i+1);
double L = (N==0) ? 0 : (ROOT::Math::gamma_quantile(alpha/2,N,1.));
double U = ROOT::Math::gamma_quantile_c(alpha/2,N+1,1) ;
dataGr->SetPoint(npois,h_SR_Prediction->GetBinCenter(i+1),h_SR_Prediction->GetBinContent(i+1));
dataGr->SetPointEYlow(npois, N-L);
dataGr->SetPointEYhigh(npois, U-N);
npois++;
}
}
}
double xG[2] = {-10,4000};
double yG[2] = {0.0,0.0};
TGraph* unityG = new TGraph(2, xG, yG);
unityG->SetLineColor(kBlue);
unityG->SetLineWidth(1);
double xPad = 0.3;
TCanvas *c_rooFit=new TCanvas("c_rooFit", "c_rooFit", 800*(1.-xPad), 600);
c_rooFit->SetFillStyle(4000);
c_rooFit->SetFrameFillColor(0);
TPad *p_1=new TPad("p_1", "p_1", 0, xPad, 1, 1);
p_1->SetFillStyle(4000);
p_1->SetFrameFillColor(0);
p_1->SetBottomMargin(0.02);
TPad* p_2 = new TPad("p_2", "p_2",0,0,1,xPad);
p_2->SetBottomMargin((1.-xPad)/xPad*0.13);
p_2->SetTopMargin(0.03);
p_2->SetFillColor(0);
p_2->SetBorderMode(0);
p_2->SetBorderSize(2);
p_2->SetFrameBorderMode(0);
p_2->SetFrameBorderMode(0);
p_1->Draw();
p_2->Draw();
p_1->cd();
int nbins = (int) (SR_hi- SR_lo)/rebin;
x.setBins(nbins);
std::cout << "chi2(data) " << aC_plot->chiSquare()<<std::endl;
//std::cout << "p-value: data under hypothesis H0: " << TMath::Prob(chi2_data->getVal(), nbins - 1) << std::endl;
aC_plot->GetXaxis()->SetRangeUser(SR_lo, SR_hi);
aC_plot->GetXaxis()->SetLabelOffset(0.02);
aC_plot->GetYaxis()->SetRangeUser(0.1, 1000.);
h_SR_Prediction->GetXaxis()->SetRangeUser(SR_lo, SR_hi);
string rebin_ = itoa(rebin);
aC_plot->GetXaxis()->SetTitle("M_{Z#gamma} [GeV] ");
aC_plot->GetYaxis()->SetTitle(("Events / "+rebin_+" GeV ").c_str());
aC_plot->SetMarkerSize(0.7);
aC_plot->GetYaxis()->SetTitleOffset(1.2);
aC_plot->Draw();
if (plotBands) {
error_curve[1]->Draw("Fsame");
error_curve[0]->Draw("Fsame");
error_curve[2]->Draw("Lsame");
dataGr->Draw("p e1 same");
}
aC_plot->SetTitle("");
TPaveText *pave = new TPaveText(0.85,0.4,0.67,0.5,"NDC");
pave->SetBorderSize(0);
pave->SetTextSize(0.05);
pave->SetTextFont(42);
pave->SetLineColor(1);
pave->SetLineStyle(1);
pave->SetLineWidth(2);
pave->SetFillColor(0);
pave->SetFillStyle(0);
char name[1000];
sprintf(name,"#chi^{2}/n = %.2f",aC_plot->chiSquare());
pave->AddText(name);
//pave->Draw();
TLegend *leg = new TLegend(0.88,0.65,0.55,0.90,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextSize(0.05);
leg->SetTextFont(42);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(0);
leg->SetFillStyle(0);
h_SR_Prediction->SetMarkerColor(kBlack);
h_SR_Prediction->SetLineColor(kBlack);
h_SR_Prediction->SetMarkerStyle(20);
h_SR_Prediction->SetMarkerSize(1.0);
//h_mMMMMa_3Tag_SR->GetXaxis()->SetTitleSize(0.09);
if (blind)
leg->AddEntry(h_SR_Prediction, "Data: sideband", "ep");
else {
if (blah == "antibtag" )
leg->AddEntry(h_SR_Prediction, "Data: anti-b-tag SR", "ep");
else
leg->AddEntry(h_SR_Prediction, "Data: b-tag SR", "ep");
}
leg->AddEntry(error_curve[2], "Fit model", "l");
leg->AddEntry(error_curve[0], "Fit #pm1#sigma", "f");
leg->AddEntry(error_curve[1], "Fit #pm2#sigma", "f");
leg->Draw();
aC_plot->Draw("axis same");
CMS_lumi( p_1, iPeriod, iPos );
p_2->cd();
RooHist* hpull;
hpull = aC_plot->pullHist();
RooPlot* frameP = x.frame() ;
frameP->SetTitle("");
frameP->GetXaxis()->SetRangeUser(SR_lo, SR_hi);
frameP->addPlotable(hpull,"P");
frameP->GetYaxis()->SetRangeUser(-7,7);
frameP->GetYaxis()->SetNdivisions(505);
frameP->GetYaxis()->SetTitle("#frac{(data-fit)}{#sigma_{stat}}");
frameP->GetYaxis()->SetTitleSize((1.-xPad)/xPad*0.06);
frameP->GetYaxis()->SetTitleOffset(1.0/((1.-xPad)/xPad));
frameP->GetXaxis()->SetTitleSize((1.-xPad)/xPad*0.06);
//frameP->GetXaxis()->SetTitleOffset(1.0);
frameP->GetXaxis()->SetLabelSize((1.-xPad)/xPad*0.05);
frameP->GetYaxis()->SetLabelSize((1.-xPad)/xPad*0.05);
frameP->Draw();
if (plotBands) {
error_curve[4]->Draw("Fsame");
error_curve[3]->Draw("Fsame");
unityG->Draw("same");
hpull->Draw("psame");
frameP->Draw("axis same");
}
c_rooFit->SaveAs((dirName+"/"+name_output+".pdf").c_str());
const int nModels = 9;
TString models[nModels] = {
"env_pdf_0_13TeV_dijet2", //0
"env_pdf_0_13TeV_exp1", //1
"env_pdf_0_13TeV_expow1", //2
"env_pdf_0_13TeV_expow2", //3 => skip
"env_pdf_0_13TeV_pow1", //4
"env_pdf_0_13TeV_lau1", //5
"env_pdf_0_13TeV_atlas1", //6
"env_pdf_0_13TeV_atlas2", //7 => skip
"env_pdf_0_13TeV_vvdijet1" //8
};
int nPars[nModels] = {
2, 1, 2, 3, 1, 1, 2, 3, 2
};
TString parNames[nModels][3] = {
"env_pdf_0_13TeV_dijet2_log1","env_pdf_0_13TeV_dijet2_log2","",
"env_pdf_0_13TeV_exp1_p1","","",
"env_pdf_0_13TeV_expow1_exp1","env_pdf_0_13TeV_expow1_pow1","",
"env_pdf_0_13TeV_expow2_exp1","env_pdf_0_13TeV_expow2_pow1","env_pdf_0_13TeV_expow2_exp2",
"env_pdf_0_13TeV_pow1_p1","","",
"env_pdf_0_13TeV_lau1_l1","","",
"env_pdf_0_13TeV_atlas1_coeff1","env_pdf_0_13TeV_atlas1_log1","",
"env_pdf_0_13TeV_atlas2_coeff1","env_pdf_0_13TeV_atlas2_log1","env_pdf_0_13TeV_atlas2_log2",
"env_pdf_0_13TeV_vvdijet1_coeff1","env_pdf_0_13TeV_vvdijet1_log1",""
}
if(bias){
//alternative model
gSystem->Load("libHiggsAnalysisCombinedLimit");
gSystem->Load("libdiphotonsUtils");
TFile *f = new TFile("antibtag_multipdf.root");
RooWorkspace* xf = (RooWorkspace*)f->Get("wtemplates");
RooWorkspace *w_alt=new RooWorkspace("Vg");
for(int i=model_number; i<=model_number; i++){
RooMultiPdf *alternative = (RooMultiPdf *)xf->pdf("model_bkg_AntiBtag");
std::cout<<"Number of pdfs "<<alternative->getNumPdfs()<<std::endl;
for (int j=0; j!=alternative->getNumPdfs(); ++j){
std::cout<<alternative->getPdf(j)->GetName()<<std::endl;
}
RooAbsPdf *alt_bg = alternative->getPdf(alternative->getCurrentIndex()+i);//->clone();
w_alt->import(*alt_bg, RooFit::RenameVariable(alt_bg->GetName(),("alt_bg_"+blah).c_str()));
w_alt->Print("V");
std::cerr<<w_alt->var("x")<<std::endl;
RooRealVar * range_ = w_alt->var("x");
range_->setRange(SR_lo,SR_hi);
char* asd = ("alt_bg_"+blah).c_str() ;
w_alt->import(nBackground2);
std::cout<<alt_bg->getVal() <<std::endl;
w_alt->pdf(asd)->fitTo(pred, RooFit::Minimizer("Minuit2"), RooFit::Range(SR_lo, SR_hi), RooFit::SumW2Error(kTRUE), RooFit::Save());
RooArgSet* altVars = w_alt->pdf(asd)->getVariables();
TIterator *it2 = altVars->createIterator();
RooRealVar* varAlt = (RooRealVar*)it2->Next();
while (varAlt) {
varAlt->setConstant(kTRUE);
varAlt = (RooRealVar*)it2->Next();
}
alt_bg->plotOn(aC_plot, RooFit::LineColor(i+1), RooFit::LineStyle(i+2));
p_1->cd();
aC_plot->GetYaxis()->SetRangeUser(0.01, maxdata*50.);
aC_plot->Draw("same");
TH1F *h=new TH1F();
h->SetLineColor(1+i);
h->SetLineStyle(i+2);
leg->AddEntry(h, alt_bg->GetName(), "l");
w_alt->SaveAs((dirName+"/w_background_alternative.root").c_str());
}
leg->Draw();
p_1->SetLogy();
c_rooFit->Update();
c_rooFit->SaveAs((dirName+"/"+name_output+blah+"_multipdf.pdf").c_str());
for (int i=0; i!=nPars[model_number]; ++i) {
std::cout<<parNames[model_number][i]<<" param "<< w_alt->var(parNames[model_number][i])->getVal()<<" "<<w_alt->var(parNames[model_number][i])->getError()<<std::endl;
}
} else {
p_1->SetLogy();
c_rooFit->Update();
c_rooFit->SaveAs((dirName+"/"+name_output+"_log.pdf").c_str());
}
RooWorkspace *w=new RooWorkspace("Vg");
w->import(bg);
w->import(nBackground);
w->SaveAs((dirName+"/w_background_GaussExp.root").c_str());
TH1F *h_mX_SR_fakeData=(TH1F*)h_mX_SR->Clone("h_mX_SR_fakeData");
h_mX_SR_fakeData->Scale(nEventsSR/h_mX_SR_fakeData->GetSumOfWeights());
RooDataHist data_obs("data_obs", "Data", RooArgList(x), h_mX_SR_fakeData);
std::cout<<" Background number of events = "<<nEventsSR<<std::endl;
RooWorkspace *w_data=new RooWorkspace("Vg");
w_data->import(data_obs);
w_data->SaveAs((dirName+"/w_data.root").c_str());
}
void plot_BSM_MCFM(int SignalOnly=0){
gROOT->ProcessLine(".x tdrstyle.cc");
double mPOLE = 125.6;
TString OUTPUT_NAME;
OUTPUT_NAME = "HtoZZ4l_MCFM_125p6_BSMPlots";
if (SignalOnly==0) OUTPUT_NAME.Append(".root");
else OUTPUT_NAME.Append("_SignalOnly.root");
TString coutput_common = user_TemplateswithTrees_dir + "../BSMReweight_GenLevel/Plots/";
gSystem->Exec("mkdir -p " + coutput_common);
TString coutput = coutput_common + OUTPUT_NAME;
TFile* foutput = new TFile(coutput, "recreate");
foutput->cd();
float ZZMass = 0;
float p0plus_VAJHU;
float p0hplus_VAJHU;
float p0minus_VAJHU;
float p0_g1prime2_VAJHU;
float p0_g1prime4_VAJHU;
float pg1g2_VAJHU;
float pg1g4_VAJHU;
float pg1g1prime2_VAJHU;
float pg1g1prime4_VAJHU;
TChain* tc[2][3][4];
TH1F* hfill[4][9];
int nbinsx = 73;
double xlimits[2]={ 160, 1620 };
if (SignalOnly==1){
xlimits[0]=100;
nbinsx = 76;
}
for (int t=0; t<4; t++){
for (int ac=0; ac<9; ac++){
hfill[t][ac]= new TH1F(Form("hSum_BSI%i_Hypo%i", t, ac), "", nbinsx, xlimits[0], xlimits[1]);
hfill[t][ac]->SetXTitle("m_{4l} (GeV)");
hfill[t][ac]->SetYTitle(Form("Events / %.0f GeV", (xlimits[1]-xlimits[0])/nbinsx));
}
}
double nCounted[2][3][9]={ { { 0 } } };
double nCountedScaled[2][3][9]={ { { 0 } } };
for (int erg_tev=7; erg_tev<9; erg_tev++){
for (int folder=0; folder<3; folder++){
TString comstring;
comstring.Form("%iTeV", erg_tev);
TString erg_dir;
erg_dir.Form("LHC_%iTeV/", erg_tev);
int EnergyIndex = 1;
if (erg_tev == 7) EnergyIndex = 0;
TString cinput_common = user_TemplateswithTrees_dir + "../BSMReweight_GenLevel/";
cinput_common.Append(+erg_dir);
cinput_common += user_folder[folder] + "/";
cout << cinput_common << endl;
for (int t=0; t<4; t++){
TString INPUT_NAME;
INPUT_NAME = "HtoZZ4l_MCFM_125p6_BSMTrees_";
INPUT_NAME += sample_suffix[t] + ".root";
INPUT_NAME.Prepend(cinput_common);
tc[EnergyIndex][folder][t] = new TChain("GenTree");
if (t!=3) tc[EnergyIndex][folder][t]->Add(INPUT_NAME);
tc[EnergyIndex][folder][t]->SetBranchAddress("ZZMass", &ZZMass);
tc[EnergyIndex][folder][t]->SetBranchAddress("p0plus_VAJHU", &p0plus_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("p0hplus_VAJHU", &p0hplus_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("p0minus_VAJHU", &p0minus_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("p0_g1prime2_VAJHU", &p0_g1prime2_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("p0_g1prime4_VAJHU", &p0_g1prime4_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("pg1g2_VAJHU", &pg1g2_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("pg1g4_VAJHU", &pg1g4_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("pg1g1prime2_VAJHU", &pg1g1prime2_VAJHU);
tc[EnergyIndex][folder][t]->SetBranchAddress("pg1g1prime4_VAJHU", &pg1g1prime4_VAJHU);
}
double nsig_counted[9] ={ 0 };
for (int ev=0; ev<tc[EnergyIndex][folder][0]->GetEntries(); ev++){
tc[EnergyIndex][folder][0]->GetEntry(ev);
if (fabs(ZZMass-mPOLE)<0.1){
nsig_counted[0] += p0plus_VAJHU;
nsig_counted[1] += p0hplus_VAJHU;
nsig_counted[2] += p0minus_VAJHU;
nsig_counted[3] += p0_g1prime2_VAJHU;
nsig_counted[4] += p0_g1prime4_VAJHU;
nsig_counted[5] += (2.*(p0plus_VAJHU+p0hplus_VAJHU)-pg1g2_VAJHU);
nsig_counted[6] += pg1g4_VAJHU;
nsig_counted[7] += pg1g1prime2_VAJHU;
nsig_counted[8] += pg1g1prime4_VAJHU;
}
}
for (int ac=0; ac<9; ac++) nCounted[EnergyIndex][folder][ac] = nsig_counted[ac];
}
}
for (int erg_tev=7; erg_tev<9; erg_tev++){
for (int folder=0; folder<3; folder++){
int EnergyIndex = 1;
if (erg_tev == 7) EnergyIndex = 0;
double nsig_SM = nSM_ScaledPeak[EnergyIndex][folder];
double SMscale = nsig_SM/nCounted[EnergyIndex][folder][0];
double scale=1;
for (int t=0; t<4; t++){
for (int ev=0; ev<tc[EnergyIndex][folder][t]->GetEntries(); ev++){
tc[EnergyIndex][folder][t]->GetEntry(ev);
if (ZZMass<xlimits[0]) continue;
if (ZZMass>=xlimits[1]) ZZMass=xlimits[1]*0.999;
scale = SMscale;
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][0] = nCounted[EnergyIndex][folder][0]*scale;
hfill[t][0]->Fill(ZZMass, p0plus_VAJHU*scale);
scale = SMscale;
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][1] = nCounted[EnergyIndex][folder][1]*scale;
hfill[t][1]->Fill(ZZMass, p0hplus_VAJHU*scale);
scale = SMscale;
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][2] = nCounted[EnergyIndex][folder][2]*scale;
hfill[t][2]->Fill(ZZMass, p0minus_VAJHU*scale);
scale = SMscale;
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][3] = nCounted[EnergyIndex][folder][3]*scale;
hfill[t][3]->Fill(ZZMass, p0_g1prime2_VAJHU*scale);
scale = SMscale;
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][4] = nCounted[EnergyIndex][folder][4]*scale;
hfill[t][4]->Fill(ZZMass, p0_g1prime4_VAJHU*scale);
scale = SMscale;
if (t>0){
hfill[t][5]->Fill(ZZMass, pg1g2_VAJHU*scale);
}
else{
scale *= (nCounted[0][2][0]+nCounted[1][2][0])/(nCounted[0][2][5]+nCounted[1][2][5]);
hfill[t][5]->Fill(ZZMass, (2.*(p0plus_VAJHU+p0hplus_VAJHU)-pg1g2_VAJHU)*scale);
}
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][5] = nCounted[EnergyIndex][folder][5]*scale;
scale = SMscale;
if (t==0) scale *= (nCounted[0][2][0]+nCounted[1][2][0])/(nCounted[0][2][6]+nCounted[1][2][6]);
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][6] = nCounted[EnergyIndex][folder][6]*scale;
hfill[t][6]->Fill(ZZMass, pg1g4_VAJHU*scale);
scale = SMscale;
if (t==0) scale *= (nCounted[0][2][0]+nCounted[1][2][0])/(nCounted[0][2][7]+nCounted[1][2][7]);
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][7] = nCounted[EnergyIndex][folder][7]*scale;
hfill[t][7]->Fill(ZZMass, pg1g1prime2_VAJHU*scale);
scale = SMscale;
// if (t==0) scale *= (nCounted[0][2][0]+nCounted[1][2][0])/(nCounted[0][2][8]+nCounted[1][2][8]);
if (t==0 && ev==0) nCountedScaled[EnergyIndex][folder][8] = nCounted[EnergyIndex][folder][8]*scale;
hfill[t][8]->Fill(ZZMass, pg1g1prime4_VAJHU*scale);
}
delete tc[EnergyIndex][folder][t];
}
}
}
for (int ac=1; ac<9; ac++){
double nTotal[2]={ 0 };
for (int erg_tev=7; erg_tev<9; erg_tev++){
for (int folder=0; folder<3; folder++){
int EnergyIndex = 1;
if (erg_tev == 7) EnergyIndex = 0;
nTotal[0] += nCountedScaled[EnergyIndex][folder][0];
nTotal[1] += nCountedScaled[EnergyIndex][folder][ac];
}
}
if (ac!=8) hfill[0][ac]->Scale(nTotal[0]/nTotal[1]);
else hfill[0][ac]->Scale(0.5);
}
double maxplot=0;
for (int t=0; t<4; t++){
for (int ac=0; ac<9; ac++){
if (SignalOnly==0 && ac<5) maxplot = max(maxplot, hfill[t][ac]->GetMaximum());
else if (SignalOnly==1 && !(ac<5 && ac>0) && t==0) maxplot = max(maxplot, hfill[t][ac]->GetMaximum());
hfill[t][ac]->SetLineWidth(2);
if (t==0 && ac>=5){
hfill[t][ac]->SetLineStyle(7);
// hfill[t][ac]->Add(hfill[1][ac]);
}
if (t==1) hfill[t][ac]->SetLineStyle(3);
if (t==3) hfill[t][ac]->SetLineStyle(9);
hfill[t][ac]->GetXaxis()->SetLabelFont(42);
hfill[t][ac]->GetXaxis()->SetLabelOffset(0.007);
hfill[t][ac]->GetXaxis()->SetLabelSize(0.04);
hfill[t][ac]->GetXaxis()->SetTitleSize(0.06);
hfill[t][ac]->GetXaxis()->SetTitleOffset(0.9);
hfill[t][ac]->GetXaxis()->SetTitleFont(42);
hfill[t][ac]->GetYaxis()->SetNdivisions(505);
hfill[t][ac]->GetYaxis()->SetLabelFont(42);
hfill[t][ac]->GetYaxis()->SetLabelOffset(0.007);
hfill[t][ac]->GetYaxis()->SetLabelSize(0.04);
hfill[t][ac]->GetYaxis()->SetTitleSize(0.06);
hfill[t][ac]->GetYaxis()->SetTitleOffset(1.1);
hfill[t][ac]->GetYaxis()->SetTitleFont(42);
}
}
TPaveText* pt = new TPaveText(0.15, 0.93, 0.85, 1, "brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetTextAlign(12);
pt->SetTextFont(42);
pt->SetTextSize(0.045);
TText* text = pt->AddText(0.025, 0.45, "#font[61]{CMS}");
text->SetTextSize(0.044);
text = pt->AddText(0.165, 0.42, "#font[52]{Simulation}");
text->SetTextSize(0.0315);
TString cErgTev = "#font[42]{19.7 fb^{-1} (8 TeV) + 5.1 fb^{-1} (7 TeV)}";
text = pt->AddText(0.537, 0.45, cErgTev);
text->SetTextSize(0.0315);
float lxmin = 0.22;
float lxwidth = 0.38;
float lymax = 0.9;
float lywidth = 0.3;
float lxmax = lxmin + lxwidth;
float lymin = lymax;
if (SignalOnly==0) lymin -= lywidth*4./5.;
else lymin -= lywidth;
float lxmin2 = 0.22+0.39;
float lymax2 = lymax;
float lxmax2 = lxmin2 + lxwidth;
float lymin2 = lymax2;
if (SignalOnly==0) lymin2 -= lywidth*2./5.;
else lymin2 -= lywidth*4./5.;
if (SignalOnly==1){
float lxmin3 = lxmin2;
float lymax3 = lymax2;
float lxmax3 = lxmax2;
float lymin3 = lymin2;
lxmin2 = lxmin;
lxmax2 = lxmax;
lymin2 = lymin;
lymax2 = lymax;
lxmin = lxmin3;
lxmax = lxmax3;
lymin = lymin3;
lymax = lymax3;
}
float pxmin = 0.756;
float pymin = 0.76;
float pxmax = 0.85;
if (SignalOnly==1){
pymin -= 0.12;
}
float pymax = pymin+0.05;
TPaveText* ptx = new TPaveText(pxmin, pymin, pxmax, pymax, "brNDC");
ptx->SetBorderSize(0);
ptx->SetTextFont(42);
ptx->SetTextSize(0.04);
ptx->SetLineColor(1);
ptx->SetLineStyle(1);
ptx->SetLineWidth(1);
ptx->SetFillColor(0);
ptx->SetFillStyle(0);
text = ptx->AddText(0.01, 0.01, "gg#rightarrow4l");
text->SetTextSize(0.04);
TString canvasname = "cCanvas_MCFMBSM_GenLevel";
if (SignalOnly==1) canvasname.Append("_SignalOnly");
TCanvas* cc = new TCanvas(canvasname, "", 8, 30, 800, 800);
cc->cd();
gStyle->SetOptStat(0);
cc->SetFillColor(0);
cc->SetBorderMode(0);
cc->SetBorderSize(2);
cc->SetTickx(1);
cc->SetTicky(1);
cc->SetLeftMargin(0.17);
cc->SetRightMargin(0.05);
cc->SetTopMargin(0.07);
cc->SetBottomMargin(0.13);
cc->SetFrameFillStyle(0);
cc->SetFrameBorderMode(0);
cc->SetFrameFillStyle(0);
cc->SetFrameBorderMode(0);
cc->SetLogy();
TLegend* ll;
TLegend* ll2;
ll = new TLegend(lxmin2, lymin2, lxmax2, lymax2);
ll2 = new TLegend(lxmin, lymin, lxmax, lymax);
ll->SetBorderSize(0);
ll->SetTextFont(42);
ll->SetTextSize(0.04);
ll->SetLineColor(1);
ll->SetLineStyle(1);
ll->SetLineWidth(1);
ll->SetFillColor(0);
ll->SetFillStyle(0);
ll2->SetBorderSize(0);
ll2->SetTextFont(42);
ll2->SetTextSize(0.04);
ll2->SetLineColor(1);
ll2->SetLineStyle(1);
ll2->SetLineWidth(1);
ll2->SetFillColor(0);
ll2->SetFillStyle(0);
TString strACtitle[9]={ "",
"f_{a2}=1", "f_{a3}=1", "f_{#Lambda1}=1", "f_{#LambdaQ}=1",
"f_{a2}=0.5, #phi_{#lower[-0.2]{a2}}=#pi", "f_{a3}=0.5", "f_{#Lambda1}=0.5", "f_{#LambdaQ}=0.5"
};
int iDraw = 2 - 2*SignalOnly;
if (SignalOnly==0) hfill[iDraw][0]->GetYaxis()->SetRangeUser(7e-3, maxplot*15.);
else{
double histmin = 7e-3;
if (hfill[iDraw][0]->GetMinimum()>0) histmin = hfill[iDraw][0]->GetMinimum();
hfill[iDraw][0]->GetYaxis()->SetRangeUser(histmin, maxplot*2000.);
}
hfill[iDraw][0]->GetXaxis()->SetRangeUser(xlimits[0], 800.);
hfill[iDraw][0]->SetLineColor(kBlack);
if (SignalOnly==0){
hfill[iDraw][0]->SetFillColor(kAzure-2);
hfill[iDraw][0]->SetFillStyle(1001);
}
hfill[iDraw][0]->Draw("hist");
hfill[iDraw][1]->SetLineColor(kBlue);
hfill[iDraw][1]->Draw("histsame");
hfill[iDraw][2]->SetLineColor(kRed);
hfill[iDraw][2]->Draw("histsame");
hfill[iDraw][3]->SetLineColor(kViolet);
hfill[iDraw][3]->Draw("histsame");
hfill[iDraw][4]->SetLineColor(kGreen+2);
hfill[iDraw][4]->Draw("histsame");
if (SignalOnly==1){
hfill[iDraw][5]->SetLineColor(kBlue);
hfill[iDraw][5]->Draw("histsame");
hfill[iDraw][6]->SetLineColor(kRed);
hfill[iDraw][6]->Draw("histsame");
hfill[iDraw][7]->SetLineColor(kViolet);
hfill[iDraw][7]->Draw("histsame");
hfill[iDraw][8]->SetLineColor(kGreen+2);
hfill[iDraw][8]->Draw("histsame");
}
if (SignalOnly==0){
hfill[1][0]->SetLineColor(kBlack);
hfill[1][0]->SetLineStyle(3);
hfill[1][0]->Draw("histsame");
}
hfill[iDraw][0]->Draw("histsame");
TLegendEntry* legendtext;
if (SignalOnly==0){
legendtext = ll->AddEntry(hfill[iDraw][0], "SM total", "f");
legendtext = ll->AddEntry(hfill[1][0], "SM bkg.", "f");
legendtext->SetFillStyle(1001);
legendtext->SetFillColor(hfill[1][0]->GetFillColor());
}
else{
legendtext = ll->AddEntry(hfill[iDraw][0], "SM signal", "f");
legendtext->SetFillStyle(3001);
}
if (SignalOnly==0){
legendtext = ll2->AddEntry(hfill[iDraw][4], Form("#Gamma_{H}=#Gamma^{SM}_{H}, %s", strACtitle[4].Data()), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][4]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][2], Form("#Gamma_{H}=#Gamma^{SM}_{H}, %s", strACtitle[2].Data()), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][2]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][1], Form("#Gamma_{H}=#Gamma^{SM}_{H}, %s", strACtitle[1].Data()), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][1]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][3], Form("#Gamma_{H}=#Gamma^{SM}_{H}, %s", strACtitle[3].Data()), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][3]->GetFillColor());
}
else{
legendtext = ll->AddEntry(hfill[iDraw][4], strACtitle[4].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][4]->GetFillColor());
legendtext = ll->AddEntry(hfill[iDraw][2], strACtitle[2].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][2]->GetFillColor());
legendtext = ll->AddEntry(hfill[iDraw][1], strACtitle[1].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][1]->GetFillColor());
legendtext = ll->AddEntry(hfill[iDraw][3], strACtitle[3].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][3]->GetFillColor());
}
if (SignalOnly==1){
legendtext = ll2->AddEntry(hfill[iDraw][8], strACtitle[8].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][8]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][5], strACtitle[5].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][5]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][7], strACtitle[7].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][7]->GetFillColor());
legendtext = ll2->AddEntry(hfill[iDraw][6], strACtitle[6].Data(), "f");
legendtext->SetFillStyle(3001);
legendtext->SetFillColor(hfill[iDraw][6]->GetFillColor());
}
ll->Draw("same");
ll2->Draw("same");
ptx->Draw();
pt->Draw();
cc->RedrawAxis();
cc->Update();
canvasname.Prepend(coutput_common);
TString canvasname_pdf = canvasname;
TString canvasname_eps = canvasname;
TString canvasname_png = canvasname;
TString canvasname_root = canvasname;
TString canvasname_c = canvasname;
canvasname_pdf.Append(".pdf");
canvasname_eps.Append(".eps");
canvasname_png.Append(".png");
canvasname_root.Append(".root");
canvasname_c.Append(".C");
cc->SaveAs(canvasname_pdf);
cc->SaveAs(canvasname_eps);
cc->SaveAs(canvasname_png);
cc->SaveAs(canvasname_root);
cc->SaveAs(canvasname_c);
foutput->WriteTObject(cc);
delete ll2;
delete ll;
cc->Close();
delete ptx;
delete pt;
for (int t=0; t<4; t++){
for (int ac=0; ac<5; ac++){
foutput->WriteTObject(hfill[t][ac]);
delete hfill[t][ac];
}
}
foutput->Close();
}
void zj_ptj_8TeV()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Thu Apr 28 10:19:34 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,800,800);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: bottomPad
TPad *bottomPad = new TPad("bottomPad", "bottomPad",0.005,0.05,0.995,0.995);
bottomPad->Draw();
bottomPad->cd();
bottomPad->Range(0,0,1,1);
bottomPad->SetFillColor(0);
bottomPad->SetFillStyle(4000);
bottomPad->SetBorderMode(0);
bottomPad->SetBorderSize(2);
bottomPad->SetTickx(1);
bottomPad->SetTicky(1);
bottomPad->SetLeftMargin(0.14);
bottomPad->SetRightMargin(0.05);
bottomPad->SetTopMargin(0.05);
bottomPad->SetBottomMargin(0.16);
bottomPad->SetFrameLineWidth(2);
bottomPad->SetFrameBorderMode(0);
bottomPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(6.172839,-0.2132695,104.9383,0.4330016);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.33);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
Double_t xAxis19[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_toterr__19 = new TH1F("z_ptj_toterr__19","z_ptj",5, xAxis19);
z_ptj_toterr__19->SetBinContent(1,0.3608738);
z_ptj_toterr__19->SetBinContent(2,0.225932);
z_ptj_toterr__19->SetBinContent(3,0.1253807);
z_ptj_toterr__19->SetBinContent(4,0.06703913);
z_ptj_toterr__19->SetBinContent(5,0.01561868);
z_ptj_toterr__19->SetBinContent(6,0.780934);
z_ptj_toterr__19->SetBinError(1,0.02073388);
z_ptj_toterr__19->SetBinError(2,0.01400757);
z_ptj_toterr__19->SetBinError(3,0.009350811);
z_ptj_toterr__19->SetBinError(4,0.004623786);
z_ptj_toterr__19->SetBinError(5,0.001395185);
z_ptj_toterr__19->SetBinError(6,0.06975923);
z_ptj_toterr__19->SetEntries(12);
z_ptj_toterr__19->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetLineColor(ci);
z_ptj_toterr__19->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr__19->SetMarkerColor(ci);
z_ptj_toterr__19->SetMarkerStyle(0);
z_ptj_toterr__19->GetXaxis()->SetNdivisions(4);
z_ptj_toterr__19->GetXaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetXaxis()->SetLabelOffset(999);
z_ptj_toterr__19->GetXaxis()->SetLabelSize(0);
z_ptj_toterr__19->GetXaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetXaxis()->SetTitleFont(132);
z_ptj_toterr__19->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_toterr__19->GetYaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetYaxis()->SetLabelSize(0.05);
z_ptj_toterr__19->GetYaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetYaxis()->SetTitleOffset(1.3);
z_ptj_toterr__19->GetYaxis()->SetTitleFont(132);
z_ptj_toterr__19->GetZaxis()->SetLabelFont(132);
z_ptj_toterr__19->GetZaxis()->SetLabelSize(0.05);
z_ptj_toterr__19->GetZaxis()->SetTitleSize(0.06);
z_ptj_toterr__19->GetZaxis()->SetTitleFont(132);
z_ptj_toterr__19->Draw("e2");
Double_t xAxis20[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_staterr__20 = new TH1F("z_ptj_staterr__20","z_ptj",5, xAxis20);
z_ptj_staterr__20->SetBinContent(1,0.3608738);
z_ptj_staterr__20->SetBinContent(2,0.225932);
z_ptj_staterr__20->SetBinContent(3,0.1253807);
z_ptj_staterr__20->SetBinContent(4,0.06703913);
z_ptj_staterr__20->SetBinContent(5,0.01561868);
z_ptj_staterr__20->SetBinContent(6,0.780934);
z_ptj_staterr__20->SetBinError(1,0.007287406);
z_ptj_staterr__20->SetBinError(2,0.005663236);
z_ptj_staterr__20->SetBinError(3,0.002951952);
z_ptj_staterr__20->SetBinError(4,0.002192817);
z_ptj_staterr__20->SetBinError(5,0.0004728113);
z_ptj_staterr__20->SetBinError(6,0.02364057);
z_ptj_staterr__20->SetEntries(12);
z_ptj_staterr__20->SetStats(0);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetLineColor(ci);
z_ptj_staterr__20->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
z_ptj_staterr__20->SetMarkerColor(ci);
z_ptj_staterr__20->SetMarkerStyle(0);
z_ptj_staterr__20->GetXaxis()->SetNdivisions(4);
z_ptj_staterr__20->GetXaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetXaxis()->SetLabelOffset(999);
z_ptj_staterr__20->GetXaxis()->SetLabelSize(0);
z_ptj_staterr__20->GetXaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetXaxis()->SetTitleFont(132);
z_ptj_staterr__20->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_staterr__20->GetYaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetYaxis()->SetLabelSize(0.05);
z_ptj_staterr__20->GetYaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetYaxis()->SetTitleFont(132);
z_ptj_staterr__20->GetZaxis()->SetLabelFont(132);
z_ptj_staterr__20->GetZaxis()->SetLabelSize(0.05);
z_ptj_staterr__20->GetZaxis()->SetTitleSize(0.06);
z_ptj_staterr__20->GetZaxis()->SetTitleFont(132);
z_ptj_staterr__20->Draw("e2same");
Double_t ptj_toterrs_fx3022[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t ptj_toterrs_fy3022[5] = {
0.3671394,
0.2404275,
0.1334798,
0.06537649,
0.01594908};
Double_t ptj_toterrs_felx3022[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fely3022[5] = {
0.02997783,
0.02293181,
0.009217184,
0.004503167,
0.001643067};
Double_t ptj_toterrs_fehx3022[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fehy3022[5] = {
0.02512764,
0.01751375,
0.008693834,
0.005222099,
0.001753864};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(5,ptj_toterrs_fx3022,ptj_toterrs_fy3022,ptj_toterrs_felx3022,ptj_toterrs_fehx3022,ptj_toterrs_fely3022,ptj_toterrs_fehy3022);
grae->SetName("ptj_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph_ptj_toterrs30153022 = new TH1F("Graph_Graph_ptj_toterrs30153022","",100,17.325,88.425);
Graph_Graph_ptj_toterrs30153022->SetMinimum(0.01287542);
Graph_Graph_ptj_toterrs30153022->SetMaximum(0.4300631);
Graph_Graph_ptj_toterrs30153022->SetDirectory(0);
Graph_Graph_ptj_toterrs30153022->SetStats(0);
Graph_Graph_ptj_toterrs30153022->SetLineWidth(2);
Graph_Graph_ptj_toterrs30153022->SetMarkerStyle(0);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetNdivisions(4);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelOffset(999);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetLabelSize(0);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetXaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetYaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30153022->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_ptj_toterrs30153022);
grae->Draw("p1p");
Double_t _fx3023[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t _fy3023[5] = {
0.3671394,
0.2404275,
0.1334798,
0.06537649,
0.01594908};
Double_t _felx3023[5] = {
0,
0,
0,
0,
0};
Double_t _fely3023[5] = {
0.01985385,
0.01348796,
0.006928003,
0.003412992,
0.001084033};
Double_t _fehx3023[5] = {
0,
0,
0,
0,
0};
Double_t _fehy3023[5] = {
0.01921597,
0.01338518,
0.006596172,
0.003531927,
0.001381032};
grae = new TGraphAsymmErrors(5,_fx3023,_fy3023,_felx3023,_fehx3023,_fely3023,_fehy3023);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph_Graph30163023 = new TH1F("Graph_Graph_Graph30163023","",100,17.325,88.425);
Graph_Graph_Graph30163023->SetMinimum(0.01337855);
Graph_Graph_Graph30163023->SetMaximum(0.4235044);
Graph_Graph_Graph30163023->SetDirectory(0);
Graph_Graph_Graph30163023->SetStats(0);
Graph_Graph_Graph30163023->SetLineWidth(2);
Graph_Graph_Graph30163023->SetMarkerStyle(0);
Graph_Graph_Graph30163023->GetXaxis()->SetNdivisions(4);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelOffset(999);
Graph_Graph_Graph30163023->GetXaxis()->SetLabelSize(0);
Graph_Graph_Graph30163023->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph30163023->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_Graph30163023->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30163023->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph30163023->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph30163023->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30163023->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30163023->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph30163023);
grae->Draw("p1p");
Double_t ptj_toterrs_fx3024[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t ptj_toterrs_fy3024[5] = {
0.389299,
0.2509263,
0.1443755,
0.07365244,
0.01874484};
Double_t ptj_toterrs_felx3024[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fely3024[5] = {
0.0283978,
0.01982285,
0.01219084,
0.00636368,
0.001684181};
Double_t ptj_toterrs_fehx3024[5] = {
0,
0,
0,
0,
0};
Double_t ptj_toterrs_fehy3024[5] = {
0.03041728,
0.02088176,
0.0127053,
0.007720362,
0.001934453};
grae = new TGraphAsymmErrors(5,ptj_toterrs_fx3024,ptj_toterrs_fy3024,ptj_toterrs_felx3024,ptj_toterrs_fehx3024,ptj_toterrs_fely3024,ptj_toterrs_fehy3024);
grae->SetName("ptj_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_ptj_toterrs30173024 = new TH1F("Graph_Graph_ptj_toterrs30173024","",100,17.175,72.075);
Graph_Graph_ptj_toterrs30173024->SetMinimum(0.01535459);
Graph_Graph_ptj_toterrs30173024->SetMaximum(0.4599818);
Graph_Graph_ptj_toterrs30173024->SetDirectory(0);
Graph_Graph_ptj_toterrs30173024->SetStats(0);
Graph_Graph_ptj_toterrs30173024->SetLineWidth(2);
Graph_Graph_ptj_toterrs30173024->SetMarkerStyle(0);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetNdivisions(4);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelOffset(999);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetLabelSize(0);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetXaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetYaxis()->SetTitleFont(132);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetLabelFont(132);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_ptj_toterrs30173024->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_ptj_toterrs30173024);
grae->Draw("p1p");
Double_t _fx3025[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t _fy3025[5] = {
0.389299,
0.2509263,
0.1443755,
0.07365244,
0.01874484};
Double_t _felx3025[5] = {
0,
0,
0,
0,
0};
Double_t _fely3025[5] = {
0.02243218,
0.01539977,
0.009090877,
0.00492242,
0.001303763};
Double_t _fehx3025[5] = {
0,
0,
0,
0,
0};
Double_t _fehy3025[5] = {
0.02411537,
0.01689699,
0.0100922,
0.005544287,
0.001554499};
grae = new TGraphAsymmErrors(5,_fx3025,_fy3025,_felx3025,_fehx3025,_fely3025,_fehy3025);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph30183025 = new TH1F("Graph_Graph_Graph30183025","",100,17.175,72.075);
Graph_Graph_Graph30183025->SetMinimum(0.01569697);
Graph_Graph_Graph30183025->SetMaximum(0.4530117);
Graph_Graph_Graph30183025->SetDirectory(0);
Graph_Graph_Graph30183025->SetStats(0);
Graph_Graph_Graph30183025->SetLineWidth(2);
Graph_Graph_Graph30183025->SetMarkerStyle(0);
Graph_Graph_Graph30183025->GetXaxis()->SetNdivisions(4);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelOffset(999);
Graph_Graph_Graph30183025->GetXaxis()->SetLabelSize(0);
Graph_Graph_Graph30183025->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph30183025->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
Graph_Graph_Graph30183025->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30183025->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph30183025->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph30183025->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph30183025->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph30183025->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph30183025);
grae->Draw("p1p");
TLegend *leg = new TLegend(0.65,0.57,0.95,0.82,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("data_stat","Data_{stat}","f");
ci = TColor::GetColor("#ff9933");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("data_tot","Data_{tot}","f");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("powheg_temp","POWHEG","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("amcatnlo_temp","aMC@NLO","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(24);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.5,0.85,0.85,0.94,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb, #sqrt{s} = 8 TeV");
pt->Draw();
Double_t xAxis21[6] = {20, 25, 30, 40, 50, 100};
TH1F *z_ptj_toterr_copy__21 = new TH1F("z_ptj_toterr_copy__21","z_ptj",5, xAxis21);
z_ptj_toterr_copy__21->SetBinContent(1,0.3608738);
z_ptj_toterr_copy__21->SetBinContent(2,0.225932);
z_ptj_toterr_copy__21->SetBinContent(3,0.1253807);
z_ptj_toterr_copy__21->SetBinContent(4,0.06703913);
z_ptj_toterr_copy__21->SetBinContent(5,0.01561868);
z_ptj_toterr_copy__21->SetBinContent(6,0.780934);
z_ptj_toterr_copy__21->SetBinError(1,0.02073388);
z_ptj_toterr_copy__21->SetBinError(2,0.01400757);
z_ptj_toterr_copy__21->SetBinError(3,0.009350811);
z_ptj_toterr_copy__21->SetBinError(4,0.004623786);
z_ptj_toterr_copy__21->SetBinError(5,0.001395185);
z_ptj_toterr_copy__21->SetBinError(6,0.06975923);
z_ptj_toterr_copy__21->SetEntries(12);
z_ptj_toterr_copy__21->SetDirectory(0);
z_ptj_toterr_copy__21->SetStats(0);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetLineColor(ci);
z_ptj_toterr_copy__21->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_ptj_toterr_copy__21->SetMarkerColor(ci);
z_ptj_toterr_copy__21->SetMarkerStyle(0);
z_ptj_toterr_copy__21->GetXaxis()->SetNdivisions(4);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelOffset(999);
z_ptj_toterr_copy__21->GetXaxis()->SetLabelSize(0);
z_ptj_toterr_copy__21->GetXaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetXaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d#it{p}_{T}^{jet} [pb/GeV]");
z_ptj_toterr_copy__21->GetYaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetYaxis()->SetLabelSize(0.05);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleOffset(1.3);
z_ptj_toterr_copy__21->GetYaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->GetZaxis()->SetLabelFont(132);
z_ptj_toterr_copy__21->GetZaxis()->SetLabelSize(0.05);
z_ptj_toterr_copy__21->GetZaxis()->SetTitleSize(0.06);
z_ptj_toterr_copy__21->GetZaxis()->SetTitleFont(132);
z_ptj_toterr_copy__21->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: lowerPad
TPad *lowerPad = new TPad("lowerPad", "lowerPad",0.005,0.05,0.995,0.995);
lowerPad->Draw();
lowerPad->cd();
lowerPad->Range(6.172839,-0.1935898,104.9383,4.832052);
lowerPad->SetFillColor(0);
lowerPad->SetFillStyle(4000);
lowerPad->SetBorderMode(0);
lowerPad->SetBorderSize(2);
lowerPad->SetTickx(1);
lowerPad->SetTicky(1);
lowerPad->SetLeftMargin(0.14);
lowerPad->SetRightMargin(0.05);
lowerPad->SetTopMargin(0.665);
lowerPad->SetBottomMargin(0.14);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
Double_t xAxis22[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp0__22 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp0__22","z_ptj",5, xAxis22);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(1,0.05745466);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(2,0.06199906);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(3,0.07457933);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(4,0.06897145);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(5,0.08932795);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetBinError(6,0.06975923);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleOffset(1.3);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0__22->Draw("e2");
Double_t xAxis23[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp1__23 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp1__23","z_ptj",5, xAxis23);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(1,0.02019378);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(2,0.02506611);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(3,0.02354391);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(4,0.03270951);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(5,0.03027217);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetBinError(6,0.02364057);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetStats(0);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp1__23->Draw("e2same");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fx3026[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fy3026[5] = {
1.017362,
1.064159,
1.064596,
0.9751989,
1.021154};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_felx3026[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fely3026[5] = {
0.08307013,
0.1014988,
0.07351357,
0.06717222,
0.1051988};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fehx3026[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp2_fehy3026[5] = {
0.06962999,
0.07751779,
0.06933948,
0.07789628,
0.1122927};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp2_fx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fy3026,zj_ptj_8TeV_400_1001_0ratiocomp2_felx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fehx3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fely3026,zj_ptj_8TeV_400_1001_0ratiocomp2_fehy3026);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp2");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026","",100,17.325,88.425);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp23026);
grae->Draw("p1");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fx3027[5] = {
23.25,
28.25,
36.5,
46.5,
82.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fy3027[5] = {
1.017362,
1.064159,
1.064596,
0.9751989,
1.021154};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_felx3027[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fely3027[5] = {
0.05501606,
0.05969921,
0.05525572,
0.05091044,
0.06940621};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fehx3027[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp3_fehy3027[5] = {
0.05324843,
0.05924428,
0.05260914,
0.05268456,
0.08842185};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp3_fx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fy3027,zj_ptj_8TeV_400_1001_0ratiocomp3_felx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fehx3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fely3027,zj_ptj_8TeV_400_1001_0ratiocomp3_fehy3027);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp3");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027","",100,17.325,88.425);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp33027);
grae->Draw("p1");
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fx3028[5] = {
21.75,
26.75,
33.5,
43.5,
67.5};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fy3028[5] = {
1.078768,
1.110628,
1.151497,
1.098648,
1.200155};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_felx3028[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fely3028[5] = {
0.07869176,
0.08773815,
0.0972306,
0.09492486,
0.1078312};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fehx3028[5] = {
0,
0,
0,
0,
0};
Double_t zj_ptj_8TeV_400_1001_0ratiocomp4_fehy3028[5] = {
0.08428785,
0.09242501,
0.1013337,
0.115162,
0.1238551};
grae = new TGraphAsymmErrors(5,zj_ptj_8TeV_400_1001_0ratiocomp4_fx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fy3028,zj_ptj_8TeV_400_1001_0ratiocomp4_felx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fehx3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fely3028,zj_ptj_8TeV_400_1001_0ratiocomp4_fehy3028);
grae->SetName("zj_ptj_8TeV_400_1001_0ratiocomp4");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028 = new TH1F("Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028","",100,17.175,72.075);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMinimum(0.51);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMaximum(1.49);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetDirectory(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetStats(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetLineWidth(2);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->SetMarkerStyle(0);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetNdivisions(4);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetXaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitle("Ratio");
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->CenterTitle(true);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetNdivisions(505);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetYaxis()->SetTitleFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetLabelFont(132);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetLabelSize(0.05);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetTitleSize(0.06);
Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_ptj_8TeV_400_1001_0ratiocomp43028);
grae->Draw("p1");
Double_t xAxis24[6] = {20, 25, 30, 40, 50, 100};
TH1F *zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24 = new TH1F("zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24","z_ptj",5, xAxis24);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(1,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(2,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(3,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(4,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(5,1);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinContent(6,0.780934);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(1,0.05745466);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(2,0.06199906);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(3,0.07457933);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(4,0.06897145);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(5,0.08932795);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetBinError(6,0.06975923);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMinimum(0.51);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMaximum(1.49);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetEntries(17);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetDirectory(0);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetLineColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMarkerColor(ci);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->SetMarkerStyle(0);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitle("#it{p}_{T}^{jet} [GeV]");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetNdivisions(4);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelOffset(0.02);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetXaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitle("Ratio");
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->CenterTitle(true);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetNdivisions(505);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTickLength(0.12325);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleOffset(1.3);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetYaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetLabelFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetLabelSize(0.05);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetTitleSize(0.06);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->GetZaxis()->SetTitleFont(132);
zj_ptj_8TeV_400_1001_0ratiocomp0_copy__24->Draw("sameaxis");
lowerPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
void forPaper_approval_2(std::string var1="h_ystar",
bool logScale=false,
std::string datafile="darko_root/goldenWithMCerror_withJESCorr_bigmatrix_corr.root",
std::string mcfile1="darko_root/bare_exclusive1Jet_zPt40_both_dressed_DYToLL_M-50_1jEnh2_2jEnh35_3jEnh40_4jEnh50_7TeV-sherpa.root",
std::string mcfile2="darko_root/bare_exclusive1Jet_zPt40_both_dressed_DYJetsToLL_TuneZ2_M-50_7TeV-madgraph-tauola.root",
std::string mcfile3="unified_angular_distributions/rebinnings/zpt40/Z_1jet_tota_cteq66_1___1___ex_m34.root",
std::string var2="",
std::string headertitle="Z#rightarrow ll + 1 jet",
std::string dataName="Data",
std::string mcName1="Sherpa",
std::string mcName2="Madgraph",
std::string mcName3="MCFM"
)
{
setTDRStyle();
gStyle->SetOptStat(0);
const double LABELSIZE = 20.0;
const int LINEWIDTH=3;
const int NHISTOS=4;
TH1F* h[NHISTOS];
char tempName[300];
if(var2 == "" )var2=var1;
std::string xtitle;
std::string output;
std::string datavar;
std::string var3;
std::string theoryName;
std::string remword3 ="h_";
std::string corrName = var1;
size_t pos3 = corrName.find(remword3);
if(pos3!= std::string::npos)
corrName.replace(pos3,remword3.length(),"");
datavar = "h_combine_" + corrName;
double ymax=1.22;
double ymin=0.78;
double xmax=3.0;
double xmin=0.0;
if(var1=="h_ystar")
{
var3 = "id4";
xtitle = "Y_{dif} = |Y_{Z}-Y_{jet}| / 2";
output = "DifYAll";
theoryName = "Ydif";
ymax = 1.45;
ymin = 0.55;
xmax = 1.7999;
}
else if(var1=="h_yB")
{
var3 = "id3";
xtitle = "Y_{sum} = |Y_{Z}+Y_{jet}| / 2";
output = "SumYAll";
theoryName = "Ysum";
ymax = 1.45;
ymin = 0.55;
xmax = 2.1999;
}
else if(var1=="h_jety")
{
var3 = "id1";
xtitle = "|Y_{jet}|";
output = "YJetAll";
datafile = "darko_root/goldenWithMCerror_withJESCorr_jety_bigmatrix.root";
theoryName = "Yjet";
xmax = 2.3999;
}
else if(var1=="h_zy")
{
var3 = "id2";
xtitle = "|Y_{Z}|";
output = "YZedAll";
theoryName = "Yzed";
xmax = 2.1999;
}
// first get the histogram files
TFile *fdata = TFile::Open(datafile.data());
cout << "Opening " << fdata->GetName() << endl;
TFile *fmc1 = TFile::Open(mcfile1.data());
cout << "Opening " << fmc1->GetName() << endl;
TFile *fmc2 = TFile::Open(mcfile2.data());
cout << "Opening " << fmc2->GetName() << endl;
TFile *fmc3 = TFile::Open(mcfile3.data());
cout << "Opening " << fmc3->GetName() << endl;
h[0] = (TH1F*)(fdata->Get(datavar.data()));
h[1] = (TH1F*)(fmc1->Get(var1.data()));
h[2] = (TH1F*)(fmc2->Get(var2.data()));
h[3] = (TH1F*)(fmc3->Get(var3.data()));
if(var1=="h_jety")
{
double value = h[3]->GetBinContent(12);
value *= 0.9;
h[3]->SetBinContent(12,value);
}
TH1D* hscale[NHISTOS];
int COLOR[NHISTOS]={1,4,2,kOrange-1};
int MARKERSTYLE[NHISTOS]={8,24,21,29};
int MARKERSIZE[NHISTOS]={1,0,0,0};
int LINESTYLE[NHISTOS]={1,1,2,6};
int FILLSTYLE[NHISTOS]={1,3345,3436,1};
for(int i=0; i < NHISTOS; i++){
hscale[i] =(TH1D*) h[0]->Clone(Form("hscale%02i",i));
hscale[i] ->SetYTitle(Form("Ratio to %s",mcName3.data()));
hscale[i] ->SetXTitle(xtitle.data());
hscale[i] ->GetXaxis()->SetDecimals();
hscale[i] ->GetYaxis()->SetDecimals();
hscale[i]->SetLineColor(COLOR[i]);
hscale[i]->SetLineWidth(LINEWIDTH);
hscale[i]->SetLineStyle(1);
hscale[i]->SetMarkerColor(COLOR[i]);
hscale[i]->SetMarkerStyle(MARKERSTYLE[i]);
hscale[i]->SetMarkerSize(MARKERSIZE[i]);
hscale[i]->SetFillColor(COLOR[i]);
hscale[i]->SetFillStyle(FILLSTYLE[i]);
hscale[i]->SetTitle("");
hscale[i]->SetMaximum(ymax);
hscale[i]->SetMinimum(ymin);
hscale[i]->SetTitleOffset(1.2,"X");
hscale[i]->SetTitleOffset(1.2,"Y");
h[i]->SetTitle("");
h[i]->SetLineStyle(LINESTYLE[i]);
h[i]->GetXaxis()->SetDecimals();
h[i]->GetYaxis()->SetDecimals();
h[i]->SetMarkerSize(1);
h[i]->SetLineColor(COLOR[i]);
h[i]->SetLineWidth(LINEWIDTH);
h[i]->SetMarkerColor(COLOR[i]);
h[i]->SetMarkerStyle(MARKERSTYLE[i]);
h[i]->SetTitleOffset(1.2,"Y");
}
h[0]->SetLineWidth(1);
hscale[0]->SetLineWidth(1);
// if normalizing to the same area, set the scale
int binLo = -1;
int binHi = -1;
int nbins = h[0]->GetNbinsX();
binLo = 1;
binHi = nbins;
double scaleFactor[NHISTOS]={1};
for(int ih=0; ih < NHISTOS; ih++){
double integral = h[ih]->Integral(binLo,binHi);
scaleFactor[ih] = integral > 0? 1.0/integral: 1;
h[ih]->Sumw2();
h[ih]->Scale(scaleFactor[ih]);
}
for(int ih=0; ih < NHISTOS; ih++)
cout << "histogram " << ih << " integral = " << h[ih]->Integral() << endl;
// get the ratio
for(int ih=0; ih < NHISTOS-1; ih++){
cout << "===================================================" << endl;
cout << "For histogram " << ih << endl;
hscale[ih]->Divide(h[ih], h[NHISTOS-1]);
hscale[ih]->SetMaximum(ymax);
hscale[ih]->SetMinimum(ymin);
} // end of loop over histograms
vector<double> maxArray;
maxArray.clear();
for(int ih=0; ih < NHISTOS; ih++){
h[ih]->GetXaxis()->SetRangeUser(xmin,xmax);
hscale[ih]->GetXaxis()->SetRangeUser(xmin,xmax);
double max_this = h[ih]->GetBinError(h[ih]->GetMaximumBin()) + h[ih]->GetMaximum();
maxArray.push_back(max_this);
}
double max = *(std::max_element(maxArray.begin(),maxArray.end()));
cout << "Max = " << max << endl;
for(int ih=0; ih < NHISTOS; ih++){
h[ih]->SetMaximum(1.1*max);
if(!logScale)
h[ih]->SetMinimum(-0.015);
else
h[ih]->SetMinimum(5e-6);
}
//////////////////////////////////////////////////////////////////////////
///
/// Making final figures and save the canvas in files
///
//////////////////////////////////////////////////////////////////////////
TCanvas* c1 = new TCanvas("c1","",700,1000);
c1->Divide(1,2,0.01,0);
c1->cd(1);
if(logScale)
gPad->SetLogy(1);
gPad->SetTopMargin(0.01);
gPad->SetBottomMargin(0);
gPad->SetRightMargin(0.04);
gStyle->SetStatFontSize(0.05);
// pad 1 and pad 2 have different sizes, need to rescale to
// make the label size the same
double temp1_pad = gPad->GetWh()*gPad->GetAbsHNDC();
cout << "pad 1 size in pixels = " << temp1_pad << endl;
double label1_size = LABELSIZE/temp1_pad;
for(int ih=0; ih < NHISTOS; ih++)
{
h[ih]->GetYaxis()->SetLabelSize(label1_size);
h[ih]->GetXaxis()->SetLabelSize(label1_size);
}
h[0]->SetYTitle("1/#sigma d#sigma/dY");
h[0]->Draw("9e1");
for(int ih=1; ih < NHISTOS; ih++)
h[ih]->Draw("9histsame");
h[0]->Draw("9e1same");
double x1NDC = 0.729812;
double y1NDC = 0.560219;
double x2NDC = 0.939633;
double y2NDC = 0.871885;
TLegend* leg = new TLegend(x1NDC,y1NDC,x2NDC,y2NDC);
leg->SetFillColor(0);
leg->SetFillStyle(0);
leg->SetTextSize(0.045);
leg->SetBorderSize(0);
leg->SetHeader(headertitle.data());
leg->AddEntry(h[0], dataName.data());
leg->AddEntry(h[1], mcName1.data(),"l");
leg->AddEntry(h[2], mcName2.data(),"l");
leg->AddEntry(h[3], mcName3.data(),"l");
leg->Draw("same");
TLatex *lar = new TLatex(0.33, 0.91, "CMS Preliminary, #sqrt{s} = 7 TeV, L_{int} = 5 fb^{-1}");
lar->SetNDC(kTRUE);
lar->SetTextSize(0.05);
lar->Draw();
TPaveText *pavetex = new TPaveText(0.17029,
0.0495665,
0.478939,
0.483501,"NDCBR");
pavetex->SetBorderSize(0);
pavetex->SetFillColor(0);
pavetex->SetFillStyle(0);
pavetex->SetLineWidth(3);
pavetex->SetTextAlign(12);
pavetex->SetTextSize(0.04);
pavetex->AddText("76 < M_{ll} < 106 GeV");
pavetex->AddText("p_{T}^{ll} > 40 GeV");
pavetex->AddText("p_{T}^{l} > 20 GeV, |#eta^{l}| < 2.1");
pavetex->AddText("N_{jet}=1, p_{T}^{jet} > 30 GeV, |#eta^{jet}| < 2.4");
pavetex->AddText("#Delta R(l,jet)>0.5");
pavetex->Draw();
cout << "pad 1 label size = " << h[0]->GetYaxis()->GetLabelSize() << endl;
//////////////////////////////////////////////////////////////////////
c1->cd(2);
gPad->SetRightMargin(0.04);
gPad->SetTopMargin(0);
gPad->SetBottomMargin(0.2);
gPad->SetTickx();
// pad 1 and pad 2 have different sizes, need to rescale to
// make the label size the same
double temp2_pad = gPad->GetWh()*gPad->GetAbsHNDC();
cout << "pad 2 size in pixels = " << temp2_pad << endl;
double label2_size = LABELSIZE/temp2_pad;
for(int ih=0; ih < NHISTOS; ih++)
{
hscale[ih]->GetYaxis()->SetLabelSize(label2_size);
hscale[ih]->GetXaxis()->SetLabelSize(label2_size);
}
hscale[0]->Draw("9e1");
for(int ih=1; ih < NHISTOS-1; ih++){
hscale[ih]->Draw("9e3same");
}
hscale[0]->Draw("9e1same");
TLine* l2 = new TLine(xmin,1.,xmax,1.);
l2->SetLineColor(kOrange-1);
l2->SetLineStyle(1);
l2->Draw("same");
TLegend* leg2 = new TLegend(0.17333,0.244648,0.392274,0.476857);
leg2->SetFillColor(0);
leg2->SetFillStyle(0);
leg2->SetTextSize(0.04);
leg2->SetBorderSize(0);
leg2->AddEntry(hscale[1], "Sherpa stat. uncertainty","f");
leg2->AddEntry(hscale[2], "Madgraph stat. uncertainty","f");
leg2->Draw("same");
gROOT->ProcessLine(".L ~/scripts/theoryErrorZed.c");
theoryErrorZed(theoryName.data());
hscale[0]->Draw("9e1same");
for(int ih=1; ih < NHISTOS-1; ih++){
hscale[ih]->Draw("9e3same");
}
hscale[0]->Draw("9e1same");
l2->Draw("same");
cout << "pad 2 label size = " << hscale[0]->GetYaxis()->GetLabelSize() << endl;
string dirName = "forPaper";
gSystem->mkdir(dirName.data());
std::string filename;
std::string psname = dirName + "/" + output;
filename = psname + ".eps";
c1->Print(filename.data());
filename = psname + ".gif";
c1->Print(filename.data());
filename = psname + ".pdf";
c1->Print(filename.data());
// c1->Close();
}
TCanvas* plotExpVersusMC(int numExpHists, TH1** expHists, int numMCHists, TH1** mcHists, TString titleString)
{
gROOT->cd();
int totalNumHists = numExpHists + numMCHists;
if(totalNumHists < 1 || (expHists == nullptr && mcHists == nullptr))
{
return nullptr;
}
const int numColors = 11;
const int colorList[numColors] = { kBlack, kBlue, kRed, kGreen, kCyan, kMagenta, kPink, kOrange, kViolet, kAzure, kYellow };
const int markerList[14] = { 20, 22, 23, 25, 24, 32, 33, 3, 5, 28, 26, 30, 29 };
//gStyle->SetTitleColor(kWhite);
TCanvas* c2 = new TCanvas(titleString, "plotManyHistCanvas", 1600, 1200);
c2->SetFillColor(kGray);
c2->cd();
c2->SetLeftMargin(.13);
c2->SetRightMargin(.06);
c2->SetGrid(1, 1);
gPad->SetTickx(1);
gPad->SetTicky(1);
gPad->SetFillColor(kWhite);
TLegend *theLegend;
if(totalNumHists >= 1)
{
double legendSize = .1 * totalNumHists;
if(legendSize < .15) legendSize = .15;
theLegend = new TLegend(0.55, .8 - legendSize, 0.87, 0.8); //Top Right
// theLegend = new TLegend(0.2,.8-legendSize,0.47,0.8); //Top Left
// theLegend = new TLegend(0.55,.4-legendSize,0.87,0.4); //Bottom Right
}
//Prep MCHists
for (int i = 0; i < numMCHists; i++)
{
mcHists[i]->SetLineColor(colorList[(i) % numColors]);
mcHists[i]->SetMarkerColor(colorList[(i) % numColors]);
mcHists[i]->SetLineStyle((i - numExpHists) % 5 + 1);
mcHists[i]->SetLineWidth(2);
mcHists[i]->GetYaxis()->SetTitleOffset(1.5);
}
//Prep MCHists
for (int i = 0; i < numExpHists; i++)
{
expHists[i]->SetLineColor(colorList[(i + numMCHists) % numColors]);
expHists[i]->SetMarkerColor(colorList[(i + numMCHists) % numColors]);
expHists[i]->SetMarkerStyle(markerList[(i) % 14]);
expHists[i]->SetLineWidth(2);
expHists[i]->SetMarkerSize(2);
expHists[i]->GetYaxis()->SetTitleOffset(1.5);
}
bool drawExpFirst = true;
for (int i = 0; i < totalNumHists; i++)
{
if(drawExpFirst && numExpHists > 0)
{
if(i < numExpHists)
{
theLegend->AddEntry(expHists[i], expHists[i]->GetTitle(), "PE");
if(i == 0)
{
expHists[i]->SetTitle("");
expHists[i]->Draw("P E1 X0");
}
else
{
expHists[i]->Draw("P E1 X0 SAME");
}
}
else
{
theLegend->AddEntry(mcHists[i - numExpHists], mcHists[i - numExpHists]->GetTitle(), "L");
mcHists[i - numExpHists]->Draw("HIST SAME");
}
}
else
{
if(i < numMCHists)
{
theLegend->AddEntry(mcHists[i], mcHists[i]->GetTitle(), "L");
if(i == 0)
{
mcHists[i]->SetTitle("");
mcHists[i]->Draw("HIST");
}
else
{
mcHists[i]->Draw("HIST SAME");
}
}
else
{
theLegend->AddEntry(expHists[i - numMCHists], expHists[i - numMCHists]->GetTitle(), "PE");
expHists[i - numMCHists]->Draw("P E1 X0 SAME");
}
}
}
if(totalNumHists >= 1)
{
theLegend->Draw();
theLegend->SetBorderSize(1);
theLegend->SetFillColor(kWhite);
}
TPaveText* textBox = new TPaveText(.14, .945, .8705314, .9825073, "NDC");
textBox->SetTextSize(.035);
textBox->AddText(titleString);
textBox->SetFillColor(kWhite);
textBox->SetLineColor(kWhite);
textBox->SetShadowColor(kWhite);
// textBox->SetFillStyle(4000);//transparent
// textBox->SetLineStyle(4000);//transparent
textBox->SetLineWidth(0);
textBox->Draw();
c2->Update();
return c2;
}
void SingleAnalysis::PostDrawMultiGraph(GraphName graphName, TMultiGraph *pMultiGraph) const
{
if(NULL == pMultiGraph)
return;
switch(graphName)
{
case E_REC:
{
pMultiGraph->GetYaxis()->SetLabelFont(42);
pMultiGraph->GetYaxis()->SetTitleSize(0.07);
pMultiGraph->GetYaxis()->SetTitleOffset(0.8);
pMultiGraph->GetYaxis()->SetTitleFont(42);
pMultiGraph->GetYaxis()->SetLabelSize(0.045);
break;
}
case E_REC_DEVIATION:
{
pMultiGraph->GetXaxis()->SetTitle("E_{beam} [GeV]");
pMultiGraph->GetXaxis()->SetLabelFont(42);
pMultiGraph->GetXaxis()->SetLabelSize(0.07);
pMultiGraph->GetXaxis()->SetTitleSize(0.1);
pMultiGraph->GetXaxis()->SetTitleOffset(1);
pMultiGraph->GetXaxis()->SetTitleFont(42);
pMultiGraph->GetYaxis()->SetLabelFont(42);
pMultiGraph->GetYaxis()->SetTitleSize(0.12);
pMultiGraph->GetYaxis()->SetTitleOffset(0.5);
pMultiGraph->GetYaxis()->SetTitleFont(42);
pMultiGraph->GetYaxis()->SetLabelSize(0.07);
break;
}
default:
{
pMultiGraph->GetXaxis()->SetTitle("E_{beam} [GeV]");
pMultiGraph->GetXaxis()->SetLabelFont(42);
pMultiGraph->GetXaxis()->SetTitleSize(0.05);
pMultiGraph->GetXaxis()->SetTitleOffset(1.);
pMultiGraph->GetXaxis()->SetTitleFont(42);
pMultiGraph->GetXaxis()->SetRangeUser(0, 90);
pMultiGraph->GetYaxis()->SetTitle("");
pMultiGraph->GetYaxis()->SetLabelFont(42);
pMultiGraph->GetYaxis()->SetTitleSize(0.045);
pMultiGraph->GetYaxis()->SetTitleOffset(1.3);
pMultiGraph->GetYaxis()->SetTitleFont(42);
pMultiGraph->GetYaxis()->SetLabelSize(0.035);
break;
}
}
if(graphName != E_REC_DEVIATION)
{
TPaveText *pt = new TPaveText(0.3, 0.2, 0.53, 0.3, "tbNDC");
pt->SetTextSize(0.05);
pt->SetTextColor(kGray+2);
pt->SetFillColor(0);
pt->SetLineWidth(0);
pt->SetBorderSize(0);
pt->AddText("CALICE SDHCAL");
pt->SetTextFont(62);
pt->Draw();
}
// plot per plot customize
switch(graphName)
{
case E_REC:
pMultiGraph->GetYaxis()->SetTitle("E_{rec}");
pMultiGraph->GetYaxis()->SetRangeUser(1, 90);
break;
case E_REC_DEVIATION:
pMultiGraph->GetYaxis()->SetTitle("#DeltaE/E_{beam}");
pMultiGraph->GetYaxis()->SetRangeUser(-0.12, 0.12);
break;
case E_RESOL:
pMultiGraph->GetYaxis()->SetTitle("#sigma_{E}/E_{rec}");
pMultiGraph->GetYaxis()->SetRangeUser(0, 0.4);
break;
case N_PFOS:
pMultiGraph->GetYaxis()->SetTitle("<N_{pfos}>");
pMultiGraph->GetYaxis()->SetRangeUser(0, 2);
break;
case N_PFOS_20GEV:
case N_PFOS_70GEV:
pMultiGraph->GetYaxis()->SetTitle("# events (normalized to unity)");
pMultiGraph->GetYaxis()->SetRangeUser(0, 1);
pMultiGraph->GetXaxis()->SetRangeUser(-0.99, 5.99);
break;
case EFFICIENCY:
pMultiGraph->GetYaxis()->SetTitle("Hit clustering efficiency #epsilon_{s}");
pMultiGraph->GetYaxis()->SetRangeUser(0.8, 1.1);
break;
default:
break;
}
}
void Overlay_NeutralEnergyDifferenceMean()
{
std::string canvasName = "canvas";
std::string canvasTitle = "Neutral energy mean";
int canvasX1 = 200;
int canvasY1 = 52;
int canvasX2 = 700;
int canvasY2 = 650;
int nPoints = 6;
int nGraphs = 5;
int neutralEnergy = 10;
double zero [] = {0, 0, 0, 0, 0, 0, 0, 0};
double energies[] = {10, 20, 30, 40, 50};
double distances[] = {5, 10, 15, 20, 25, 30};
std::string graphXTitle = "Distance between showers [cm]";
std::string graphYTitle = "<E_{n,rec} - E_{n,meas}>";
double rangeUserXMin = 0;
double rangeUserXMax = 35;
double rangeUserYMin = -10;
double rangeUserYMax = 10;
TCanvas *pCanvas = new TCanvas(canvasName.c_str(), canvasTitle.c_str(), canvasX1, canvasY1, canvasX2, canvasY2);
pCanvas->SetFillColor(0);
pCanvas->SetBorderMode(0);
pCanvas->SetBorderSize(2);
pCanvas->SetTickx(1);
pCanvas->SetTicky(1);
pCanvas->SetLeftMargin(0.15);
pCanvas->SetRightMargin(0.03);
pCanvas->SetTopMargin(0.05);
pCanvas->SetBottomMargin(0.14);
pCanvas->SetFrameBorderMode(0);
pCanvas->SetFrameBorderMode(0);
TH1F *pH = new TH1F("pouette", "", rangeUserXMax-rangeUserXMin, rangeUserXMin, rangeUserXMax);
pH->SetMinimum(rangeUserYMin);
pH->SetMaximum(rangeUserYMax);
pH->Draw();
gStyle->SetOptStat(0);
pH->GetYaxis()->SetTitle(graphYTitle.c_str());
pH->GetXaxis()->SetTitle(graphXTitle.c_str());
pH->GetYaxis()->SetLabelFont(42);
pH->GetYaxis()->SetTitleSize(0.05);
pH->GetYaxis()->SetTitleOffset(1);
pH->GetYaxis()->SetTitleFont(42);
pH->GetYaxis()->SetLabelSize(0.045);
TLegend *leg = new TLegend(0.3,0.55,0.8,0.8,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetLineColor(0);
leg->SetLineStyle(1);
leg->SetLineWidth(1);
leg->SetFillColor(0);
leg->SetTextSize(0.03);
int currentColor = 1;
for(unsigned int e=0 ; e<nGraphs ; e++)
{
int energy = energies[e];
double *data = new double[nPoints];
double *dataError = zero;//new double[nPoints];
std::stringstream fileName;
fileName << "data_neutralERecDifference_ArborPFA_Test\ Beam_" << energy << "Gev.txt";
readData(fileName.str(), data, nPoints);
// fileName.str("");
// fileName << "fichier_data_error_energy_" << energy << "GeV.txt";
// readData(fileName.str(), dataError, nPoints);
TGraphErrors *gre = 0;
std::stringstream graphName;
graphName << "Charged particle energy = " << energy << " GeV";
gre = new TGraphErrors(nPoints, distances, data, zero, dataError);
gre->SetName(graphName.str().c_str());
gre->SetTitle(graphName.str().c_str());
gre->SetLineColor(currentColor);
gre->SetMarkerColor(currentColor);
gre->SetMarkerStyle(23);
gre->SetMarkerSize(1.2);
gre->Draw("lp same");
TLegendEntry *entry = 0;
entry=leg->AddEntry(gre, graphName.str().c_str(), "lp");
entry->SetFillStyle(1001);
entry->SetMarkerColor(currentColor);
entry->SetMarkerStyle(23);
entry->SetMarkerSize(1);
currentColor++;
}
TPaveText *pt = new TPaveText(0.55, 0.2, 0.93, 0.3, "tbNDC");
pt->SetTextSize(0.05);
pt->SetTextColor(kGray+2);
pt->SetFillColor(0);
pt->SetLineWidth(0);
pt->SetBorderSize(0);
pt->AddText("CALICE Preliminary");
pt->SetTextFont(62);
pt->Draw();
leg->Draw();
pCanvas->Modified();
pCanvas->cd();
pCanvas->SetSelected(pCanvas);
}
void compareThreeVersions(TimeInfoAllSteps TIAS_1,TimeInfoAllSteps TIAS_2,TimeInfoAllSteps TIAS_3,TString file_postfix) {
gStyle->SetOptStat(0);
TPaveText* labelcms = new TPaveText(0.20,0.79,0.54,0.92,"NDCBR");
labelcms->SetTextAlign(12);
labelcms->SetTextSize(0.04);
labelcms->SetFillColor(kWhite);
labelcms->AddText("CMS Preliminary Simulation");
//labelcms->AddText("CMS Preliminary");
//labelcms->AddText("Simulation");
labelcms->AddText("#sqrt{s} = 8 TeV, t#bar{t}+PU");
labelcms->SetBorderSize(0);
labelcms->SetTextFont(42);
labelcms->SetLineWidth(2);
TH1F *htime_iter_1 = new TH1F("time_iter_1"+file_postfix,"time_iter_1",8,0,8);
timingTestPerIter(TIAS_1,htime_iter_1);
htime_iter_1->SetLineStyle(TIAS_1.style());
htime_iter_1->SetLineColor(TIAS_1.color());
TH1F *htime_iter_2 = new TH1F("time_iter_2"+file_postfix,"time_iter_2",8,0,8);
timingTestPerIter(TIAS_2,htime_iter_2);
htime_iter_2->SetLineStyle(TIAS_2.style());
htime_iter_2->SetLineColor(TIAS_2.color());
TH1F *htime_iter_3 = new TH1F("time_iter_3"+file_postfix,"time_iter_3",8,0,8);
timingTestPerIter(TIAS_3,htime_iter_3);
htime_iter_3->SetLineStyle(TIAS_3.style());
htime_iter_3->SetLineColor(TIAS_3.color());
float iter_Iter0PU20 = htime_iter_1->GetBinContent(1);
htime_iter_1->Scale(1./iter_Iter0PU20);
htime_iter_2->Scale(1./iter_Iter0PU20);
htime_iter_3->Scale(1./iter_Iter0PU20);
htime_iter_1->GetYaxis()->SetRangeUser(0,1.1*TMath::Max(htime_iter_1->GetBinContent(htime_iter_1->GetMaximumBin()),TMath::Max(htime_iter_2->GetBinContent(htime_iter_2->GetMaximumBin()),htime_iter_3->GetBinContent(htime_iter_3->GetMaximumBin()))));
TLegend* leg = new TLegend(0.2,0.59,0.45,0.79);
leg->SetNColumns(1);
leg->SetFillColor(kWhite);
leg->SetLineColor(kWhite);
leg->SetBorderSize(0);
leg->SetTextSize(0.04);
leg->SetTextFont(42);
leg->AddEntry(htime_iter_1,TIAS_1.legend(),"L");
leg->AddEntry(htime_iter_2,TIAS_2.legend(),"L");
leg->AddEntry(htime_iter_3,TIAS_3.legend(),"L");
htime_iter_1->GetYaxis()->SetTitleOffset(1.2);
htime_iter_1->GetXaxis()->SetRangeUser(0,7);
htime_iter_1->GetYaxis()->SetTitle("time [1/iter0@<PU>=20]");
htime_iter_1->GetYaxis()->SetTitleSize(0.04);
htime_iter_1->GetXaxis()->SetTitleSize(0.04);
htime_iter_1->GetYaxis()->SetLabelSize(0.04);
TCanvas c0;
c0.SetTicks(1,1);
htime_iter_1->Draw("H");
htime_iter_2->Draw("H,same");
htime_iter_3->Draw("H,same");
leg->Draw();
labelcms->Draw();
c0.SetGridy();
c0.SaveAs("timingNew_iter_"+file_postfix+".png");
TH1F *htime_step_1 = new TH1F("time_step_1"+file_postfix,"time_step_1",5,0,5);
timingTestPerStep(TIAS_1,htime_step_1);
htime_step_1->SetLineStyle(TIAS_1.style());
htime_step_1->SetLineColor(TIAS_1.color());
htime_step_1->GetYaxis()->SetTitle("time [1/building@<PU>=20]");
TH1F *htime_step_2 = new TH1F("time_step_2"+file_postfix,"time_step_2",5,0,5);
timingTestPerStep(TIAS_2,htime_step_2);
htime_step_2->SetLineStyle(TIAS_2.style());
htime_step_2->SetLineColor(TIAS_2.color());
TH1F *htime_step_3 = new TH1F("time_step_3"+file_postfix,"time_step_3",5,0,5);
timingTestPerStep(TIAS_3,htime_step_3);
htime_step_3->SetLineStyle(TIAS_3.style());
htime_step_3->SetLineColor(TIAS_3.color());
float step_BuildPU20 = htime_step_1->GetBinContent(3);
htime_step_1->Scale(1./step_BuildPU20);
htime_step_2->Scale(1./step_BuildPU20);
htime_step_3->Scale(1./step_BuildPU20);
htime_step_1->GetYaxis()->SetTitleOffset(1.2);
htime_step_1->GetXaxis()->SetRangeUser(0,7);
htime_step_1->GetYaxis()->SetTitleSize(0.04);
htime_step_1->GetXaxis()->SetTitleSize(0.04);
htime_step_1->GetYaxis()->SetLabelSize(0.04);
TCanvas c1;
c1.SetTicks(1,1);
//htime_step_1->Scale(1./htime_step_1->Integral());
//htime_step_2->Scale(1./htime_step_2->Integral());
//htime_step_3->Scale(1./htime_step_3->Integral());
htime_step_1->GetYaxis()->SetRangeUser(0,1.1*TMath::Max(htime_step_1->GetBinContent(htime_step_1->GetMaximumBin()),TMath::Max(htime_step_2->GetBinContent(htime_step_2->GetMaximumBin()),htime_step_3->GetBinContent(htime_step_3->GetMaximumBin()))));
htime_step_1->GetXaxis()->SetRangeUser(1,5);
htime_step_1->Draw("H");
htime_step_2->Draw("H,same");
htime_step_3->Draw("H,same");
leg->SetX1NDC(0.60);
leg->SetX2NDC(0.85);
labelcms->SetX1NDC(0.60);
labelcms->SetX2NDC(0.94);
leg->Draw();
labelcms->Draw();
c1.SetGridy();
c1.SaveAs("timingNew_step_"+file_postfix+".png");
TH1F *htime_track_1 = (TH1F*) htime_iter_1->Clone("time_track_1"+file_postfix);
htime_track_1->Reset();
makeTimePerTrackPlot(htime_track_1,htime_iter_1,TIAS_1.mtvfile());
TH1F *htime_track_2 = (TH1F*) htime_iter_2->Clone("time_track_2"+file_postfix);
htime_track_2->Reset();
makeTimePerTrackPlot(htime_track_2,htime_iter_2,TIAS_2.mtvfile());
TH1F *htime_track_3 = (TH1F*) htime_iter_3->Clone("time_track_3"+file_postfix);
htime_track_3->Reset();
makeTimePerTrackPlot(htime_track_3,htime_iter_3,TIAS_3.mtvfile());
float track_Iter0PU20 = htime_track_1->GetBinContent(1);
htime_track_1->Scale(1./track_Iter0PU20);
htime_track_2->Scale(1./track_Iter0PU20);
htime_track_3->Scale(1./track_Iter0PU20);
htime_track_1->GetYaxis()->SetRangeUser(0,1.1*TMath::Max(htime_track_1->GetBinContent(htime_track_1->GetMaximumBin()),TMath::Max(htime_track_2->GetBinContent(htime_track_2->GetMaximumBin()),htime_track_3->GetBinContent(htime_track_3->GetMaximumBin()))));
TCanvas c2;
c2.SetTicks(1,1);
htime_track_1->GetYaxis()->SetTitleOffset(1.2);
htime_track_1->GetXaxis()->SetRangeUser(0,7);
htime_track_1->GetYaxis()->SetTitle("time/track [1/iter0@<PU>=20]");
htime_track_1->GetYaxis()->SetTitleSize(0.04);
htime_track_1->GetXaxis()->SetTitleSize(0.04);
htime_track_1->GetYaxis()->SetLabelSize(0.04);
htime_track_1->Draw("H");
htime_track_2->Draw("H,same");
htime_track_3->Draw("H,same");
leg->SetX1NDC(0.2);
leg->SetX2NDC(0.45);
labelcms->SetX1NDC(0.2);
labelcms->SetX2NDC(0.54);
leg->Draw();
labelcms->Draw();
c2.SetGridy();
c2.SaveAs("timingNew_track_"+file_postfix+".png");
}
void wmj_fit_stack()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Tue May 3 07:35:30 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,700,500);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(-0.1730123,-5358.989,1.06279,32919.5);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.14);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3005[20] = {
0.025,
0.075,
0.125,
0.175,
0.225,
0.275,
0.325,
0.375,
0.425,
0.475,
0.525,
0.575,
0.625,
0.675,
0.725,
0.775,
0.825,
0.875,
0.925,
0.975};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3005[20] = {
0,
0,
9,
26,
118,
315,
881,
1812,
3086,
4575,
5856,
6537,
6824,
6348,
5662,
4769,
4941,
6726,
13179,
28019};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3005[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3005[20] = {
0,
0,
2.943461,
5.066015,
10.37428,
17.25528,
29.18586,
42.07053,
55.05403,
67.14059,
76.02614,
80.35326,
82.10902,
79.17591,
74.74792,
68.55976,
69.79403,
81.51372,
114.3009,
166.8895};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3005[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3005[20] = {
1.147874,
1.147874,
4.110204,
6.164324,
11.37428,
18.25528,
30.18586,
43.07053,
56.05403,
68.14059,
77.02614,
81.35326,
83.10902,
80.17591,
75.74792,
69.55976,
70.79403,
82.51372,
115.3009,
167.8895};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(20,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3005,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3005);
grae->SetName("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]");
grae->SetTitle("Histogram of combData_plot__muminus_ISO_");
grae->SetFillColor(1);
grae->SetFillStyle(0);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005 = new TH1F("Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMinimum(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetDirectory(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetStats(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetLineWidth(2);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]56300130023005);
grae->Draw("e1pa");
Double_t h5_fx7[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t h5_fy7[47] = {
0,
0,
0,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
8.630173e-06,
3.789762,
3.789762,
71.12063,
71.12063,
211.4394,
211.4394,
857.7007,
857.7007,
1798.534,
1798.534,
3099.609,
3099.609,
4425.901,
4425.901,
5601.584,
5601.584,
6533.687,
6533.687,
6910.208,
6910.208,
6392.171,
6392.171,
5777.193,
5777.193,
5020.919,
5020.919,
5207.834,
5207.834,
6798.159,
6798.159,
13033.98,
13033.98,
27939.27,
27939.27,
0,
0,
0,
0};
TGraph *graph = new TGraph(47,h5_fx7,h5_fy7);
graph->SetName("h5");
graph->SetTitle("Projection of ");
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#3333ff");
graph->SetFillColor(ci);
ci = TColor::GetColor("#3333ff");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#3333ff");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_h57 = new TH1F("Graph_h57","Projection of ",100,-0.1292683,1.129268);
Graph_h57->SetMinimum(0);
Graph_h57->SetMaximum(30733.2);
Graph_h57->SetDirectory(0);
Graph_h57->SetStats(0);
Graph_h57->SetLineWidth(2);
Graph_h57->SetMarkerStyle(0);
Graph_h57->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_h57->GetXaxis()->SetRange(11,90);
Graph_h57->GetXaxis()->SetNdivisions(1005);
Graph_h57->GetXaxis()->SetLabelFont(132);
Graph_h57->GetXaxis()->SetLabelOffset(0.02);
Graph_h57->GetXaxis()->SetLabelSize(0.05);
Graph_h57->GetXaxis()->SetTitleSize(0.06);
Graph_h57->GetXaxis()->SetTitleFont(132);
Graph_h57->GetYaxis()->SetTitle("Events / (0.05)");
Graph_h57->GetYaxis()->SetLabelFont(132);
Graph_h57->GetYaxis()->SetLabelSize(0.05);
Graph_h57->GetYaxis()->SetTitleSize(0.06);
Graph_h57->GetYaxis()->SetTitleFont(132);
Graph_h57->GetZaxis()->SetLabelFont(132);
Graph_h57->GetZaxis()->SetLabelSize(0.05);
Graph_h57->GetZaxis()->SetTitleSize(0.06);
Graph_h57->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_h57);
graph->Draw("f");
Double_t sumBkgPdf_totbkg_fitW_fx8[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t sumBkgPdf_totbkg_fitW_fy8[47] = {
0,
0,
0,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
2.878569e-06,
3.789757,
3.789757,
71.12063,
71.12063,
211.4394,
211.4394,
855.0132,
855.0132,
1794.513,
1794.513,
3088.88,
3088.88,
4399.772,
4399.772,
5568.272,
5568.272,
6460.972,
6460.972,
6790.725,
6790.725,
6142.168,
6142.168,
5333.208,
5333.208,
4062.953,
4062.953,
3103.771,
3103.771,
2394.858,
2394.858,
1734.507,
1734.507,
1650.883,
1650.883,
0,
0,
0,
0};
graph = new TGraph(47,sumBkgPdf_totbkg_fitW_fx8,sumBkgPdf_totbkg_fitW_fy8);
graph->SetName("sumBkgPdf_totbkg_fitW");
graph->SetTitle("Projection of ");
ci = TColor::GetColor("#ffff66");
graph->SetFillColor(ci);
ci = TColor::GetColor("#ffff66");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#ffff66");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618 = new TH1F("Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618","Projection of ",100,-0.1292683,1.129268);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMinimum(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMaximum(7469.797);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetDirectory(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetStats(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetLineWidth(2);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetRange(11,90);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_Graph_Graph_Graph_sumBkgPdf_totbkg_fitW131618);
graph->Draw("f");
Double_t sumPdf_Norm[muminus_ISO_]_fx9[47] = {
-0.02439024,
-0.02439024,
0,
1.0201e-08,
0.04999999,
0.05000001,
0.09999999,
0.1,
0.15,
0.15,
0.2,
0.2,
0.25,
0.25,
0.3,
0.3,
0.35,
0.35,
0.4,
0.4,
0.45,
0.45,
0.5,
0.5,
0.55,
0.55,
0.6,
0.6,
0.65,
0.65,
0.7,
0.7,
0.75,
0.75,
0.8,
0.8,
0.85,
0.85,
0.9,
0.9,
0.95,
0.95,
1,
1,
1,
1.02439,
1.02439};
Double_t sumPdf_Norm[muminus_ISO_]_fy9[47] = {
0,
0,
0,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
8.630165e-06,
3.789756,
3.789756,
71.12054,
71.12054,
211.4391,
211.4391,
857.6998,
857.6998,
1798.532,
1798.532,
3099.606,
3099.606,
4425.897,
4425.897,
5601.578,
5601.578,
6533.68,
6533.68,
6910.201,
6910.201,
6392.164,
6392.164,
5777.187,
5777.187,
5020.914,
5020.914,
5207.829,
5207.829,
6798.152,
6798.152,
13033.97,
13033.97,
27939.24,
27939.24,
0,
0,
0,
0};
graph = new TGraph(47,sumPdf_Norm[muminus_ISO_]_fx9,sumPdf_Norm[muminus_ISO_]_fy9);
graph->SetName("sumPdf_Norm[muminus_ISO_]");
graph->SetTitle("Projection of ");
ci = TColor::GetColor("#ff0000");
graph->SetFillColor(ci);
graph->SetFillStyle(0);
ci = TColor::GetColor("#ff0000");
graph->SetLineColor(ci);
graph->SetLineWidth(2);
ci = TColor::GetColor("#ff0000");
graph->SetMarkerColor(ci);
graph->SetMarkerStyle(0);
TH1F *Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839 = new TH1F("Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839","Projection of ",100,-0.1292683,1.129268);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMinimum(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMaximum(30733.17);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetDirectory(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetStats(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetLineWidth(2);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetRange(11,90);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839->GetZaxis()->SetTitleFont(132);
graph->SetHistogram(Graph_Graph_Graph_Graph_sumPdf_Norm[muminus_ISO_]151839);
graph->Draw("");
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3006[20] = {
0.025,
0.075,
0.125,
0.175,
0.225,
0.275,
0.325,
0.375,
0.425,
0.475,
0.525,
0.575,
0.625,
0.675,
0.725,
0.775,
0.825,
0.875,
0.925,
0.975};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3006[20] = {
0,
0,
9,
26,
118,
315,
881,
1812,
3086,
4575,
5856,
6537,
6824,
6348,
5662,
4769,
4941,
6726,
13179,
28019};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3006[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3006[20] = {
0,
0,
2.943461,
5.066015,
10.37428,
17.25528,
29.18586,
42.07053,
55.05403,
67.14059,
76.02614,
80.35326,
82.10902,
79.17591,
74.74792,
68.55976,
69.79403,
81.51372,
114.3009,
166.8895};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3006[20] = {
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025,
0.025};
Double_t h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3006[20] = {
1.147874,
1.147874,
4.110204,
6.164324,
11.37428,
18.25528,
30.18586,
43.07053,
56.05403,
68.14059,
77.02614,
81.35326,
83.10902,
80.17591,
75.74792,
69.55976,
70.79403,
82.51372,
115.3009,
167.8895};
grae = new TGraphAsymmErrors(20,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fy3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_felx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehx3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fely3006,h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]_fehy3006);
grae->SetName("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]");
grae->SetTitle("Histogram of combData_plot__muminus_ISO_");
grae->SetFillColor(1);
grae->SetFillStyle(0);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006 = new TH1F("Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMinimum(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetDirectory(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetStats(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetLineWidth(2);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph_Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]563001300230053006);
grae->Draw("e1p");
TLegend *leg = new TLegend(0.2,0.53,0.6,0.78,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]","Data, #sqrt{s}=8 TeV","lep");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("h5","Electroweak","f");
ci = TColor::GetColor("#3333ff");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
ci = TColor::GetColor("#3333ff");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("sumBkgPdf_totbkg_fitW","QCD","f");
ci = TColor::GetColor("#ffff66");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
ci = TColor::GetColor("#ffff66");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("sumPdf_Norm[muminus_ISO_]","Fit","l");
ci = TColor::GetColor("#ff0000");
entry->SetLineColor(ci);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.2,0.8,0.4,0.87,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb");
pt->Draw();
TH1F *Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy = new TH1F("Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy","Histogram of combData_plot__muminus_ISO_",100,0,1.1);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMinimum(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMaximum(31005.58);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetDirectory(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetStats(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetLineWidth(2);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->SetMarkerStyle(0);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitle("#it{p}_{T}^{#it{#mu}} / #it{p}_{T}^{#it{#mu}-jet} ");
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetRange(1,91);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetNdivisions(1005);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelOffset(0.02);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetXaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitle("Events / (0.05)");
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleOffset(1.1);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetYaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetLabelFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetLabelSize(0.05);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetTitleSize(0.06);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->GetZaxis()->SetTitleFont(132);
Graph_Graph_Graph_Graph_h_combData_Cut[samples == samples::fitW_etam_0 || samples == samples::fitW_etam_1 || samples == samples::fitW_etam_2 || samples == samples::fitW_etam_3]5630013002_copy->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
//___________________________________________________________________________
Double_t* IfitBin(TH1F* dataInput, TH1F* sigTemplate, TH1F* bkgTemplate,
int fit_data=1)
{
cout << "Input files are " << dataInput->GetName() << "\t" << sigTemplate->GetName() << "\t" << bkgTemplate->GetName() << endl;
TCanvas *c1 = new TCanvas("HF1", "Histos1", 0, 0, 600, 600);
dataCollBin.clear();
sigCollBin.clear();
bkgCollBin.clear();
Double_t* fitted = new Double_t[8];
fitted[0] = fitted[1] = fitted[2] = fitted[3] = fitted[4] = fitted[5] = fitted[6] = fitted[7] = 0.0;
TH1F *hsum = new TH1F();
float ntemplate = 1.;
float sigfrac = 0.1;
TH1F *hsum_norm = new TH1F();
TH1F *hdata;
TH1F *hsig = (TH1F*)sigTemplate->Clone();
hsig->SetName("hsig");
hsig->Rebin(6);
TH1F *hbkg = (TH1F*)bkgTemplate->Clone();
hbkg->SetName("hbkg");
hbkg->Rebin(6);
float ndata=0;
if ( fit_data>0 ) {
hdata = (TH1F*)dataInput->Clone();
hdata -> SetName("hdata");
ndata = hdata->Integral();
}else {
hsum = (TH1F*)hsig->Clone();
hsum->Add(hbkg,1);
cout << "For histogram " << sigTemplate->GetName() << " and " << bkgTemplate->GetName() << " sum = " <<
hsum->Integral() << endl;
if (hsum->Integral()>1.) ntemplate = hsum->Integral();
sigfrac = hsig->Integral()/ntemplate;
hsum_norm = (TH1F*)hsum->Clone();
hsum_norm->Scale(1./hsum->Integral());
hdata = (TH1F*)hsum_norm->Clone();
hdata -> SetName("hdata");
ndata=ntemplate;
hdata->FillRandom(hsum_norm, ndata);
}
if(ndata==0) {
printf(" --- no events in the fit \n");
return fitted;
}
printf(" --------- before the fit ------------- \n");
printf("Nsig %2.3f, Nbg %2.3f, Ntemplate %3.3f \n", hsig->Integral(), hbkg->Integral(), ntemplate);
// printf("Purity %2.3f, init size %4.3f, test sample size %4d\n", hsig->Integral()/hsum->Integral(), hsum->Integral(), ndata);
printf(" -------------------------------------- \n");
hdata->Rebin(6);
int nbins = hdata->GetNbinsX();
hsig->Scale(1./hsig->Integral());
hbkg->Scale(1./hbkg->Integral());
for (int ibin=1; ibin<=nbins; ibin++) {
dataCollBin.push_back(hdata->GetBinContent(ibin));
sigCollBin.push_back(hsig->GetBinContent(ibin));
bkgCollBin.push_back(hbkg->GetBinContent(ibin));
}
printf( " ----- Got %d, %d, %d events for fit ----- \n ", dataCollBin.size(),
sigCollBin.size(), bkgCollBin.size() );
if ( dataCollBin.size() != sigCollBin.size() || sigCollBin.size()!=bkgCollBin.size() ) {
printf(" error ... inconsistent hit collection size \n");
return fitted;
}
//--------------------------------------------------
//init parameters for fit
Double_t vstart[10] = {1., 1.};
vstart[0] = sigfrac*ndata;
vstart[1] = (1-sigfrac)*ndata;
TMinuit *gMinuit = new TMinuit(NPARBIN);
gMinuit->Command("SET STR 1");
gMinuit->SetFCN(fcnBin);
Double_t arglist[10];
Int_t ierflg = 0;
arglist[0] = 1;
gMinuit->mnexcm("SET ERR", arglist ,1,ierflg);
arglist[0] = 1;
gMinuit->mnexcm("SET PRINT", arglist ,1,ierflg);
Double_t step[] = { 0.1, 0.1,};
gMinuit->mnparm(0, "Signal yield" , vstart[0], step[0], 0., ndata*2. , ierflg);
gMinuit->mnparm(1, "background yield" , vstart[1], step[1], 0., ndata*2. , ierflg);
printf(" --------------------------------------------------------- \n");
printf(" Now ready for minimization step \n --------------------------------------------------------- \n");
arglist[0] = 2000; // number of iteration
arglist[1] = 1.;
gMinuit->mnexcm("MIGRAD", arglist ,2,ierflg);
printf (" -------------------------------------------- \n");
printf("Finished. ierr = %d \n", ierflg);
infoBin.clear();
infoBin_err.clear();
double para[NPARBIN+1],errpara[NPARBIN+1];
if ( ierflg == 0 )
{
for(int j=0; j<=NPARBIN-1;j++) {
gMinuit->GetParameter(j, para[j],errpara[j]);
para[NPARBIN] = dataCollBin.size();
infoBin.push_back(para[j]);
infoBin_err.push_back(errpara[j]);
printf("Parameter (yeild) %d = %f +- %f\n",j,para[j],errpara[j]);
}
printf(" fitted yield %2.3f \n", (para[0]+para[1])/ndata );
infoBin.push_back(sigCollBin.size());
}
else {
printf(" *********** Fit failed! ************\n");
gMinuit->GetParameter(0, para[0],errpara[0]);
gMinuit->GetParameter(1, para[1],errpara[1]);
para[0]=0.; errpara[0]=0.;
}
// Print results
Double_t amin,edm,errdef;
Int_t nvpar,nparx,icstat;
gMinuit->mnstat(amin,edm,errdef,nvpar,nparx,icstat);
gMinuit->mnprin(1,amin);
gMinuit->mnmatu(1);
printf(" ========= happy ending !? =========================== \n");
printf("FCN = %3.3f \n", amin);
double yerr[20];
for(int i=0;i<20;i++){
yerr[i] = 0.;
}
hsig->Scale(para[0]);
hbkg->Scale(para[1]);
TH1F *hfit = (TH1F*)hsig->Clone();
hfit->Add(hbkg);
hsig->SetLineColor(1);
hsig->SetFillColor(5);
hsig->SetFillStyle(3001);
hbkg->SetLineWidth(2);
// plot
c1->Draw();
//gPad->SetLogy();
hdata->SetLineColor(1);
hdata->SetNdivisions(505,"XY");
hdata->SetXTitle("Iso_{ECAL}+Iso_{HCAL}+Iso_{TRK} (GeV)");
hdata->SetYTitle("Entries");
hdata->SetTitle("");
hdata->SetMarkerStyle(8);
hdata->SetMinimum(0.);
hdata->SetMaximum(hdata->GetMaximum()*1.5);
hdata->Draw("p e");
hsig->Draw("hist same");
hbkg->SetMarkerStyle(0);
hbkg->SetFillColor(8);
hbkg->SetLineWidth(1);
hbkg->SetFillStyle(3013);
hbkg->SetError(yerr);
hbkg->Draw("hist same");
hfit->SetMarkerStyle(0);
hfit->SetLineColor(1);
hfit->SetLineWidth(2);
hfit->SetError(yerr);
hfit->Draw("hist same");
double chi2ForThisBin=0;
int nbinForThisBin=0;
chi2NbinsHisto(hfit, hdata, chi2ForThisBin, nbinForThisBin);
TPaveText *pavetex = new TPaveText(0.43, 0.87, 0.90, 0.92,"NDCBR");
pavetex->SetBorderSize(0);
pavetex->SetFillColor(0);
pavetex->SetFillStyle(0);
pavetex->SetLineWidth(3);
pavetex->SetTextAlign(12);
pavetex->SetTextSize(0.03);
pavetex->AddText(Form("#chi^{2}/NDF=%.1f/%d",chi2ForThisBin, nbinForThisBin));
pavetex->Draw();
char text[1000];
TLegend *tleg = new TLegend(0.43, 0.60, 0.90, 0.87);
tleg->SetHeader(dataInput->GetTitle());
tleg->SetTextSize(0.03);
tleg->SetFillColor(0);
tleg->SetShadowColor(0);
tleg->SetBorderSize(0);
sprintf(text,"Data %5.1f events",hdata->Integral());
tleg->AddEntry(hdata,text,"pl");
sprintf(text,"Fitted %5.1f events",hfit->Integral());
tleg->AddEntry(hfit,text,"l");
sprintf(text,"SIG %5.1f #pm %5.1f events",para[0], errpara[0]);
tleg->AddEntry(hsig,text,"f");
sprintf(text,"BKG %5.1f #pm %5.1f events",para[1], errpara[1]);
tleg->AddEntry(hbkg,text,"f");
tleg->Draw();
gPad->RedrawAxis();
cout << dataInput->GetName() << endl;
char fname[300];
sprintf(fname,"plots/Ifit_%s.eps",dataInput->GetName());
c1->SaveAs(fname);
sprintf(fname,"plots/Ifit_%s.gif",dataInput->GetName());
c1->SaveAs(fname);
printf("----- fit results with signal projection ----------- \n");
// ftemplate->Close();
int purityMaxBin = hsig->FindBin(5.0)-1;
Double_t scale_signal = hsig->Integral(1,purityMaxBin)/hsig->Integral();
Double_t scale_background = hbkg->Integral(1,purityMaxBin)/hbkg->Integral();
fitted[0] = para[0];
fitted[1] = errpara[0];
fitted[2] = para[1];
fitted[3] = errpara[1];
// for integral up to 5 GeV
fitted[4] = para[0]*scale_signal;
fitted[5] = errpara[0]*scale_signal;
fitted[6] = para[1]*scale_background;
fitted[7] = errpara[1]*scale_background;
return fitted;
}
void zj_dphi_8TeV_log()
{
//=========Macro generated from canvas: default_Canvas/defaultCanvas
//========= (Thu Apr 28 10:19:34 2016) by ROOT version6.04/10
TCanvas *default_Canvas = new TCanvas("default_Canvas", "defaultCanvas",0,0,800,800);
gStyle->SetOptStat(0);
gStyle->SetOptTitle(0);
default_Canvas->SetHighLightColor(2);
default_Canvas->Range(0,0,1,1);
default_Canvas->SetFillColor(0);
default_Canvas->SetBorderMode(0);
default_Canvas->SetBorderSize(2);
default_Canvas->SetTickx(1);
default_Canvas->SetTicky(1);
default_Canvas->SetLeftMargin(0.14);
default_Canvas->SetRightMargin(0.05);
default_Canvas->SetTopMargin(0.05);
default_Canvas->SetBottomMargin(0.16);
default_Canvas->SetFrameLineWidth(2);
default_Canvas->SetFrameBorderMode(0);
// ------------>Primitives in pad: bottomPad
TPad *bottomPad = new TPad("bottomPad", "bottomPad",0.005,0.05,0.995,0.995);
bottomPad->Draw();
bottomPad->cd();
bottomPad->Range(0,0,1,1);
bottomPad->SetFillColor(0);
bottomPad->SetFillStyle(4000);
bottomPad->SetBorderMode(0);
bottomPad->SetBorderSize(2);
bottomPad->SetTickx(1);
bottomPad->SetTicky(1);
bottomPad->SetLeftMargin(0.14);
bottomPad->SetRightMargin(0.05);
bottomPad->SetTopMargin(0.05);
bottomPad->SetBottomMargin(0.16);
bottomPad->SetFrameLineWidth(2);
bottomPad->SetFrameBorderMode(0);
bottomPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: upperPad
TPad *upperPad = new TPad("upperPad", "upperPad",0.005,0.05,0.995,0.995);
upperPad->Draw();
upperPad->cd();
upperPad->Range(-0.5429913,-2.062865,3.335518,1.327192);
upperPad->SetFillColor(0);
upperPad->SetFillStyle(4000);
upperPad->SetBorderMode(0);
upperPad->SetBorderSize(2);
upperPad->SetLogy();
upperPad->SetTickx(1);
upperPad->SetTicky(1);
upperPad->SetLeftMargin(0.14);
upperPad->SetRightMargin(0.05);
upperPad->SetTopMargin(0.05);
upperPad->SetBottomMargin(0.33);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
upperPad->SetFrameLineWidth(2);
upperPad->SetFrameBorderMode(0);
TH1F *z_dphi_toterr__37 = new TH1F("z_dphi_toterr__37","z_dphi",6,0,3.141593);
z_dphi_toterr__37->SetBinContent(1,0.261202);
z_dphi_toterr__37->SetBinContent(2,0.3162337);
z_dphi_toterr__37->SetBinContent(3,0.3852335);
z_dphi_toterr__37->SetBinContent(4,0.7150939);
z_dphi_toterr__37->SetBinContent(5,1.796506);
z_dphi_toterr__37->SetBinContent(6,7.258657);
z_dphi_toterr__37->SetBinContent(7,3.800624);
z_dphi_toterr__37->SetBinError(1,0.03375302);
z_dphi_toterr__37->SetBinError(2,0.03634611);
z_dphi_toterr__37->SetBinError(3,0.0427511);
z_dphi_toterr__37->SetBinError(4,0.05747346);
z_dphi_toterr__37->SetBinError(5,0.1278198);
z_dphi_toterr__37->SetBinError(6,0.3295685);
z_dphi_toterr__37->SetBinError(7,0.1725617);
z_dphi_toterr__37->SetEntries(14);
z_dphi_toterr__37->SetStats(0);
Int_t ci; // for color index setting
TColor *color; // for color definition with alpha
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetLineColor(ci);
z_dphi_toterr__37->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr__37->SetMarkerColor(ci);
z_dphi_toterr__37->SetMarkerStyle(0);
z_dphi_toterr__37->GetXaxis()->SetNdivisions(4);
z_dphi_toterr__37->GetXaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetXaxis()->SetLabelOffset(999);
z_dphi_toterr__37->GetXaxis()->SetLabelSize(0);
z_dphi_toterr__37->GetXaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetXaxis()->SetTitleFont(132);
z_dphi_toterr__37->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_toterr__37->GetYaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetYaxis()->SetLabelSize(0.05);
z_dphi_toterr__37->GetYaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetYaxis()->SetTitleOffset(1.3);
z_dphi_toterr__37->GetYaxis()->SetTitleFont(132);
z_dphi_toterr__37->GetZaxis()->SetLabelFont(132);
z_dphi_toterr__37->GetZaxis()->SetLabelSize(0.05);
z_dphi_toterr__37->GetZaxis()->SetTitleSize(0.06);
z_dphi_toterr__37->GetZaxis()->SetTitleFont(132);
z_dphi_toterr__37->Draw("e2");
TH1F *z_dphi_staterr__38 = new TH1F("z_dphi_staterr__38","z_dphi",6,0,3.141593);
z_dphi_staterr__38->SetBinContent(1,0.261202);
z_dphi_staterr__38->SetBinContent(2,0.3162337);
z_dphi_staterr__38->SetBinContent(3,0.3852335);
z_dphi_staterr__38->SetBinContent(4,0.7150939);
z_dphi_staterr__38->SetBinContent(5,1.796506);
z_dphi_staterr__38->SetBinContent(6,7.258657);
z_dphi_staterr__38->SetBinContent(7,3.800624);
z_dphi_staterr__38->SetBinError(1,0.01764622);
z_dphi_staterr__38->SetBinError(2,0.01951211);
z_dphi_staterr__38->SetBinError(3,0.02108379);
z_dphi_staterr__38->SetBinError(4,0.03021233);
z_dphi_staterr__38->SetBinError(5,0.04899558);
z_dphi_staterr__38->SetBinError(6,0.1002775);
z_dphi_staterr__38->SetBinError(7,0.05250516);
z_dphi_staterr__38->SetEntries(14);
z_dphi_staterr__38->SetStats(0);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetLineColor(ci);
z_dphi_staterr__38->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
z_dphi_staterr__38->SetMarkerColor(ci);
z_dphi_staterr__38->SetMarkerStyle(0);
z_dphi_staterr__38->GetXaxis()->SetNdivisions(4);
z_dphi_staterr__38->GetXaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetXaxis()->SetLabelOffset(999);
z_dphi_staterr__38->GetXaxis()->SetLabelSize(0);
z_dphi_staterr__38->GetXaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetXaxis()->SetTitleFont(132);
z_dphi_staterr__38->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_staterr__38->GetYaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetYaxis()->SetLabelSize(0.05);
z_dphi_staterr__38->GetYaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetYaxis()->SetTitleFont(132);
z_dphi_staterr__38->GetZaxis()->SetLabelFont(132);
z_dphi_staterr__38->GetZaxis()->SetLabelSize(0.05);
z_dphi_staterr__38->GetZaxis()->SetTitleSize(0.06);
z_dphi_staterr__38->GetZaxis()->SetTitleFont(132);
z_dphi_staterr__38->Draw("e2same");
Double_t dphi_toterrs_fx3043[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t dphi_toterrs_fy3043[6] = {
0.1926519,
0.2207154,
0.3546338,
0.6864677,
2.01515,
7.779723};
Double_t dphi_toterrs_felx3043[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fely3043[6] = {
0.03415733,
0.03577312,
0.05682396,
0.09598937,
0.1897309,
0.5023379};
Double_t dphi_toterrs_fehx3043[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fehy3043[6] = {
0.04673612,
0.05211068,
0.07287587,
0.1207667,
0.143423,
0.3949306};
TGraphAsymmErrors *grae = new TGraphAsymmErrors(6,dphi_toterrs_fx3043,dphi_toterrs_fy3043,dphi_toterrs_felx3043,dphi_toterrs_fehx3043,dphi_toterrs_fely3043,dphi_toterrs_fehy3043);
grae->SetName("dphi_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_dphi_toterrs3043 = new TH1F("Graph_dphi_toterrs3043","",100,0.07853982,3.220132);
Graph_dphi_toterrs3043->SetMinimum(0.1426451);
Graph_dphi_toterrs3043->SetMaximum(8.976269);
Graph_dphi_toterrs3043->SetDirectory(0);
Graph_dphi_toterrs3043->SetStats(0);
Graph_dphi_toterrs3043->SetLineWidth(2);
Graph_dphi_toterrs3043->SetMarkerStyle(0);
Graph_dphi_toterrs3043->GetXaxis()->SetNdivisions(4);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelOffset(999);
Graph_dphi_toterrs3043->GetXaxis()->SetLabelSize(0);
Graph_dphi_toterrs3043->GetXaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetXaxis()->SetTitleFont(132);
Graph_dphi_toterrs3043->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_dphi_toterrs3043->GetYaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetYaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3043->GetYaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetYaxis()->SetTitleFont(132);
Graph_dphi_toterrs3043->GetZaxis()->SetLabelFont(132);
Graph_dphi_toterrs3043->GetZaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3043->GetZaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3043->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_dphi_toterrs3043);
grae->Draw("p1p");
Double_t _fx3044[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t _fy3044[6] = {
0.1926519,
0.2207154,
0.3546338,
0.6864677,
2.01515,
7.779723};
Double_t _felx3044[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fely3044[6] = {
0.02873042,
0.03223617,
0.04948444,
0.07740101,
0.1117696,
0.3251484};
Double_t _fehx3044[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fehy3044[6] = {
0.03948655,
0.04420709,
0.06698785,
0.1028379,
0.1104173,
0.2636435};
grae = new TGraphAsymmErrors(6,_fx3044,_fy3044,_felx3044,_fehx3044,_fely3044,_fehy3044);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_Graph3044 = new TH1F("Graph_Graph3044","",100,0.07853982,3.220132);
Graph_Graph3044->SetMinimum(0.1475293);
Graph_Graph3044->SetMaximum(8.831311);
Graph_Graph3044->SetDirectory(0);
Graph_Graph3044->SetStats(0);
Graph_Graph3044->SetLineWidth(2);
Graph_Graph3044->SetMarkerStyle(0);
Graph_Graph3044->GetXaxis()->SetNdivisions(4);
Graph_Graph3044->GetXaxis()->SetLabelFont(132);
Graph_Graph3044->GetXaxis()->SetLabelOffset(999);
Graph_Graph3044->GetXaxis()->SetLabelSize(0);
Graph_Graph3044->GetXaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetXaxis()->SetTitleFont(132);
Graph_Graph3044->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_Graph3044->GetYaxis()->SetLabelFont(132);
Graph_Graph3044->GetYaxis()->SetLabelSize(0.05);
Graph_Graph3044->GetYaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetYaxis()->SetTitleFont(132);
Graph_Graph3044->GetZaxis()->SetLabelFont(132);
Graph_Graph3044->GetZaxis()->SetLabelSize(0.05);
Graph_Graph3044->GetZaxis()->SetTitleSize(0.06);
Graph_Graph3044->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph3044);
grae->Draw("p1p");
Double_t dphi_toterrs_fx3045[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t dphi_toterrs_fy3045[6] = {
0.3071206,
0.3788296,
0.5235754,
0.9026085,
2.09773,
8.042891};
Double_t dphi_toterrs_felx3045[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fely3045[6] = {
0.03223837,
0.05049376,
0.0513308,
0.07859791,
0.1671821,
0.6048482};
Double_t dphi_toterrs_fehx3045[6] = {
0,
0,
0,
0,
0,
0};
Double_t dphi_toterrs_fehy3045[6] = {
0.05069997,
0.04670219,
0.05741781,
0.0972528,
0.2140117,
0.6377428};
grae = new TGraphAsymmErrors(6,dphi_toterrs_fx3045,dphi_toterrs_fy3045,dphi_toterrs_felx3045,dphi_toterrs_fehx3045,dphi_toterrs_fely3045,dphi_toterrs_fehy3045);
grae->SetName("dphi_toterrs");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_dphi_toterrs3045 = new TH1F("Graph_dphi_toterrs3045","",100,0,3.063053);
Graph_dphi_toterrs3045->SetMinimum(0.247394);
Graph_dphi_toterrs3045->SetMaximum(9.521209);
Graph_dphi_toterrs3045->SetDirectory(0);
Graph_dphi_toterrs3045->SetStats(0);
Graph_dphi_toterrs3045->SetLineWidth(2);
Graph_dphi_toterrs3045->SetMarkerStyle(0);
Graph_dphi_toterrs3045->GetXaxis()->SetNdivisions(4);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelOffset(999);
Graph_dphi_toterrs3045->GetXaxis()->SetLabelSize(0);
Graph_dphi_toterrs3045->GetXaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetXaxis()->SetTitleFont(132);
Graph_dphi_toterrs3045->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_dphi_toterrs3045->GetYaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetYaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3045->GetYaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetYaxis()->SetTitleFont(132);
Graph_dphi_toterrs3045->GetZaxis()->SetLabelFont(132);
Graph_dphi_toterrs3045->GetZaxis()->SetLabelSize(0.05);
Graph_dphi_toterrs3045->GetZaxis()->SetTitleSize(0.06);
Graph_dphi_toterrs3045->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_dphi_toterrs3045);
grae->Draw("p1p");
Double_t _fx3046[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t _fy3046[6] = {
0.3071206,
0.3788296,
0.5235754,
0.9026085,
2.09773,
8.042891};
Double_t _felx3046[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fely3046[6] = {
0.02710968,
0.03295216,
0.04124188,
0.06402796,
0.1338949,
0.4694553};
Double_t _fehx3046[6] = {
0,
0,
0,
0,
0,
0};
Double_t _fehy3046[6] = {
0.03469016,
0.04020747,
0.048986,
0.07362466,
0.1488417,
0.5027875};
grae = new TGraphAsymmErrors(6,_fx3046,_fy3046,_felx3046,_fehx3046,_fely3046,_fehy3046);
grae->SetName("");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_Graph3046 = new TH1F("Graph_Graph3046","",100,0,3.063053);
Graph_Graph3046->SetMinimum(0.2520099);
Graph_Graph3046->SetMaximum(9.372245);
Graph_Graph3046->SetDirectory(0);
Graph_Graph3046->SetStats(0);
Graph_Graph3046->SetLineWidth(2);
Graph_Graph3046->SetMarkerStyle(0);
Graph_Graph3046->GetXaxis()->SetNdivisions(4);
Graph_Graph3046->GetXaxis()->SetLabelFont(132);
Graph_Graph3046->GetXaxis()->SetLabelOffset(999);
Graph_Graph3046->GetXaxis()->SetLabelSize(0);
Graph_Graph3046->GetXaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetXaxis()->SetTitleFont(132);
Graph_Graph3046->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
Graph_Graph3046->GetYaxis()->SetLabelFont(132);
Graph_Graph3046->GetYaxis()->SetLabelSize(0.05);
Graph_Graph3046->GetYaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetYaxis()->SetTitleFont(132);
Graph_Graph3046->GetZaxis()->SetLabelFont(132);
Graph_Graph3046->GetZaxis()->SetLabelSize(0.05);
Graph_Graph3046->GetZaxis()->SetTitleSize(0.06);
Graph_Graph3046->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_Graph3046);
grae->Draw("p1p");
TLegend *leg = new TLegend(0.2,0.57,0.5,0.82,NULL,"brNDC");
leg->SetBorderSize(0);
leg->SetTextFont(132);
leg->SetLineColor(1);
leg->SetLineStyle(1);
leg->SetLineWidth(2);
leg->SetFillColor(10);
leg->SetFillStyle(0);
TLegendEntry *entry=leg->AddEntry("data_stat","Data_{stat}","f");
ci = TColor::GetColor("#ff9933");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("data_tot","Data_{tot}","f");
ci = TColor::GetColor("#ffff00");
entry->SetFillColor(ci);
entry->SetFillStyle(1001);
entry->SetLineStyle(1);
entry->SetLineWidth(2);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(21);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("powheg_temp","POWHEG","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(20);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
entry=leg->AddEntry("amcatnlo_temp","aMC@NLO","p");
entry->SetLineColor(1);
entry->SetLineStyle(1);
entry->SetLineWidth(1);
entry->SetMarkerColor(1);
entry->SetMarkerStyle(24);
entry->SetMarkerSize(1);
entry->SetTextFont(132);
leg->Draw();
TPaveText *pt = new TPaveText(0.2,0.8,0.55,0.95,"brNDC");
pt->SetBorderSize(0);
pt->SetFillStyle(0);
pt->SetLineWidth(2);
pt->SetTextFont(132);
TText *AText = pt->AddText("LHCb, #sqrt{s} = 8 TeV");
pt->Draw();
TH1F *z_dphi_toterr_copy__39 = new TH1F("z_dphi_toterr_copy__39","z_dphi",6,0,3.141593);
z_dphi_toterr_copy__39->SetBinContent(1,0.261202);
z_dphi_toterr_copy__39->SetBinContent(2,0.3162337);
z_dphi_toterr_copy__39->SetBinContent(3,0.3852335);
z_dphi_toterr_copy__39->SetBinContent(4,0.7150939);
z_dphi_toterr_copy__39->SetBinContent(5,1.796506);
z_dphi_toterr_copy__39->SetBinContent(6,7.258657);
z_dphi_toterr_copy__39->SetBinContent(7,3.800624);
z_dphi_toterr_copy__39->SetBinError(1,0.03375302);
z_dphi_toterr_copy__39->SetBinError(2,0.03634611);
z_dphi_toterr_copy__39->SetBinError(3,0.0427511);
z_dphi_toterr_copy__39->SetBinError(4,0.05747346);
z_dphi_toterr_copy__39->SetBinError(5,0.1278198);
z_dphi_toterr_copy__39->SetBinError(6,0.3295685);
z_dphi_toterr_copy__39->SetBinError(7,0.1725617);
z_dphi_toterr_copy__39->SetEntries(14);
z_dphi_toterr_copy__39->SetDirectory(0);
z_dphi_toterr_copy__39->SetStats(0);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetLineColor(ci);
z_dphi_toterr_copy__39->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
z_dphi_toterr_copy__39->SetMarkerColor(ci);
z_dphi_toterr_copy__39->SetMarkerStyle(0);
z_dphi_toterr_copy__39->GetXaxis()->SetNdivisions(4);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelOffset(999);
z_dphi_toterr_copy__39->GetXaxis()->SetLabelSize(0);
z_dphi_toterr_copy__39->GetXaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetXaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->GetYaxis()->SetTitle("d#sigma(#it{Zj})/d|#Delta#it{#phi} [pb/rad]");
z_dphi_toterr_copy__39->GetYaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetYaxis()->SetLabelSize(0.05);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleOffset(1.3);
z_dphi_toterr_copy__39->GetYaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->GetZaxis()->SetLabelFont(132);
z_dphi_toterr_copy__39->GetZaxis()->SetLabelSize(0.05);
z_dphi_toterr_copy__39->GetZaxis()->SetTitleSize(0.06);
z_dphi_toterr_copy__39->GetZaxis()->SetTitleFont(132);
z_dphi_toterr_copy__39->Draw("sameaxis");
upperPad->Modified();
default_Canvas->cd();
// ------------>Primitives in pad: lowerPad
TPad *lowerPad = new TPad("lowerPad", "lowerPad",0.005,0.05,0.995,0.995);
lowerPad->Draw();
lowerPad->cd();
lowerPad->Range(-0.5429913,-0.4371795,3.335518,5.614103);
lowerPad->SetFillColor(0);
lowerPad->SetFillStyle(4000);
lowerPad->SetBorderMode(0);
lowerPad->SetBorderSize(2);
lowerPad->SetTickx(1);
lowerPad->SetTicky(1);
lowerPad->SetLeftMargin(0.14);
lowerPad->SetRightMargin(0.05);
lowerPad->SetTopMargin(0.665);
lowerPad->SetBottomMargin(0.14);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
lowerPad->SetFrameLineWidth(2);
lowerPad->SetFrameBorderMode(0);
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp0__40 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp0__40","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(1,0.1292219);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(2,0.1149343);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(3,0.1109745);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(4,0.08037191);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(5,0.07114914);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(6,0.04540351);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetBinError(7,0.1725617);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleOffset(1.3);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0__40->Draw("e2");
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp1__41 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp1__41","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(1,0.06755777);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(2,0.06170157);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(3,0.0547299);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(4,0.04224946);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(5,0.02727271);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(6,0.01381488);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetBinError(7,0.05250516);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetStats(0);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetFillColor(ci);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetLineWidth(2);
ci = TColor::GetColor("#ff9933");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp1__41->Draw("e2same");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fx3047[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fy3047[6] = {
0.7375589,
0.6979502,
0.9205685,
0.9599686,
1.121705,
1.071786};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_felx3047[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fely3047[6] = {
0.1307698,
0.1131224,
0.1475052,
0.1342332,
0.1056111,
0.06920535};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehx3047[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehy3047[6] = {
0.1789271,
0.1647854,
0.1891732,
0.1688822,
0.07983441,
0.05440822};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fy3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_felx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehx3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fely3047,zj_dphi_8TeV_log_400_1001_0ratiocomp2_fehy3047);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp2");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047","",100,0.07853982,3.220132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp23047);
grae->Draw("p1");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fx3048[6] = {
0.3403392,
0.863938,
1.387537,
1.911136,
2.434734,
2.958333};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fy3048[6] = {
0.7375589,
0.6979502,
0.9205685,
0.9599686,
1.121705,
1.071786};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_felx3048[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fely3048[6] = {
0.1099931,
0.1019378,
0.1284531,
0.1082389,
0.06221498,
0.04479457};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehx3048[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehy3048[6] = {
0.1511724,
0.1397925,
0.1738889,
0.1438103,
0.06146226,
0.03632125};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fy3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_felx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehx3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fely3048,zj_dphi_8TeV_log_400_1001_0ratiocomp3_fehy3048);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp3");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(24);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048","",100,0.07853982,3.220132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp33048);
grae->Draw("p1");
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fx3049[6] = {
0.1832596,
0.7068583,
1.230457,
1.754056,
2.277655,
2.801253};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fy3049[6] = {
1.175797,
1.197942,
1.359112,
1.262224,
1.167673,
1.108041};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_felx3049[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fely3049[6] = {
0.1234231,
0.1596723,
0.1332459,
0.1099127,
0.09305962,
0.08332785};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehx3049[6] = {
0,
0,
0,
0,
0,
0};
Double_t zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehy3049[6] = {
0.1941025,
0.1476825,
0.1490468,
0.136,
0.1191266,
0.08785962};
grae = new TGraphAsymmErrors(6,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fy3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_felx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehx3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fely3049,zj_dphi_8TeV_log_400_1001_0ratiocomp4_fehy3049);
grae->SetName("zj_dphi_8TeV_log_400_1001_0ratiocomp4");
grae->SetTitle("");
grae->SetFillColor(1);
grae->SetLineWidth(2);
grae->SetMarkerStyle(20);
TH1F *Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049 = new TH1F("Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049","",100,0,3.063053);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMinimum(0.41);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMaximum(1.59);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetDirectory(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetStats(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetLineWidth(2);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->SetMarkerStyle(0);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetNdivisions(4);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelOffset(0.02);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetXaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitle("Ratio");
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->CenterTitle(true);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetNdivisions(505);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetYaxis()->SetTitleFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetLabelFont(132);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetLabelSize(0.05);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetTitleSize(0.06);
Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049->GetZaxis()->SetTitleFont(132);
grae->SetHistogram(Graph_zj_dphi_8TeV_log_400_1001_0ratiocomp43049);
grae->Draw("p1");
TH1F *zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42 = new TH1F("zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42","z_dphi",6,0,3.141593);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(1,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(2,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(3,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(4,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(5,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(6,1);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinContent(7,3.800624);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(1,0.1292219);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(2,0.1149343);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(3,0.1109745);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(4,0.08037191);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(5,0.07114914);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(6,0.04540351);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetBinError(7,0.1725617);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMinimum(0.41);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMaximum(1.59);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetEntries(20);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetDirectory(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetStats(0);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetFillColor(ci);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetLineColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetLineWidth(2);
ci = TColor::GetColor("#ffff00");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMarkerColor(ci);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->SetMarkerStyle(0);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitle("|#Delta#it{#phi}| [rad]");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetNdivisions(4);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelOffset(0.02);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetXaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitle("Ratio");
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->CenterTitle(true);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetNdivisions(505);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTickLength(0.12325);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleOffset(1.3);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetYaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetLabelFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetLabelSize(0.05);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetTitleSize(0.06);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->GetZaxis()->SetTitleFont(132);
zj_dphi_8TeV_log_400_1001_0ratiocomp0_copy__42->Draw("sameaxis");
lowerPad->Modified();
default_Canvas->cd();
default_Canvas->Modified();
default_Canvas->cd();
default_Canvas->SetSelected(default_Canvas);
}
