void compare_samples(TString sample="Singlemu", TString mode="rate") {
TString ifile1, ifile2, data;
if (sample=="Minbias") {
ifile1="l1Ntuple_ZeroBias_RunD_oldPMTs_HF-.root";
ifile2="l1Ntuple_ZeroBias_RunD_upgradedPMTs_HF-.root";
data="MinBias Run2012D";
} else if (sample=="Singlemu") {
// ifile1="l1Ntuple_IsoMu_RunC_1509.root";
ifile1="l1Ntuple_SingleMu_RunD_oldPMTs.root";
ifile2="l1Ntuple_SingleMu_RunD_upgradedPMTs.root";
data="singleMu 2012D";
}
TString dirname="l1ExtraJets/";
TH1D *ptJet_e, *ptJet_c, *ptJet_l;
TH1D *ptJetHF_e, *ptJetHF_c, *ptJetHF_l;
// Compute cummulative histograms
if (mode=="rate") {
// // Central + Tau jets
ptJet_e=hTemp(ifile1,dirname,"ptJet_early");
ptJet_c=hTemp(ifile1,dirname,"ptJet_central");
ptJet_l=hTemp(ifile1,dirname,"ptJet_late");
// HFjets
ptJetHF_e=hTemp(ifile1,dirname,"ptJetHF_early");
ptJetHF_c=hTemp(ifile1,dirname,"ptJetHF_central");
ptJetHF_l=hTemp(ifile1,dirname,"ptJetHF_late");
ptJetHF2_e=hTemp(ifile2,dirname,"ptJetHF_early");
ptJetHF2_c=hTemp(ifile2,dirname,"ptJetHF_central");
ptJetHF2_l=hTemp(ifile2,dirname,"ptJetHF_late");
} else if (mode=="norate") {
ptJet_e=readHist(dirname+"ptJet_early",ifile1,4);
ptJet_c=readHist(dirname+"ptJet_central",ifile1,4);
ptJet_l=readHist(dirname+"ptJet_late",ifile1,4);
ptJetHF_e=readHist(dirname+"ptJetHF_early",ifile1,4);
ptJetHF_c=readHist(dirname+"ptJetHF_central",ifile1,4);
ptJetHF_l=readHist(dirname+"ptJetHF_late",ifile1,4);
ptJetHF2_e=readHist(dirname+"ptJetHF_early",ifile2,4);
ptJetHF2_c=readHist(dirname+"ptJetHF_central",ifile2,4);
ptJetHF2_l=readHist(dirname+"ptJetHF_late",ifile2,4);
}
TH1D *etaJet_e=readHist(dirname+"etaJet_early",ifile1,2);
TH1D *etaJet_c=readHist(dirname+"etaJet_central",ifile1,2);
TH1D *etaJet_l=readHist(dirname+"etaJet_late",ifile1,2);
TH1D *etaJetHF_e=readHist(dirname+"etaJetHF_early",ifile1,2);
TH1D *etaJetHF_c=readHist(dirname+"etaJetHF_central",ifile1,2);
TH1D *etaJetHF_l=readHist(dirname+"etaJetHF_late",ifile1,2);
TH1D *phiJet_e=readHist(dirname+"phiJet_early",ifile1,0);
TH1D *phiJet_c=readHist(dirname+"phiJet_central",ifile1,0);
TH1D *phiJet_l=readHist(dirname+"phiJet_late",ifile1,0);
setTDRStyle();
// Central+Tau Jets
TCanvas *c1=getaCanvas("c1");
c1->Divide(1,2);
TPad *pad1 = new TPad("pad1", "The pad with the function",0.05,0.4,0.95,0.95);
pad1->Draw();
// pad1->Range(-112.6742,-73.17708,1143.438,551.3021);
pad1->SetFillColor(0);
pad1->SetBorderMode(0);
pad1->SetBorderSize(2);
// pad1->SetGridx();
// pad1->SetGridy();
pad1->SetLeftMargin(0.1271439);
pad1->SetRightMargin(0.07307979);
pad1->SetTopMargin(0.08215179);
pad1->SetBottomMargin(0.117181);
pad1->SetFrameBorderMode(0);
pad1->SetFrameBorderMode(0);
TPad *pad2 = new TPad("pad2", "The pad with the histogram",0.05,0.05,0.95,0.4);
pad2->Draw();
// pad2->Range(-153.3652,-2.142584,1185.427,5.367464);
pad2->SetFillColor(0);
pad2->SetBorderMode(0);
pad2->SetBorderSize(2);
pad2->SetGridx();
pad2->SetGridy();
pad2->SetLeftMargin(0.1215511);
pad2->SetRightMargin(0.07867263);
pad2->SetTopMargin(0.04892967);
pad2->SetBottomMargin(0.3521407);
pad2->SetFrameBorderMode(0);
pad2->SetFrameBorderMode(0);
pad1->cd();
gPad->SetLogy();
// gPad->SetGridx(); gPad->SetGridy();
TLegend *leg = legend();
leg->SetHeader(data);
// leg->AddEntry(ptJet_e,"BX=-1","LF");
leg->AddEntry(ptJet_c,"BX=0","LF");
leg->AddEntry(ptJet_e,"BX=-1","LF");
// leg->AddEntry(ptJet_l,"BX=+1","LF");
ptJet_c->Draw("EHIST");
ptJet_c->SetLineColor(1); ptJet_c->SetLineWidth(2);
ptJet_c->GetXaxis()->SetRangeUser(0.,260.);
ptJet_c->GetYaxis()->SetRangeUser(1.01,100000000.);
if (mode=="rate") {
ptJet_c->SetTitle(";L1Jet E_{T} Threshold (GeV);a.u.");
} else {
ptJet_c->SetTitle(";L1Jet E_{T} (GeV);a.u.");
}
ptJet_e->Draw("EHISTSAMES");
ptJet_e->SetLineColor(1); ptJet_e->SetLineColor(2);
ptJet_e->SetLineStyle(1);ptJet_e->SetLineWidth(2);
// ptJet_l->Draw("EHISTSAMES");
// ptJet_l->SetLineColor(1); ptJet_l->SetLineColor(4);
// ptJet_l->SetLineStyle(3);
leg->Draw("SAME");
// TPad *p2=getaPad_dn("p2"); p2->Draw();
// p2->cd();
//c1->cd(2);
pad2->cd();
gPad->SetGridx(); gPad->SetGridy();
TH1D *r1_e = ptJet_e->Clone();
r1_e->Divide(ptJet_c);
TH1D *r1_l = ptJet_l->Clone();
r1_l->Divide(ptJet_c);
TLegend *leg = legend();
leg->SetHeader(data);
leg->AddEntry(r1_e,"BX=-1","LF");
//leg->AddEntry(r1_l,"BX=+1","LF");
r1_e->Draw("EHIST");r1_e->SetLineColor(4);
// r1_l->Draw("EHISTSAME"); r1_l->SetLineColor(4);
r1_e->GetYaxis()->SetRangeUser(-0.01,1.0);
r1_e->GetXaxis()->SetRangeUser(0.,260.);
r1_e->SetLineWidth(2); r1_l->SetLineWidth(2);
if (mode=="rate") {
r1_e->SetTitle(";L1Jet E_{T} Threshold (GeV);Ratio");
} else {
r1_e->SetTitle(";L1Jet E_{T} (GeV);Ratio");
}
// r1_e->SetTitle(";L1Jet E_{T} (GeV);ratio");
r1_e->GetXaxis()->SetLabelSize(.1);
r1_e->GetYaxis()->SetLabelSize(.1);
r1_e->GetXaxis()->SetTitleSize(.09);
r1_e->GetYaxis()->SetTitleSize(.095);
r1_e->GetXaxis()->SetTitleOffset(1.0);
r1_e->GetYaxis()->SetTitleOffset(0.6);
//leg->Draw("SAME");
//return;
// Eta Jets
TCanvas *c2=getaCanvas("c2");
c2->Divide(1,2);
// TPad *p1=getaPad_up("p1"); p1->Draw();
// p1->cd();
TPad *pad3 = new TPad("pad3", "The pad with the function",0.05,0.4,0.95,0.95);
pad3->Draw();
// pad3->Range(-112.6742,-73.17708,1143.438,551.3021);
pad3->SetFillColor(0);
pad3->SetBorderMode(0);
pad3->SetBorderSize(2);
// pad3->SetGridx();
// pad3->SetGridy();
pad3->SetLeftMargin(0.1271439);
pad3->SetRightMargin(0.07307979);
pad3->SetTopMargin(0.08215179);
pad3->SetBottomMargin(0.117181);
pad3->SetFrameBorderMode(0);
pad3->SetFrameBorderMode(0);
TPad *pad4 = new TPad("pad4", "The pad with the histogram",0.05,0.1,0.95,0.4);
pad4->Draw();
// pad4->Range(-153.3652,-2.142584,1185.427,5.367464);
pad4->SetFillColor(0);
pad4->SetBorderMode(0);
pad4->SetBorderSize(2);
pad4->SetGridx();
pad4->SetGridy();
pad4->SetLeftMargin(0.1215511);
pad4->SetRightMargin(0.07867263);
pad4->SetTopMargin(0.04892967);
pad4->SetBottomMargin(0.3521407);
pad4->SetFrameBorderMode(0);
pad4->SetFrameBorderMode(0);
pad3->cd();
// c2->cd(1);
// gPad->SetGridx(); gPad->SetGridy();
TLegend *leg = legend();
leg->SetHeader(data);
// leg->AddEntry(etaJet_e,"BX=-1","LF");
leg->AddEntry(etaJet_c,"BX=0","LF");
// leg->AddEntry(etaJet_l,"BX=+1","LF");
leg->AddEntry(etaJet_e,"BX=-1","LF");
etaJet_c->Draw("EHIST");
etaJet_c->SetLineColor(1); etaJet_c->SetLineWidth(2);
etaJet_c->SetTitle(";L1Jet #eta;a.u.");
// etaJet_c->GetYaxis()->SetRangeUser(1.01,100000000.);
etaJet_e->Draw("EHISTSAMES");
etaJet_e->SetLineColor(1); etaJet_e->SetLineColor(2); etaJet_e->SetLineStyle(1);etaJet_e->SetLineWidth(2);
// etaJet_l->Draw("EHISTSAMES");
// etaJet_l->SetLineColor(1); etaJet_l->SetLineColor(4); etaJet_l->SetLineStyle(3);
leg->Draw("SAME");
pad4->cd();
// c2->cd(2);
gPad->SetGridx(); gPad->SetGridy();
TH1D *r2_e = etaJet_e->Clone();
r2_e->Divide(etaJet_c);
TH1D *r2_l = etaJet_l->Clone();
r2_l->Divide(etaJet_c);
TLegend *leg = legend();
leg->SetHeader(data);
leg->AddEntry(r2_e,"BX=-1","LF");
// leg->AddEntry(etaJet_c,"BX=0","LF");
// leg->AddEntry(r2_l,"BX=+1","LF");
r2_e->Draw("HIST");r2_e->SetLineColor(4);r2_e->SetLineWidth(2);
// r2_l->Draw("HISTSAME"); r2_l->SetLineColor(4);
r2_e->GetYaxis()->SetRangeUser(-0.01,1.0);
//r2_e->GetXaxis()->SetRangeUser(0.,300.);
r2_e->SetLineWidth(2); r2_l->SetLineWidth(2);
r2_e->SetTitle(";L1Jet #eta;Ratio ");
r2_e->GetXaxis()->SetTitleSize(.09);
r2_e->GetYaxis()->SetTitleSize(.095);
r2_e->GetXaxis()->SetTitleOffset(1.0);
r2_e->GetYaxis()->SetTitleOffset(0.6);
r2_e->GetXaxis()->SetLabelSize(.1);
r2_e->GetYaxis()->SetLabelSize(.1);
// leg->Draw("SAME");
// HF jets
TCanvas *c3=getaCanvas("c3");
c3->Divide(1,2);
TPad *pad5 = new TPad("pad5", "The pad with the function",0.05,0.4,0.95,0.95);
pad5->Draw();
// pad5->Range(-112.6742,-73.17708,1143.438,551.3021);
pad5->SetFillColor(0);
pad5->SetBorderMode(0);
pad5->SetBorderSize(2);
// pad5->SetGridx();
// pad5->SetGridy();
pad5->SetLeftMargin(0.1271439);
pad5->SetRightMargin(0.07307979);
pad5->SetTopMargin(0.08215179);
pad5->SetBottomMargin(0.117181);
pad5->SetFrameBorderMode(0);
pad5->SetFrameBorderMode(0);
TPad *pad6 = new TPad("pad6", "The pad with the histogram",0.05,0.1,0.95,0.4);
pad6->Draw();
// pad6->Range(-153.3652,-2.142584,1185.427,5.367464);
pad6->SetFillColor(0);
pad6->SetBorderMode(0);
pad6->SetBorderSize(2);
pad6->SetGridx();
pad6->SetGridy();
pad6->SetLeftMargin(0.1215511);
pad6->SetRightMargin(0.07867263);
pad6->SetTopMargin(0.04892967);
pad6->SetBottomMargin(0.3521407);
pad6->SetFrameBorderMode(0);
pad6->SetFrameBorderMode(0);
pad5->cd();
// c3->cd(1);
gPad->SetLogy();
// gPad->SetGridx(); gPad->SetGridy();
TLegend *leg = legend();
leg->SetHeader(data);
// leg->AddEntry(ptJetHF_e,"BX=-1","LF");
leg->AddEntry(ptJetHF_c,"old PMTs region: BX=0 Jet","LF");
leg->AddEntry(ptJetHF2_c,"upgraded PMTs region: BX=0 Jet","LF");
leg->AddEntry(ptJetHF_e,"old PMTs region: BX=-1 Jet","LF");
leg->AddEntry(ptJetHF2_e,"upgraded PMTs region: BX=-1 Jet","LF");
ptJetHF_c->Rebin(2); ptJetHF2_c->Rebin(2);
ptJetHF_e->Rebin(2); ptJetHF2_e->Rebin(2);
ptJetHF_c->Draw("EHIST");
ptJetHF2_c->Draw("EHISTSAME");
ptJetHF2_c->SetLineColor(1); //ptJetHF2_c->SetLineWidth(2);
// ptJetHF2_c->SetFillStyle(3002); ptJetHF2_c->SetFillColor(1);
ptJetHF_c->SetLineColor(1); ptJetHF_c->SetLineWidth(2);
ptJetHF_c->GetXaxis()->SetRangeUser(0.,260.);
ptJetHF_c->GetYaxis()->SetRangeUser(1.01,100000000.);
// ptJetHF_c->SetTitle(";L1 FwdJet E_{T} (GeV);a.u.");
if (mode=="rate") {
ptJetHF_c->SetTitle(";L1 FwdJet E_{T} Threshold (GeV);a.u.");
} else {
ptJetHF_c->SetTitle(";L1 FwdJet E_{T} (GeV);a.u.");
}
ptJetHF_e->Draw("EHISTSAME");
ptJetHF2_e->Draw("EHISTSAME");
ptJetHF2_e->SetLineColor(2);// ptJetHF2_e->SetLineWidth(2);
// ptJetHF2_e->SetFillStyle(3002);ptJetHF2_e->SetFillColor(2);
ptJetHF_e->SetLineColor(2);
ptJetHF_e->SetLineStyle(1); ptJetHF_e->SetLineWidth(2);
// ptJetHF_l->Draw("EHISTSAMES");
// ptJetHF_l->SetLineColor(1); ptJetHF_l->SetLineColor(4);
// ptJetHF_l->SetLineStyle(3);
leg->Draw("SAME");
// TPad *p2=getaPad_dn("p2"); p2->Draw();
// p2->cd();
// c3->cd(2);
pad6->cd();
gPad->SetGridx(); gPad->SetGridy();
TH1D *r1HF_e = ptJetHF_e->Clone();
r1HF_e->Divide(ptJetHF_c);
TH1D *r1HF2_e = ptJetHF2_e->Clone();
r1HF2_e->Divide(ptJetHF2_c);
TLegend *leg = legend();
leg->AddEntry(r1HF_e,"old PMTs region","LF");
leg->AddEntry(r1HF2_e,"upgraded PMTs region","LF");
// leg->AddEntry(ptJetHF_c,"BX=0","LF");
// leg->AddEntry(r1HF_l,"BX=+1","LF");
r1HF_e->Draw("EHIST");r1HF_e->SetLineColor(4);
// r1HF_l->Draw("EHISTSAME"); r1HF_l->SetLineColor(4);
r1HF2_e->Draw("EHISTSAME"); //r1HF2_e->SetLineColor(kBlack-6);
//r1HF2_e->SetFillColor(4);
r1HF2_e->SetLineColor(4); r1HF2_e->SetLineWidth(0.5);
r1HF_e->GetYaxis()->SetRangeUser(-0.01,1.0);
r1HF_e->GetXaxis()->SetRangeUser(0.,260.);
r1HF_e->SetLineWidth(2); //r1HF2_e->SetLineWidth(1);
// r1HF_e->SetTitle(";L1 FwdJet E_{T} (GeV);ratio");
if (mode=="rate") {
r1HF_e->SetTitle(";L1 FwdJet E_{T} Threshold (GeV);Ratio");
} else {
r1HF_e->SetTitle(";L1 FwdJet E_{T} (GeV);Ratio");
}
r1HF_e->GetXaxis()->SetTitleSize(.09);
r1HF_e->GetYaxis()->SetTitleSize(.095);
r1HF_e->GetXaxis()->SetTitleOffset(1.0);
r1HF_e->GetYaxis()->SetTitleOffset(0.6);
r1HF_e->GetXaxis()->SetLabelSize(.1);
r1HF_e->GetYaxis()->SetLabelSize(.1);
leg->Draw("SAME");
TCanvas *cf=getanotherCanvas("cf");
gPad->SetGridx(); gPad->SetGridy();
TLegend *leg = legend();
leg->AddEntry(r1HF_e,"old PMTs region","LF");
leg->AddEntry(r1HF2_e,"upgraded PMTs region","LF");
r1HF_e->Draw("EHIST");r1HF_e->SetLineColor(4);
// r1HF2_e->SetLineColor(1);
r1HF2_e->Draw("EHISTSAME"); //r1HF2_e->SetLineColor(kBlack-6);
r1HF2_e->SetLineWidth(2); r1HF2_e->SetLineStyle(2);
r1HF_e->GetXaxis()->SetRangeUser(0.,170.);
r1HF_e->GetYaxis()->SetRangeUser(0.001,0.3);
r1HF_e->GetXaxis()->SetTitleSize(.05);
r1HF_e->GetYaxis()->SetTitleSize(.05);
r1HF_e->GetXaxis()->SetTitleOffset(1.0);
r1HF_e->GetYaxis()->SetTitleOffset(0.6);
r1HF_e->GetXaxis()->SetLabelSize(.05);
r1HF_e->GetYaxis()->SetLabelSize(.05);
leg->Draw("SAME");
}
// -----------------------------------------------------------------------------
//
TCanvas* aDrawBkgdPlots( TString path,
TString canvas_name,
TString name,
TString dirmame,
int rebin,
bool norm,
bool log,
TDirectory* file )
{
// SetSomeStyles();
// Create legend
TLegend* legend = new TLegend( 0.75, 0.65, 0.92, 0.92, NULL, "brNDC" );
legend->SetFillColor(0);
legend->SetLineColor(0);
// Create canvas
TCanvas* aCanvas = getaCanvas( canvas_name, file, log );
// Create histograms
TH1D* qcd = readHist( path, name, "IC5Calo_QCD_Pythia_Merged.root", dirmame, rebin );
TH1D* tt_jets = readHist( path, name, "IC5Calo_TTbarJets.root", dirmame, rebin );
TH1D* w_jets = readHist( path, name, "IC5Calo_WJets.root", dirmame, rebin );
TH1D* z_inv = readHist( path, name, "IC5Calo_Zinv.root", dirmame, rebin );
TH1D* z_jets = readHist( path, name, "IC5Calo_ZJets.root", dirmame, rebin );
TH1D* lm0 = readHist( path, name, "IC5Calo_LM0.root", dirmame, rebin );
TH1D* lm1 = readHist( path, name, "IC5Calo_LM1.root", dirmame, rebin );
// Combine Z+jets and Z->inv
TH1D* z_all = z_inv->Clone();
z_all->Add(z_jets,1);
// TH1D* z_all = 0;
// if ( z_inv && z_jets ) {
// z_all = z_inv->Clone();
// z_all->Add(z_jets,1);
// } else if ( z_inv ) {
// z_all = z_inv->Clone();
// } else if ( z_jets ) {
// z_all = z_jets->Clone();
// }
// Select Z+jets and Z->inv separate or not
bool combine = true;
// Line colour and fill
if ( qcd ) qcd->SetLineColor(kGreen+2);
if ( qcd ) qcd->SetFillColor(kGreen+2);
if ( qcd ) qcd->SetFillStyle(3003);
if ( tt_jets ) tt_jets->SetLineColor(kBlue);
if ( tt_jets ) tt_jets->SetLineStyle(1);
if ( tt_jets ) tt_jets->SetLineWidth(1);
w_jets->SetLineColor(kBlue);
w_jets->SetLineStyle(3);
w_jets->SetLineWidth(1);
if ( combine ) {
z_all->SetLineColor(kBlack);
z_all->SetLineStyle(3);
z_all->SetLineWidth(1);
} else {
z_inv->SetLineColor(kBlack);
z_inv->SetLineStyle(1);
z_inv->SetLineWidth(1);
z_jets->SetLineColor(kBlack);
z_jets->SetLineStyle(3);
z_jets->SetLineWidth(1);
}
lm0->SetLineColor(kRed);
lm0->SetLineStyle(1);
lm0->SetLineWidth(2);
lm1->SetLineColor(kRed);
lm1->SetLineStyle(3);
lm1->SetLineWidth(2);
// Populate legend
legend->AddEntry( qcd, " QCD", "f" );
legend->AddEntry( lm0, " SUSY LM0", "L" );
legend->AddEntry( lm1, " SUSY LM1", "L" );
legend->AddEntry( tt_jets, " t#bar{t}+jets", "L" );
legend->AddEntry( w_jets, " W+jets", "L" );
if ( combine ) {
legend->AddEntry( z_all, " Z", "L" );
} else {
legend->AddEntry( z_jets, " Z+jets", "L" );
legend->AddEntry( z_inv, " Z#rightarrow#nu#nu", "L" );
}
// Calc maximum number of entries
double aMax = 0.;
if ( qcd->GetMaximum() > aMax ) { aMax = qcd->GetMaximum(); }
if ( lm0->GetMaximum() > aMax ) { aMax = lm0->GetMaximum(); }
if ( lm1->GetMaximum() > aMax ) { aMax = lm1->GetMaximum(); }
if ( tt_jets->GetMaximum() > aMax ) { aMax = tt_jets->GetMaximum(); }
if ( w_jets->GetMaximum() > aMax ) { aMax = w_jets->GetMaximum(); }
if ( combine ) {
if ( z_all->GetMaximum() > aMax ) { aMax = z_all->GetMaximum(); }
} else {
if ( z_inv->GetMaximum() > aMax ) { aMax = z_inv->GetMaximum(); }
if ( z_jets->GetMaximum() > aMax ) { aMax = z_jets->GetMaximum(); }
}
// Calc minimum number of entries
double aMin = 1.e12;
if ( qcd->GetMinimum(1.e-12) < aMin ) { aMin = qcd->GetMinimum(1.e-12); }
if ( lm0->GetMinimum(1.e-12) < aMin ) { aMin = lm0->GetMinimum(1.e-12); }
if ( lm1->GetMinimum(1.e-12) < aMin ) { aMin = lm1->GetMinimum(1.e-12); }
if ( tt_jets->GetMinimum(1.e-12) < aMin ) { aMin = tt_jets->GetMinimum(1.e-12); }
if ( w_jets->GetMinimum(1.e-12) < aMin ) { aMin = w_jets->GetMinimum(1.e-12); }
if ( combine ) {
if ( z_all->GetMinimum(1.e-12) < aMin ) { aMin = z_all->GetMinimum(1.e-12); }
} else {
if ( z_inv->GetMinimum(1.e-12) < aMin ) { aMin = z_inv->GetMinimum(1.e-12); }
if ( z_jets->GetMinimum(1.e-12) < aMin ) { aMin = z_jets->GetMinimum(1.e-12); }
}
if ( qcd ) qcd->GetYaxis()->SetTitleOffset(1.43);
if ( qcd ) qcd->GetYaxis()->SetTitleSize(0.06);
if ( qcd ) qcd->GetXaxis()->SetTitleSize(0.06);
if ( qcd ) qcd->GetXaxis()->SetTitleOffset(0.9);
if ( log ) {
if ( qcd ) qcd->SetMaximum( aMax * 10. );
if ( qcd ) qcd->SetMinimum( aMin * 0.1 );
} else {
if ( qcd ) qcd->SetMaximum( aMax * 1.1 );
if ( qcd ) qcd->SetMinimum( aMin * 0.9 );
}
if ( norm ) {
if ( qcd ) qcd->DrawNormalized("Ehist");
if ( lm0->GetEntries() > 0. ) { lm0->DrawNormalized("hsame"); }
if ( lm1->GetEntries() > 0. ) { lm1->DrawNormalized("hsame"); }
if ( tt_jets->GetEntries() > 0. ) { tt_jets->DrawNormalized("hsame"); }
if ( w_jets->GetEntries() > 0. ) { w_jets->DrawNormalized("hsame"); }
if ( combine ) {
if ( z_all->GetEntries() > 0. ) { z_all->DrawNormalized("hsame"); }
} else {
if ( z_inv->GetEntries() > 0. ) { z_inv->DrawNormalized("hsame"); }
if ( z_jets->GetEntries() > 0. ) { z_jets->DrawNormalized("hsame"); }
}
} else {
if ( qcd ) qcd->Draw("h");
lm0->Draw("sameH");
lm1->Draw("sameH");
if ( tt_jets ) tt_jets->Draw("sameh");
w_jets->Draw("sameH");
if ( combine ) {
z_all->Draw("sameH");
} else {
z_inv->Draw("sameH");
z_jets->Draw("sameH");
}
}
file->cd();
legend->Draw("same");
aCanvas->Write();
return aCanvas;
}
