void plotsNoFitETau(bool mine=false,TString name="collMass",TString file="LFV_vbf_collMass_fakeRate_zjetsEmbed_newSignal.root", TString dir="preselectionMetFix_Nov2", TString dirInternal="vbfmutau",double blindA=100, double blindB=160, bool blind=false, TString Xaxis="M_{#mu,#tau}_{coll} [GeV]", TString Yaxis="Events / 20 GeV", TString Tex="#mu#tau_{h} 0 Jet", double xmin=0, double xmax=300, int rebinning=1, bool setLogY=false, double legx1=0.6, double legy1=0.9, double legx2=0.9, double legy2=0.5, double MAX=-1, double MIN=-1, double CORRFR=1., double scaleSignal=1.0, double ymin=-0.95, double ymax=0.95){
//gROOT->LoadMacro("tdrstyle.C");
setTDRStyle();
//gROOT->LoadMacro("CMS_lumi.C");
writeExtraText = false;
int iPeriod = 2;
// second parameter in example_plot is iPos, which drives the position of the CMS logo in the plot
int iPos=11;// : top-left, left-aligned
// iPos=33 : top-right, right-aligned
// iPos=22 : center, centered
// mode generally :
// iPos = 10*(alignement 1/2/3) + position (1/2/3 = left/center/right)
// example_plot( iPeriod, 11 ); // left-aligned
// example_plot( iPeriod, 33 ); // right-aligned
// example_plot( iPeriod, 0 ); // out of frame (in exceptional cases)
// example_plot( iPeriod, 11 ); // default: left-aligned
// example_plot( iPeriod, 22 ); // centered
// example_plot( iPeriod, 33 ); // right-aligned
double kForPlotting=1;
double branchingratioTauTau=0.063;
double branchingratioTauMu=0.1;
double Lumi=19717;
// Get Plots
// Get Plots
TFile *_file0= new TFile(dir+"/"+file);
TString find="DYrightarrowlljets";
if(!mine) find="zjetsother";
TH1F *hDY=_file0->Get(dirInternal+find); hDY->SetName("DY");
find="Diboson";
if(!mine) find="ww";
TH1F *hWW=_file0->Get(dirInternal+find); hWW->SetName("WW");
find="SingleTop";
if(!mine) find="singlet";
TH1F *hTOP=_file0->Get(dirInternal+find); hTOP->SetName("TOP");
find="tbart";
if(!mine) find="ttbar";
TH1F *hTT=_file0->Get(dirInternal+find); hTT->SetName("TT");
find="Zrightarrowtautauembedded";
if(!mine) find="ztautau";
TH1F *hZTauTau=_file0->Get(dirInternal+find); hZTauTau->SetName("ZTauTau");
find="Fakes";
if(!mine) find="fakes";
TH1F *hFAKES=_file0->Get(dirInternal+find); hFAKES->SetName("fakes");
hFAKES->Scale(CORRFR);
TH1F *hLFVVBF126=_file0->Get(dirInternal+"LFVvbfHiggsBR=1"); hLFVVBF126->SetName("LFVVBF126");
TH1F *hSMVBF126=_file0->Get(dirInternal+"SMvbfHiggs"); hSMVBF126->SetName("SMVBF126");
TH1F *hLFVGG126=_file0->Get(dirInternal+"LFVggHiggsBR=1"); hLFVGG126->SetName("LFVGG126");
TH1F *hSMGG126=_file0->Get(dirInternal+"SMHiggs"); hSMGG126->SetName("SMGG126");
find="Observed";
if(!mine) find="data_obs";
TH1F* hdata_obsNoErrFix=_file0->Get(dirInternal+find);
TH1F* hdata_obs = hdata_obsNoErrFix->Clone();
std::cout << hdata_obsNoErrFix->GetEntries() << std::endl;
for (i = 1; i <= hdata_obs->GetNbinsX(); i++){
if (hdata_obs->GetBinContent(i) == 0){
hdata_obs->SetBinContent(i,0.0);
hdata_obs->SetBinError(i,1.8);
}
}
hdata_obs->SetName("data_obs");
TH1F* hSMHIGGS=hSMGG126->Clone(); hSMHIGGS->SetName("HIGGSSM");
hSMHIGGS->Add(hSMVBF126);
// For the Control plots only we want signal to be scaled to 100% Br
hLFVVBF126->Scale(scaleSignal);
hLFVGG126->Scale(scaleSignal);
// From the fit
/* hDY->Scale(1.122280);
hWW->Scale(1.254172);
hTOP->Scale(1.348162);
hTT->Scale(1.082784);
hZTauTau->Scale(1.157596);
hFAKES->Scale(0.901000);
*/
// Daniel's Colors
/*
WGammaStar=kCyan
ZTauTauEmbedded=kOrange-4
ZLL_residual=kAzure+3
TTBar=40
SingleTop=kGreen-2
EWKDiBoson=kRed+2
WJets/QCD Multijets=kMagenta-10
SMHToTauTau=kMagenta
*/
hFAKES->SetFillColor(kMagenta-10); hFAKES->SetLineColor(kMagenta+4); hFAKES->SetLineWidth(1);
hZTauTau->SetFillColor(kOrange-4); hZTauTau->SetLineColor(kOrange+4); hZTauTau->SetLineWidth(1);
hDY->SetFillColor(kAzure+3); hDY->SetLineColor(kAzure+4); hDY->SetLineWidth(1);
hWW->SetFillColor(kAzure+3); hWW->SetLineColor(kAzure+3); hWW->SetLineWidth(1);
hTOP->SetFillColor(kGreen-2); hTOP->SetLineColor(kGreen+4); hTOP->SetLineWidth(1);
hTT->SetFillColor(kGreen-2); hTT->SetLineColor(kGreen-2); hTT->SetLineWidth(1);
hLFVGG126->SetLineColor(kBlue+1); hLFVGG126->SetLineWidth(3);
hLFVVBF126->SetLineColor(kBlue+1); hLFVVBF126->SetLineWidth(3); hLFVVBF126->SetLineStyle(kDashed);
hSMVBF126->SetLineColor(kMagenta); hSMVBF126->SetLineWidth(3); hSMVBF126->SetLineStyle(kDashed);
hSMGG126->SetLineColor(kMagenta); hSMGG126->SetLineWidth(3);
hSMHIGGS->SetFillColor(kMagenta); hSMHIGGS->SetLineColor(kMagenta+1); hSMHIGGS->SetLineWidth(1);
hdata_obs->SetMarkerSize(1); // Closer to Daniel's
// Rebin
hFAKES->Rebin(rebinning);
hZTauTau->Rebin(rebinning);
hDY->Rebin(rebinning);
hTOP->Rebin(rebinning);
hTT->Rebin(rebinning);
hWW->Rebin(rebinning);
hLFVGG126->Rebin(rebinning);
hLFVVBF126->Rebin(rebinning);
hSMVBF126->Rebin(rebinning);
hSMGG126->Rebin(rebinning);
hdata_obs->Rebin(rebinning);
hSMHIGGS->Rebin(rebinning);
// PLOT
TCanvas *c1 = new TCanvas("canvas_"+name);
TPad *Pad1= new TPad("pad1","",0,0.2,1,1); Pad1->Draw(); Pad1->cd();;
Pad1->SetLeftMargin(0.2147651);
Pad1->SetRightMargin(0.06543624);
Pad1->SetTopMargin(0.07);
Pad1->SetBottomMargin(0.04);
for (int i=0; i<hFAKES->GetNbinsX()+1; i++){
double content=hFAKES->GetBinContent(i);
double contentErr=0;//sqrt(content);//hFAKES->GetBinError(i);
double err=sqrt(0.3*0.3*content*content+contentErr*contentErr);
hFAKES->SetBinError(i,err);
}
TH1F* fullMC2=hZTauTau->Clone(); fullMC2->Add(hFAKES); fullMC2->Add(hTT); fullMC2->Add(hWW);
fullMC2->Add(hDY); fullMC2->Add(hTOP);
fullMC2->SetFillColorAlpha(kGray+2, 0.35);
//fullMC2->SetFillStyle(3002);
fullMC2->SetMarkerSize(0);
fullMC2->Draw("hist");
fullMC2->GetXaxis()->SetTitle("");
fullMC2->GetYaxis()->SetTitle(Yaxis);
fullMC2->GetXaxis()->SetRangeUser(xmin,xmax);
fullMC2->GetYaxis()->SetTitleOffset(1.2);
fullMC2->GetYaxis()->SetTitleSize(0.05);
fullMC2->GetXaxis()->SetNdivisions(505);
fullMC2->GetYaxis()->SetLabelSize(0.04);
fullMC2->GetXaxis()->SetLabelSize(0);
THStack* stack = new THStack("stack","");
stack->Add(hFAKES);
stack->Add(hWW);
stack->Add(hDY);
stack->Add(hTT);
stack->Add(hTOP);
stack->Add(hZTauTau);
stack->Add(hSMHIGGS);
cout<<"Yields"<<endl;
cout<<"VV " <<hWW->Integral()<<endl;
cout<<"TOP " <<hTOP->Integral()<<endl;
cout<<"TT " <<hTT->Integral()<<endl;
cout<<"DY " <<hDY->Integral()<<endl;
cout<<"ZTauTau " <<hZTauTau->Integral()<<endl;
cout<<"FAKES " <<hFAKES->Integral()<<endl;
cout<<"LFVGG126 " <<hLFVGG126->Integral()<<endl;
cout<<"LFVVBF126 " <<hLFVVBF126->Integral()<<endl;
stack->Draw("samehist");
fullMC2->Draw("sames,E2");
hLFVVBF126->Draw("sameshist");
hLFVGG126->Draw("sameshist");
double maxData=hdata_obs->GetMaximum();
double maxMC=stack->GetMaximum()*1.2;
double maxLFV=hLFVGG126->GetMaximum();
double maxLFV2=hLFVVBF126->GetMaximum();
double minMC=stack->GetMinimum();
if(maxData>maxMC) {maxMC=1.2*maxData;}
if(maxLFV>maxMC) {maxMC=1.2*maxLFV;}
if(maxLFV2>maxMC) {maxMC=1.2*maxLFV2;}
if(MAX!=-1) {maxMC=MAX;}
stack->SetMaximum(maxMC);
stack->GetYaxis()->SetRangeUser(minMC,maxMC);
fullMC2->SetMaximum(maxMC);
fullMC2->GetYaxis()->SetRangeUser(minMC,maxMC);
if(setLogY){fullMC2->SetMinimum(0.1); fullMC2->GetYaxis()->SetRangeUser(1,maxMC*1000);}
TH1F* histoDataUnblindedV4=hdata_obs->Clone(); histoDataUnblindedV4->SetName("data_unblinded");
if(!blind) hdata_obs->Draw("sames");
else {
int findBinA=hdata_obs->FindBin(blindA);
int findBinB=hdata_obs->FindBin(blindB);
for (int i=findBinA; i<findBinB; i++) hdata_obs->SetBinContent(i,-1000);
hdata_obs->Draw("sames");
TPave *pave = new TPave(blindA,minMC,blindB,maxMC,4,"br");
//pave->SetFillColorAlpha(kGray+3,0.35);
pave->SetFillColor(1);
pave->SetFillStyle(3003);
pave->SetDrawOption(0);
pave->SetBorderSize(0);
pave->Draw();
}
TLegend *leg = new TLegend(legx1,legy1,legx2,legy2,NULL,"brNDC");
leg->SetFillColor(0);
leg->SetBorderSize(0);
leg->SetFillStyle(0);
TLegendEntry *entry;
entry=leg->AddEntry(hdata_obs,"Data, e#tau_{h}","p");
entry=leg->AddEntry(fullMC2,"Bkgd. uncertainty","f");
entry=leg->AddEntry(hSMHIGGS,"SM Higgs","f");
eblindBy=leg->AddEntry(hZTauTau,"Z#rightarrow#tau#tau","f");
entry=leg->AddEntry(hTOP,"t#bar{t}, t, #bar{t}","f");
entry=leg->AddEntry(hDY,"Other","f");
entry=leg->AddEntry(hFAKES,"MisID'd #tau","f");
entry=leg->AddEntry(hLFVGG126,"LFV GG Higgs (B=100%)","l");
entry=leg->AddEntry(hLFVVBF126,"LFV VBF Higgs (B=100%)","l");
leg->Draw();
CMS_lumi( Pad1, iPeriod, iPos );
TLatex latex;
latex.SetNDC();
latex.SetTextAngle(0);
latex.SetTextColor(kBlack);
latex.SetTextFont(42);
latex.SetTextAlign(31);
latex.SetTextSize(.05);
latex.DrawLatex(0.39,0.8,Tex);
//cmsPrelim(Lumi);
Pad1->SetLogy(setLogY);
gPad->RedrawAxis();
c1->cd(); TPad *Pad2= new TPad("pad2","",0,0,1,0.23); Pad2->Draw(); Pad2->cd(); Pad2->SetGridy();
Pad2->SetLeftMargin(0.2147651);
Pad2->SetRightMargin(0.06543624);
Pad2->SetTopMargin(0.0);
Pad2->SetBottomMargin(0.37);
Pad2->SetFillStyle(0);
TH1F * Ratio=fullMC2->Clone(); Ratio->SetName("Ratio");
Ratio->Scale(-1);
Ratio->Add(hdata_obs,1);
Ratio->Divide(fullMC2);
Ratio->SetLineColor(kBlack);
Ratio->GetXaxis()->SetLabelFont(42);
Ratio->GetXaxis()->SetTitleFont(42);
Ratio->GetYaxis()->SetNdivisions(505);
Ratio->GetYaxis()->SetLabelFont(42);
Ratio->GetYaxis()->SetLabelSize(0.122);
Ratio->GetYaxis()->SetRangeUser(ymin,ymax);
Ratio->GetXaxis()->SetRangeUser(xmin,xmax);
Ratio->GetXaxis()->SetLabelSize(0.12);
Ratio->GetXaxis()->SetLabelFont(42);
Ratio->SetYTitle("#frac{Data-Bkgd}{Bkgd}");
Ratio->SetXTitle(Xaxis);
Ratio->GetXaxis()->SetNdivisions(505);
Ratio->GetYaxis()->CenterTitle(true);
Ratio->GetYaxis()->SetTitleOffset(0.4);
Ratio->GetYaxis()->SetTitleSize(0.11);
Ratio->GetXaxis()->SetTitleOffset(0.75);
Ratio->GetXaxis()->SetTitleSize(0.20);
Ratio->SetMarkerSize(1.);
TH1F* RatioError = Ratio->Clone(); RatioError->SetName("RatioError");
for (int i=0; i<RatioError->GetNbinsX()+1; i++){
double error=fullMC2->GetBinError(i)*hdata_obs->GetBinContent(i)/fullMC2->GetBinContent(i)/fullMC2->GetBinContent(i);
RatioError->SetBinContent(i,0);
RatioError->SetBinError(i,error);
}
//RatioError->SetFillStyle(3002);
RatioError->SetFillColorAlpha(kGray+2,0.35); RatioError->SetMarkerSize(0);
RatioError->Draw("E2");
//Ratio->Draw("samesE0"); //errorbarswithoutdatapoints
Ratio->Draw("sames");
if(!blind) Ratio->Draw("sames");
else {
int findBinA=Ratio->FindBin(blindA);
int findBinB=Ratio->FindBin(blindB);
for (int i=findBinA; i<findBinB; i++) Ratio->SetBinContent(i,-100);
Ratio->Draw("sames");
TPave *pave = new TPave(blindA,-1.0,blindB,1.0,4,"br");
pave->SetFillColor(1);
pave->SetFillStyle(3003);
//pave->SetFillColorAlpha(kGray+3,0.35);
pave->SetDrawOption(0);
pave->SetBorderSize(0);
pave->Draw();
}
if(!setLogY){
// c1->SaveAs(name+"ErrorBarsWithoutDataPoints.png");
// c1->SaveAs(name+"ErrorBarsWithoutDataPoints.pdf");
c1->SaveAs(name+".png");
c1->SaveAs(name+".pdf");
}
else {
c1->SaveAs(name+"_log.png");
c1->SaveAs(name+"_log.pdf");
}
}
void HistoData::drawResult(TH1 *Summary, bool Vertical, bool SetBinLabel) {
// add label to the summary if desired
if (SetBinLabel) {
Summary->GetXaxis()->SetBinLabel(bin,getRefHisto()->GetTitle());
//Summary->GetXaxis()->SetBinLabel(bin,name.c_str());
}
else
Summary->GetXaxis()->SetBinLabel(bin,name.c_str());
double minimum = Summary->GetMinimum();
// determine where to draw the result (score axis)
// 1: solid bar starts
// 2: solid bar ends, hatched bar starts
// 3: hatched bar ends
double score1 = minimum;
double score2 = (lowScore == 10. || lowScore < minimum) ? minimum : lowScore;
double score3 = (lowScore == 10.) ? 1 : ((highScore < minimum) ? minimum : highScore);
// determine where to draw the result (binning axis)
double binCenter = Summary->GetBinCenter(bin);
double binWidth = Summary->GetBinWidth(bin);
double bin1 = binCenter - binWidth / 3;
double bin2 = binCenter + binWidth / 3;
// set coordinates of display based on plot alignment
double solidX1, solidY1, solidX2, solidY2;
double hatchedX1, hatchedY1, hatchedX2, hatchedY2;
double axisX1, axisY1, axisX2, axisY2;
if (Vertical) {
solidX1 = bin1; solidX2 = bin2;
solidY1 = score1; solidY2 = score2;
hatchedX1 = bin1; hatchedX2 = bin2;
hatchedY1 = score2; hatchedY2 = score3;
axisX1 = bin1; axisX2 = bin2;
axisY1 = minimum; axisY2 = minimum;
} else {
solidX1 = score1; solidX2 = score2;
solidY1 = bin1; solidY2 = bin2;
hatchedX1 = score2; hatchedX2 = score3;
hatchedY1 = bin1; hatchedY2 = bin2;
axisX1 = minimum; axisX2 = minimum;
axisY1 = bin1; axisY2 = bin2;
}
// a solid bar is drawn from zero to the lowest score
if (lowScore > minimum && lowScore != 10.) {
TPave *solidBar = new TPave(solidX1,solidY1,solidX2,solidY2,1,"");
solidBar->SetBit(kCanDelete);
solidBar->SetLineColor(lineUseFillColor ? solidFillColor : solidLineColor);
solidBar->SetFillColor(solidFillColor);
solidBar->SetFillStyle(solidFillStyle);
solidBar->Draw();
}
// a hatched bar is drawn from the lowest score to the highest score
if ((lowScore != highScore && highScore > minimum) || lowScore == 10.) {
TPave *hatchedBar = new TPave(hatchedX1,hatchedY1,hatchedX2,hatchedY2,1,"");
hatchedBar->SetBit(kCanDelete);
hatchedBar->SetLineColor(lineUseFillColor ? shadedFillColor : shadedLineColor);
hatchedBar->SetFillColor(shadedFillColor);
hatchedBar->SetFillStyle(3004);
hatchedBar->Draw();
}
// paste a line over the base axis to fix border color artifacts
TLine *axisLine = new TLine(axisX1,axisY1,axisX2,axisY2);
if (Vertical) axisLine->SetLineColor(Summary->GetAxisColor("X"));
else axisLine->SetLineColor(Summary->GetAxisColor("Y"));
axisLine->SetBit(kCanDelete);
axisLine->Draw("SAME");
//paste a line before (in this proceeding that means after) Barrel and Endcap Plots
if ( name == "ERPt" || name == "BRPt" ) {
axisY1=axisY2=binCenter+binWidth/2;
axisX2=Summary->GetMaximum();
axisX1=Summary->GetMinimum()-200;
TLine *regionLine = new TLine(axisX1,axisY1,axisX2,axisY2);
regionLine->SetLineColor(Summary->GetAxisColor("X"));
regionLine->SetBit(kCanDelete);
regionLine->SetLineWidth(3);
regionLine->Draw();
}
}
