| 359 |
|
oneline->SetLineWidth(1); |
| 360 |
|
oneline->Draw("same"); |
| 361 |
|
} |
| 362 |
+ |
|
| 363 |
|
float make_one_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, float ymax, bool logscale, |
| 364 |
|
string xlabel, string filename, bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) { |
| 365 |
|
|
| 437 |
|
|
| 438 |
|
} |
| 439 |
|
|
| 440 |
< |
void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale, string xlabel, string filename, |
| 441 |
< |
bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) { |
| 440 |
> |
// Compare data to data |
| 441 |
> |
float make_data_comparison_plot(string variable, TCut cut, int nbins, float min, float max, float ymax, bool logscale, |
| 442 |
> |
string xlabel, string filename, float legendPosition=0.55) { |
| 443 |
> |
|
| 444 |
> |
TCut ibasiccut=basiccut&&cut; |
| 445 |
> |
|
| 446 |
> |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 447 |
> |
TPad* rcan = new TPad("rcan","rcan",0,0,1,1); |
| 448 |
> |
rcan->SetLogy(logscale); |
| 449 |
> |
rcan->cd(); |
| 450 |
> |
|
| 451 |
> |
std::cout << "Data comparison: variable = " << variable << ", cut = " << ibasiccut.GetTitle() << std::endl; |
| 452 |
> |
|
| 453 |
> |
TH1F *data1 = allsamples.Draw("data1",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity); |
| 454 |
> |
TH1F *data2 = comparesamples.Draw("data2",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity); |
| 455 |
> |
|
| 456 |
> |
data1->SetMarkerSize(DataMarkerSize); |
| 457 |
> |
data2->SetLineColor(kRed); |
| 458 |
> |
|
| 459 |
> |
if ( !logscale ) { |
| 460 |
> |
data1->SetMinimum(0); // Defaults |
| 461 |
> |
} else { |
| 462 |
> |
data1->SetMinimum(0.5); |
| 463 |
> |
} |
| 464 |
> |
if (ymax<0) { |
| 465 |
> |
if ( logscale ) data1->SetMaximum(5.3*data1->GetMaximum()); |
| 466 |
> |
else data1->SetMaximum(1.5*data1->GetMaximum()); |
| 467 |
> |
} else { |
| 468 |
> |
data1->SetMaximum(ymax); |
| 469 |
> |
} |
| 470 |
> |
|
| 471 |
> |
float ymaxSet = data1->GetMaximum(); |
| 472 |
> |
|
| 473 |
> |
data1->GetXaxis()->SetTitle(xlabel.c_str()); |
| 474 |
> |
data1->GetYaxis()->SetTitle("Events"); |
| 475 |
> |
data1->GetXaxis()->CenterTitle(); |
| 476 |
> |
data1->GetYaxis()->CenterTitle(); |
| 477 |
> |
|
| 478 |
> |
TLegend *mleg = make_legend("",legendPosition,0.7,false,legendPosition+0.2); |
| 479 |
> |
mleg->AddEntry(data1, "Old 5/fb", "PL"); |
| 480 |
> |
mleg->AddEntry(data2, "New x/fb", "L"); |
| 481 |
> |
|
| 482 |
> |
data1->Draw("E1"); |
| 483 |
> |
data2->Draw("HIST,SAMES"); |
| 484 |
> |
mleg->Draw(); |
| 485 |
> |
DrawPrelim(); |
| 486 |
> |
save_with_ratio( data1, data2, rcan, "compareData/" + filename, false, false, "old/new" ); |
| 487 |
> |
|
| 488 |
> |
data1->Delete(); |
| 489 |
> |
data2->Delete(); |
| 490 |
> |
delete mleg; |
| 491 |
> |
delete ckin; |
| 492 |
> |
|
| 493 |
> |
} |
| 494 |
> |
|
| 495 |
> |
void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale, |
| 496 |
> |
string xlabel, string filename, float legendPosition=0.55) { |
| 497 |
|
|
| 498 |
|
string mllcuts[] = { "mll>20","mll>20&&mll<70", "mll>75&&mll<105", "mll>120" }; |
| 499 |
|
string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "75 < m_{ll} < 105 GeV", "m_{ll} > 120 GeV" }; |
| 507 |
|
if ( addcut != "" ) mllcuts[i] += "&&"+addcut; |
| 508 |
|
if ( i==start ) { |
| 509 |
|
ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel, |
| 510 |
< |
filename+plotname[i], plotratio, loadlastminmax, legendPosition ); |
| 510 |
> |
filename+plotname[i], legendPosition ); |
| 511 |
|
} else { |
| 512 |
|
make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 513 |
< |
filename+plotname[i], plotratio, loadlastminmax, legendPosition ); |
| 513 |
> |
filename+plotname[i], legendPosition ); |
| 514 |
|
} |
| 515 |
|
make_one_OFSF_plot(variable, "id1==1&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 516 |
< |
filename+plotname[i]+"_mm", plotratio, loadlastminmax, legendPosition ); |
| 516 |
> |
filename+plotname[i]+"_mm", legendPosition ); |
| 517 |
|
make_one_OFSF_plot(variable, "id1==0&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 518 |
< |
filename+plotname[i]+"_ee", plotratio, loadlastminmax, legendPosition ); |
| 518 |
> |
filename+plotname[i]+"_ee", legendPosition ); |
| 519 |
|
} |
| 520 |
|
|
| 521 |
|
} |
| 733 |
|
bool nolog=false; |
| 734 |
|
|
| 735 |
|
bool doOFSF = false; |
| 736 |
< |
bool doKin = true; |
| 736 |
> |
bool doKin = false; |
| 737 |
> |
bool doDataComp = true; |
| 738 |
> |
|
| 739 |
|
|
| 740 |
|
if(doPF) write_warning(__FUNCTION__,"Please use caution when trying to produce PF plots; not all versions of the JZB trees have these variables!"); |
| 741 |
|
float mll_low=50; |
| 746 |
|
} |
| 747 |
|
|
| 748 |
|
if ( doOFSF ) { |
| 749 |
< |
make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll", false, false); |
| 750 |
< |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets", false, false); |
| 751 |
< |
make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets", false, false); |
| 752 |
< |
make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1", false, false); |
| 753 |
< |
make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2", false, false); |
| 754 |
< |
// make_OFSF_plots("pt1", "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1", false, false); |
| 755 |
< |
// make_OFSF_plots("pt2", "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2", false, false); |
| 756 |
< |
make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", false, false, 0.15); |
| 757 |
< |
make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", false, false, 0.15); |
| 758 |
< |
// make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", false, false, 0.2); |
| 759 |
< |
// make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", false, false, 0.2); |
| 760 |
< |
// make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", false, false, 0.2); |
| 761 |
< |
make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", false, false, 0.2); |
| 762 |
< |
make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", false, false, 0.2); |
| 763 |
< |
// make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", false, false, 0.2); |
| 764 |
< |
// make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", false, false, 0.2); |
| 765 |
< |
make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilc", false, false); |
| 766 |
< |
make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijc", false, false); |
| 767 |
< |
make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", "dphilco", false, false); |
| 768 |
< |
make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(jet,MET)|)", "dphijco", false, false); |
| 749 |
> |
make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll"); |
| 750 |
> |
|
| 751 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets"); |
| 752 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto"); |
| 753 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30"); |
| 754 |
> |
|
| 755 |
> |
make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets"); |
| 756 |
> |
make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1"); |
| 757 |
> |
make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2"); |
| 758 |
> |
// make_OFSF_plots("pt1", "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1"); |
| 759 |
> |
// make_OFSF_plots("pt2", "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2"); |
| 760 |
> |
make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", 0.15); |
| 761 |
> |
make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", 0.15); |
| 762 |
> |
// make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", 0.2); |
| 763 |
> |
// make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", 0.2); |
| 764 |
> |
// make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", 0.2); |
| 765 |
> |
make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", 0.2); |
| 766 |
> |
make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", 0.2); |
| 767 |
> |
// make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", 0.2); |
| 768 |
> |
// make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", 0.2); |
| 769 |
> |
make_OFSF_plots("TMath::Min(TMath::Abs(dphiMet1), TMath::Abs(dphiMet2))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(l,MET)|)", "dphilc"); |
| 770 |
> |
make_OFSF_plots("TMath::Min(TMath::Abs(pfJetDphiMet[0]), TMath::Min(TMath::Abs(pfJetDphiMet[1]), TMath::Abs(pfJetDphiMet[2])))", "met[4]>100", 16, 0, 3.2, false, "Min(|#Delta#phi(jet,MET)|)", "dphijc"); |
| 771 |
> |
make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(dphiMet1)), (TMath::Pi() - TMath::Abs(dphiMet2)))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(l,MET)|)", "dphilco"); |
| 772 |
> |
make_OFSF_plots("TMath::Min((TMath::Pi()-TMath::Abs(pfJetDphiMet[0])), TMath::Min( (TMath::Pi()-TMath::Abs(pfJetDphiMet[1])), (TMath::Pi() - TMath::Abs(pfJetDphiMet[2]))))", "met[4]>100", 16, 0, 3.2, false, "Min(#pi - |#Delta#phi(jet,MET)|)", "dphijco"); |
| 773 |
> |
} |
| 774 |
> |
|
| 775 |
> |
if ( doDataComp && !PlottingSetup::openBox ) { |
| 776 |
> |
TCut mllCut(""); |
| 777 |
> |
if ( !PlottingSetup::openBox ) mllCut = "mll>120"; |
| 778 |
> |
|
| 779 |
> |
TCut cutSignal = cutmass&&cutnJets&&"met[4]>100"; |
| 780 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut, 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig"); |
| 781 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee"); |
| 782 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm"); |
| 783 |
> |
|
| 784 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto"); |
| 785 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_btagVeto"); |
| 786 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_btagVeto"); |
| 787 |
> |
|
| 788 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet"); |
| 789 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_AtLeastOneBJet"); |
| 790 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_AtLeastOneBJet"); |
| 791 |
> |
|
| 792 |
> |
make_data_comparison_plot("mll", cutOSOF&&cutSignal, 60, 20., 320.,-1., false, "m_{ll}", "mll_OF_sig"); |
| 793 |
> |
make_data_comparison_plot("mll", cutmass&&"met[4]>100&&met[4]<150&&pfJetGoodNum40<3", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_CR"); |
| 794 |
|
} |
| 795 |
|
|
| 796 |
< |
string mllCut(""); |
| 714 |
< |
if ( !PlottingSetup::openBox ) mllCut = "&&mll>120"; |
| 796 |
> |
|
| 797 |
|
if ( doKin ) { |
| 798 |
+ |
string mllCut(""); |
| 799 |
+ |
if ( !PlottingSetup::openBox ) mllCut = "&&mll>120"; |
| 800 |
+ |
|
| 801 |
|
// Plots in signal region |
| 802 |
|
make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true); |
| 803 |
|
make_kin_plot("mll","mll>20&&met[4]>100"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true); |
