| 33 |
|
|
| 34 |
|
|
| 35 |
|
|
| 36 |
< |
void find_one_peak_combination(TCut specialcut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="") |
| 36 |
> |
void find_one_peak_combination(TCut specialcut, bool SwitchOffNJetsCut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &MCSigmaError, float &DataSigma, float& DataSigmaError, stringstream &result, bool doPUreweighting = true, string saveas="") |
| 37 |
|
{ |
| 38 |
|
// Temporarily switch off PU reweighting, if asked |
| 39 |
|
TCut weightbackup=cutWeight; |
| 42 |
|
int nbins=100; |
| 43 |
|
if(PlottingSetup::DoBTag) nbins=25; |
| 44 |
|
|
| 45 |
+ |
TCut nJetsCut(cutnJets); |
| 46 |
+ |
if(SwitchOffNJetsCut) nJetsCut=specialcut; |
| 47 |
+ |
|
| 48 |
|
TCanvas *tempcan = new TCanvas("tempcan","Temporary canvas for peak finding preparations"); |
| 49 |
< |
TH1F *rawJZBeemmMC = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,mc, luminosity); |
| 50 |
< |
TH1F *rawJZBeemmData = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&specialcut,data, luminosity); |
| 51 |
< |
TH1F *rawJZBemMC = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,mc, luminosity); |
| 52 |
< |
TH1F *rawJZBemData = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&specialcut,data, luminosity); |
| 49 |
> |
TH1F *rawJZBeemmMC = allsamples.Draw("rawJZBeemmMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,mc, luminosity); |
| 50 |
> |
TH1F *rawJZBeemmData = allsamples.Draw("rawJZBeemmData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSSF&&nJetsCut&&specialcut,data, luminosity); |
| 51 |
> |
TH1F *rawJZBemMC = allsamples.Draw("rawJZBemMC",jzbvariablemc,nbins,-50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,mc, luminosity); |
| 52 |
> |
TH1F *rawJZBemData = allsamples.Draw("rawJZBemData",jzbvariabledata,nbins, -50,50, "JZB [GeV]", "events", cutmass&&cutOSOF&&nJetsCut&&specialcut,data, luminosity); |
| 53 |
|
TH1F *rawttbarjzbeemmMC; |
| 54 |
|
|
| 55 |
|
if(method==doKM) { |
| 56 |
|
//we only need this histo for the KM fitting... |
| 57 |
< |
rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&&specialcut,mc,luminosity,allsamples.FindSample("TTJet")); |
| 57 |
> |
rawttbarjzbeemmMC = allsamples.Draw("rawttbarjzbeemmMC",jzbvariablemc,nbins, -50,50, "JZB [GeV]", "events",cutmass&&cutOSSF&&nJetsCut&&specialcut,mc,luminosity,allsamples.FindSample("TTJet")); |
| 58 |
|
MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas); |
| 59 |
|
DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas); |
| 60 |
|
delete rawttbarjzbeemmMC; |
| 87 |
|
delete tempcan; |
| 88 |
|
} |
| 89 |
|
|
| 90 |
< |
void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb) |
| 90 |
> |
void find_peaks(float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, stringstream &result, bool doPUreweighting, stringstream &datajzb, stringstream &mcjzb, string sSpecialCut="", bool SwitchOffNJetsCut=false) |
| 91 |
|
{ |
| 92 |
|
switch_overunderflow(false); |
| 93 |
< |
|
| 93 |
> |
|
| 94 |
> |
TCut SpecialCut("mll>=0"); |
| 95 |
> |
if(sSpecialCut!="") SpecialCut=TCut(sSpecialCut.c_str()); |
| 96 |
> |
|
| 97 |
|
bool DoInvidualeemmPeaks=false; |
| 98 |
|
|
| 99 |
|
float mcpeak, datapeak; |
| 108 |
|
float mcSigma,mcSigmaError, dataSigma, dataSigmaError; |
| 109 |
|
|
| 110 |
|
dout << "Finding global peak : " << endl; |
| 111 |
< |
find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,""); |
| 111 |
> |
find_one_peak_combination(SpecialCut,SwitchOffNJetsCut,mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,""); |
| 112 |
|
|
| 113 |
|
if(DoInvidualeemmPeaks) { |
| 114 |
|
dout << "Finding peak for electrons : " << endl; |
| 115 |
< |
find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele"); |
| 115 |
> |
find_one_peak_combination(SpecialCut&&TCut("id1==0"),SwitchOffNJetsCut,mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele"); |
| 116 |
|
dout << "Finding peak for muons : " << endl; |
| 117 |
< |
find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu"); |
| 117 |
> |
find_one_peak_combination(SpecialCut&&TCut("id1==1"),SwitchOffNJetsCut,mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu"); |
| 118 |
|
|
| 119 |
|
datajzb << "(" << jzbvariabledata; |
| 120 |
|
mcjzb << "(" << jzbvariablemc; |
| 373 |
|
} |
| 374 |
|
|
| 375 |
|
float make_one_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, float ymax, bool logscale, |
| 376 |
< |
string xlabel, string filename, bool plotratio=true, bool loadlastminmax=false, float legendPosition=0.55) { |
| 376 |
> |
string xlabel, string filename, float legendPosition=0.55) { |
| 377 |
|
|
| 378 |
|
TCut ibasiccut=basiccut; |
| 379 |
|
bool draw_separation_lines=false; |
| 391 |
|
rcan->SetLogy(logscale); |
| 392 |
|
rcan->cd(); |
| 393 |
|
|
| 394 |
+ |
std::cout << "OF/SF comparison: variable = " << variable << ", cut = " << cutSF.GetTitle() << std::endl; |
| 395 |
|
TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cutSF,data,luminosity); |
| 396 |
|
TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cutOF,data,luminosity); |
| 397 |
< |
string signal("LM3"); |
| 398 |
< |
TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max); |
| 399 |
< |
int idx = signalsamples.FindSample(signal)[0]; |
| 400 |
< |
(signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF); |
| 401 |
< |
signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity); |
| 402 |
< |
signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor); |
| 403 |
< |
signalhisto->SetLineStyle(2); |
| 397 |
> |
// string signal("LM3"); |
| 398 |
> |
// TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max); |
| 399 |
> |
// int idx = signalsamples.FindSample(signal)[0]; |
| 400 |
> |
// (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF); |
| 401 |
> |
// signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity); |
| 402 |
> |
// signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor); |
| 403 |
> |
// signalhisto->SetLineStyle(2); |
| 404 |
|
datahistoSF->SetMarkerSize(DataMarkerSize); |
| 405 |
|
datahistoOF->SetLineColor(kRed); |
| 406 |
|
|
| 430 |
|
} else { |
| 431 |
|
mleg->AddEntry((TObject*)0, "", ""); |
| 432 |
|
} |
| 433 |
< |
mleg->AddEntry(signalhisto, "LM3", "L"); |
| 433 |
> |
//mleg->AddEntry(signalhisto, "LM3", "L"); |
| 434 |
|
|
| 435 |
|
datahistoSF->Draw("E1"); |
| 436 |
|
if (datahistoOF->Integral()>0) datahistoOF->Draw("HIST,SAMES"); |
| 437 |
< |
signalhisto->Draw("HIST,SAMES"); |
| 437 |
> |
//signalhisto->Draw("HIST,SAMES"); |
| 438 |
|
mleg->Draw(); |
| 439 |
|
DrawPrelim(); |
| 440 |
< |
if (datahistoOF->Integral()>0) save_with_ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" ); |
| 441 |
< |
else CompleteSave(ckin, "SFOF/" + filename); |
| 440 |
> |
if (datahistoOF->Integral()>0) { |
| 441 |
> |
save_with_ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" ); |
| 442 |
> |
} else { |
| 443 |
> |
CompleteSave(rcan, "SFOF/" + filename); |
| 444 |
> |
delete rcan; |
| 445 |
> |
} |
| 446 |
|
|
| 447 |
|
datahistoSF->Delete(); |
| 448 |
|
datahistoOF->Delete(); |
| 449 |
< |
signalhisto->Delete(); |
| 449 |
> |
//signalhisto->Delete(); |
| 450 |
|
delete mleg; |
| 451 |
|
delete ckin; |
| 452 |
|
|
| 493 |
|
data1->GetYaxis()->CenterTitle(); |
| 494 |
|
|
| 495 |
|
TLegend *mleg = make_legend("",legendPosition,0.7,false,legendPosition+0.2); |
| 496 |
< |
mleg->AddEntry(data1, "Old 5/fb", "PL"); |
| 497 |
< |
mleg->AddEntry(data2, "New x/fb", "L"); |
| 496 |
> |
mleg->AddEntry(data1, "New 3.8/fb", "PL"); |
| 497 |
> |
mleg->AddEntry(data2, "Old 5.1/fb", "L"); |
| 498 |
|
|
| 499 |
|
data1->Draw("E1"); |
| 500 |
|
data2->Draw("HIST,SAMES"); |
| 512 |
|
void make_OFSF_plots(string variable, string addcut, int nbins, float min, float max, bool logscale, |
| 513 |
|
string xlabel, string filename, float legendPosition=0.55) { |
| 514 |
|
|
| 515 |
< |
string mllcuts[] = { "mll>20","mll>20&&mll<70", "mll>75&&mll<105", "mll>120" }; |
| 516 |
< |
string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "75 < m_{ll} < 105 GeV", "m_{ll} > 120 GeV" }; |
| 515 |
> |
string mllcuts[] = { "mll>15","mll>15&&mll<70", "mll>75&&mll<105", "mll>120" }; |
| 516 |
> |
string mllcutname[] = { "m_{ll} > 15 GeV", "15 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 120 GeV" }; |
| 517 |
|
string plotname[] = {"_all","_low","_peak","_high"}; |
| 518 |
|
float ymax; |
| 519 |
|
|
| 523 |
|
for ( int i=start; i<4; ++i ) { |
| 524 |
|
if ( addcut != "" ) mllcuts[i] += "&&"+addcut; |
| 525 |
|
if ( i==start ) { |
| 526 |
< |
ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel, |
| 526 |
> |
ymax = make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, -1, logscale, xlabel, |
| 527 |
|
filename+plotname[i], legendPosition ); |
| 528 |
|
} else { |
| 529 |
< |
make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 529 |
> |
make_one_OFSF_plot(variable, mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 530 |
|
filename+plotname[i], legendPosition ); |
| 531 |
|
} |
| 532 |
|
make_one_OFSF_plot(variable, "id1==1&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 765 |
|
if ( doOFSF ) { |
| 766 |
|
make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll"); |
| 767 |
|
|
| 768 |
< |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets"); |
| 769 |
< |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto"); |
| 770 |
< |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30"); |
| 768 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets"); |
| 769 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto"); |
| 770 |
> |
make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30"); |
| 771 |
|
|
| 772 |
|
make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets"); |
| 773 |
+ |
make_OFSF_plots("pfJetGoodPtBtag[0]", "met[4]>100&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", "ptb1"); |
| 774 |
+ |
|
| 775 |
|
make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1"); |
| 776 |
|
make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2"); |
| 777 |
|
// make_OFSF_plots("pt1", "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1"); |
| 793 |
|
|
| 794 |
|
if ( doDataComp && !PlottingSetup::openBox ) { |
| 795 |
|
TCut mllCut(""); |
| 796 |
< |
if ( !PlottingSetup::openBox ) mllCut = "mll>120"; |
| 796 |
> |
float massmin = 15.; |
| 797 |
> |
float massmax = 315; |
| 798 |
> |
int massnbins = 60; |
| 799 |
> |
if ( !PlottingSetup::openBox ) { |
| 800 |
> |
mllCut = "mll>120"; |
| 801 |
> |
massmin = 120; |
| 802 |
> |
massmax = 360; |
| 803 |
> |
massnbins = 14; |
| 804 |
> |
} |
| 805 |
|
|
| 806 |
|
TCut cutSignal = cutmass&&cutnJets&&"met[4]>100"; |
| 807 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut, 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig"); |
| 808 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee"); |
| 788 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm"); |
| 789 |
< |
|
| 790 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto"); |
| 791 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_btagVeto"); |
| 792 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30==0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_btagVeto"); |
| 793 |
< |
|
| 794 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet"); |
| 795 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==0&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_ee_AtLeastOneBJet"); |
| 796 |
< |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"id1==1&&pfJetGoodNumBtag30>0", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_sig_mm_AtLeastOneBJet"); |
| 807 |
> |
make_data_comparison_plot("mll", cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive"); |
| 808 |
> |
make_data_comparison_plot("mll", cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_inclusive"); |
| 809 |
|
|
| 810 |
< |
make_data_comparison_plot("mll", cutOSOF&&cutSignal, 60, 20., 320.,-1., false, "m_{ll}", "mll_OF_sig"); |
| 811 |
< |
make_data_comparison_plot("mll", cutmass&&"met[4]>100&&met[4]<150&&pfJetGoodNum40<3", 60, 20., 320.,-1., false, "m_{ll}", "mll_SF_CR"); |
| 810 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig"); |
| 811 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto"); |
| 812 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet"); |
| 813 |
> |
|
| 814 |
> |
make_data_comparison_plot("mll", mllCut&&cutOSOF&&cutSignal, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig"); |
| 815 |
> |
make_data_comparison_plot("mll", cutmass&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_CR"); |
| 816 |
> |
|
| 817 |
> |
make_data_comparison_plot("pfJetGoodNum40", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_SF_sig"); |
| 818 |
> |
make_data_comparison_plot("pfJetGoodNum40", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_OF_sig"); |
| 819 |
> |
make_data_comparison_plot("pfJetGoodNumBtag30", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_SF_sig"); |
| 820 |
> |
make_data_comparison_plot("pfJetGoodNumBtag30", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_OF_sig"); |
| 821 |
|
} |
| 822 |
|
|
| 823 |
|
|
