| 33 |
|
|
| 34 |
|
|
| 35 |
|
|
| 36 |
< |
void find_one_peak_combination(TCut specialcut, float &MCPeak,float &MCPeakError, float &DataPeak, float &DataPeakError, float &MCSigma, float &DataSigma, 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")); |
| 58 |
< |
MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method,saveas); |
| 59 |
< |
DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method,saveas); |
| 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; |
| 61 |
|
} |
| 62 |
|
else { |
| 65 |
|
reducedMC->Add(rawJZBemMC,-1); |
| 66 |
|
reducedData->Add(rawJZBemData,-1); |
| 67 |
|
//this is Kostas' way of doing it - we subtract em to get rid of some of the ttbar contribution (in reality, of flavor-symmetric contribution) |
| 68 |
< |
MCPeak=find_peak(reducedMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method,saveas); |
| 69 |
< |
DataPeak=find_peak(reducedData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method,saveas); |
| 68 |
> |
MCPeak=find_peak(reducedMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,MCSigmaError,method,saveas); |
| 69 |
> |
DataPeak=find_peak(reducedData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,DataSigmaError,method,saveas); |
| 70 |
|
delete reducedMC; |
| 71 |
|
delete reducedData; |
| 72 |
|
} |
| 76 |
|
|
| 77 |
|
// MCPeak=find_peak(rawJZBeemmMC, rawttbarjzbeemmMC, -40, 40, mc, MCPeakError,MCSigma,method); |
| 78 |
|
// DataPeak=find_peak(rawJZBeemmData, rawJZBeemmData, -40, 40, data, DataPeakError,DataSigma,method); |
| 79 |
< |
dout << " We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- ?? (not impl.)" << endl; |
| 80 |
< |
result << " We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- ?? (not impl.)" << endl; |
| 81 |
< |
dout << " We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- ?? (not impl.)" << endl; |
| 82 |
< |
result << " We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- ?? (not impl.)" << endl; |
| 79 |
> |
dout << " We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- " << DataSigmaError << endl; |
| 80 |
> |
result << " We have found the peak in Data at " << DataPeak << " +/- " << DataPeakError << " with sigma=" << DataSigma << " +/- " << DataSigmaError << endl; |
| 81 |
> |
dout << " We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- " << MCSigmaError << endl; |
| 82 |
> |
result << " We have found the peak in MC at " << MCPeak << " +/- " << MCPeakError << " with sigma=" << MCSigma << " +/- " << MCSigmaError << endl; |
| 83 |
|
delete rawJZBeemmData; |
| 84 |
|
delete rawJZBeemmMC; |
| 85 |
|
delete rawJZBemData; |
| 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 |
+ |
bool overunderflowstatus=addoverunderflowbins; |
| 93 |
+ |
switch_overunderflow(false); |
| 94 |
+ |
|
| 95 |
+ |
TCut SpecialCut("mll>=0"); |
| 96 |
+ |
if(sSpecialCut!="") SpecialCut=TCut(sSpecialCut.c_str()); |
| 97 |
|
|
| 98 |
|
bool DoInvidualeemmPeaks=false; |
| 99 |
|
|
| 106 |
|
float datammpeak,dataeepeak; |
| 107 |
|
float datammpeakerr,dataeepeakerr; |
| 108 |
|
|
| 109 |
< |
float mcSigma, dataSigma; |
| 109 |
> |
float mcSigma,mcSigmaError, dataSigma, dataSigmaError; |
| 110 |
|
|
| 111 |
|
dout << "Finding global peak : " << endl; |
| 112 |
< |
find_one_peak_combination(TCut(""),mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma, dataSigma,result,doPUreweighting,""); |
| 112 |
> |
find_one_peak_combination(SpecialCut,SwitchOffNJetsCut,mcpeak,mcpeakerr, datapeak,datapeakerr,mcSigma,mcSigmaError, dataSigma,dataSigmaError,result,doPUreweighting,""); |
| 113 |
|
|
| 114 |
|
if(DoInvidualeemmPeaks) { |
| 115 |
|
dout << "Finding peak for electrons : " << endl; |
| 116 |
< |
find_one_peak_combination(TCut("id1==0"),mceepeak,mceepeakerr, dataeepeak,dataeepeakerr,mcSigma, dataSigma,result,doPUreweighting,"_ele"); |
| 116 |
> |
find_one_peak_combination(SpecialCut&&TCut("id1==0"),SwitchOffNJetsCut,mceepeak,mceepeakerr,dataeepeak,dataeepeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_ele"); |
| 117 |
|
dout << "Finding peak for muons : " << endl; |
| 118 |
< |
find_one_peak_combination(TCut("id1==1"),mcmmpeak,mcmmpeakerr, datammpeak,datammpeakerr,mcSigma, dataSigma,result,doPUreweighting,"_mu"); |
| 118 |
> |
find_one_peak_combination(SpecialCut&&TCut("id1==1"),SwitchOffNJetsCut,mcmmpeak,mcmmpeakerr,datammpeak,datammpeakerr,mcSigma,mcSigmaError,dataSigma,dataSigmaError,result,doPUreweighting,"_mu"); |
| 119 |
|
|
| 120 |
|
datajzb << "(" << jzbvariabledata; |
| 121 |
|
mcjzb << "(" << jzbvariablemc; |
| 141 |
|
else mcjzb << "+ (id1!=id2)*" << TMath::Abs(mcpeak) << " "; |
| 142 |
|
|
| 143 |
|
mcjzb << ")"; |
| 144 |
+ |
|
| 145 |
|
} else { |
| 146 |
|
datajzb << "(" << jzbvariabledata; |
| 147 |
|
mcjzb << "(" << jzbvariablemc; |
| 156 |
|
|
| 157 |
|
mcjzb << ")"; |
| 158 |
|
} |
| 150 |
– |
|
| 159 |
|
|
| 160 |
+ |
MCPeak=mcpeak; |
| 161 |
+ |
MCPeakError=mcpeakerr; |
| 162 |
+ |
DataPeak=datapeak; |
| 163 |
+ |
DataPeakError=datapeakerr; |
| 164 |
+ |
switch_overunderflow(overunderflowstatus); |
| 165 |
|
} |
| 166 |
|
|
| 167 |
< |
void make_special_obs_pred_mll_plot(string mcjzb, float jzbthreshold) { |
| 167 |
> |
void make_special_obs_pred_mll_plot(string datajzb, string mcjzb, float jzbthreshold, float binWidth = 5.0) { |
| 168 |
|
float min=70.0; |
| 169 |
|
float max=115.0; |
| 170 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
| 171 |
< |
min=10; |
| 172 |
< |
max=150; |
| 171 |
> |
min=20; |
| 172 |
> |
max=300; |
| 173 |
|
} |
| 174 |
< |
int nbins=int((max-min)/5); |
| 174 |
> |
int nbins=int((max-min)/binWidth); |
| 175 |
|
|
| 176 |
|
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 177 |
|
|
| 178 |
< |
/* stringstream largerzeroS; |
| 179 |
< |
largerzeroS << "(" << mcjzb << ">" << jzbthreshold << ")"; |
| 180 |
< |
TCut largerzero(largerzeroS.str().c_str()); |
| 178 |
> |
stringstream largerzeroS; |
| 179 |
> |
largerzeroS << "(" << datajzb << ">" << jzbthreshold << ")"; |
| 180 |
> |
TCut largerzeroD(largerzeroS.str().c_str()); |
| 181 |
|
|
| 182 |
|
stringstream smallerzeroS; |
| 183 |
< |
smallerzeroS << "(" << mcjzb << "<-" << jzbthreshold << ")"; |
| 184 |
< |
TCut smallerzero(smallerzeroS.str().c_str());*/ |
| 185 |
< |
|
| 186 |
< |
TH1F *RcorrJZBeemm = allsamples.Draw("RcorrJZBeemm", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data,luminosity); |
| 187 |
< |
THStack mcRcorrJZBeemm = allsamples.DrawStack("mcRcorrJZBeemm","mll",nbins,min,max, "m_{ll} [GeV}", "events", cutmass&&cutOSSF&&cutnJets,mc,luminosity); |
| 188 |
< |
TH1F *LcorrJZBeemm = allsamples.Draw("LcorrJZBeemm", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSSF&&cutnJets,data,luminosity); |
| 189 |
< |
TH1F *RcorrJZBem = allsamples.Draw("RcorrJZBem", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSOF&&cutnJets,data,luminosity); |
| 190 |
< |
TH1F *LcorrJZBem = allsamples.Draw("LcorrJZBem", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSOF&&cutnJets,data,luminosity); |
| 183 |
> |
smallerzeroS << "(" << datajzb << "<-" << jzbthreshold << ")"; |
| 184 |
> |
TCut smallerzeroD(smallerzeroS.str().c_str()); |
| 185 |
> |
|
| 186 |
> |
|
| 187 |
> |
stringstream largerzeroMS; |
| 188 |
> |
largerzeroMS << "(" << mcjzb << ">" << jzbthreshold << ")"; |
| 189 |
> |
TCut largerzeroM(largerzeroMS.str().c_str()); |
| 190 |
> |
|
| 191 |
> |
TH1F *RcorrJZBeemm = allsamples.Draw("RcorrJZBeemm", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&largerzeroD,data,luminosity); |
| 192 |
> |
THStack mcRcorrJZBeemm = allsamples.DrawStack("mcRcorrJZBeemm","mll",nbins,min,max, "m_{ll} [GeV}", "events", cutmass&&cutOSSF&&cutnJets&&largerzeroM,mc,luminosity); |
| 193 |
> |
TH1F *LcorrJZBeemm = allsamples.Draw("LcorrJZBeemm", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&smallerzeroD,data,luminosity); |
| 194 |
> |
TH1F *RcorrJZBem = allsamples.Draw("RcorrJZBem", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&largerzeroD,data,luminosity); |
| 195 |
> |
TH1F *LcorrJZBem = allsamples.Draw("LcorrJZBem", "mll",nbins,min,max, "m_{ll} [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&smallerzeroD,data,luminosity); |
| 196 |
|
|
| 197 |
|
TH1F *RcorrJZBSBem; |
| 198 |
|
TH1F *LcorrJZBSBem; |
| 201 |
|
|
| 202 |
|
// TH1F *RcorrJZBeemmNoS; |
| 203 |
|
|
| 204 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 205 |
< |
RcorrJZBSBem = allsamples.Draw("RcorrJZBSBem", "mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,data, luminosity); |
| 206 |
< |
LcorrJZBSBem = allsamples.Draw("LcorrJZBSBem", "mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,data, luminosity); |
| 207 |
< |
RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm","mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,data, luminosity); |
| 208 |
< |
LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm","mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,data, luminosity); |
| 209 |
< |
} |
| 210 |
< |
|
| 211 |
< |
TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred"); |
| 212 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 213 |
< |
Bpred->Add(RcorrJZBem,1.0/3.); |
| 214 |
< |
Bpred->Add(LcorrJZBem,-1.0/3.); |
| 215 |
< |
Bpred->Add(RcorrJZBSBem,1.0/3.); |
| 216 |
< |
Bpred->Add(LcorrJZBSBem,-1.0/3.); |
| 217 |
< |
Bpred->Add(RcorrJZBSBeemm,1.0/3.); |
| 218 |
< |
Bpred->Add(LcorrJZBSBeemm,-1.0/3.); |
| 219 |
< |
} else { |
| 220 |
< |
Bpred->Add(RcorrJZBem,1.0); |
| 221 |
< |
Bpred->Add(LcorrJZBem,-1.0); |
| 222 |
< |
} |
| 223 |
< |
|
| 204 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 205 |
> |
RcorrJZBSBem = allsamples.Draw("RcorrJZBSBem", "mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets&&largerzeroD,data, luminosity); |
| 206 |
> |
LcorrJZBSBem = allsamples.Draw("LcorrJZBSBem", "mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets&&smallerzeroD,data, luminosity); |
| 207 |
> |
RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm","mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets&&largerzeroD,data, luminosity); |
| 208 |
> |
LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm","mll",nbins,min,max, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets&&smallerzeroD,data, luminosity); |
| 209 |
> |
} |
| 210 |
> |
|
| 211 |
> |
// Separate predictions |
| 212 |
> |
TH1F* SFN = (TH1F*)LcorrJZBeemm->Clone("SFN"); |
| 213 |
> |
TH1F* OFP = (TH1F*)RcorrJZBem->Clone("OFP"); |
| 214 |
> |
TH1F* OFN = (TH1F*)LcorrJZBem->Clone("OFN"); |
| 215 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 216 |
> |
OFP->Scale(1.0/3.0); |
| 217 |
> |
OFP->Add(RcorrJZBSBem,1.0/3.); |
| 218 |
> |
OFP->Add(RcorrJZBSBeemm,1.0/3.); |
| 219 |
> |
OFN->Scale(1.0/3.0); |
| 220 |
> |
OFN->Add(LcorrJZBSBem,1.0/3.); |
| 221 |
> |
OFN->Add(LcorrJZBSBeemm,1.0/3.); |
| 222 |
> |
} |
| 223 |
> |
|
| 224 |
> |
TH1F* Bpred = (TH1F*)SFN->Clone("Bpred"); |
| 225 |
> |
Bpred->Add(OFP); |
| 226 |
> |
Bpred->Add(OFN,-1); |
| 227 |
|
Bpred->SetLineColor(kRed); |
| 228 |
|
|
| 229 |
+ |
RcorrJZBeemm->SetTitleOffset(1.3,"y"); |
| 230 |
|
RcorrJZBeemm->Draw(); |
| 231 |
< |
mcRcorrJZBeemm.Draw("same"); |
| 231 |
> |
mcRcorrJZBeemm.Draw("histo,same"); |
| 232 |
|
Bpred->Draw("histo,same"); |
| 233 |
|
RcorrJZBeemm->Draw("same"); |
| 234 |
|
|
| 241 |
|
stringstream saveas; |
| 242 |
|
saveas << "kin/Mll_After_Cut/Cut_At" << jzbthreshold; |
| 243 |
|
CompleteSave(ckin,saveas.str()); |
| 244 |
< |
|
| 244 |
> |
|
| 245 |
> |
// Draw all predictions overlayed |
| 246 |
> |
unsigned int w = gStyle->GetHistLineWidth()+1; // Make line a bit wider, since we dash it |
| 247 |
> |
SFN->SetLineColor(kGreen+2); |
| 248 |
> |
SFN->SetLineStyle(2); |
| 249 |
> |
SFN->SetLineWidth(w); |
| 250 |
> |
OFP->SetLineColor(kBlue+2); |
| 251 |
> |
OFP->SetLineStyle(2); |
| 252 |
> |
OFP->SetLineWidth(w); |
| 253 |
> |
OFN->SetLineColor(kMagenta+2); |
| 254 |
> |
OFN->SetLineStyle(3); |
| 255 |
> |
OFN->SetLineWidth(w); |
| 256 |
> |
|
| 257 |
> |
RcorrJZBeemm->Draw(); |
| 258 |
> |
SFN->Draw("histo,same"); |
| 259 |
> |
OFP->Draw("histo,same"); |
| 260 |
> |
OFN->Draw("histo,same"); |
| 261 |
> |
Bpred->Draw("histo,same"); |
| 262 |
> |
RcorrJZBeemm->Draw("same"); |
| 263 |
> |
|
| 264 |
> |
TLegend *leg2 = make_legend("",0.52,0.7); |
| 265 |
> |
leg2->AddEntry(RcorrJZBeemm,"observed (data)","lp"); |
| 266 |
> |
leg2->AddEntry(Bpred,"predicted (data)","l"); |
| 267 |
> |
leg2->AddEntry(SFN, " SF JZB<0","l"); |
| 268 |
> |
leg2->AddEntry(OFN, " OF JZB<0","l"); |
| 269 |
> |
leg2->AddEntry(OFP, " OF JZB>0","l"); |
| 270 |
> |
leg2->Draw("same"); |
| 271 |
> |
|
| 272 |
> |
saveas.str(""); |
| 273 |
> |
saveas << "kin/Mll_After_Cut/Cut_At" << jzbthreshold << "_nomc"; |
| 274 |
> |
CompleteSave(ckin,saveas.str()); |
| 275 |
> |
|
| 276 |
|
delete RcorrJZBeemm; |
| 277 |
|
delete LcorrJZBeemm; |
| 278 |
|
delete RcorrJZBem; |
| 279 |
|
delete LcorrJZBem; |
| 280 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 280 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 281 |
|
delete RcorrJZBSBeemm; |
| 282 |
|
delete LcorrJZBSBeemm; |
| 283 |
|
delete RcorrJZBSBem; |
| 305 |
|
datahistoOSSF->SetMarkerSize(DataMarkerSize); |
| 306 |
|
datahistoOSSF->Draw(); |
| 307 |
|
|
| 308 |
< |
mcstackOSSF.Draw("same"); |
| 308 |
> |
mcstackOSSF.Draw("histo,same"); |
| 309 |
|
datahistoOSSF->Draw("same"); |
| 310 |
|
|
| 311 |
|
datahistoOSOF->SetMarkerColor(TColor::GetColor("#FE642E")); |
| 371 |
|
oneline->Draw("same"); |
| 372 |
|
} |
| 373 |
|
|
| 374 |
+ |
float make_one_OFSF_plot(string variable, string addcut, string legendTitle, int nbins, float min, float max, float ymax, bool logscale, |
| 375 |
+ |
string xlabel, string filename, float legendPosition=0.55) { |
| 376 |
+ |
|
| 377 |
+ |
TCut ibasiccut=basiccut; |
| 378 |
+ |
bool draw_separation_lines=false; |
| 379 |
+ |
|
| 380 |
+ |
if(addcut != "") ibasiccut = ibasiccut && addcut.c_str(); |
| 381 |
+ |
|
| 382 |
+ |
TCut cutSF; |
| 383 |
+ |
TCut cutOF; |
| 384 |
+ |
|
| 385 |
+ |
cutOF = cutOSOF&&cutnJets&&ibasiccut; |
| 386 |
+ |
cutSF = cutOSSF&&cutnJets&&ibasiccut; |
| 387 |
+ |
|
| 388 |
+ |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 389 |
+ |
TPad* rcan = new TPad("rcan","rcan",0,0,1,1); |
| 390 |
+ |
rcan->SetLogy(logscale); |
| 391 |
+ |
rcan->cd(); |
| 392 |
+ |
|
| 393 |
+ |
std::cout << "OF/SF comparison: variable = " << variable << ", cut = " << cutSF.GetTitle() << std::endl; |
| 394 |
+ |
TH1F *datahistoSF = allsamples.Draw("datahistoSF",variable,nbins,min,max, xlabel, "events",cutSF,data,luminosity); |
| 395 |
+ |
TH1F *datahistoOF = allsamples.Draw("datahistoOF",variable,nbins,min,max, xlabel, "events",cutOF,data,luminosity); |
| 396 |
+ |
// string signal("LM3"); |
| 397 |
+ |
// TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max); |
| 398 |
+ |
// int idx = signalsamples.FindSample(signal)[0]; |
| 399 |
+ |
// (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cutSF); |
| 400 |
+ |
// signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity); |
| 401 |
+ |
// signalhisto->SetLineColor((signalsamples.collection)[idx].samplecolor); |
| 402 |
+ |
// signalhisto->SetLineStyle(2); |
| 403 |
+ |
datahistoSF->SetMarkerSize(DataMarkerSize); |
| 404 |
+ |
datahistoOF->SetLineColor(kRed); |
| 405 |
+ |
|
| 406 |
+ |
if ( !logscale ) { |
| 407 |
+ |
datahistoSF->SetMinimum(0); // Defaults |
| 408 |
+ |
} else { |
| 409 |
+ |
datahistoSF->SetMinimum(0.5); |
| 410 |
+ |
} |
| 411 |
+ |
if (ymax<0) { |
| 412 |
+ |
if ( logscale ) datahistoSF->SetMaximum(5.3*datahistoSF->GetMaximum()); |
| 413 |
+ |
else datahistoSF->SetMaximum(1.8*datahistoSF->GetMaximum()); |
| 414 |
+ |
} else { |
| 415 |
+ |
datahistoSF->SetMaximum(ymax); |
| 416 |
+ |
} |
| 417 |
+ |
|
| 418 |
+ |
float ymaxSet = datahistoSF->GetMaximum(); |
| 419 |
+ |
|
| 420 |
+ |
datahistoSF->GetXaxis()->SetTitle(xlabel.c_str()); |
| 421 |
+ |
datahistoSF->GetYaxis()->SetTitle("Events"); |
| 422 |
+ |
datahistoSF->GetXaxis()->CenterTitle(); |
| 423 |
+ |
datahistoSF->GetYaxis()->CenterTitle(); |
| 424 |
+ |
|
| 425 |
+ |
TLegend *mleg = make_legend(legendTitle.c_str(),legendPosition,0.7,false,legendPosition+0.2); |
| 426 |
+ |
mleg->AddEntry(datahistoSF, "Same-flavor", "PL"); |
| 427 |
+ |
if (datahistoOF->Integral()>0) { |
| 428 |
+ |
mleg->AddEntry(datahistoOF, "Opposite-flavor", "L"); |
| 429 |
+ |
} else { |
| 430 |
+ |
mleg->AddEntry((TObject*)0, "", ""); |
| 431 |
+ |
} |
| 432 |
+ |
//mleg->AddEntry(signalhisto, "LM3", "L"); |
| 433 |
+ |
|
| 434 |
+ |
datahistoSF->Draw("E1"); |
| 435 |
+ |
if (datahistoOF->Integral()>0) datahistoOF->Draw("HIST,SAMES"); |
| 436 |
+ |
//signalhisto->Draw("HIST,SAMES"); |
| 437 |
+ |
mleg->Draw(); |
| 438 |
+ |
DrawPrelim(); |
| 439 |
+ |
if (datahistoOF->Integral()>0) { |
| 440 |
+ |
save_with_ratio( datahistoSF, datahistoOF, rcan, "SFOF/" + filename, false, false, "SF/OF" ); |
| 441 |
+ |
} else { |
| 442 |
+ |
CompleteSave(rcan, "SFOF/" + filename); |
| 443 |
+ |
delete rcan; |
| 444 |
+ |
} |
| 445 |
+ |
|
| 446 |
+ |
datahistoSF->Delete(); |
| 447 |
+ |
datahistoOF->Delete(); |
| 448 |
+ |
//signalhisto->Delete(); |
| 449 |
+ |
delete mleg; |
| 450 |
+ |
delete ckin; |
| 451 |
+ |
|
| 452 |
+ |
return ymaxSet; |
| 453 |
+ |
|
| 454 |
+ |
} |
| 455 |
+ |
|
| 456 |
+ |
// Compare data to data |
| 457 |
+ |
float make_data_comparison_plot(string variable, TCut cut, int nbins, float min, float max, float ymax, bool logscale, |
| 458 |
+ |
string xlabel, string filename, float legendPosition=0.55) { |
| 459 |
+ |
|
| 460 |
+ |
TCut ibasiccut=basiccut&&cut; |
| 461 |
+ |
|
| 462 |
+ |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 463 |
+ |
TPad* rcan = new TPad("rcan","rcan",0,0,1,1); |
| 464 |
+ |
rcan->SetLogy(logscale); |
| 465 |
+ |
rcan->cd(); |
| 466 |
+ |
|
| 467 |
+ |
std::cout << "Data comparison: variable = " << variable << ", cut = " << ibasiccut.GetTitle() << std::endl; |
| 468 |
+ |
|
| 469 |
+ |
TH1F *data1 = allsamples.Draw("data1",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity); |
| 470 |
+ |
TH1F *data2 = comparesamples.Draw("data2",variable,nbins,min,max, xlabel, "events",ibasiccut,data,luminosity); |
| 471 |
+ |
data1->Scale(5.0/9.2); |
| 472 |
+ |
|
| 473 |
+ |
data1->SetMarkerSize(DataMarkerSize); |
| 474 |
+ |
data2->SetLineColor(kRed); |
| 475 |
+ |
|
| 476 |
+ |
if ( !logscale ) { |
| 477 |
+ |
data1->SetMinimum(0); // Defaults |
| 478 |
+ |
} else { |
| 479 |
+ |
data1->SetMinimum(0.5); |
| 480 |
+ |
} |
| 481 |
+ |
if (ymax<0) { |
| 482 |
+ |
if ( logscale ) data1->SetMaximum(5.3*data1->GetMaximum()); |
| 483 |
+ |
else data1->SetMaximum(1.5*data1->GetMaximum()); |
| 484 |
+ |
} else { |
| 485 |
+ |
data1->SetMaximum(ymax); |
| 486 |
+ |
} |
| 487 |
+ |
|
| 488 |
+ |
float ymaxSet = data1->GetMaximum(); |
| 489 |
+ |
|
| 490 |
+ |
data1->GetXaxis()->SetTitle(xlabel.c_str()); |
| 491 |
+ |
data1->GetYaxis()->SetTitle("Events"); |
| 492 |
+ |
data1->GetXaxis()->CenterTitle(); |
| 493 |
+ |
data1->GetYaxis()->CenterTitle(); |
| 494 |
+ |
|
| 495 |
+ |
TLegend *mleg = make_legend("",legendPosition,0.7,false,legendPosition+0.2); |
| 496 |
+ |
mleg->AddEntry(data1, "2012 data (scaled)", "PL"); |
| 497 |
+ |
mleg->AddEntry(data2, "2011 data", "L"); |
| 498 |
+ |
|
| 499 |
+ |
data1->Draw("E1"); |
| 500 |
+ |
data2->Draw("HIST,SAMES"); |
| 501 |
+ |
mleg->Draw(); |
| 502 |
+ |
DrawPrelim(); |
| 503 |
+ |
save_with_ratio( data1, data2, rcan, "compareData/" + filename, false, false, "new/old" ); |
| 504 |
+ |
|
| 505 |
+ |
data1->Delete(); |
| 506 |
+ |
data2->Delete(); |
| 507 |
+ |
delete mleg; |
| 508 |
+ |
delete ckin; |
| 509 |
+ |
|
| 510 |
+ |
} |
| 511 |
+ |
|
| 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>15&&mll<70", "mll>75&&mll<105", "mll>120" }; |
| 516 |
+ |
string mllcutname[] = { "m_{ll} > 20 GeV", "20 < m_{ll} < 70 GeV", "70 < m_{ll} < 110 GeV", "m_{ll} > 120 GeV" }; |
| 517 |
+ |
string plotname[] = {"_all","_low","_peak","_high"}; |
| 518 |
+ |
float ymax; |
| 519 |
+ |
|
| 520 |
+ |
int start = 0; |
| 521 |
+ |
if ( !PlottingSetup::openBox ) start = 3; |
| 522 |
+ |
|
| 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, |
| 527 |
+ |
filename+plotname[i], legendPosition ); |
| 528 |
+ |
} else { |
| 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, |
| 533 |
+ |
filename+plotname[i]+"_mm", legendPosition ); |
| 534 |
+ |
make_one_OFSF_plot(variable, "id1==0&&id1==id2&&"+mllcuts[i], mllcutname[i], nbins, min, max, ymax, logscale, xlabel, |
| 535 |
+ |
filename+plotname[i]+"_ee", legendPosition ); |
| 536 |
+ |
} |
| 537 |
+ |
|
| 538 |
+ |
} |
| 539 |
+ |
|
| 540 |
+ |
TCut ReplaceInCut(TCut cut, string replacethis, string withthis) { |
| 541 |
+ |
string scut=(string)(const char*)cut; |
| 542 |
+ |
string acut=ReplaceAll(scut,replacethis,withthis); |
| 543 |
+ |
return TCut(acut.c_str()); |
| 544 |
+ |
} |
| 545 |
+ |
|
| 546 |
+ |
void DoMCSystPlot(string datavariable, string mcvariable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) { |
| 547 |
+ |
TCut cut=addcut&&basiccut; |
| 548 |
+ |
|
| 549 |
+ |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 550 |
+ |
ckin->SetLogy(dolog); |
| 551 |
+ |
cout << "\r Creating " << filename << " : data (1/6)" << std::flush; |
| 552 |
+ |
TH1F *datahisto = allsamples.Draw("datahisto",datavariable,nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 553 |
+ |
datahisto->SetMarkerSize(DataMarkerSize); |
| 554 |
+ |
if ( !dolog ) datahisto->SetMinimum(0); // Defaults |
| 555 |
+ |
else datahisto->SetMinimum(0.5); |
| 556 |
+ |
if (dolog) datahisto->SetMaximum(5.3*datahisto->GetMaximum()); |
| 557 |
+ |
else datahisto->SetMaximum(1.3*datahisto->GetMaximum()); |
| 558 |
+ |
cout << "\r Creating " << filename << " : MC central (2/6)" << std::flush; |
| 559 |
+ |
THStack mcstack = allsamples.DrawStack("mcstack", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 560 |
+ |
TH1F *MCcentral = CollapseStack(mcstack); |
| 561 |
+ |
|
| 562 |
+ |
TCut bkpcut = cut; |
| 563 |
+ |
cut = ReplaceInCut(cut,"pfJetGoodNum40","pfJetGoodNum40p1sigma"); |
| 564 |
+ |
cout << "\r Creating " << filename << " : MC JES up (3/6)" << std::flush; |
| 565 |
+ |
TH1F *MCJESup = allsamples.Draw("MCJESup", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 566 |
+ |
|
| 567 |
+ |
cut = ReplaceInCut(cut,"pfJetGoodNum40p1sigma","pfJetGoodNum40n1sigma"); |
| 568 |
+ |
cout << "\r Creating " << filename << " : MC JES down (4/6)" << std::flush; |
| 569 |
+ |
TH1F *MCJESdn = allsamples.Draw("MCJESdn", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 570 |
+ |
|
| 571 |
+ |
cut=bkpcut; |
| 572 |
+ |
TCut bkpweight = cutWeight; |
| 573 |
+ |
cutWeight= ReplaceInCut(cutWeight,"weight","PUweightUP"); |
| 574 |
+ |
cout << "\r Creating " << filename << " : MC PU up (5/6)" << std::flush; |
| 575 |
+ |
TH1F *MCPUup = allsamples.Draw("MCPUup", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 576 |
+ |
|
| 577 |
+ |
cutWeight= ReplaceInCut(cutWeight,"PUweightUP","PUweightDOWN"); |
| 578 |
+ |
cout << "\r Creating " << filename << " : MC PU down (6/6)" << std::flush; |
| 579 |
+ |
TH1F *MCPUdn = allsamples.Draw("MCPUdn", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 580 |
+ |
cutWeight=bkpweight; |
| 581 |
+ |
cout << "\r Creating " << filename << " : Processing systematics ..." << std::flush; |
| 582 |
+ |
TH1F *Systematic = (TH1F*)MCPUdn->Clone("Systematic"); |
| 583 |
+ |
for(int i=1;i<=Systematic->GetNbinsX();i++) { |
| 584 |
+ |
float jesdn = abs(MCcentral->GetBinContent(i)-MCJESdn->GetBinContent(i)); |
| 585 |
+ |
float jesup = abs(MCcentral->GetBinContent(i)-MCJESup->GetBinContent(i)); |
| 586 |
+ |
float jes = jesdn>jesup?jesdn:jesup; |
| 587 |
+ |
|
| 588 |
+ |
float PUup = abs(MCcentral->GetBinContent(i)-MCPUup->GetBinContent(i)); |
| 589 |
+ |
float PUdn = abs(MCcentral->GetBinContent(i)-MCPUdn->GetBinContent(i)); |
| 590 |
+ |
float pu = PUdn>PUup?PUdn:PUup; |
| 591 |
+ |
|
| 592 |
+ |
float sys=sqrt(jes*jes+pu*pu); |
| 593 |
+ |
sys/=MCcentral->GetBinContent(i); |
| 594 |
+ |
if(MCcentral->GetBinContent(i)==0) sys=0; |
| 595 |
+ |
Systematic->SetBinContent(i,sys); |
| 596 |
+ |
} |
| 597 |
+ |
|
| 598 |
+ |
datahisto->Draw("e1"); |
| 599 |
+ |
ckin->Update(); |
| 600 |
+ |
mcstack.Draw("histo,same"); |
| 601 |
+ |
|
| 602 |
+ |
datahisto->Draw("same,e1"); |
| 603 |
+ |
TLegend *kinleg = allsamples.allbglegend(); |
| 604 |
+ |
kinleg->Draw(); |
| 605 |
+ |
|
| 606 |
+ |
TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1); |
| 607 |
+ |
kinpad->SetLogy(dolog); |
| 608 |
+ |
kinpad->cd(); |
| 609 |
+ |
datahisto->Draw("e1"); |
| 610 |
+ |
mcstack.Draw("histo,same"); |
| 611 |
+ |
datahisto->Draw("same,e1"); |
| 612 |
+ |
datahisto->Draw("same,axis"); |
| 613 |
+ |
|
| 614 |
+ |
kinleg->Draw("same"); |
| 615 |
+ |
DrawPrelim(); |
| 616 |
+ |
string saveas="kin/"+filename; |
| 617 |
+ |
|
| 618 |
+ |
cout << "Passing filename to be saved : " << filename << " as " << saveas << endl; |
| 619 |
+ |
save_with_ratio_and_sys_band( datahisto, CollapseStack(mcstack), kinpad->cd(), saveas, false, false, "data/mc",Systematic ); |
| 620 |
+ |
|
| 621 |
+ |
ckin->cd(); |
| 622 |
+ |
MCcentral->SetFillColor(kWhite); |
| 623 |
+ |
MCcentral->SetLineColor(kBlack); |
| 624 |
+ |
MCcentral->SetLineWidth(2); |
| 625 |
+ |
MCPUdn->SetLineColor(kRed); |
| 626 |
+ |
MCPUup->SetLineColor(kYellow); |
| 627 |
+ |
MCJESdn->SetLineColor(kBlue); |
| 628 |
+ |
MCJESup->SetLineColor(kGreen); |
| 629 |
+ |
|
| 630 |
+ |
TLegend *legsmall = new TLegend(0.0,0.85,1.0,1.00); |
| 631 |
+ |
legsmall->SetNColumns(5); |
| 632 |
+ |
legsmall->AddEntry(MCcentral,"central","L"); |
| 633 |
+ |
legsmall->AddEntry(MCPUdn,"PU down","L"); |
| 634 |
+ |
legsmall->AddEntry(MCPUup,"PU up","L"); |
| 635 |
+ |
legsmall->AddEntry(MCJESdn,"JES down","L"); |
| 636 |
+ |
legsmall->AddEntry(MCJESup,"JES up","L"); |
| 637 |
+ |
legsmall->SetFillColor(kWhite); |
| 638 |
+ |
legsmall->SetBorderSize(0); |
| 639 |
+ |
|
| 640 |
+ |
datahisto->Draw("e1"); |
| 641 |
+ |
MCcentral->Draw("histo,same"); |
| 642 |
+ |
MCPUdn->Draw("histo,same"); |
| 643 |
+ |
MCPUup->Draw("histo,same"); |
| 644 |
+ |
MCJESdn->Draw("histo,same"); |
| 645 |
+ |
MCJESup->Draw("histo,same"); |
| 646 |
+ |
datahisto->Draw("e1,same"); |
| 647 |
+ |
legsmall->Draw(); |
| 648 |
+ |
CompleteSave(ckin,saveas+"___AllLines"); |
| 649 |
+ |
|
| 650 |
+ |
delete datahisto; |
| 651 |
+ |
delete MCcentral; |
| 652 |
+ |
delete MCPUdn; |
| 653 |
+ |
delete MCPUup; |
| 654 |
+ |
delete MCJESdn; |
| 655 |
+ |
delete MCJESup; |
| 656 |
+ |
delete ckin; |
| 657 |
+ |
} |
| 658 |
+ |
|
| 659 |
+ |
void DoMCSystPlot(string variable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) { |
| 660 |
+ |
DoMCSystPlot(variable, variable, nbins, min, max, addcut, dolog, xlabel, filename); |
| 661 |
+ |
} |
| 662 |
+ |
|
| 663 |
|
float lastrange_min=0; |
| 664 |
|
float lastrange_max=0; |
| 665 |
|
|
| 668 |
|
// TCut basiccut("(pfJetGoodNum>=2&&pfJetGoodID[0])&&(pfJetGoodNum>=2&&pfJetGoodID[1])&&(passed_triggers||!is_data)"); |
| 669 |
|
TCut ibasiccut=basiccut; |
| 670 |
|
bool draw_separation_lines=false; |
| 671 |
+ |
bool drawsignal = false; |
| 672 |
|
|
| 673 |
|
if(isPF) ibasiccut=basiccut&&"pfjzb[0]>-998"; |
| 674 |
|
|
| 690 |
|
//Step 2: Refine the cut |
| 691 |
|
TCut cut; |
| 692 |
|
cut=cutmass&&cutOSSF&&cutnJets&&ibasiccut; |
| 693 |
< |
if(filename=="nJets") cut=cutmass&&cutOSSF&&ibasiccut; |
| 694 |
< |
if(filename=="nJets_osof") cut=cutmass&&cutOSOF&&ibasiccut; |
| 693 |
> |
if(filename=="nJets" || filename=="nJets_inclusive" || filename=="nJets_met100" || filename=="nJets_inclusive_met100") cut=cutmass&&cutOSSF&&ibasiccut; |
| 694 |
> |
if(filename=="nJets_osof" || filename=="nJets_osof_inclusive" || filename=="nJets_osof_met100" || filename=="nJets_osof_inclusive_met100") cut=cutmass&&cutOSOF&&ibasiccut; |
| 695 |
|
if(filename=="nJets_nocuts_except_mll_ossf") cut=cutmass&&cutOSSF; |
| 696 |
< |
if(filename=="mll") { |
| 696 |
> |
if(filename=="mll"||filename=="mll_met100") { |
| 697 |
|
cut=cutOSSF&&cutnJets&&ibasiccut; |
| 698 |
|
draw_separation_lines=true; |
| 699 |
|
} |
| 700 |
< |
if(filename=="mll_ee") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0"; |
| 701 |
< |
if(filename=="mll_osof") { |
| 700 |
> |
if(filename=="mll_ee"||filename=="mll_ee_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==0"; |
| 701 |
> |
if(filename=="mll_mm"||filename=="mll_mm_met100") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1"; |
| 702 |
> |
if(filename=="mll_osof"||filename=="mll_osof_met100") { |
| 703 |
|
cut=cutOSOF&&cutnJets&&ibasiccut; |
| 704 |
|
draw_separation_lines=true; |
| 705 |
|
} |
| 706 |
< |
if(filename=="mll_mm") cut=cutOSSF&&cutnJets&&ibasiccut&&"id1==1"; |
| 707 |
< |
if(Contains(filename,"mll_aboveJZB100")) cut=cutOSSF&&cutnJets&&ibasiccut; |
| 706 |
> |
if(Contains(filename,"aboveJZB")) cut=cutOSSF&&cutnJets&&ibasiccut; |
| 707 |
> |
if(Contains(filename,"mll_ee_above")) cut=cut&&"id1==0"; |
| 708 |
> |
if(Contains(filename,"mll_mm_above")) cut=cut&&"id1==1"; |
| 709 |
|
if(Contains(filename,"mll_osof_aboveJZB")) cut=cutOSOF&&cutnJets&&ibasiccut; |
| 710 |
|
if(filename=="mll_inclusive"||filename=="mll_inclusive_highrange") cut=cutOSSF; |
| 711 |
|
if(filename=="mll_inclusive_osof") cut=cutOSOF; |
| 713 |
|
if(filename=="mll_inclusive_mm") cut=cutOSSF&&"id1==1"; |
| 714 |
|
if(filename=="pfJetGoodEta_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets; |
| 715 |
|
if(filename=="pfJetGoodPt_0") cut=cutOSSF&&cutmass&&ibasiccut&&cutnJets; |
| 716 |
+ |
if(filename=="numVtx") cut=cutmass&&ibasiccut; |
| 717 |
|
|
| 718 |
|
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 719 |
|
ckin->SetLogy(logscale); |
| 720 |
|
TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 721 |
|
datahisto->SetMarkerSize(DataMarkerSize); |
| 722 |
< |
THStack mcstack = allsamples.DrawStack("mcstack",variable,nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 722 |
> |
if ( !logscale ) datahisto->SetMinimum(0); // Defaults |
| 723 |
> |
else datahisto->SetMinimum(0.5); |
| 724 |
> |
// Custom max. |
| 725 |
> |
if(variable=="TMath::Abs(pfJetDphiMet[0])") datahisto->SetMaximum(1.5*datahisto->GetMaximum()); |
| 726 |
|
if(variable=="pfJetGoodPt[0]") datahisto->SetMaximum(10*datahisto->GetMaximum()); |
| 727 |
|
if(variable=="pt") datahisto->SetMaximum(10*datahisto->GetMaximum()); |
| 728 |
|
if(filename=="mll_inclusive"||filename=="mll_inclusive_mm"||filename=="mll_inclusive_ee") datahisto->SetMinimum(1); |
| 729 |
|
if(filename=="mll_osof") datahisto->SetMaximum(10*datahisto->GetMaximum()); |
| 730 |
|
if(filename=="mll_osof") datahisto->SetMinimum(9); |
| 731 |
< |
datahisto->SetMaximum(5.3*datahisto->GetMaximum()); |
| 732 |
< |
datahisto->Draw("e1"); |
| 733 |
< |
ckin->Update(); |
| 731 |
> |
if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum()); |
| 732 |
> |
else datahisto->SetMaximum(1.3*datahisto->GetMaximum()); |
| 733 |
> |
|
| 734 |
> |
cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl; |
| 735 |
|
if(loadlastminmax) { |
| 736 |
|
datahisto->SetMinimum(lastrange_min); |
| 737 |
|
datahisto->SetMaximum(lastrange_max); |
| 741 |
|
} |
| 742 |
|
} |
| 743 |
|
|
| 744 |
< |
mcstack.Draw("same"); |
| 744 |
> |
// Draw signal by hand (for some reason I don't manage to use the sample class: it adds it to the stack!) |
| 745 |
> |
string signal("LM3"); |
| 746 |
> |
TH1F* signalhisto = new TH1F("signalhisto",signal.c_str(),nbins,min,max); |
| 747 |
> |
int idx = signalsamples.FindSample(signal)[0]; |
| 748 |
> |
if(drawsignal) (signalsamples.collection)[idx].events->Project("signalhisto",variable.c_str(),cut); |
| 749 |
> |
if(drawsignal) signalhisto->Scale((signalsamples.collection)[idx].weight*luminosity); |
| 750 |
> |
if(drawsignal) signalhisto->SetLineColor(kOrange); |
| 751 |
> |
|
| 752 |
> |
THStack mcstack = allsamples.DrawStack("mcstack", variable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 753 |
> |
datahisto->Draw("e1"); |
| 754 |
> |
ckin->Update(); |
| 755 |
> |
mcstack.Draw("histo,same"); |
| 756 |
> |
|
| 757 |
|
datahisto->Draw("same,e1"); |
| 758 |
|
TLegend *kinleg = allsamples.allbglegend(); |
| 759 |
|
kinleg->Draw(); |
| 795 |
|
kinpad->cd(); |
| 796 |
|
kinpad->SetLogy(logscale); |
| 797 |
|
datahisto->Draw("e1"); |
| 798 |
< |
mcstack.Draw("same"); |
| 798 |
> |
mcstack.Draw("histo,same"); |
| 799 |
> |
if(drawsignal) signalhisto->Draw("same"); |
| 800 |
|
datahisto->Draw("same,e1"); |
| 801 |
|
datahisto->Draw("same,axis"); |
| 802 |
|
if(RestrictToMassPeak&&draw_separation_lines) { |
| 804 |
|
upperboundary->Draw("same"); |
| 805 |
|
} |
| 806 |
|
|
| 807 |
+ |
if(drawsignal) kinleg->AddEntry("signalhisto",signal.c_str(),"l"); |
| 808 |
|
kinleg->Draw(); |
| 809 |
|
write_cut->Draw(); |
| 810 |
|
DrawPrelim(); |
| 819 |
|
else CompleteSave(ckin,"kin/"+filename); |
| 820 |
|
} |
| 821 |
|
datahisto->Delete(); |
| 822 |
+ |
delete signalhisto; |
| 823 |
|
delete ckin; |
| 824 |
|
} |
| 825 |
|
|
| 826 |
< |
void make_JES_plot() { |
| 826 |
> |
|
| 827 |
> |
void make_plain_kin_plot(string variable, string addcut, int nbins, float min, float max, bool logscale, |
| 828 |
> |
string xlabel, string filename, bool isPF=true, bool plotratio=true, bool loadlastminmax=false ) { |
| 829 |
> |
TCut cut=TCut(addcut.c_str())&&basiccut; |
| 830 |
> |
|
| 831 |
> |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 832 |
> |
ckin->SetLogy(logscale); |
| 833 |
> |
TH1F *datahisto = allsamples.Draw("datahisto",variable,nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 834 |
> |
datahisto->SetMarkerSize(DataMarkerSize); |
| 835 |
> |
if ( !logscale ) datahisto->SetMinimum(0); // Defaults |
| 836 |
> |
else datahisto->SetMinimum(0.5); |
| 837 |
> |
// Custom max. |
| 838 |
> |
if (logscale) datahisto->SetMaximum(5.3*datahisto->GetMaximum()); |
| 839 |
> |
else datahisto->SetMaximum(1.3*datahisto->GetMaximum()); |
| 840 |
> |
|
| 841 |
> |
cout << "******** Cut used : " << (const char*) cut << " for plot " << filename << endl; |
| 842 |
> |
if(loadlastminmax) { |
| 843 |
> |
datahisto->SetMinimum(lastrange_min); |
| 844 |
> |
datahisto->SetMaximum(lastrange_max); |
| 845 |
> |
if(logscale) { |
| 846 |
> |
datahisto->SetMinimum(pow(10,lastrange_min)); |
| 847 |
> |
datahisto->SetMaximum(pow(10,lastrange_max)); |
| 848 |
> |
} |
| 849 |
> |
} |
| 850 |
> |
|
| 851 |
> |
THStack mcstack = allsamples.DrawStack("mcstack", variable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 852 |
> |
datahisto->Draw("e1"); |
| 853 |
> |
ckin->Update(); |
| 854 |
> |
mcstack.Draw("histo,same"); |
| 855 |
> |
|
| 856 |
> |
datahisto->Draw("same,e1"); |
| 857 |
> |
TLegend *kinleg = allsamples.allbglegend(); |
| 858 |
> |
kinleg->Draw(); |
| 859 |
> |
|
| 860 |
> |
lastrange_min=ckin->GetUymin(); |
| 861 |
> |
lastrange_max=ckin->GetUymax(); |
| 862 |
> |
|
| 863 |
> |
if ( plotratio ) { |
| 864 |
> |
TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1); |
| 865 |
> |
kinpad->cd(); |
| 866 |
> |
kinpad->SetLogy(logscale); |
| 867 |
> |
datahisto->Draw("e1"); |
| 868 |
> |
mcstack.Draw("histo,same"); |
| 869 |
> |
datahisto->Draw("same,e1"); |
| 870 |
> |
datahisto->Draw("same,axis"); |
| 871 |
> |
|
| 872 |
> |
kinleg->Draw(); |
| 873 |
> |
DrawPrelim(); |
| 874 |
> |
string saveas="kin/"+filename; |
| 875 |
> |
if(isPF) saveas="kin/"+filename+"__PF"; |
| 876 |
> |
save_with_ratio(datahisto,mcstack,kinpad->cd(),saveas); |
| 877 |
> |
} else { |
| 878 |
> |
if(isPF) CompleteSave(ckin,"kin/"+filename+"__PF"); |
| 879 |
> |
else CompleteSave(ckin,"kin/"+filename); |
| 880 |
> |
} |
| 881 |
> |
datahisto->Delete(); |
| 882 |
> |
delete ckin; |
| 883 |
> |
} |
| 884 |
> |
|
| 885 |
> |
|
| 886 |
> |
void make_JES_plot(TCut cut, string name) { |
| 887 |
|
|
| 888 |
|
int nbins=10; |
| 889 |
|
float min=-0.5; |
| 890 |
|
float max = 9.5; |
| 891 |
|
bool logscale=true; |
| 892 |
|
string xlabel="nJets"; |
| 468 |
– |
TCut cut=cutmass&&cutOSSF&&basiccut; |
| 893 |
|
|
| 894 |
|
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 895 |
|
ckin->SetLogy(logscale); |
| 896 |
< |
TH1F *datahisto = allsamples.Draw("datahisto","pfJetGoodNum",nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 896 |
> |
TH1F *datahisto = allsamples.Draw("datahisto","pfJetGoodNum40",nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 897 |
|
datahisto->SetMarkerSize(DataMarkerSize); |
| 898 |
< |
THStack mcstack = allsamples.DrawStack("mcstack","pfJetGoodNum",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 899 |
< |
TH1F *JESup = allsamples.Draw("JESup","pfJetGoodNum315",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 900 |
< |
TH1F *JESdn = allsamples.Draw("JESdn","pfJetGoodNum285",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 898 |
> |
THStack mcstack = allsamples.DrawStack("mcstack","pfJetGoodNum40",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 899 |
> |
TH1F *JESup = allsamples.Draw("JESup","pfJetGoodNum40p1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 900 |
> |
TH1F *JESdn = allsamples.Draw("JESdn","pfJetGoodNum40n1sigma",nbins,min,max, xlabel, "events",cut,mc,luminosity); |
| 901 |
|
|
| 902 |
|
datahisto->SetMinimum(1); |
| 903 |
|
datahisto->SetMaximum(5.3*datahisto->GetMaximum()); // in line with kinematic plots style |
| 916 |
|
JESunc->SetFillColor(TColor::GetColor("#00ADE1")); |
| 917 |
|
JESunc->SetFillStyle(3002); |
| 918 |
|
datahisto->Draw("e1"); |
| 919 |
< |
mcstack.Draw("same"); |
| 919 |
> |
mcstack.Draw("histo,same"); |
| 920 |
|
JESunc->Draw("2"); |
| 921 |
|
datahisto->Draw("same,e1"); |
| 922 |
|
TLegend *kinleg = allsamples.allbglegend(); |
| 923 |
|
kinleg->AddEntry(JESunc,"JES uncertainty","f"); |
| 924 |
|
kinleg->Draw(); |
| 925 |
< |
CompleteSave(ckin,"Systematics/JES"); |
| 925 |
> |
CompleteSave(ckin,"Systematics/JES"+name); |
| 926 |
|
datahisto->Delete(); |
| 927 |
|
delete ckin; |
| 928 |
|
|
| 929 |
|
} |
| 930 |
|
|
| 931 |
+ |
string Cut2Str(TCut cut) { |
| 932 |
+ |
return ((string)(const char*)cut); |
| 933 |
+ |
} |
| 934 |
+ |
|
| 935 |
+ |
void MakeElegantTwoThreeComparisons() { |
| 936 |
+ |
|
| 937 |
+ |
TCut signal("met[4]>100&&met[4]<150&&mll>20"&&PlottingSetup::basicqualitycut); |
| 938 |
+ |
TCut TwoJets("pfJetGoodNum40==2"); |
| 939 |
+ |
TCut ThreeJets("pfJetGoodNum40>=3"); |
| 940 |
+ |
|
| 941 |
+ |
float min=20; |
| 942 |
+ |
float max=300; |
| 943 |
+ |
int nbins=int((max-min)/5.0); |
| 944 |
+ |
string xlabel="m_{ll} [GeV]"; |
| 945 |
+ |
|
| 946 |
+ |
TCanvas *can2 = new TCanvas("can2","can2"); |
| 947 |
+ |
TH1F *TwoOF = allsamples.Draw("TwoOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_")); |
| 948 |
+ |
TH1F *TwoSF = allsamples.Draw("TwoSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&TwoJets,mc,luminosity,allsamples.FindSample("TTJets_")); |
| 949 |
+ |
TH1F *ThreeOF = allsamples.Draw("ThreeOF","mll",nbins,min,max, xlabel, "events",cutOSOF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_")); |
| 950 |
+ |
TH1F *ThreeSF = allsamples.Draw("ThreeSF","mll",nbins,min,max, xlabel, "events",cutOSSF&&signal&&ThreeJets,mc,luminosity,allsamples.FindSample("TTJets_")); |
| 951 |
+ |
|
| 952 |
+ |
can2->cd(1); |
| 953 |
+ |
|
| 954 |
+ |
TwoSF->SetFillColor(kWhite); |
| 955 |
+ |
TwoOF->SetFillColor(kWhite); |
| 956 |
+ |
ThreeSF->SetFillColor(kWhite); |
| 957 |
+ |
ThreeOF->SetFillColor(kWhite); |
| 958 |
+ |
|
| 959 |
+ |
TPad *kin3 = new TPad("kin3","kin3",0,0,1,1); |
| 960 |
+ |
kin3->cd(); |
| 961 |
+ |
TwoOF->Draw("histo"); |
| 962 |
+ |
TwoSF->Draw("e1,same"); |
| 963 |
+ |
TLegend *leg_SFOF_Two = make_legend(); |
| 964 |
+ |
leg_SFOF_Two->SetHeader("t#bar{t} MC:"); |
| 965 |
+ |
leg_SFOF_Two->AddEntry(TwoOF,"OF, nJets==2","l"); |
| 966 |
+ |
leg_SFOF_Two->AddEntry(TwoSF,"SF, nJets==2","lp"); |
| 967 |
+ |
leg_SFOF_Two->Draw(); |
| 968 |
+ |
save_with_ratio(TwoSF,TwoOF,kin3->cd(),"DoubleRatios/SFOF_Two"); |
| 969 |
+ |
|
| 970 |
+ |
TPad *kin4 = new TPad("kin4","kin4",0,0,1,1); |
| 971 |
+ |
kin4->cd(); |
| 972 |
+ |
ThreeOF->Draw("histo"); |
| 973 |
+ |
ThreeSF->Draw("e1,same"); |
| 974 |
+ |
TLegend *leg_SFOF_Three = make_legend(); |
| 975 |
+ |
leg_SFOF_Three->SetHeader("t#bar{t} MC:"); |
| 976 |
+ |
leg_SFOF_Three->AddEntry(TwoOF,"OF, nJet#geq3","l"); |
| 977 |
+ |
leg_SFOF_Three->AddEntry(TwoSF,"SF, nJets#geq3","lp"); |
| 978 |
+ |
leg_SFOF_Three->Draw(); |
| 979 |
+ |
save_with_ratio(ThreeSF,ThreeOF,kin4->cd(),"DoubleRatios/SFOF_Three"); |
| 980 |
+ |
|
| 981 |
+ |
|
| 982 |
+ |
TwoSF->Scale(ThreeSF->Integral()/TwoSF->Integral()); |
| 983 |
+ |
TPad *kin = new TPad("kin","kin",0,0,1,1); |
| 984 |
+ |
kin->cd(); |
| 985 |
+ |
TwoSF->Draw("histo"); |
| 986 |
+ |
ThreeSF->Draw("e1,same"); |
| 987 |
+ |
TLegend *leg_SF_Two_Three = make_legend(); |
| 988 |
+ |
leg_SF_Two_Three->SetHeader("t#bar{t} MC:"); |
| 989 |
+ |
leg_SF_Two_Three->AddEntry(TwoSF,"#splitline{SF, nJets==2}{ (scaled)}","l"); |
| 990 |
+ |
leg_SF_Two_Three->AddEntry(ThreeSF,"SF, nJets#geq3","lp"); |
| 991 |
+ |
leg_SF_Two_Three->Draw(); |
| 992 |
+ |
save_with_ratio(TwoSF,ThreeSF,kin->cd(),"DoubleRatios/SF_Two_Three"); |
| 993 |
+ |
|
| 994 |
+ |
TwoOF->Scale(ThreeOF->Integral()/TwoOF->Integral()); |
| 995 |
+ |
TPad *kin2 = new TPad("kin2","kin2",0,0,1,1); |
| 996 |
+ |
kin2->cd(); |
| 997 |
+ |
TwoOF->Draw("histo"); |
| 998 |
+ |
ThreeOF->Draw("e1,same"); |
| 999 |
+ |
TLegend *leg_OF_Two_Three = make_legend(); |
| 1000 |
+ |
leg_OF_Two_Three->SetHeader("t#bar{t} MC:"); |
| 1001 |
+ |
leg_OF_Two_Three->AddEntry(TwoSF,"#splitline{OF, nJets==2}{ (scaled)}","l"); |
| 1002 |
+ |
leg_OF_Two_Three->AddEntry(ThreeSF,"OF, nJets#geq3","lp"); |
| 1003 |
+ |
leg_OF_Two_Three->Draw(); |
| 1004 |
+ |
save_with_ratio(TwoOF,ThreeOF,kin2->cd(),"DoubleRatios/OF_Two_Three"); |
| 1005 |
+ |
|
| 1006 |
+ |
} |
| 1007 |
+ |
|
| 1008 |
+ |
void PresentRange(TH1F *Down, TH1F *Up, TH1F *central, TVirtualPad *pad, string title) { |
| 1009 |
+ |
pad->cd(); |
| 1010 |
+ |
TGraphErrors *gr = new TGraphErrors(Down->GetNbinsX()); |
| 1011 |
+ |
for(int i=1;i<=Down->GetNbinsX();i++) { |
| 1012 |
+ |
float average=0.5*(Down->GetBinContent(i)+Up->GetBinContent(i)); |
| 1013 |
+ |
float error=abs(Down->GetBinContent(i)-average); |
| 1014 |
+ |
gr->SetPoint(i-1,Down->GetBinCenter(i),average); |
| 1015 |
+ |
gr->SetPointError(i-1,Down->GetBinWidth(i)/2.0,error); |
| 1016 |
+ |
} |
| 1017 |
+ |
gr->SetFillColor(TColor::GetColor("#F79F81")); |
| 1018 |
+ |
gr->GetYaxis()->SetRangeUser(0,2); |
| 1019 |
+ |
gr->GetXaxis()->SetTitle(central->GetXaxis()->GetTitle()); |
| 1020 |
+ |
gr->GetYaxis()->SetTitle(central->GetYaxis()->GetTitle()); |
| 1021 |
+ |
gr->GetXaxis()->CenterTitle(); |
| 1022 |
+ |
gr->GetYaxis()->CenterTitle(); |
| 1023 |
+ |
|
| 1024 |
+ |
gr->Draw("A2"); |
| 1025 |
+ |
central->Draw("e1,same"); |
| 1026 |
+ |
DrawMCPrelim(); |
| 1027 |
+ |
TText *rtitle = write_text(0.8,0.8,title); |
| 1028 |
+ |
TLine *line = new TLine(Down->GetBinLowEdge(1),1.0,Down->GetBinLowEdge(Down->GetNbinsX())+Down->GetBinWidth(Down->GetNbinsX()),1.0); |
| 1029 |
+ |
line->SetLineColor(kBlue); |
| 1030 |
+ |
line->SetLineStyle(2); |
| 1031 |
+ |
line->Draw("same"); |
| 1032 |
+ |
rtitle->Draw("same"); |
| 1033 |
+ |
} |
| 1034 |
+ |
|
| 1035 |
+ |
TH1F* ProduceTTbarRatio(TH1F *datahisto, TH1F *fullmc, TH1F *originalttbar, TH1F *httbar) { |
| 1036 |
+ |
TH1F *basehisto = (TH1F*)fullmc->Clone("basehisto"); |
| 1037 |
+ |
basehisto->Add(originalttbar,-1); |
| 1038 |
+ |
basehisto->Add(httbar); |
| 1039 |
+ |
|
| 1040 |
+ |
string sname=httbar->GetName(); |
| 1041 |
+ |
sname="R"+sname.substr(1,sname.length()); |
| 1042 |
+ |
TH1F *ratio = (TH1F*)datahisto->Clone(sname.c_str()); |
| 1043 |
+ |
ratio->Divide(basehisto); |
| 1044 |
+ |
delete basehisto; |
| 1045 |
+ |
return ratio; |
| 1046 |
+ |
} |
| 1047 |
+ |
|
| 1048 |
+ |
void make_ttbar_comparison(string datavariable, string mcvariable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) { |
| 1049 |
+ |
|
| 1050 |
+ |
TCut cut=addcut&&basiccut; |
| 1051 |
+ |
|
| 1052 |
+ |
TCanvas *ckin = new TCanvas("ckin","Kinematic Plots (in the making)",600,600); |
| 1053 |
+ |
ckin->SetLogy(dolog); |
| 1054 |
+ |
cout << "\r Creating " << filename << " : data (1/11)" << std::flush; |
| 1055 |
+ |
TH1F *datahisto = allsamples.Draw("datahisto",datavariable,nbins,min,max, xlabel, "events",cut,data,luminosity); |
| 1056 |
+ |
datahisto->SetMarkerSize(DataMarkerSize); |
| 1057 |
+ |
if ( !dolog ) datahisto->SetMinimum(0); // Defaults |
| 1058 |
+ |
else datahisto->SetMinimum(0.5); |
| 1059 |
+ |
if (dolog) datahisto->SetMaximum(5.3*datahisto->GetMaximum()); |
| 1060 |
+ |
else datahisto->SetMaximum(1.3*datahisto->GetMaximum()); |
| 1061 |
+ |
cout << "\r Creating " << filename << " : MC central (stack) (2/11)" << std::flush; |
| 1062 |
+ |
THStack mcstack = allsamples.DrawStack("mcstack", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity); |
| 1063 |
+ |
TH1F *MCcentral = CollapseStack(mcstack); |
| 1064 |
+ |
|
| 1065 |
+ |
cout << "\r Creating " << filename << " : MC central (histo) (3/11)" << std::flush; |
| 1066 |
+ |
TH1F *TCentral = allsamples.Draw ("TCentral", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,allsamples.FindSample("TTJets_")); |
| 1067 |
+ |
TH1F *RCentral = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TCentral); |
| 1068 |
+ |
|
| 1069 |
+ |
cout << "\r Creating " << filename << " : MC Matching up (4/11)" << std::flush; |
| 1070 |
+ |
TH1F *TMatchingUp = systsamples.Draw("TMatchingUp", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_matchingup")); |
| 1071 |
+ |
TH1F *RMatchingUp = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TMatchingUp); |
| 1072 |
+ |
|
| 1073 |
+ |
cout << "\r Creating " << filename << " : MC Matching down (5/11)" << std::flush; |
| 1074 |
+ |
TH1F *TMatchingDown = systsamples.Draw("TMatchingDown",mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_matchingdown")); |
| 1075 |
+ |
TH1F *RMatchingDown = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TMatchingDown); |
| 1076 |
+ |
|
| 1077 |
+ |
cout << "\r Creating " << filename << " : MC Scale up (6/11)" << std::flush; |
| 1078 |
+ |
TH1F *TScaleUp = systsamples.Draw("TScaleUp", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_scaleup")); |
| 1079 |
+ |
TH1F *RScaleUp = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TScaleUp); |
| 1080 |
+ |
|
| 1081 |
+ |
cout << "\r Creating " << filename << " : MC Scale down (7/11)" << std::flush; |
| 1082 |
+ |
TH1F *TScaleDown = systsamples.Draw("TScaleDown", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_scaledown")); |
| 1083 |
+ |
TH1F *RScaleDown = ProduceTTbarRatio(datahisto,MCcentral,TCentral,TScaleDown); |
| 1084 |
+ |
|
| 1085 |
+ |
cout << "\r Creating " << filename << " : MC mass=166 (8/11)" << std::flush; |
| 1086 |
+ |
TH1F *Ttmass166 = systsamples.Draw("Ttmass166", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_mass166")); |
| 1087 |
+ |
TH1F *Rtmass166 = ProduceTTbarRatio(datahisto,MCcentral,TCentral,Ttmass166); |
| 1088 |
+ |
|
| 1089 |
+ |
cout << "\r Creating " << filename << " : MC mass=169 (9/11)" << std::flush; |
| 1090 |
+ |
TH1F *Ttmass169 = systsamples.Draw("Ttmass169", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_mass169")); |
| 1091 |
+ |
TH1F *Rtmass169 = ProduceTTbarRatio(datahisto,MCcentral,TCentral,Ttmass169); |
| 1092 |
+ |
|
| 1093 |
+ |
cout << "\r Creating " << filename << " : MC mass=175 (10/11)" << std::flush; |
| 1094 |
+ |
TH1F *Ttmass175 = systsamples.Draw("Ttmass175", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_mass175")); |
| 1095 |
+ |
TH1F *Rtmass175 = ProduceTTbarRatio(datahisto,MCcentral,TCentral,Ttmass175); |
| 1096 |
+ |
|
| 1097 |
+ |
cout << "\r Creating " << filename << " : MC mass=178 (11/11)" << std::flush; |
| 1098 |
+ |
TH1F *Ttmass178 = systsamples.Draw("Ttmass178", mcvariable,nbins,min,max,xlabel,"events",cut,mc,luminosity,systsamples.FindSample("TTJets_mass178")); |
| 1099 |
+ |
TH1F *Rtmass178 = ProduceTTbarRatio(datahisto,MCcentral,TCentral,Ttmass178); |
| 1100 |
+ |
|
| 1101 |
+ |
|
| 1102 |
+ |
datahisto->Draw("e1"); |
| 1103 |
+ |
ckin->Update(); |
| 1104 |
+ |
mcstack.Draw("histo,same"); |
| 1105 |
+ |
|
| 1106 |
+ |
datahisto->Draw("same,e1"); |
| 1107 |
+ |
TLegend *kinleg = allsamples.allbglegend(); |
| 1108 |
+ |
kinleg->Draw(); |
| 1109 |
+ |
|
| 1110 |
+ |
TPad *kinpad = new TPad("kinpad","kinpad",0,0,1,1); |
| 1111 |
+ |
kinpad->SetLogy(dolog); |
| 1112 |
+ |
kinpad->cd(); |
| 1113 |
+ |
datahisto->Draw("e1"); |
| 1114 |
+ |
mcstack.Draw("histo,same"); |
| 1115 |
+ |
datahisto->Draw("same,e1"); |
| 1116 |
+ |
datahisto->Draw("same,axis"); |
| 1117 |
+ |
|
| 1118 |
+ |
kinleg->Draw("same"); |
| 1119 |
+ |
DrawPrelim(); |
| 1120 |
+ |
string saveas="kin/"+filename; |
| 1121 |
+ |
|
| 1122 |
+ |
cout << "Passing filename to be saved : " << filename << " as " << saveas << endl; |
| 1123 |
+ |
save_with_ratio( datahisto, CollapseStack(mcstack), kinpad->cd(), saveas, false, false, "data/mc"); |
| 1124 |
+ |
|
| 1125 |
+ |
ckin->cd(); |
| 1126 |
+ |
TCentral->SetLineColor(kBlack); |
| 1127 |
+ |
|
| 1128 |
+ |
TMatchingUp->SetLineColor(TColor::GetColor("#2c17b1"));//blue |
| 1129 |
+ |
TMatchingUp->SetLineStyle(2);//dashed |
| 1130 |
+ |
TMatchingDown->SetLineColor(TColor::GetColor("#2c17b1")); |
| 1131 |
+ |
TMatchingDown->SetLineStyle(3);//dotted |
| 1132 |
+ |
|
| 1133 |
+ |
TScaleUp->SetLineColor(TColor::GetColor("#FF8500"));//orange |
| 1134 |
+ |
TScaleUp->SetLineStyle(2);//dashed |
| 1135 |
+ |
TScaleDown->SetLineColor(TColor::GetColor("#FF8500")); |
| 1136 |
+ |
TScaleDown->SetLineStyle(3);//dotted |
| 1137 |
+ |
|
| 1138 |
+ |
Ttmass166->SetLineColor(TColor::GetColor("#0acf00"));//green |
| 1139 |
+ |
Ttmass166->SetLineStyle(3);//dotted |
| 1140 |
+ |
Ttmass178->SetLineColor(TColor::GetColor("#0acf00")); |
| 1141 |
+ |
Ttmass178->SetLineStyle(2);//dashed |
| 1142 |
+ |
|
| 1143 |
+ |
|
| 1144 |
+ |
Ttmass169->SetLineColor(TColor::GetColor("#FF52F1"));//purple |
| 1145 |
+ |
Ttmass169->SetLineStyle(3);//dotted |
| 1146 |
+ |
Ttmass175->SetLineColor(TColor::GetColor("#FF52F1")); |
| 1147 |
+ |
Ttmass175->SetLineStyle(2);//dashed |
| 1148 |
+ |
|
| 1149 |
+ |
|
| 1150 |
+ |
|
| 1151 |
+ |
|
| 1152 |
+ |
MCcentral->SetLineStyle(3); |
| 1153 |
+ |
MCcentral->SetLineColor(kRed); |
| 1154 |
+ |
|
| 1155 |
+ |
TCentral->Draw("e1,same"); |
| 1156 |
+ |
TCentral->Draw("histo,same"); |
| 1157 |
+ |
TMatchingUp->Draw("histo,same"); |
| 1158 |
+ |
TMatchingDown->Draw("histo,same"); |
| 1159 |
+ |
TScaleUp->Draw("histo,same"); |
| 1160 |
+ |
TScaleDown->Draw("histo,same"); |
| 1161 |
+ |
Ttmass166->Draw("histo,same"); |
| 1162 |
+ |
Ttmass169->Draw("histo,same"); |
| 1163 |
+ |
Ttmass175->Draw("histo,same"); |
| 1164 |
+ |
Ttmass178->Draw("histo,same"); |
| 1165 |
+ |
|
| 1166 |
+ |
TLegend *sleg = make_legend(); |
| 1167 |
+ |
sleg->AddEntry(TCentral,"Central","L"); |
| 1168 |
+ |
sleg->AddEntry(TMatchingUp,"matching up","L"); |
| 1169 |
+ |
sleg->AddEntry(TMatchingDown,"matching down","L"); |
| 1170 |
+ |
sleg->AddEntry(TScaleUp,"scale up","L"); |
| 1171 |
+ |
sleg->AddEntry(TScaleDown,"scale down","L"); |
| 1172 |
+ |
sleg->AddEntry(Ttmass166,"m_{t}=166","L"); |
| 1173 |
+ |
sleg->AddEntry(Ttmass169,"m_{t}=169","L"); |
| 1174 |
+ |
sleg->AddEntry(Ttmass175,"m_{t}=175","L"); |
| 1175 |
+ |
sleg->AddEntry(Ttmass178,"m_{t}=178","L"); |
| 1176 |
+ |
sleg->Draw(); |
| 1177 |
+ |
|
| 1178 |
+ |
CompleteSave(ckin,saveas+"___AllLines"); |
| 1179 |
+ |
ckin->cd(); |
| 1180 |
+ |
ckin->SetLogy(0); |
| 1181 |
+ |
|
| 1182 |
+ |
RCentral->SetLineColor(kBlack); |
| 1183 |
+ |
|
| 1184 |
+ |
RMatchingUp->SetLineColor(TColor::GetColor("#2c17b1"));//blue |
| 1185 |
+ |
RMatchingUp->SetLineStyle(2);//dashed |
| 1186 |
+ |
RMatchingDown->SetLineColor(TColor::GetColor("#2c17b1")); |
| 1187 |
+ |
RMatchingDown->SetLineStyle(3);//dotted |
| 1188 |
+ |
|
| 1189 |
+ |
RScaleUp->SetLineColor(TColor::GetColor("#FF8500"));//orange |
| 1190 |
+ |
RScaleUp->SetLineStyle(2);//dashed |
| 1191 |
+ |
RScaleDown->SetLineColor(TColor::GetColor("#FF8500")); |
| 1192 |
+ |
RScaleDown->SetLineStyle(3);//dotted |
| 1193 |
+ |
|
| 1194 |
+ |
Rtmass166->SetLineColor(TColor::GetColor("#0acf00"));//green |
| 1195 |
+ |
Rtmass166->SetLineStyle(3);//dotted |
| 1196 |
+ |
Rtmass178->SetLineColor(TColor::GetColor("#0acf00")); |
| 1197 |
+ |
Rtmass178->SetLineStyle(2);//dashed |
| 1198 |
+ |
|
| 1199 |
+ |
|
| 1200 |
+ |
Rtmass169->SetLineColor(TColor::GetColor("#FF52F1"));//purple |
| 1201 |
+ |
Rtmass169->SetLineStyle(3);//dotted |
| 1202 |
+ |
Rtmass175->SetLineColor(TColor::GetColor("#FF52F1")); |
| 1203 |
+ |
Rtmass175->SetLineStyle(2);//dashed |
| 1204 |
+ |
|
| 1205 |
+ |
|
| 1206 |
+ |
RCentral->GetYaxis()->SetRangeUser(0,2); |
| 1207 |
+ |
RCentral->Draw("e1"); |
| 1208 |
+ |
RMatchingUp->Draw("histo,same"); |
| 1209 |
+ |
RMatchingDown->Draw("histo,same"); |
| 1210 |
+ |
RScaleUp->Draw("histo,same"); |
| 1211 |
+ |
RScaleDown->Draw("histo,same"); |
| 1212 |
+ |
Rtmass166->Draw("histo,same"); |
| 1213 |
+ |
Rtmass169->Draw("histo,same"); |
| 1214 |
+ |
Rtmass175->Draw("histo,same"); |
| 1215 |
+ |
Rtmass178->Draw("histo,same"); |
| 1216 |
+ |
|
| 1217 |
+ |
CompleteSave(ckin,saveas+"___AllRatios"); |
| 1218 |
+ |
|
| 1219 |
+ |
TCanvas *multicanvas = new TCanvas("multicanvas","multicanvas",1000,1000); |
| 1220 |
+ |
multicanvas->Divide(2,2); |
| 1221 |
+ |
PresentRange(RMatchingUp,RMatchingDown,RCentral,multicanvas->cd(1),"Matching"); |
| 1222 |
+ |
PresentRange(RScaleUp,RScaleDown,RCentral,multicanvas->cd(2),"Scale"); |
| 1223 |
+ |
PresentRange(Rtmass166,Rtmass178,RCentral,multicanvas->cd(3),"m_{t} (2#sigma)"); |
| 1224 |
+ |
PresentRange(Rtmass169,Rtmass175,RCentral,multicanvas->cd(4),"m_{t} (1#sigma)"); |
| 1225 |
+ |
CompleteSave(multicanvas,saveas+"___RangeIllustration"); |
| 1226 |
+ |
|
| 1227 |
+ |
|
| 1228 |
+ |
delete datahisto; |
| 1229 |
+ |
delete TCentral; |
| 1230 |
+ |
delete TMatchingUp; |
| 1231 |
+ |
delete TMatchingDown; |
| 1232 |
+ |
delete TScaleUp; |
| 1233 |
+ |
delete TScaleDown; |
| 1234 |
+ |
delete Ttmass166; |
| 1235 |
+ |
delete Ttmass169; |
| 1236 |
+ |
delete Ttmass175; |
| 1237 |
+ |
delete Ttmass178; |
| 1238 |
+ |
|
| 1239 |
+ |
delete ckin; |
| 1240 |
+ |
|
| 1241 |
+ |
} |
| 1242 |
+ |
|
| 1243 |
+ |
void make_ttbar_comparison(string variable, int nbins, float min, float max, TCut addcut, bool dolog, string xlabel, string filename) { |
| 1244 |
+ |
make_ttbar_comparison(variable, variable, nbins, min, max, addcut, dolog, xlabel, filename); |
| 1245 |
+ |
} |
| 1246 |
+ |
|
| 1247 |
|
void do_kinematic_plots(string mcjzb, string datajzb, bool doPF=false) |
| 1248 |
|
{ |
| 1249 |
+ |
// switch_overunderflow(true); |
| 1250 |
|
bool dolog=true; |
| 1251 |
|
bool nolog=false; |
| 1252 |
+ |
|
| 1253 |
+ |
bool doOFSF = false; |
| 1254 |
+ |
bool doKin = true; |
| 1255 |
+ |
bool doDataComp = false; |
| 1256 |
+ |
bool MakeTwoThreeJetComparison = false; |
| 1257 |
+ |
|
| 1258 |
+ |
|
| 1259 |
|
if(doPF) write_warning(__FUNCTION__,"Please use caution when trying to produce PF plots; not all versions of the JZB trees have these variables!"); |
| 1260 |
|
float mll_low=50; |
| 1261 |
|
float mll_hi=160; |
| 1262 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
| 1263 |
|
mll_low=20; |
| 1264 |
< |
mll_hi=210; |
| 1264 |
> |
mll_hi=320; |
| 1265 |
|
} |
| 1266 |
|
|
| 1267 |
< |
make_kin_plot("met[4]","",40,0,200,dolog,"MET [GeV]","met",doPF,true); |
| 1268 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true); |
| 1269 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true); |
| 1270 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true); |
| 1271 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true); |
| 1272 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true); |
| 1273 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true); |
| 1274 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low)),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true); |
| 1275 |
< |
make_kin_plot("mll","",(int)((mll_hi-mll_low))/2,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true); |
| 1276 |
< |
make_kin_plot("mll","",(int)((350-mll_low)),mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF); |
| 1277 |
< |
make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF); |
| 1278 |
< |
// make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF); |
| 1279 |
< |
// make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF); |
| 1280 |
< |
make_kin_plot("pt","",50,0,400,dolog,"Z p_{T} [GeV]","Zpt",doPF); |
| 1281 |
< |
make_kin_plot("pt1","",50,0,100,nolog,"p_{T} [GeV]","pt1",doPF); |
| 1282 |
< |
make_kin_plot("pt2","",50,0,100,nolog,"p_{T} [GeV]","pt2",doPF); |
| 1283 |
< |
make_kin_plot("eta1","",40,-5,5,nolog,"#eta_{l}","eta",doPF); |
| 1284 |
< |
make_kin_plot("jzb[1]","",100,-150,150,dolog,"JZB [GeV]","jzb_ossf",doPF); |
| 1285 |
< |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets",doPF); |
| 1286 |
< |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF); |
| 1287 |
< |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF); |
| 1288 |
< |
if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]"); |
| 1289 |
< |
stringstream jzbcut; |
| 1290 |
< |
jzbcut << "((is_data&&("<<datajzb<<")>100)||(!is_data&&("<<mcjzb<<")>100))"; |
| 1291 |
< |
make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true); |
| 1292 |
< |
make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_aboveJZB100",doPF,true); |
| 1293 |
< |
make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true); |
| 1294 |
< |
make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_aboveJZB100",doPF,true); |
| 1295 |
< |
stringstream jzbcut2; |
| 1296 |
< |
jzbcut2 << "((is_data&&("<<datajzb<<")>150)||(!is_data&&("<<mcjzb<<")>150))"; |
| 1297 |
< |
make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_aboveJZB150",doPF,true); |
| 1298 |
< |
make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_aboveJZB150",doPF,true); |
| 1299 |
< |
stringstream jzbcut3; |
| 1300 |
< |
jzbcut3 << "((is_data&&("<<datajzb<<")>50)||(!is_data&&("<<mcjzb<<")>50))"; |
| 1301 |
< |
make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_aboveJZB50",doPF,true); |
| 1302 |
< |
make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_aboveJZB50",doPF,true,true); |
| 1303 |
< |
|
| 1304 |
< |
make_special_obs_pred_mll_plot(mcjzb,0); |
| 1305 |
< |
make_special_obs_pred_mll_plot(mcjzb,50); |
| 1306 |
< |
make_special_obs_pred_mll_plot(mcjzb,100); |
| 1307 |
< |
make_special_obs_pred_mll_plot(mcjzb,150); |
| 1308 |
< |
make_special_obs_pred_mll_plot(mcjzb,200); |
| 1309 |
< |
make_special_obs_pred_mll_plot(mcjzb,250); |
| 1267 |
> |
if(MakeTwoThreeJetComparison) MakeElegantTwoThreeComparisons(); |
| 1268 |
> |
|
| 1269 |
> |
if ( doOFSF ) { |
| 1270 |
> |
make_OFSF_plots("mll", "met[4]>100", 60, 20., 320., false, "m_{ll}", "mll"); |
| 1271 |
> |
|
| 1272 |
> |
// make_OFSF_plots("pfJetGoodNum40", "met[4]>100", 7, 3, 10, true, "#(jets)", "njets"); |
| 1273 |
> |
// make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30==0", 7, 3, 10, true, "#(jets)", "njets_btagVeto"); |
| 1274 |
> |
// make_OFSF_plots("pfJetGoodNum40", "met[4]>100&&pfJetGoodNumBtag30>0", 7, 3, 10, true, "#(jets)", "njets_AtLeastOneBJet30"); |
| 1275 |
> |
|
| 1276 |
> |
// make_OFSF_plots("pfJetGoodNumBtag30", "met[4]>100", 5, 0, 5, true, "#(b-jets)", "nbjets"); |
| 1277 |
> |
// make_OFSF_plots("pfJetGoodPtBtag[0]", "met[4]>100&&pfJetGoodNumBtag30>0", 20, 0, 400, true, "p_{T}(leading b-jet)", "ptb1"); |
| 1278 |
> |
|
| 1279 |
> |
// make_OFSF_plots("iso1", "met[4]>100", 20, 0, 0.3, true, "lepton 1 isolation", "iso1"); |
| 1280 |
> |
// make_OFSF_plots("iso2", "met[4]>100", 20, 0, 0.3, true, "lepton 2 isolation", "iso2"); |
| 1281 |
> |
// make_OFSF_plots("pt1", "met[4]>100", 30, 0., 300., true, "p_{T,1}", "pt1"); |
| 1282 |
> |
// make_OFSF_plots("pt2", "met[4]>100", 22, 0., 220., true, "p_{T,2}", "pt2"); |
| 1283 |
> |
// make_OFSF_plots("eta1", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{1}", "eta1", 0.15); |
| 1284 |
> |
// make_OFSF_plots("eta2", "met[4]>100", 10, -2.5, 2.5, false, "#eta_{2}", "eta2", 0.15); |
| 1285 |
> |
// make_OFSF_plots("phi1", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{1}", "phi1", 0.2); |
| 1286 |
> |
// make_OFSF_plots("phi2", "met[4]>100", 10, -TMath::Pi(), TMath::Pi(), false, "#phi_{2}", "phi2", 0.2); |
| 1287 |
> |
// make_OFSF_plots("pfJetGoodPt[0]/pfJetGoodPt[1]", "met[4]>100", 20, 1, 10, true, "pt_{j}^{1}/pt_{j}^{2}", "jpt1pt2", 0.2); |
| 1288 |
> |
make_OFSF_plots("TMath::Abs(pfJetDphiMet[0])", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(jet1,MET)|", "dphij1met", 0.2); |
| 1289 |
> |
// make_OFSF_plots("TMath::Abs(dphi)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,l2)|", "dphi", 0.2); |
| 1290 |
> |
// make_OFSF_plots("TMath::Abs(dphiMet1)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l1,MET)|", "dphiMet1", 0.2); |
| 1291 |
> |
// make_OFSF_plots("TMath::Abs(dphiMet2)", "met[4]>100", 16, 0, 3.2, false, "|#Delta#phi(l2,MET)|", "dphiMet2", 0.2); |
| 1292 |
> |
// 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"); |
| 1293 |
> |
// 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"); |
| 1294 |
> |
// 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"); |
| 1295 |
> |
// 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"); |
| 1296 |
> |
} |
| 1297 |
> |
|
| 1298 |
> |
if ( doDataComp) { |
| 1299 |
> |
TCut mllCut(""); |
| 1300 |
> |
float massmin = 15.; |
| 1301 |
> |
float massmax = 315; |
| 1302 |
> |
int massnbins = 60; |
| 1303 |
> |
/* if ( !PlottingSetup::openBox ) { |
| 1304 |
> |
mllCut = "mll>70"; |
| 1305 |
> |
massmin = 120; |
| 1306 |
> |
massmax = 360; |
| 1307 |
> |
massnbins = 14; |
| 1308 |
> |
}*/ |
| 1309 |
> |
|
| 1310 |
> |
TCut cutSignal = cutmass&&cutnJets&&"met[4]>100"; |
| 1311 |
> |
|
| 1312 |
> |
/* make_data_comparison_plot("numVtx", "",40, -0.5, 39.5,-1., true, "N(Vertices)", "numVtx"); |
| 1313 |
> |
make_data_comparison_plot("pfJetGoodEta[0]", "",50, -3.0, 3.0,-1., true, "N(Jets)", "pfJetGoodEta0"); |
| 1314 |
> |
|
| 1315 |
> |
|
| 1316 |
> |
make_data_comparison_plot("pfJetGoodNum40", cutOSOF,10, -0.5, 9.5,-1., true, "N(Jets)", "njets_OF"); |
| 1317 |
> |
make_data_comparison_plot("met[4]", cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_inclusive"); |
| 1318 |
> |
make_data_comparison_plot("met[1]", cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_inclusive"); |
| 1319 |
> |
make_data_comparison_plot("met[4]", cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_inclusive"); |
| 1320 |
> |
make_data_comparison_plot("met[1]", cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_inclusive"); |
| 1321 |
> |
|
| 1322 |
> |
make_data_comparison_plot("met[4]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "pfMET", "met_SF_exclusive"); |
| 1323 |
> |
// make_data_comparison_plot("met[1]", cutnJets&&cutOSSF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_SF_exclusive"); |
| 1324 |
> |
make_data_comparison_plot("met[4]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "pfMET", "met_OF_exclusive"); |
| 1325 |
> |
// make_data_comparison_plot("met[1]", cutnJets&&cutOSOF,60, 0., 300.,-1., true, "type 1 corrected pfMET", "mett1_OF_exclusive"); |
| 1326 |
> |
*/ |
| 1327 |
> |
|
| 1328 |
> |
string basiccutsaved = (const char*)basiccut; |
| 1329 |
> |
|
| 1330 |
> |
string newbasiccut=basiccutsaved; |
| 1331 |
> |
|
| 1332 |
> |
size_t found; |
| 1333 |
> |
found = newbasiccut.find( "pt2>20" ); |
| 1334 |
> |
while (found!=string::npos){ |
| 1335 |
> |
newbasiccut.replace( found, string( "pt2>20" ).length(), "pt2>10" ); |
| 1336 |
> |
found = newbasiccut.find( "pt2>20" ); |
| 1337 |
> |
} |
| 1338 |
> |
// basiccut=TCut(newbasiccut.c_str()); |
| 1339 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSSF,60, 0., 300.,-1., false, "mll", "mll_SF_Aachen_pt2010_met"); |
| 1340 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[4]>150")&&cutOSOF,60, 0., 300.,-1., false, "mll", "mll_OF_Aachen_pt2010_met"); |
| 1341 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_Aachen_pt2010_t1cpfmet"); |
| 1342 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=2&&met[1]>150")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_Aachen_pt2010_t1cpfmet"); |
| 1343 |
> |
// basiccut=TCut(basiccutsaved.c_str()); |
| 1344 |
> |
|
| 1345 |
> |
// make_data_comparison_plot("pfJetGoodNum40", cutOSSF, 8, 0., 8.,-1., true, "#(jets)", "njets_SF"); |
| 1346 |
> |
// make_data_comparison_plot("pfJetGoodNum40", cutOSOF, 8, 0., 8.,-1., true, "#(jets)", "njets_OF"); |
| 1347 |
> |
|
| 1348 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_met"); |
| 1349 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[4]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_met"); |
| 1350 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSSF,60, 0., 300.,-1., true, "mll", "mll_SF_ETH_t1cpfmet"); |
| 1351 |
> |
// make_data_comparison_plot("mll", TCut("pfJetGoodNum40>=3&&met[1]>100")&&cutOSOF,60, 0., 300.,-1., true, "mll", "mll_OF_ETH_t1cpfmet"); |
| 1352 |
> |
|
| 1353 |
> |
|
| 1354 |
> |
/* make_data_comparison_plot("mll", cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive"); |
| 1355 |
> |
make_data_comparison_plot("mll", cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_ee"); |
| 1356 |
> |
make_data_comparison_plot("mll", cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_inclusive_mm"); |
| 1357 |
> |
make_data_comparison_plot("mll", cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_inclusive"); |
| 1358 |
> |
make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==0"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_ee"); |
| 1359 |
> |
make_data_comparison_plot("mll", cutnJets&&cutOSSF&&TCut("id1==1"),60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive_mm"); |
| 1360 |
> |
make_data_comparison_plot("mll", cutnJets&&cutOSSF,60, 15., 315.,-1., true, "m_{ll}", "mll_SF_exclusive"); |
| 1361 |
> |
make_data_comparison_plot("mll", cutnJets&&cutOSOF,60, 15., 315.,-1., true, "m_{ll}", "mll_OF_exclusive"); |
| 1362 |
> |
*/ |
| 1363 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig"); |
| 1364 |
> |
make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig"); |
| 1365 |
> |
|
| 1366 |
> |
make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSSF&&cutSignal&&mllCut, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_central"); |
| 1367 |
> |
make_data_comparison_plot("mll", TCut("abs(eta1)<1.4&&abs(eta2)<1.4") && cutOSOF&&cutSignal&&mllCut, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_central"); |
| 1368 |
> |
|
| 1369 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_btagVeto"); |
| 1370 |
> |
make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30==0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_btagVeto"); |
| 1371 |
> |
|
| 1372 |
> |
make_data_comparison_plot("mll", cutOSSF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_sig_AtLeastOneBJet"); |
| 1373 |
> |
make_data_comparison_plot("mll", cutOSOF&&cutSignal&&mllCut&&"pfJetGoodNumBtag30>0", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig_AtLeastOneBJet"); |
| 1374 |
> |
|
| 1375 |
> |
/* |
| 1376 |
> |
make_data_comparison_plot("mll", mllCut&&cutOSOF&&cutSignal, massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_sig"); |
| 1377 |
> |
make_data_comparison_plot("mll", cutmass&&cutOSSF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_SF_CR"); |
| 1378 |
> |
make_data_comparison_plot("mll", cutmass&&cutOSOF&&"met[4]>100&&met[4]<150&&pfJetGoodNum40==2", massnbins, 15., 315.,-1., false, "m_{ll}", "mll_OF_CR"); |
| 1379 |
> |
|
| 1380 |
> |
make_data_comparison_plot("pfJetGoodNum40", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_SF_sig"); |
| 1381 |
> |
make_data_comparison_plot("pfJetGoodNum40", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(jets)", "njets_OF_sig"); |
| 1382 |
> |
make_data_comparison_plot("pfJetGoodNumBtag30", cutOSSF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_SF_sig"); |
| 1383 |
> |
make_data_comparison_plot("pfJetGoodNumBtag30", cutOSOF&&cutSignal&&mllCut, 8, 0., 8.,-1., false, "#(b-jets)", "nbjets_OF_sig"); |
| 1384 |
> |
*/ |
| 1385 |
> |
} |
| 1386 |
> |
|
| 1387 |
> |
|
| 1388 |
> |
if ( doKin ) { |
| 1389 |
> |
string mllCut(""); |
| 1390 |
> |
|
| 1391 |
> |
|
| 1392 |
> |
/* |
| 1393 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_0j",doPF,true); |
| 1394 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_1j",doPF,true); |
| 1395 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_2j",doPF,true); |
| 1396 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_30to70_3j",doPF,true); |
| 1397 |
> |
|
| 1398 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_0j",doPF,true); |
| 1399 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_1j",doPF,true); |
| 1400 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_2j",doPF,true); |
| 1401 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/MET_OF_Z_3j",doPF,true); |
| 1402 |
> |
|
| 1403 |
> |
|
| 1404 |
> |
make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets),120,0,1200,dolog,"H_{T}","HT_3jets_OF",doPF,true); |
| 1405 |
> |
make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF",doPF,true); |
| 1406 |
> |
make_plain_kin_plot("pfTightHT",Cut2Str(cutOSOF&&cutnJets&&TCut("met[4]>100")),120,0,1200,dolog,"H_{T}","HT_3jets_OF_MET100",doPF,true); |
| 1407 |
> |
make_plain_kin_plot("pfTightHT",Cut2Str(cutOSSF&&cutnJets&&TCut("met[4]>100")) ,120,0,1200,dolog,"H_{T}","HT_3jets_SF_MET100",doPF,true); |
| 1408 |
> |
|
| 1409 |
> |
|
| 1410 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_0j",doPF,true); |
| 1411 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_1j",doPF,true); |
| 1412 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_2j",doPF,true); |
| 1413 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_30to70_3j",doPF,true); |
| 1414 |
> |
|
| 1415 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_0j",doPF,true); |
| 1416 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_1j",doPF,true); |
| 1417 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_2j",doPF,true); |
| 1418 |
> |
make_plain_kin_plot("met[1]",Cut2Str(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),40,0,400,dolog,"PFMET","Debunking/T1MET_OF_Z_3j",doPF,true); |
| 1419 |
> |
|
| 1420 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70",doPF,true); |
| 1421 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70",doPF,true); |
| 1422 |
> |
|
| 1423 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120",doPF,true); |
| 1424 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&cutnJets&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120",doPF,true);*/ |
| 1425 |
> |
|
| 1426 |
> |
/* make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_SF_20to70___2jetsonly",doPF,true); |
| 1427 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>20&&mll<70")),50,0,300,dolog,"PFMET","MET_OF_20to70___2jetsonly",doPF,true); |
| 1428 |
> |
|
| 1429 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSSF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_SF_above120___2jetsonly",doPF,true); |
| 1430 |
> |
make_plain_kin_plot("met[4]",Cut2Str(cutOSOF&&TCut("pfJetGoodNum40==2")&&TCut("mll>120")),50,0,300,dolog,"PFMET","MET_OF_above120___2jetsonly",doPF,true); |
| 1431 |
> |
*/ |
| 1432 |
> |
|
| 1433 |
> |
|
| 1434 |
> |
|
| 1435 |
> |
|
| 1436 |
> |
|
| 1437 |
> |
|
| 1438 |
> |
|
| 1439 |
> |
// make_kin_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","",45,1,10,dolog,"pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",doPF,true); |
| 1440 |
> |
// make_kin_plot("TMath::Abs(pfJetDphiMet[0])","",32,0,3.2,nolog,"|#Delta#phi(jet1,MET)|","dphiJ1MET",doPF,true); |
| 1441 |
> |
// Plots in signal region |
| 1442 |
> |
// make_kin_plot("met[4]","",70,0,350,dolog,"MET [GeV]","met",doPF,true); |
| 1443 |
> |
|
| 1444 |
> |
|
| 1445 |
> |
|
| 1446 |
> |
|
| 1447 |
> |
|
| 1448 |
> |
make_kin_plot("MetFactor","",20,0,2,nolog,"MetFactor","MetFactor",doPF,true); |
| 1449 |
> |
|
| 1450 |
> |
make_ttbar_comparison("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j__ALLMCSAMPLES"); |
| 1451 |
> |
|
| 1452 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_OF_30to70_0j"); |
| 1453 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_OF_30to70_1j"); |
| 1454 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_30to70_2j"); |
| 1455 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>30&&mll<70&&pfTightHT>100&&pfJetGoodNum40==3")),dolog,"PFMET","Debunking/MET_OF_30to70_3j"); |
| 1456 |
> |
|
| 1457 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==0")),dolog,"PFMET","Debunking/MET_OF_Z_0j"); |
| 1458 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfJetGoodNum40==1")),dolog,"PFMET","Debunking/MET_OF_Z_1j"); |
| 1459 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==2")),dolog,"PFMET","Debunking/MET_OF_Z_2j"); |
| 1460 |
> |
DoMCSystPlot("met[4]",40,0,400,(cutOSOF&&TCut("mll>80&&mll<100&&pfTightHT>100&&pfJetGoodNum40==3")),dolog,"PFMET","Debunking/MET_OF_Z_3j"); |
| 1461 |
> |
|
| 1462 |
> |
|
| 1463 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF,dolog,"m_{ll} [GeV]","kinSys/mll"); |
| 1464 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF&&TCut("id1==0"),dolog,"m_{ll} [GeV]","kinSys/mll_ee"); |
| 1465 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSSF&&TCut("id1==1"),dolog,"m_{ll} [GeV]","kinSys/mll_mm"); |
| 1466 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutnJets&&cutOSOF,dolog,"m_{ll} [GeV]","kinSys/mll_osof"); |
| 1467 |
> |
|
| 1468 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutOSSF,dolog,"m_{ll} [GeV]","kinSys/mll_inclusive"); |
| 1469 |
> |
DoMCSystPlot("mll", (int)((mll_hi-mll_low)/5),mll_low,mll_hi,cutOSOF,dolog,"m_{ll} [GeV]","kinSys/mll_inclusive_osof"); |
| 1470 |
> |
|
| 1471 |
> |
DoMCSystPlot("pfJetGoodNum40", 9,-0.5,8.5,cutOSSF,dolog,"nJets","kinSys/nJets_ossf"); |
| 1472 |
> |
DoMCSystPlot("pfJetGoodNum40", 9,-0.5,8.5,cutOSOF,dolog,"nJets","kinSys/nJets_osof"); |
| 1473 |
> |
|
| 1474 |
> |
DoMCSystPlot("pfJetGoodPt[0]/pfJetGoodPt[1]",45,1,10,cutnJets&&cutOSSF,dolog,"p_{T}^{J1} / p_{T}^{J2}","kinSys/j1j2ratio"); |
| 1475 |
> |
DoMCSystPlot("TMath::Abs(pfJetDphiMet[0])",32,0,3.2,cutnJets&&cutOSSF,dolog,"|#Delta#phi(jet1,MET)|","kinSys/dphiJ1MET"); |
| 1476 |
> |
|
| 1477 |
> |
DoMCSystPlot(datajzb,mcjzb, (int)((PlottingSetup::jzbHigh+110)/10),-110,PlottingSetup::jzbHigh,cutnJets&&cutOSSF&&Restrmasscut,dolog,"m_{ll} [GeV]","kinSys/jzb_OS_SFZP"); |
| 1478 |
> |
DoMCSystPlot(datajzb,mcjzb, (int)((PlottingSetup::jzbHigh+110)/10),-110,PlottingSetup::jzbHigh,cutnJets&&cutOSOF&&Restrmasscut,dolog,"m_{ll} [GeV]","kinSys/jzb_OS_OFZP"); |
| 1479 |
> |
|
| 1480 |
> |
|
| 1481 |
> |
|
| 1482 |
> |
make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll",doPF,true); |
| 1483 |
> |
make_kin_plot("mll","mll>20",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof",doPF,true,true); |
| 1484 |
> |
make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_ee",doPF,true); |
| 1485 |
> |
make_kin_plot("mll","mll>20"+mllCut,(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_mm",doPF,true); |
| 1486 |
> |
|
| 1487 |
> |
make_kin_plot("pfJetGoodNum40",mllCut,9,-0.5,8.5,dolog,"nJets","nJets",doPF); |
| 1488 |
> |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof",doPF); |
| 1489 |
> |
|
| 1490 |
> |
// 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_met100",doPF,true); |
| 1491 |
> |
// make_kin_plot("mll","mll>20&&met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_osof_met100",doPF,true,true); |
| 1492 |
> |
// 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_ee_met100",doPF,true); |
| 1493 |
> |
// 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_mm_met100",doPF,true); |
| 1494 |
> |
|
| 1495 |
> |
|
| 1496 |
> |
// make_kin_plot("pfJetGoodNum40","met[4]>100"+mllCut,9,-0.5,8.5,dolog,"nJets","nJets_met100",doPF); |
| 1497 |
> |
// make_kin_plot("pfJetGoodNum40","met[4]>100",9,-0.5,8.5,dolog,"nJets","nJets_osof_met100",doPF); |
| 1498 |
> |
|
| 1499 |
> |
// Further inclusive invariant mass plots |
| 1500 |
> |
make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive",doPF,true); |
| 1501 |
> |
// make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_ee",doPF,true); |
| 1502 |
> |
// make_kin_plot("mll","",(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_mm",doPF,true); |
| 1503 |
> |
make_kin_plot("mll",mllCut,(int)((mll_hi-mll_low))/5,mll_low,mll_hi,dolog,"m_{ll} [GeV]","mll_inclusive_osof",doPF,true); |
| 1504 |
> |
|
| 1505 |
> |
//make_kin_plot("mll","",(int)((350-mll_low))/5,mll_low,350,dolog,"m_{ll} [GeV]","mll_inclusive_highrange",doPF); |
| 1506 |
> |
//if(!doPF) make_special_mll_plot((int)((mll_hi-mll_low)/5),mll_low,mll_hi,dolog,"m_{ll} [GeV]"); |
| 1507 |
> |
|
| 1508 |
> |
|
| 1509 |
> |
// Number of jets |
| 1510 |
> |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_inclusive",doPF); |
| 1511 |
> |
make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_osof_inclusive",doPF); |
| 1512 |
> |
//make_kin_plot("pfJetGoodNum40","",9,-0.5,8.5,dolog,"nJets","nJets_nocuts_except_mll_ossf",doPF); |
| 1513 |
> |
|
| 1514 |
> |
// Others |
| 1515 |
> |
make_kin_plot("numVtx","",(int)(30.5-(-0.5)),-0.5,30.5,nolog,"N(Vtx)","numVtx",doPF); |
| 1516 |
> |
// make_kin_plot("jetpt[0]","",40,0,200,dolog,"leading jet p_{T} [GeV]","pfJetGoodPt_0",doPF); |
| 1517 |
> |
// make_kin_plot("jeteta[0]","",40,-5,5,nolog,"leading jet #eta","pfJetGoodEta_0",doPF); |
| 1518 |
> |
make_kin_plot("pt","",50,0,500,dolog,"Z p_{T} [GeV]","Zpt",doPF); |
| 1519 |
> |
make_kin_plot("pt1","",50,0,200,nolog,"p_{T} [GeV]","pt1",doPF); |
| 1520 |
> |
make_kin_plot("pt2","",50,0,200,nolog,"p_{T} [GeV]","pt2",doPF); |
| 1521 |
> |
make_kin_plot("eta1","",40,-3,3,nolog,"#eta_{l}","eta",doPF); |
| 1522 |
> |
make_kin_plot("jzb[1]","",100,-150,200,dolog,"JZB [GeV]","jzb_ossf",doPF); |
| 1523 |
> |
// stringstream jzbcut; |
| 1524 |
> |
// jzbcut << "((is_data&&("<<datajzb<<")>100)||(!is_data&&("<<mcjzb<<")>100))"; |
| 1525 |
> |
// make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB100",doPF,true); |
| 1526 |
> |
// make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB100",doPF,true); |
| 1527 |
> |
// make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB100",doPF,true); |
| 1528 |
> |
// make_kin_plot("mll",jzbcut.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB100",doPF,true); |
| 1529 |
> |
// stringstream jzbcut2; |
| 1530 |
> |
// jzbcut2 << "((is_data&&("<<datajzb<<")>150)||(!is_data&&("<<mcjzb<<")>150))"; |
| 1531 |
> |
// make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB150",doPF,true); |
| 1532 |
> |
// make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB150",doPF,true); |
| 1533 |
> |
// make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB150",doPF,true); |
| 1534 |
> |
// make_kin_plot("mll",jzbcut2.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB150",doPF,true); |
| 1535 |
> |
// stringstream jzbcut3; |
| 1536 |
> |
// jzbcut3 << "((is_data&&("<<datajzb<<")>50)||(!is_data&&("<<mcjzb<<")>50))"; |
| 1537 |
> |
// make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_aboveJZB50",doPF,true); |
| 1538 |
> |
// make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_osof_aboveJZB50",doPF,true,true); |
| 1539 |
> |
// make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_ee_aboveJZB50",doPF,true); |
| 1540 |
> |
// make_kin_plot("mll",jzbcut3.str(),(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV]","mll_mm_aboveJZB50",doPF,true); |
| 1541 |
> |
|
| 1542 |
> |
// make_kin_plot("mll","met[4]>100",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ll} [GeV] (MET>100GeV)","mll_met100_ll",doPF,true); |
| 1543 |
> |
//make_kin_plot("mll","met[4]>150&&id1==0",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{ee} [GeV] (MET>150GeV)","mll_met150_ee",doPF,true); |
| 1544 |
> |
//make_kin_plot("mll","met[4]>150&&id1==1",(int)((mll_hi-mll_low)/5),mll_low,mll_hi,nolog,"m_{#mu#mu} [GeV] (MET>150GeV)","mll_met150_mm",doPF,true); |
| 1545 |
> |
} |
| 1546 |
> |
|
| 1547 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,0); |
| 1548 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,50); |
| 1549 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,80); |
| 1550 |
> |
make_special_obs_pred_mll_plot(datajzb,mcjzb,100); |
| 1551 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,150); |
| 1552 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,200); |
| 1553 |
> |
// make_special_obs_pred_mll_plot(datajzb,mcjzb,250); |
| 1554 |
> |
|
| 1555 |
> |
// make_JES_plot(cutmass&&cutOSSF&&basiccut,"_ossf"); |
| 1556 |
> |
// make_JES_plot(cutmass&&cutOSOF&&basiccut,"_osof"); |
| 1557 |
> |
|
| 1558 |
> |
switch_overunderflow(false); |
| 1559 |
|
} |
| 1560 |
|
|
| 1561 |
|
void make_comp_plot( string var, string xlabel, string filename, float jzbcut, string mcjzb, string datajzb, |
| 1562 |
|
int nbins, float xmin, float xmax, bool log, |
| 1563 |
|
float ymin=0, float ymax=0, bool leftJustified=false ) { |
| 1564 |
< |
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- WATCH OUT: the argument in the function changed! |
| 1564 |
> |
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- WATCH OUT: the argument in the function changed! |
| 1565 |
|
|
| 1566 |
|
TCut weightbackup=cutWeight;//backing up the correct weight (restoring below!) |
| 1567 |
|
if(weightbackup==TCut("1.0")||weightbackup==TCut("1")) write_warning(__FUNCTION__,"WATCH OUT THE WEIGHT HAS POSSIBLY NOT BEEN RESET!!!! PLEASE CHANGE LINE "+any2string(__LINE__)); |
| 1568 |
< |
if(var=="numVtx") cutWeight=TCut("1.0"); |
| 1568 |
> |
//if(var=="numVtx") cutWeight=TCut("1.0"); |
| 1569 |
|
TCut jzbData[]= { TCut(TString(datajzb+">"+any2string(jzbcut))),TCut(TString(datajzb+"<-"+any2string(jzbcut))) }; |
| 1570 |
|
TCut jzbMC[] = { TCut(TString(mcjzb+">"+any2string(jzbcut))),TCut(TString(mcjzb+"<-"+any2string(jzbcut))) }; |
| 1571 |
|
|
| 1572 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- below: the next ~20 lines changed! |
| 1573 |
|
int nRegions=4; |
| 1574 |
< |
if(!PlottingSetup::RestrictToMassPeak) { |
| 1574 |
> |
if(!PlottingSetup::RestrictToMassPeak||!PlottingSetup::UseSidebandsForcJZB) { |
| 1575 |
|
nRegions=2; |
| 1576 |
|
} |
| 1577 |
|
|
| 1589 |
|
TCanvas *ccomp = new TCanvas("ccomp","Comparison plot",600,400); |
| 1590 |
|
ccomp->SetLogy(log); |
| 1591 |
|
TH1F *datahisto = allsamples.Draw("datahisto", var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbData[ijzb],data,luminosity); |
| 1592 |
< |
TH1F *lm4histo = signalsamples.Draw("lm4histo", var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbData[ijzb],data,luminosity,signalsamples.FindSample("LM4")); |
| 1592 |
> |
TH1F *lm3histo = signalsamples.Draw("lm3histo", var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbMC[ijzb],data,luminosity,signalsamples.FindSample("LM3")); |
| 1593 |
|
THStack mcstack = allsamples.DrawStack("mcstack",var,nbins,xmin,xmax,xlabel,"events",kRegions[iregion]&&jzbMC[ijzb], mc, luminosity); |
| 1594 |
|
datahisto->SetMarkerSize(DataMarkerSize); |
| 1595 |
|
if (ymax>ymin) datahisto->SetMaximum(ymax); |
| 1596 |
< |
lm4histo->SetLineStyle(2); |
| 1596 |
> |
lm3histo->SetLineStyle(2); |
| 1597 |
|
datahisto->Draw("e1"); |
| 1598 |
< |
mcstack.Draw("same"); |
| 1598 |
> |
mcstack.Draw("histo,same"); |
| 1599 |
|
datahisto->Draw("same,e1"); |
| 1600 |
< |
// lm4histo->Draw("hist,same"); |
| 1600 |
> |
lm3histo->Draw("hist,same"); |
| 1601 |
|
TLegend *kinleg = allsamples.allbglegend((sRegions[iregion]+(ijzb?"neg":"pos")).c_str()); |
| 1602 |
|
if ( leftJustified ) { |
| 1603 |
|
Float_t w = kinleg->GetX2()-kinleg->GetX1(); |
| 1604 |
|
kinleg->SetX1(0.2); |
| 1605 |
|
kinleg->SetX2(0.2+w); |
| 1606 |
|
} |
| 1607 |
< |
// kinleg->AddEntry(lm4histo,"LM4","l"); |
| 1607 |
> |
kinleg->AddEntry(lm3histo,"LM3","l"); |
| 1608 |
|
kinleg->Draw(); |
| 1609 |
|
TText* write_variable = write_text(0.99,0.01,var); |
| 1610 |
|
write_variable->SetTextAlign(31); |
| 1613 |
|
CompleteSave(ccomp,"compare/JZBcut_at_"+any2string(jzbcut)+"/"+filename+"/"+filename+sRegions[iregion]+(ijzb?"neg":"pos")); |
| 1614 |
|
delete datahisto; |
| 1615 |
|
delete ccomp; |
| 1616 |
< |
delete lm4histo; |
| 1616 |
> |
delete lm3histo; |
| 1617 |
|
} |
| 1618 |
|
cutWeight=weightbackup; |
| 1619 |
|
} |
| 1625 |
|
|
| 1626 |
|
// switch_overunderflow(true); // switching overflow/underflow bins on |
| 1627 |
|
|
| 1628 |
< |
|
| 1628 |
> |
switch_overunderflow(true); |
| 1629 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- the arguments changed |
| 1630 |
|
for(int ijzb=0;ijzb<(int)jzb_cuts.size();ijzb++) { |
| 1631 |
|
float jzbcut=jzb_cuts[ijzb]; // Comparison plots are done for this JZB cut |
| 1632 |
|
float mll_low=50;float mll_high=170; |
| 1633 |
|
if(!PlottingSetup::RestrictToMassPeak) { |
| 1634 |
< |
mll_high=200; |
| 1635 |
< |
mll_low=50; |
| 1634 |
> |
mll_high=300; |
| 1635 |
> |
mll_low=20; |
| 1636 |
|
} |
| 1637 |
< |
make_comp_plot("mll","m_{ll} [GeV]","mll",jzbcut,mcjzb,datajzb,30,mll_low,mll_high,false,0,16.); |
| 1637 |
> |
make_comp_plot("pfJetGoodPt[0]/pfJetGoodPt[1]","pt_{j}^{1}/pt_{j}^{2}","j1j2ratio",jzbcut,mcjzb,datajzb,100,0,10,true); |
| 1638 |
> |
make_comp_plot("TMath::Abs(pfJetDphiMet[0])","|#Delta#phi(jet1,MET)|","dphiJ1MET",jzbcut,mcjzb,datajzb,32,0,3.2,false,0,0,true); |
| 1639 |
> |
|
| 1640 |
> |
make_comp_plot("mll","m_{ll} [GeV]","mll",jzbcut,mcjzb,datajzb,56,mll_low,mll_high,false,0,16.); |
| 1641 |
|
make_comp_plot("met[4]","pfMET [GeV]","pfmet",jzbcut,mcjzb,datajzb,18,0,360,false,0,16.); |
| 1642 |
< |
make_comp_plot("pfJetGoodNum","#(jets)","njets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.); |
| 1642 |
> |
make_comp_plot("pfJetGoodNum40","#(jets)","njets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.); |
| 1643 |
> |
make_comp_plot("pfJetGoodNumBtag","#(b-jets)","nBjets",jzbcut,mcjzb,datajzb,10,0,10, false,0,35.); |
| 1644 |
|
make_comp_plot("pt","Z p_{T} [GeV]","Zpt",jzbcut,mcjzb,datajzb,26,0,525,false,0.,21.); |
| 1645 |
< |
make_comp_plot("numVtx","#(prim. vertices)","nvtx",jzbcut,mcjzb,datajzb,20,0.,20.,false,0,16.); |
| 1645 |
> |
make_comp_plot("numVtx","#(prim. vertices)","nvtx",jzbcut,mcjzb,datajzb,40,0.,40.,false,0,16.); |
| 1646 |
|
make_comp_plot("TMath::Abs(dphi)","#Delta#phi(leptons)","dphilep",jzbcut,mcjzb,datajzb,10,0.,3.1415,false,0,16.,true); |
| 1647 |
|
make_comp_plot("TMath::Abs(dphi_sumJetVSZ[1])","#Delta#phi(Z,jets)","dphiZjets",jzbcut,mcjzb,datajzb,10,0.,3.1415,false,0,16.,true); |
| 1648 |
+ |
make_comp_plot("eta1","#eta_1","eta1",jzbcut,mcjzb,datajzb,10,0.,2.5,false,0,16.); |
| 1649 |
+ |
make_comp_plot("eta2","#eta_2","eta2",jzbcut,mcjzb,datajzb,10,0.,2.5,false,0,16.); |
| 1650 |
|
} |
| 1651 |
|
|
| 1652 |
|
switch_overunderflow(false); // switching overflow/underflow bins off |
| 1845 |
|
// the line above is not necessary anymore as we're now looking at a prediction without Z+Jets, and not multiplied with (1.0/3) |
| 1846 |
|
TF1 *ttbarlogpar = do_logpar_fit_to_plot(Tprediction); |
| 1847 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 1848 |
< |
if(PlottingSetup::RestrictToMassPeak) ttbarlogpar->SetParameter(0,1.0/3*ttbarlogpar->GetParameter(0));//correcting for the fact that we didn't multiply with (1.0/3); |
| 1848 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) ttbarlogpar->SetParameter(0,1.0/3*ttbarlogpar->GetParameter(0));//correcting for the fact that we didn't multiply with (1.0/3); |
| 1849 |
|
|
| 1850 |
|
|
| 1851 |
|
TF1 *ttbarlogparP = new TF1("ttbarlogparP",LogParabolaP,0,(prediction->GetXaxis()->GetBinLowEdge(prediction->GetNbinsX())+prediction->GetXaxis()->GetBinWidth(prediction->GetNbinsX())),8); |
| 1942 |
|
|
| 1943 |
|
void do_prediction_plot(string jzb, TCanvas *globalcanvas, float high, int use_data, bool overlay_signal = false,string subdir="" ) |
| 1944 |
|
{ |
| 1945 |
< |
// switch_overunderflow(true); |
| 1945 |
> |
|
| 1946 |
|
bool is_data=false; |
| 1947 |
|
bool use_signal=false; |
| 1948 |
|
if(use_data==1) is_data=true; |
| 1949 |
|
if(use_data==2) use_signal=true; |
| 1950 |
< |
int nbins=50;//100; |
| 1951 |
< |
if(is_data) nbins=50; |
| 1950 |
> |
int nbins=int(high/10);//100; |
| 1951 |
> |
if(is_data) nbins=int(high/10); |
| 1952 |
|
float low=0; |
| 1953 |
< |
float hi=500; |
| 1953 |
> |
float hi=high; |
| 1954 |
> |
|
| 1955 |
> |
stringstream cutpositiveS; |
| 1956 |
> |
cutpositiveS << "(" << jzb << ">0)"; |
| 1957 |
> |
TCut cutpositive(cutpositiveS.str().c_str()); |
| 1958 |
> |
stringstream cutnegativeS; |
| 1959 |
> |
cutnegativeS << "(" << jzb << "<0)"; |
| 1960 |
> |
TCut cutnegative(cutnegativeS.str().c_str()); |
| 1961 |
> |
|
| 1962 |
|
|
| 1963 |
|
TH1F *blankback = new TH1F("blankback","blankback",int(high/10),0,high); |
| 1964 |
< |
TH1F *RcorrJZBeemm = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1965 |
< |
TH1F *LcorrJZBeemm = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1966 |
< |
TH1F *RcorrJZBem = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1967 |
< |
TH1F *LcorrJZBem = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1964 |
> |
TH1F *RcorrJZBeemm = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1965 |
> |
TH1F *LcorrJZBeemm = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1966 |
> |
TH1F *RcorrJZBem = allsamples.Draw("RcorrJZBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1967 |
> |
TH1F *LcorrJZBem = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1968 |
|
|
| 1969 |
|
blankback->GetXaxis()->SetTitle(RcorrJZBeemm->GetXaxis()->GetTitle()); |
| 1970 |
|
blankback->GetYaxis()->SetTitle(RcorrJZBeemm->GetYaxis()->GetTitle()); |
| 1979 |
|
|
| 1980 |
|
TH1F *RcorrJZBeemmNoS; |
| 1981 |
|
|
| 1982 |
< |
//these are for the ratio |
| 1983 |
< |
|
| 1984 |
< |
TH1F *JRcorrJZBeemm = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1985 |
< |
TH1F *JLcorrJZBeemm = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1986 |
< |
TH1F *JRcorrJZBem = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 976 |
< |
TH1F *JLcorrJZBem = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1982 |
> |
//these are for the ratio |
| 1983 |
> |
TH1F *JRcorrJZBeemm = allsamples.Draw("JRcorrJZBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1984 |
> |
TH1F *JLcorrJZBeemm = allsamples.Draw("JLcorrJZBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1985 |
> |
TH1F *JRcorrJZBem = allsamples.Draw("JRcorrJZBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1986 |
> |
TH1F *JLcorrJZBem = allsamples.Draw("JLcorrJZBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&cutmass&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1987 |
|
|
| 1988 |
|
TH1F *JRcorrJZBSBem; |
| 1989 |
|
TH1F *JLcorrJZBSBem; |
| 1990 |
|
TH1F *JRcorrJZBSBeemm; |
| 1991 |
|
TH1F *JLcorrJZBSBeemm; |
| 1992 |
|
|
| 1993 |
< |
if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false); |
| 1993 |
> |
if(use_data==2 || overlay_signal) RcorrJZBeemmNoS = allsamples.Draw("RcorrJZBeemmNoS",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,false); |
| 1994 |
|
|
| 1995 |
|
|
| 1996 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 1997 |
< |
RcorrJZBSBem = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1998 |
< |
LcorrJZBSBem = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1999 |
< |
RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2000 |
< |
LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 1996 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 1997 |
> |
RcorrJZBSBem = allsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1998 |
> |
LcorrJZBSBem = allsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 1999 |
> |
RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2000 |
> |
LcorrJZBSBeemm = allsamples.Draw("LcorrJZBSBeemm",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2001 |
|
|
| 2002 |
|
//these are for the ratio |
| 2003 |
< |
JRcorrJZBSBem = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 2004 |
< |
JLcorrJZBSBem = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 2005 |
< |
JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2006 |
< |
JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2003 |
> |
JRcorrJZBSBem = allsamples.Draw("JRcorrJZBSBem",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 2004 |
> |
JLcorrJZBSBem = allsamples.Draw("JLcorrJZBSBem",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 2005 |
> |
JRcorrJZBSBeemm = allsamples.Draw("JRcorrJZBSBeemm",jzb.c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutpositive&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2006 |
> |
JLcorrJZBSBeemm = allsamples.Draw("JLcorrJZBSBeemm",("-"+jzb).c_str(),PlottingSetup::global_ratio_binning, "JZB [GeV]", "events", cutnegative&&sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 2007 |
|
} |
| 2008 |
|
|
| 2009 |
|
TH1F *lm4RcorrJZBeemm; |
| 2010 |
< |
if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM")); |
| 2010 |
> |
if(overlay_signal || use_data == 2 || use_data == 1) lm4RcorrJZBeemm = allsamples.Draw("lm4RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutpositive&&cutmass&&cutOSSF&&cutnJets,is_data, luminosity,allsamples.FindSample("LM")); |
| 2011 |
|
|
| 2012 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed. |
| 2013 |
|
|
| 2020 |
|
TH1F *BpredSys = new TH1F("Bpredsys","Bpredsys",PlottingSetup::global_ratio_binning.size()-1,&PlottingSetup::global_ratio_binning[0]); |
| 2021 |
|
ClearHisto(BpredSys); |
| 2022 |
|
|
| 2023 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 2023 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 2024 |
|
Bpred->Add(RcorrJZBem,1.0/3.); |
| 2025 |
|
Bpred->Add(LcorrJZBem,-1.0/3.); |
| 2026 |
|
Bpred->Add(RcorrJZBSBem,1.0/3.); |
| 2071 |
|
SQRT(BpredSys); |
| 2072 |
|
BpredSys->Divide(JBpred); |
| 2073 |
|
|
| 1064 |
– |
|
| 2074 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak ---- prediction changed |
| 2075 |
|
TH1F *Tpred = (TH1F*)RcorrJZBem->Clone("Bpred"); |
| 2076 |
|
Tpred->Add(LcorrJZBem,-1.0); |
| 2077 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 2077 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 2078 |
|
Tpred->Add(RcorrJZBSBem,1.0); |
| 2079 |
|
Tpred->Add(LcorrJZBSBem,-1.0); |
| 2080 |
|
Tpred->Add(RcorrJZBSBeemm,1.0); |
| 2123 |
|
Bpred->Draw("hist,same"); |
| 2124 |
|
//analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same"); |
| 2125 |
|
RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP! |
| 2126 |
< |
lm4RcorrJZBeemm->Draw("hist,same"); |
| 2126 |
> |
//lm4RcorrJZBeemm->Draw("hist,same"); |
| 2127 |
|
legBpred->AddEntry(RcorrJZBeemm,"observed","p"); |
| 2128 |
|
legBpred->AddEntry(Bpred,"predicted","l"); |
| 2129 |
|
// legBpred->AddEntry(analytical_function[1],"predicted fit","l"); |
| 2130 |
|
// legBpred->AddEntry(analytical_function[2],"stat. uncert.","l"); |
| 2131 |
< |
legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2131 |
> |
// legBpred->AddEntry(lm4RcorrJZBeemm,(allsamples.collection[allsamples.FindSample("LM")[0]].samplename).c_str(),"l"); |
| 2132 |
|
legBpred->Draw(); |
| 2133 |
|
DrawPrelim(); |
| 2134 |
|
|
| 2138 |
|
Bpred->SetLineWidth(2); |
| 2139 |
|
predcomppad->cd(); |
| 2140 |
|
predcomppad->SetLogy(1); |
| 1132 |
– |
|
| 2141 |
|
TH1F *jzbnegative = (TH1F*)LcorrJZBeemm->Clone("jzbnegative"); |
| 2142 |
|
TH1F *sidebandsemu = (TH1F*)Bpred->Clone("sidebandsemu"); |
| 2143 |
|
sidebandsemu->Add(jzbnegative,-1); |
| 2159 |
|
// lm4RcorrJZBeemm->SetLineColor(kOrange+1); |
| 2160 |
|
lm4RcorrJZBeemm->SetLineWidth(2); |
| 2161 |
|
//lm4RcorrJZBeemm->SetLineWidth(2); // paper style. overruled. |
| 2162 |
< |
lm4RcorrJZBeemm->Draw("histo,same"); |
| 2162 |
> |
// lm4RcorrJZBeemm->Draw("histo,same"); |
| 2163 |
|
DrawPrelim(); |
| 2164 |
|
TLegend *speciallegBpred = make_legend("",0.45,0.55); |
| 2165 |
|
//TLegend *speciallegBpred = make_legend("",0.35,0.55); // paper style. overruled. |
| 2166 |
|
speciallegBpred->AddEntry(RcorrJZBeemm,"Data","pl"); |
| 2167 |
|
speciallegBpred->AddEntry(Bpred,"Total background","l"); |
| 2168 |
|
speciallegBpred->AddEntry(jzbnegative,"JZB<0 (data)","f"); |
| 2169 |
< |
if(PlottingSetup::RestrictToMassPeak) speciallegBpred->AddEntry(sidebandsemu,"Sidebands/e#mu (data)","f"); |
| 2169 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) speciallegBpred->AddEntry(sidebandsemu,"Sidebands/e#mu (data)","f"); |
| 2170 |
|
else speciallegBpred->AddEntry(sidebandsemu,"e#mu (data)","f"); |
| 2171 |
|
// speciallegBpred->AddEntry(lm4RcorrJZBeemmC,"LM4","l"); |
| 2172 |
< |
speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2172 |
> |
// speciallegBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2173 |
|
speciallegBpred->Draw(); |
| 2174 |
|
save_with_ratio(JRcorrJZBeemm,JBpred,predcomppad,subdir+"Bpred_Data_____PredictionComposition",true,true,"data/pred",BpredSys); |
| 2175 |
|
|
| 2178 |
|
// THStack kostack = allsamples.DrawStack("RcorrJZBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,!is_data, luminosity,use_signal); |
| 2179 |
|
blankback->Draw(); |
| 2180 |
|
// kostack.Draw("same"); |
| 2181 |
< |
predcomposition.Draw(); |
| 2181 |
> |
predcomposition.Draw("hist"); |
| 2182 |
|
Bpred->Draw("hist,same"); |
| 2183 |
|
//analytical_function[0]->Draw("same"); analytical_function[1]->Draw("same");analytical_function[2]->Draw("same"); |
| 2184 |
|
RcorrJZBeemm->Draw("e1x0,same");//HAVE IT ON TOP! |
| 2187 |
|
CompleteSave(specialcanv,subdir+"Bpred_Data_____PredictionCompositioninMC"); |
| 2188 |
|
Bpred->SetLineWidth((int)CurrentBpredLineWidth); |
| 2189 |
|
|
| 2190 |
+ |
|
| 2191 |
+ |
//for(int i=1;i<=Bpred->GetNbinsX();i++) cout << Bpred->GetBinLowEdge(i) << ";" << Bpred->GetBinLowEdge(i)+Bpred->GetBinWidth(i) << ";;" << RcorrJZBeemm->GetBinContent(i) << ";" << LcorrJZBeemm->GetBinContent(i) << ";" << RcorrJZBem->GetBinContent(i) << ";" << LcorrJZBem->GetBinContent(i) << endl; |
| 2192 |
+ |
|
| 2193 |
|
delete speciallegBpred; |
| 2194 |
|
delete Zjetspred; |
| 2195 |
|
delete TTbarpred; |
| 2205 |
|
legBpred->AddEntry(Bpred,"MC predicted","l"); |
| 2206 |
|
if(versok) legBpred->AddEntry((TObject*)0,"",""); // Just for alignment // causes seg fault on root v5.18 |
| 2207 |
|
if(versok) legBpred->AddEntry((TObject*)0,"",""); // causes seg fault on root v5.18 |
| 2208 |
< |
if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2208 |
> |
// if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2209 |
|
legBpred->Draw(); |
| 2210 |
|
DrawMCPrelim(); |
| 2211 |
|
Bpredsaveas="Bpred_MC"; |
| 2237 |
|
|
| 2238 |
|
RcorrJZBeemmNoS->SetLineStyle(2); |
| 2239 |
|
legBpred2->AddEntry(RcorrJZBeemmNoS,"MC B","l"); |
| 2240 |
< |
legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l"); |
| 2240 |
> |
// legBpred2->AddEntry(lm4RcorrJZBeemm,"MC S","l"); |
| 2241 |
|
legBpred2->Draw(); |
| 2242 |
|
RcorrJZBeemmNoS->SetLineColor(TColor::GetColor("#61210B")); |
| 2243 |
|
RcorrJZBeemmNoS->Draw("histo,same"); |
| 2254 |
|
string ytitle("ratio"); |
| 2255 |
|
if ( use_data==1 ) ytitle = "data/pred"; |
| 2256 |
|
//save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,Bpredsaveas,true,use_data!=1,ytitle); |
| 2257 |
< |
save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,true,ytitle,BpredSys);//not extending the y range anymore up to 4 |
| 2257 |
> |
save_with_ratio(JRcorrJZBeemm,JBpred,kinpad,subdir+Bpredsaveas,true,false,ytitle,BpredSys);//not extending the y range anymore up to 4 |
| 2258 |
> |
|
| 2259 |
|
|
| 2260 |
|
|
| 2261 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 2262 |
|
// The part below is meaningless for the offpeak analysis (it's a comparison of the different estimates but there is but one estimate!) |
| 2263 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 2263 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 2264 |
|
TH1F *Bpredem = (TH1F*)LcorrJZBeemm->Clone("Bpredem"); |
| 2265 |
|
Bpredem->Add(RcorrJZBem); |
| 2266 |
|
Bpredem->Add(LcorrJZBem,-1); |
| 2315 |
|
legBpredc->AddEntry(Bpredem,"MC OFZP","l"); |
| 2316 |
|
legBpredc->AddEntry(BpredSBem,"MC OFSB","l"); |
| 2317 |
|
legBpredc->AddEntry(BpredSBeemm,"MC SFSB","l"); |
| 2318 |
< |
if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2318 |
> |
// if ( overlay_signal ) legBpred->AddEntry(lm4RcorrJZBeemm,"LM4","l"); |
| 2319 |
|
legBpredc->Draw(); |
| 2320 |
|
CompleteSave(globalcanvas,subdir+"Bpred_MCwithS_comparison"); |
| 2321 |
|
} |
| 2339 |
|
delete blankback; |
| 2340 |
|
|
| 2341 |
|
delete BpredSys; |
| 2342 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 2342 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 2343 |
|
delete RcorrJZBSBem; |
| 2344 |
|
delete LcorrJZBSBem; |
| 2345 |
|
delete RcorrJZBSBeemm; |
| 2351 |
|
delete JLcorrJZBSBeemm; |
| 2352 |
|
} |
| 2353 |
|
if(overlay_signal || use_data==2) delete lm4RcorrJZBeemm; |
| 1342 |
– |
switch_overunderflow(false); |
| 2354 |
|
} |
| 2355 |
|
|
| 2356 |
|
void do_prediction_plots(string mcjzb, string datajzb, float DataSigma, float MCSigma, bool overlay_signal ) { |
| 2357 |
+ |
switch_overunderflow(true); |
| 2358 |
|
TCanvas *globalcanvas = new TCanvas("globalcanvas","Prediction Canvas"); |
| 2359 |
|
do_prediction_plot(datajzb,globalcanvas,jzbHigh ,data,overlay_signal); |
| 2360 |
|
if ( !PlottingSetup::Approved ) { |
| 2363 |
|
} else { |
| 2364 |
|
write_info(__FUNCTION__,"You set approved to true, therefore not producing prediction/observation plots for MC with and without signal."); |
| 2365 |
|
} |
| 2366 |
< |
} |
| 2367 |
< |
|
| 1356 |
< |
void do_ratio_plot(int is_data,vector<float> binning, string jzb, TCanvas *can, float high=-9999) { |
| 1357 |
< |
bool do_data=0; |
| 1358 |
< |
bool dosignal=0; |
| 1359 |
< |
if(is_data==1) do_data=1; |
| 1360 |
< |
if(is_data==2) dosignal=1; |
| 1361 |
< |
TH1F *RcorrJZBeemm = allsamples.Draw("RcorrJZBeemm",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal); |
| 1362 |
< |
TH1F *LcorrJZBeemm = allsamples.Draw("LcorrJZBeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSSF&&cutnJets,do_data, luminosity,dosignal); |
| 1363 |
< |
TH1F *RcorrJZBem = allsamples.Draw("RcorrJZBem",jzb.c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal); |
| 1364 |
< |
TH1F *LcorrJZBem = allsamples.Draw("LcorrJZBem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", cutmass&&cutOSOF&&cutnJets,do_data, luminosity,dosignal); |
| 1365 |
< |
|
| 1366 |
< |
TH1F *RcorrJZBSBem; |
| 1367 |
< |
TH1F *LcorrJZBSBem; |
| 1368 |
< |
TH1F *RcorrJZBSBeemm; |
| 1369 |
< |
TH1F *LcorrJZBSbeemm; |
| 1370 |
< |
|
| 1371 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 1372 |
< |
RcorrJZBSBem = allsamples.Draw("RcorrJZBSbem",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal); |
| 1373 |
< |
LcorrJZBSBem = allsamples.Draw("LcorrJZBSbem",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,do_data, luminosity,dosignal); |
| 1374 |
< |
RcorrJZBSBeemm = allsamples.Draw("RcorrJZBSbeemm",jzb.c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal); |
| 1375 |
< |
LcorrJZBSbeemm = allsamples.Draw("LcorrJZBSbeemm",("-"+jzb).c_str(),binning, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,do_data, luminosity,dosignal); |
| 1376 |
< |
} |
| 1377 |
< |
|
| 1378 |
< |
|
| 1379 |
< |
|
| 1380 |
< |
|
| 1381 |
< |
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 1382 |
< |
TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred"); |
| 1383 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 1384 |
< |
Bpred->Add(RcorrJZBem,1.0/3); |
| 1385 |
< |
Bpred->Add(LcorrJZBem,-1.0/3); |
| 1386 |
< |
Bpred->Add(RcorrJZBSBem,1.0/3); |
| 1387 |
< |
Bpred->Add(LcorrJZBSBem,-1.0/3); |
| 1388 |
< |
Bpred->Add(RcorrJZBSBeemm,1.0/3); |
| 1389 |
< |
Bpred->Add(LcorrJZBSbeemm,-1.0/3); |
| 1390 |
< |
} else { |
| 1391 |
< |
Bpred->Add(RcorrJZBem,1.0); |
| 1392 |
< |
Bpred->Add(LcorrJZBem,-1.0); |
| 1393 |
< |
} |
| 1394 |
< |
|
| 1395 |
< |
can->cd(); |
| 1396 |
< |
can->SetLogy(0); |
| 1397 |
< |
Bpred->SetLineColor(kRed); |
| 1398 |
< |
Bpred->SetStats(0); |
| 1399 |
< |
if(high>0) Bpred->GetXaxis()->SetRangeUser(0,high); |
| 1400 |
< |
TH1F *JZBratioforfitting=(TH1F*)RcorrJZBeemm->Clone("JZBratioforfitting"); |
| 1401 |
< |
JZBratioforfitting->Divide(Bpred); |
| 1402 |
< |
TGraphAsymmErrors *JZBratio = histRatio(RcorrJZBeemm,Bpred,is_data,binning,false); |
| 1403 |
< |
|
| 1404 |
< |
|
| 1405 |
< |
JZBratio->SetTitle(""); |
| 1406 |
< |
JZBratio->GetYaxis()->SetRangeUser(0.0,9.0); |
| 1407 |
< |
// if(is_data==1) JZBratio->GetXaxis()->SetRangeUser(0,jzbHigh ); |
| 1408 |
< |
|
| 1409 |
< |
TF1 *pol0 = new TF1("pol0","[0]",0,1000); |
| 1410 |
< |
TF1 *pol0d = new TF1("pol0","[0]",0,1000); |
| 1411 |
< |
// straightline_fit->SetParameter(0,1); |
| 1412 |
< |
JZBratioforfitting->Fit(pol0,"Q0R","",0,30); |
| 1413 |
< |
pol0d->SetParameter(0,pol0->GetParameter(0)); |
| 1414 |
< |
|
| 1415 |
< |
JZBratio->GetYaxis()->SetTitle("Observed/Predicted"); |
| 1416 |
< |
JZBratio->GetXaxis()->SetTitle("JZB [GeV]"); |
| 1417 |
< |
if ( high>0 ) JZBratio->GetXaxis()->SetRangeUser(0.0,high); |
| 1418 |
< |
JZBratio->GetYaxis()->SetNdivisions(519); |
| 1419 |
< |
JZBratio->GetYaxis()->SetRangeUser(0.0,4.0); |
| 1420 |
< |
JZBratio->GetYaxis()->CenterTitle(); |
| 1421 |
< |
JZBratio->GetXaxis()->CenterTitle(); |
| 1422 |
< |
JZBratio->SetMarkerSize(DataMarkerSize); |
| 1423 |
< |
JZBratio->Draw("AP"); |
| 1424 |
< |
/////---------------------------- |
| 1425 |
< |
TPaveText *writeresult = new TPaveText(0.15,0.78,0.49,0.91,"blNDC"); |
| 1426 |
< |
writeresult->SetFillStyle(4000); |
| 1427 |
< |
writeresult->SetFillColor(kWhite); |
| 1428 |
< |
writeresult->SetTextFont(42); |
| 1429 |
< |
writeresult->SetTextSize(0.03); |
| 1430 |
< |
writeresult->SetTextAlign(12); |
| 1431 |
< |
ostringstream jzb_agreement_data_text; |
| 1432 |
< |
jzb_agreement_data_text<< setprecision(2) << "mean =" << pol0->GetParameter(0) << " #pm " << setprecision(1) << pol0->GetParError(0); |
| 1433 |
< |
if(is_data==1) fitresultconstdata=pol0->GetParameter(0);// data |
| 1434 |
< |
if(is_data==0) fitresultconstmc=pol0->GetParameter(0); // monte carlo, no signal |
| 1435 |
< |
/* if(is_data) writeresult->AddText("Data closure test"); |
| 1436 |
< |
else writeresult->AddText("MC closure test"); |
| 1437 |
< |
*/ |
| 1438 |
< |
writeresult->AddText(jzb_agreement_data_text.str().c_str()); |
| 1439 |
< |
// writeresult->Draw("same"); |
| 1440 |
< |
// pol0d->Draw("same"); |
| 1441 |
< |
TF1 *topline = new TF1("","1.5",0,1000); |
| 1442 |
< |
TF1 *bottomline = new TF1("","0.5",0,1000); |
| 1443 |
< |
topline->SetLineColor(kBlue); |
| 1444 |
< |
topline->SetLineStyle(2); |
| 1445 |
< |
bottomline->SetLineColor(kBlue); |
| 1446 |
< |
bottomline->SetLineStyle(2); |
| 1447 |
< |
// topline->Draw("same"); |
| 1448 |
< |
// bottomline->Draw("same"); |
| 1449 |
< |
TF1 *oneline = new TF1("","1.0",0,1000); |
| 1450 |
< |
oneline->SetLineColor(kBlue); |
| 1451 |
< |
oneline->SetLineStyle(1); |
| 1452 |
< |
oneline->Draw("same"); |
| 1453 |
< |
TLegend *phony_leg = make_legend("ratio",0.6,0.55,false);//this line is just to have the default CMS Preliminary (...) on the canvas as well. |
| 1454 |
< |
if(is_data==1) DrawPrelim(); |
| 1455 |
< |
else DrawMCPrelim(); |
| 1456 |
< |
TLegend *leg = new TLegend(0.55,0.75,0.89,0.89); |
| 1457 |
< |
leg->SetTextFont(42); |
| 1458 |
< |
leg->SetTextSize(0.04); |
| 1459 |
< |
// if(is_data==1) leg->SetHeader("Ratio (data)"); |
| 1460 |
< |
// else leg->SetHeader("Ratio (MC)"); |
| 1461 |
< |
|
| 1462 |
< |
TString MCtitle("MC "); |
| 1463 |
< |
if (is_data==1) MCtitle = ""; |
| 1464 |
< |
|
| 1465 |
< |
leg->SetFillStyle(4000); |
| 1466 |
< |
leg->SetFillColor(kWhite); |
| 1467 |
< |
leg->SetTextFont(42); |
| 1468 |
< |
// leg->AddEntry(topline,"+20\% sys envelope","l"); |
| 1469 |
< |
leg->AddEntry(JZBratio,MCtitle+"obs / "+MCtitle+"pred","p"); |
| 1470 |
< |
leg->AddEntry(oneline,"ratio = 1","l"); |
| 1471 |
< |
// leg->AddEntry(pol0d,"fit in [0,30] GeV","l"); |
| 1472 |
< |
// leg->AddEntry(bottomline,"#pm50% envelope","l"); |
| 1473 |
< |
|
| 1474 |
< |
|
| 1475 |
< |
//leg->Draw("same"); // no longer drawing legend |
| 1476 |
< |
|
| 1477 |
< |
if(is_data==1) CompleteSave(can, "jzb_ratio_data"); |
| 1478 |
< |
if(is_data==0) CompleteSave(can, "jzb_ratio_mc"); |
| 1479 |
< |
if(is_data==2) CompleteSave(can, "jzb_ratio_mc_BandS");//special case, MC with signal! |
| 1480 |
< |
|
| 1481 |
< |
delete RcorrJZBeemm; |
| 1482 |
< |
delete LcorrJZBeemm; |
| 1483 |
< |
delete RcorrJZBem; |
| 1484 |
< |
delete LcorrJZBem; |
| 1485 |
< |
|
| 1486 |
< |
delete RcorrJZBSBem; |
| 1487 |
< |
delete LcorrJZBSBem; |
| 1488 |
< |
delete RcorrJZBSBeemm; |
| 1489 |
< |
delete LcorrJZBSbeemm; |
| 1490 |
< |
} |
| 1491 |
< |
|
| 1492 |
< |
void do_ratio_plots(string mcjzb,string datajzb,vector<float> ratio_binning) { |
| 1493 |
< |
TCanvas *globalc = new TCanvas("globalc","Ratio Plot Canvas"); |
| 1494 |
< |
globalc->SetLogy(0); |
| 1495 |
< |
|
| 1496 |
< |
do_ratio_plot(mc,ratio_binning,mcjzb,globalc, jzbHigh ); |
| 1497 |
< |
do_ratio_plot(data,ratio_binning,datajzb,globalc, jzbHigh ); |
| 1498 |
< |
do_ratio_plot(mcwithsignal,ratio_binning,mcjzb,globalc, jzbHigh ); |
| 2366 |
> |
delete globalcanvas; |
| 2367 |
> |
switch_overunderflow(false); |
| 2368 |
|
} |
| 2369 |
|
|
| 2370 |
|
string give_jzb_expression(float peak, int type) { |
| 2552 |
|
datahisto->SetMinimum(0.1); |
| 2553 |
|
datahisto->SetMarkerSize(DataMarkerSize); |
| 2554 |
|
datahisto->Draw("e1"); |
| 2555 |
< |
mcstack.Draw("same"); |
| 2555 |
> |
mcstack.Draw("histo,same"); |
| 2556 |
|
datahisto->Draw("same,e1"); |
| 2557 |
|
|
| 2558 |
|
TLegend *leg; |
| 2559 |
< |
if(is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("OFZP",datahisto); |
| 2560 |
< |
else if(!is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("SFZP",datahisto); |
| 2561 |
< |
else if( is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("OFSB",datahisto); |
| 2562 |
< |
else if(!is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("SFSB",datahisto); |
| 2563 |
< |
else { |
| 2564 |
< |
std::cerr << "Unable to decode cut: " << cut.GetTitle() << std::endl; |
| 2565 |
< |
exit(-1); |
| 2566 |
< |
} |
| 2559 |
> |
if (!PlottingSetup::RestrictToMassPeak||!PlottingSetup::UseSidebandsForcJZB) { |
| 2560 |
> |
if(is_OF(cut)) leg = allsamples.allbglegend("Opposite flavor",datahisto); |
| 2561 |
> |
else leg = allsamples.allbglegend("Same flavor",datahisto); |
| 2562 |
> |
} else { |
| 2563 |
> |
if(is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("OFZP",datahisto); |
| 2564 |
> |
else if(!is_OF(cut) && is_ZP(cut)) leg = allsamples.allbglegend("SFZP",datahisto); |
| 2565 |
> |
else if( is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("OFSB",datahisto); |
| 2566 |
> |
else if(!is_OF(cut) && !is_ZP(cut)) leg = allsamples.allbglegend("SFSB",datahisto); |
| 2567 |
> |
else { |
| 2568 |
> |
std::cerr << "Unable to decode cut: " << cut.GetTitle() << std::endl; |
| 2569 |
> |
exit(-1); |
| 2570 |
> |
} |
| 2571 |
> |
} |
| 2572 |
|
leg->Draw(); |
| 2573 |
|
string write_cut = decipher_cut(cut,""); |
| 2574 |
|
TText *writeline1 = write_cut_on_canvas(write_cut.c_str()); |
| 2583 |
|
datahisto->SetMinimum(0.1); |
| 2584 |
|
datahisto->SetMarkerSize(DataMarkerSize); |
| 2585 |
|
datahisto->Draw("e1"); |
| 2586 |
< |
mcstack.Draw("same"); |
| 2586 |
> |
mcstack.Draw("histo,same"); |
| 2587 |
|
datahisto->Draw("same,e1"); |
| 2588 |
|
leg->SetHeader(""); |
| 2589 |
|
leg->Draw(); |
| 2619 |
|
draw_pure_jzb_histo(cutOSSF&&cutnJets&&cutmass&&"id1==1",datajzb,mcjzb,"Dibosons/mm/jzb_OS_SFZP",dican,binning); |
| 2620 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",datajzb,mcjzb,"Dibosons/ee/jzb_OS_OFZP",dican,binning); |
| 2621 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,mcjzb,"Dibosons/mm/jzb_OS_OFZP",dican,binning); |
| 2622 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_SFSB",dican,binning); |
| 2623 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_OFSB",dican,binning); |
| 2622 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_SFSB",dican,binning); |
| 2623 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_OFSB",dican,binning); |
| 2624 |
|
|
| 2625 |
|
draw_pure_jzb_histo(cutOSSF&&cutnJets&&cutmass,datajzb,mcjzb,"Dibosons/jzb_OS_SFZP_coarse",dican,coarse_binning); |
| 2626 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass,datajzb,mcjzb,"Dibosons/jzb_OS_OFZP_coarse",dican,coarse_binning); |
| 2627 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_SFSB_coarse",dican,coarse_binning); |
| 2628 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_OFSB_coarse",dican,coarse_binning); |
| 2627 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_SFSB_coarse",dican,coarse_binning); |
| 2628 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"Dibosons/jzb_OS_OFSB_coarse",dican,coarse_binning); |
| 2629 |
|
|
| 2630 |
|
delete dican; |
| 2631 |
|
} |
| 2632 |
|
|
| 2633 |
|
|
| 2634 |
|
void diboson_plots(string mcjzb, string datajzb,vector<float> ratio_binning) { |
| 2635 |
+ |
switch_overunderflow(true); |
| 2636 |
|
vector<int> SamplesToBeModified = allsamples.FindSampleBySampleName("WW/WZ/ZZ"); |
| 2637 |
|
|
| 2638 |
|
if(SamplesToBeModified.size()==0 || SamplesToBeModified[0]==-1) { |
| 2672 |
|
dout << " Reset xs for sample " << (allsamples.collection)[SamplesToBeModified[i]].samplename << " from " << Upxs << " to " << (allsamples.collection)[SamplesToBeModified[i]].xs << " (by a factor of " << stretchfactor << ") and reset the correct name (from " << Upname << ")" << endl; |
| 2673 |
|
|
| 2674 |
|
} |
| 2675 |
+ |
// switch_overunderflow(false); |
| 2676 |
+ |
} |
| 2677 |
+ |
|
| 2678 |
+ |
|
| 2679 |
+ |
void draw_normalized_data_vs_data_histo(TCut cut1, TCut cut2, string variable, string legentry1, string legentry2, string savename, TCanvas *can,vector<float> binning) { |
| 2680 |
+ |
TPad *jzbpad = new TPad("jzbpad","jzbpad",0,0,1,1); |
| 2681 |
+ |
jzbpad->cd(); |
| 2682 |
+ |
jzbpad->SetLogy(1); |
| 2683 |
+ |
string xlabel="JZB [GeV]"; |
| 2684 |
+ |
|
| 2685 |
+ |
TH1F *datahisto1 = allsamples.Draw("datahisto1",variable,binning, xlabel, "events",cut1,data,luminosity); |
| 2686 |
+ |
datahisto1->SetLineColor(kRed); |
| 2687 |
+ |
datahisto1->SetMarkerColor(kRed); |
| 2688 |
+ |
TH1F *datahisto2 = allsamples.Draw("datahisto2",variable,binning, xlabel, "events",cut2,data,luminosity); |
| 2689 |
+ |
datahisto2->SetLineColor(kBlue); |
| 2690 |
+ |
datahisto2->SetMarkerColor(kBlue); |
| 2691 |
+ |
|
| 2692 |
+ |
datahisto2->SetMarkerSize(DataMarkerSize); |
| 2693 |
+ |
datahisto1->DrawNormalized("e1"); |
| 2694 |
+ |
datahisto2->DrawNormalized("histo,same"); |
| 2695 |
+ |
datahisto1->DrawNormalized("same,e1"); |
| 2696 |
+ |
|
| 2697 |
+ |
TLegend *leg = make_legend(); |
| 2698 |
+ |
leg->AddEntry(datahisto1,legentry1.c_str()); |
| 2699 |
+ |
leg->AddEntry(datahisto2,legentry2.c_str()); |
| 2700 |
+ |
leg->Draw(); |
| 2701 |
+ |
|
| 2702 |
+ |
save_with_ratio(datahisto1,datahisto2,jzbpad->cd(),("jzb/"+savename)); |
| 2703 |
|
|
| 2704 |
+ |
datahisto1->Delete(); |
| 2705 |
+ |
datahisto2->Delete(); |
| 2706 |
|
} |
| 2707 |
+ |
|
| 2708 |
+ |
|
| 2709 |
|
void jzb_plots(string mcjzb, string datajzb,vector<float> ratio_binning) { |
| 2710 |
+ |
switch_overunderflow(true); |
| 2711 |
|
TCanvas *can = new TCanvas("can","JZB Plots Canvas"); |
| 2712 |
|
float max=jzbHigh ; |
| 2713 |
|
float min=-120; |
| 2728 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",datajzb,mcjzb,"ee/jzb_OS_OFZP",can,binning); |
| 2729 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,mcjzb,"mm/jzb_OS_OFZP",can,binning); |
| 2730 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 2731 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_SFSB",can,binning); |
| 2732 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB",can,binning); |
| 2731 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_SFSB",can,binning); |
| 2732 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB",can,binning); |
| 2733 |
> |
draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm",can,binning); |
| 2734 |
> |
draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm_coarse",can,coarse_binning); |
| 2735 |
> |
draw_normalized_data_vs_data_histo(cutOSOF&&cutnJets&&cutmass&&"id1==0",cutOSOF&&cutnJets&&cutmass&&"id1==1",datajzb,"ee","mm","jzb_ee_vs_mm_coarsest",can,coarsest_binning); |
| 2736 |
> |
|
| 2737 |
|
} |
| 2738 |
|
|
| 2739 |
|
draw_pure_jzb_histo(cutOSSF&&cutnJets&&cutmass,datajzb,mcjzb,"jzb_OS_SFZP_coarse",can,coarse_binning); |
| 2740 |
|
if ( !PlottingSetup::Approved ) { |
| 2741 |
|
draw_pure_jzb_histo(cutOSOF&&cutnJets&&cutmass,datajzb,mcjzb,"jzb_OS_OFZP_coarse",can,coarse_binning); |
| 2742 |
|
flag_this_change(__FUNCTION__,__LINE__,true);//PlottingSetup::RestrictToMassPeak |
| 2743 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_SFSB_coarse",can,coarse_binning); |
| 2744 |
< |
if(PlottingSetup::RestrictToMassPeak) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB_coarse",can,coarse_binning); |
| 2743 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSSF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_SFSB_coarse",can,coarse_binning); |
| 2744 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) draw_pure_jzb_histo(cutOSOF&&cutnJets&&sidebandcut,datajzb,mcjzb,"jzb_OS_OFSB_coarse",can,coarse_binning); |
| 2745 |
|
} |
| 2746 |
|
delete can; |
| 2747 |
+ |
switch_overunderflow(false); |
| 2748 |
|
} |
| 2749 |
|
|
| 2750 |
|
|
| 2794 |
|
TCanvas *yica = new TCanvas("yica","yield canvas"); |
| 2795 |
|
|
| 2796 |
|
int nRegions=4; |
| 2797 |
< |
if(!PlottingSetup::RestrictToMassPeak) nRegions=2; |
| 2797 |
> |
if(!PlottingSetup::RestrictToMassPeak||!PlottingSetup::UseSidebandsForcJZB) nRegions=2; |
| 2798 |
|
string tsRegions[] = {"SFZP","OFZP","SFSB","OFSB"}; |
| 2799 |
|
string posneg[] = {"pos","neg"}; |
| 2800 |
|
TCut tkRegions[] = {cutOSSF&&cutnJets&&cutmass,cutOSOF&&cutnJets&&cutmass,cutOSSF&&cutnJets&&sidebandcut,cutOSOF&&cutnJets&&sidebandcut}; |
| 3085 |
|
TCanvas *cancorr = new TCanvas("cancorr","Canvas for Response Correction"); |
| 3086 |
|
cancorr->SetLogz(); |
| 3087 |
|
cancorr->SetRightMargin(0.13); |
| 3088 |
< |
TCut zptforresponsepresentation("pt<600"&&cutmass&&cutOSSF&&"((sumJetPt[1]/pt)<5.0)"&&SpecialCut&&passtrig); |
| 3088 |
> |
TCut zptforresponsepresentation(Restrmasscut&&SpecialCut&&passtrig); |
| 3089 |
> |
|
| 3090 |
|
if(PlottingSetup::DoBTag) zptforresponsepresentation=zptforresponsepresentation&&PlottingSetup::bTagRequirement; |
| 3091 |
|
TH2F *niceresponseplotd = new TH2F("niceresponseplotd","",100,0,600,100,0,5); |
| 3092 |
+ |
niceresponseplotd->Sumw2(); |
| 3093 |
+ |
TH2F* emuResponse = (TH2F*)niceresponseplotd->Clone("emuResponse"); |
| 3094 |
|
vector<int> SampleIndices=allsamples.FindSample(FindKeyword); |
| 3095 |
|
for(int iSample=0;iSample<(int)SampleIndices.size();iSample++) { |
| 3096 |
|
if((allsamples.collection)[SampleIndices[iSample]].is_data && !dodata) continue; |
| 3097 |
|
if((allsamples.collection)[SampleIndices[iSample]].is_data ==false && dodata) continue; |
| 3098 |
|
|
| 3099 |
|
dout << " Response correction : Using sample " << (allsamples.collection)[SampleIndices[iSample]].filename << " for " << FindKeyword << endl; |
| 3100 |
< |
(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+niceresponseplotd",zptforresponsepresentation*cutWeight); |
| 3100 |
> |
(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+niceresponseplotd",(zptforresponsepresentation&&cutOSSF)*cutWeight); |
| 3101 |
> |
(allsamples.collection)[SampleIndices[iSample]].events->Draw("sumJetPt[1]/pt:pt>>+emuResponse",(zptforresponsepresentation&&cutOSOF)*cutWeight); |
| 3102 |
|
} |
| 3103 |
< |
|
| 3103 |
> |
niceresponseplotd->Add(emuResponse,-1); |
| 3104 |
> |
|
| 3105 |
|
niceresponseplotd->SetStats(0); |
| 3106 |
|
niceresponseplotd->GetXaxis()->SetTitle("Z p_{T} [GeV]"); |
| 3107 |
|
niceresponseplotd->GetYaxis()->SetTitle("Response"); |
| 3109 |
|
niceresponseplotd->GetYaxis()->CenterTitle(); |
| 3110 |
|
niceresponseplotd->Draw("COLZ"); |
| 3111 |
|
TProfile * profd = (TProfile*)niceresponseplotd->ProfileX(); |
| 3112 |
+ |
profd->Rebin(2); |
| 3113 |
|
profd->SetMarkerSize(DataMarkerSize); |
| 3114 |
|
profd->Fit("pol0","","same,e1",100,400); |
| 3115 |
|
DrawPrelim(); |
| 3129 |
|
CompleteSave(cancorr,"ResponseCorrection/Response_Correction_Illustration_New_"+SaveAs); |
| 3130 |
|
delete cancorr; |
| 3131 |
|
delete niceresponseplotd; |
| 3132 |
+ |
delete profd; |
| 3133 |
|
return correction; |
| 3134 |
|
} |
| 3135 |
|
|
| 3242 |
|
configfile<<"\n\n"; |
| 3243 |
|
configfile<<"luminosity="<<luminosity<<";\n"; |
| 3244 |
|
configfile<<"RestrictToMassPeak="<<RestrictToMassPeak<<";//defines the type of analysis we're running\n"; |
| 3245 |
+ |
configfile<<"UseSidebandsForcJZB="<<UseSidebandsForcJZB<<";//tells us whether to use the sidebands or not\n"; |
| 3246 |
|
|
| 3247 |
|
configfile<<"\n\ncout << \"Configuration successfully loaded!\" << endl; \n \n } \n \n"; |
| 3248 |
|
|
| 3285 |
|
//in the case of the on peak analysis, we compare the 3 control regions to the real value |
| 3286 |
|
//in the case of the OFF peak analysis, we compare our control region to the real value |
| 3287 |
|
TCut weightbackup=cutWeight; |
| 3288 |
< |
// cutWeight="1.0"; |
| 3288 |
> |
switch_overunderflow(true); |
| 3289 |
> |
|
| 3290 |
> |
bool doPURW=false; |
| 3291 |
> |
|
| 3292 |
> |
|
| 3293 |
> |
if(!doPURW) { |
| 3294 |
> |
dout << "Not doing PU reweighting for ttbar closure test" << endl; |
| 3295 |
> |
cutWeight="1.0"; |
| 3296 |
> |
// Do it without PU re-weighting |
| 3297 |
> |
float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0; |
| 3298 |
> |
stringstream resultsNoPU; |
| 3299 |
> |
stringstream noPUdatajzb; |
| 3300 |
> |
stringstream noPUmcjzb; |
| 3301 |
> |
|
| 3302 |
> |
stringstream mcjzbnoPU; |
| 3303 |
> |
find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,resultsNoPU,true,noPUdatajzb,noPUmcjzb); |
| 3304 |
> |
mcjzb = noPUmcjzb.str(); |
| 3305 |
> |
} |
| 3306 |
> |
|
| 3307 |
> |
|
| 3308 |
|
float simulatedlumi = luminosity; //in pb please - adjust to your likings |
| 3309 |
|
|
| 3310 |
< |
TH1F *TZem = systsamples.Draw("TZem",mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3311 |
< |
TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3310 |
> |
TH1F *TZem = systsamples.Draw("TZem", mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3311 |
> |
TH1F *nTZem = systsamples.Draw("nTZem","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3312 |
|
TH1F *TSem; |
| 3313 |
|
TH1F *nTSem; |
| 3314 |
< |
TH1F *TZeemm = systsamples.Draw("TZeemm",mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3315 |
< |
TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3314 |
> |
TH1F *TZeemm = systsamples.Draw("TZeemm", mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3315 |
> |
TH1F *nTZeemm = systsamples.Draw("nTZeemm","-"+mcjzb,binning,"JZB [GeV]","events",cutmass&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3316 |
|
TH1F *TSeemm; |
| 3317 |
|
TH1F *nTSeemm; |
| 3318 |
|
|
| 3319 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3320 |
< |
TSem = systsamples.Draw("TSem",mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3321 |
< |
nTSem = systsamples.Draw("nTSem","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3322 |
< |
TSeemm = systsamples.Draw("TSeemm",mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3323 |
< |
nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTJets")); |
| 3319 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3320 |
> |
TSem = systsamples.Draw("TSem", mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3321 |
> |
nTSem = systsamples.Draw("nTSem", "-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSOF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3322 |
> |
TSeemm = systsamples.Draw("TSeemm", mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3323 |
> |
nTSeemm = systsamples.Draw("nTSeemm","-"+mcjzb,binning,"JZB [GeV]","events",sidebandcut&&cutOSSF&&cutnJets,mc,simulatedlumi,systsamples.FindSample("TTT")); |
| 3324 |
|
} |
| 3325 |
|
|
| 3326 |
|
TCanvas *tcan = new TCanvas("tcan","tcan"); |
| 2387 |
– |
|
| 3327 |
|
tcan->SetLogy(1); |
| 3328 |
|
|
| 3329 |
|
TZeemm->SetLineColor(kBlack); |
| 3332 |
|
TZem->SetMarkerColor(kRed); |
| 3333 |
|
|
| 3334 |
|
|
| 3335 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3335 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3336 |
|
TSem->SetLineColor(TColor::GetColor("#00A616")); |
| 3337 |
|
TSeemm->SetLineColor(kMagenta+2); |
| 3338 |
|
TSem->SetMarkerColor(TColor::GetColor("#00A616")); |
| 3346 |
|
TZeemm->GetXaxis()->SetRangeUser(-100,binning[binning.size()-1]); |
| 3347 |
|
histos.push_back(TZem); |
| 3348 |
|
histos.push_back(TZeemm); |
| 3349 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3349 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3350 |
|
TSeemm->GetXaxis()->SetRangeUser(-100,binning[binning.size()-1]); |
| 3351 |
|
TSem->GetXaxis()->SetRangeUser(-100,binning[binning.size()-1]); |
| 3352 |
|
histos.push_back(TSem); |
| 3353 |
|
histos.push_back(TSeemm); |
| 3354 |
|
} |
| 3355 |
< |
draw_all_ttbar_histos(tcan,histos,"histo",0.04); |
| 3355 |
> |
draw_all_ttbar_histos(tcan,histos,"histo",8); |
| 3356 |
|
|
| 3357 |
|
TLegend *leg = make_legend("MC t#bar{t}",0.6,0.65,false); |
| 3358 |
|
leg->AddEntry(TZeemm,"SFZP","l"); |
| 3359 |
|
leg->AddEntry(TZem,"OFZP","l"); |
| 3360 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3360 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3361 |
|
leg->AddEntry(TSeemm,"SFSB","l"); |
| 3362 |
|
leg->AddEntry(TSem,"OFSB","l"); |
| 3363 |
|
} |
| 3364 |
|
leg->Draw("same"); |
| 3365 |
|
DrawMCPrelim(simulatedlumi); |
| 3366 |
|
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison"); |
| 3367 |
+ |
|
| 3368 |
|
TH1F *TZemcopy = (TH1F*)TZem->Clone("TZemcopy"); |
| 3369 |
|
TH1F *TZeemmcopy = (TH1F*)TZeemm->Clone("TZeemmcopy"); |
| 3370 |
|
TH1F *TSeemmcopy; |
| 3371 |
|
TH1F *TSemcopy; |
| 3372 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3372 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3373 |
|
TSeemmcopy = (TH1F*)TSeemm->Clone("TSeemmcopy"); |
| 3374 |
|
TSemcopy = (TH1F*)TSem->Clone("TSemcopy"); |
| 3375 |
|
} |
| 3376 |
|
|
| 3377 |
|
TZem->Divide(TZeemm); |
| 3378 |
|
TZem->SetMarkerStyle(21); |
| 3379 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3379 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3380 |
|
TSem->Divide(TZeemm); |
| 3381 |
|
TSeemm->Divide(TZeemm); |
| 3382 |
|
TSem->SetMarkerStyle(24); |
| 3383 |
|
TSeemm->SetMarkerStyle(32); |
| 3384 |
< |
} |
| 3384 |
> |
} |
| 3385 |
|
|
| 3386 |
|
tcan->SetLogy(0); |
| 3387 |
|
TZem->GetYaxis()->SetRangeUser(0,2.5); |
| 3388 |
|
TZem->GetYaxis()->SetTitle("ratio"); |
| 3389 |
|
TZem->Draw(); |
| 3390 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3390 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3391 |
|
TSem->Draw("same"); |
| 3392 |
|
TSeemm->Draw("same"); |
| 3393 |
|
} |
| 3417 |
|
|
| 3418 |
|
TLegend *leg2 = make_legend("MC t#bar{t}",0.55,0.75,false); |
| 3419 |
|
leg2->AddEntry(TZem,"OFZP / SFZP","ple"); |
| 3420 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3420 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3421 |
|
leg2->AddEntry(TSeemm,"SFSB / SFZP","ple"); |
| 3422 |
|
leg2->AddEntry(TSem,"OFSB / SFZP","ple"); |
| 3423 |
|
} |
| 3430 |
|
DrawMCPrelim(simulatedlumi); |
| 3431 |
|
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_shape_comparison_ratio"); |
| 3432 |
|
|
| 3433 |
+ |
|
| 3434 |
+ |
if (0) { // Turn this off: we don't need this |
| 3435 |
|
|
| 3436 |
< |
///-------------- second part: only look at the quantity we actually care about! |
| 3437 |
< |
TH1F *leftsfzp = (TH1F*)nTZeemm->Clone("leftsfzp"); |
| 3438 |
< |
TH1F *rightsfzp = (TH1F*)TZeemmcopy->Clone("rightsfzp"); |
| 3439 |
< |
rightsfzp->Add(leftsfzp,-1); |
| 3440 |
< |
TH1F *leftofzp = (TH1F*)nTZem->Clone("leftofzp"); |
| 3441 |
< |
TH1F *rightofzp = (TH1F*)TZemcopy->Clone("rightofzp"); |
| 3442 |
< |
rightofzp->Add(leftofzp,-1); |
| 3443 |
< |
TH1F *leftofsb; |
| 3444 |
< |
TH1F *rightofsb; |
| 3445 |
< |
TH1F *leftsfsb; |
| 3446 |
< |
TH1F *rightsfsb; |
| 3447 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3448 |
< |
leftofsb = (TH1F*)nTSem->Clone("leftofsb"); |
| 3449 |
< |
rightofsb = (TH1F*)TSemcopy->Clone("rightofsb"); |
| 3450 |
< |
rightofsb->Add(leftofsb,-1); |
| 3451 |
< |
leftsfsb = (TH1F*)nTSeemm->Clone("leftsfsb"); |
| 3452 |
< |
rightsfsb = (TH1F*)TSeemmcopy->Clone("rightsfsb"); |
| 3453 |
< |
rightsfsb->Add(leftsfsb,-1); |
| 3454 |
< |
} |
| 3436 |
> |
///-------------- second part: only look at the quantity we actually care about! |
| 3437 |
> |
TH1F *leftsfzp = (TH1F*)nTZeemm->Clone("leftsfzp"); |
| 3438 |
> |
TH1F *rightsfzp = (TH1F*)TZeemmcopy->Clone("rightsfzp"); |
| 3439 |
> |
rightsfzp->Add(leftsfzp,-1); |
| 3440 |
> |
TH1F *leftofzp = (TH1F*)nTZem->Clone("leftofzp"); |
| 3441 |
> |
TH1F *rightofzp = (TH1F*)TZemcopy->Clone("rightofzp"); |
| 3442 |
> |
rightofzp->Add(leftofzp,-1); |
| 3443 |
> |
TH1F *leftofsb; |
| 3444 |
> |
TH1F *rightofsb; |
| 3445 |
> |
TH1F *leftsfsb; |
| 3446 |
> |
TH1F *rightsfsb; |
| 3447 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3448 |
> |
leftofsb = (TH1F*)nTSem->Clone("leftofsb"); |
| 3449 |
> |
rightofsb = (TH1F*)TSemcopy->Clone("rightofsb"); |
| 3450 |
> |
rightofsb->Add(leftofsb,-1); |
| 3451 |
> |
leftsfsb = (TH1F*)nTSeemm->Clone("leftsfsb"); |
| 3452 |
> |
rightsfsb = (TH1F*)TSeemmcopy->Clone("rightsfsb"); |
| 3453 |
> |
rightsfsb->Add(leftsfsb,-1); |
| 3454 |
> |
} |
| 3455 |
|
|
| 3456 |
< |
tcan->SetLogy(1); |
| 3457 |
< |
rightsfzp->GetXaxis()->SetRangeUser(0,binning[binning.size()-1]); |
| 3458 |
< |
rightsfzp->GetYaxis()->SetTitle("#deltaJZB / 25 GeV"); |
| 3459 |
< |
rightsfzp->Draw("histo"); |
| 3460 |
< |
rightofzp->Draw("histo,same"); |
| 3461 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3462 |
< |
rightofsb->Draw("histo,same"); |
| 3463 |
< |
rightsfsb->Draw("histo,same"); |
| 3464 |
< |
} |
| 3465 |
< |
DrawMCPrelim(simulatedlumi); |
| 3456 |
> |
tcan->SetLogy(1); |
| 3457 |
> |
rightsfzp->GetXaxis()->SetRangeUser(0,binning[binning.size()-1]); |
| 3458 |
> |
rightsfzp->GetYaxis()->SetTitle("#deltaJZB / 25 GeV"); |
| 3459 |
> |
rightsfzp->Draw("histo"); |
| 3460 |
> |
rightofzp->Draw("histo,same"); |
| 3461 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3462 |
> |
rightofsb->Draw("histo,same"); |
| 3463 |
> |
rightsfsb->Draw("histo,same"); |
| 3464 |
> |
} |
| 3465 |
> |
DrawMCPrelim(simulatedlumi); |
| 3466 |
|
|
| 3467 |
< |
TLegend *legA = make_legend("MC t#bar{t}",0.6,0.65,false); |
| 3468 |
< |
legA->AddEntry(rightsfzp,"SFZP","l"); |
| 3469 |
< |
legA->AddEntry(rightofzp,"OFZP","l"); |
| 3470 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3471 |
< |
legA->AddEntry(rightofsb,"SFSB","l"); |
| 3472 |
< |
legA->AddEntry(rightsfsb,"OFSB","l"); |
| 3473 |
< |
} |
| 3474 |
< |
legA->Draw(); |
| 3467 |
> |
TLegend *legA = make_legend("MC t#bar{t}",0.6,0.65,false); |
| 3468 |
> |
legA->AddEntry(rightsfzp,"SFZP","l"); |
| 3469 |
> |
legA->AddEntry(rightofzp,"OFZP","l"); |
| 3470 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3471 |
> |
legA->AddEntry(rightofsb,"SFSB","l"); |
| 3472 |
> |
legA->AddEntry(rightsfsb,"OFSB","l"); |
| 3473 |
> |
} |
| 3474 |
> |
legA->Draw(); |
| 3475 |
|
|
| 3476 |
< |
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison"); |
| 3476 |
> |
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison"); |
| 3477 |
|
|
| 3478 |
< |
tcan->SetLogy(0); |
| 3479 |
< |
rightofzp->Divide(rightsfzp); |
| 3480 |
< |
rightofzp->GetXaxis()->SetRangeUser(0.0,binning[binning.size()-1]); |
| 3481 |
< |
rightofzp->GetYaxis()->SetRangeUser(0.0,2.5); |
| 3482 |
< |
rightofzp->GetYaxis()->SetTitle("#deltaJZB ratio"); |
| 3483 |
< |
rightofzp->Draw(); |
| 3484 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3485 |
< |
rightofsb->Divide(rightsfzp); |
| 3486 |
< |
rightsfsb->Divide(rightsfzp); |
| 3487 |
< |
rightofsb->Draw("same"); |
| 3488 |
< |
rightsfsb->Draw("same"); |
| 3489 |
< |
} |
| 3490 |
< |
|
| 3491 |
< |
TLine *top2 = new TLine(0.0,1.0+linepos,binning[binning.size()-1],1.0+linepos); |
| 3492 |
< |
TLine *center2 = new TLine(0.0,1.0,binning[binning.size()-1],1.0); |
| 3493 |
< |
TLine *bottom2 = new TLine(0.0,1.0-linepos,binning[binning.size()-1],1.0-linepos); |
| 3494 |
< |
|
| 3495 |
< |
top2->SetLineColor(kBlue);top2->SetLineStyle(2); |
| 3496 |
< |
bottom2->SetLineColor(kBlue);bottom2->SetLineStyle(2); |
| 3497 |
< |
center2->SetLineColor(kBlue); |
| 3498 |
< |
|
| 3499 |
< |
top2->Draw("same"); |
| 3500 |
< |
center2->Draw("same"); |
| 3501 |
< |
bottom2->Draw("same"); |
| 3502 |
< |
|
| 3503 |
< |
rightofzp->SetMarkerStyle(21); |
| 3504 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3505 |
< |
rightofsb->SetMarkerStyle(24); |
| 3506 |
< |
rightsfsb->SetMarkerStyle(32); |
| 3507 |
< |
} |
| 3508 |
< |
|
| 3509 |
< |
TLegend *leg3 = make_legend("MC t#bar{t}",0.55,0.75,false); |
| 3510 |
< |
leg3->AddEntry(rightofzp,"OFZP / SFZP","ple"); |
| 3511 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3512 |
< |
leg3->AddEntry(rightsfsb,"SFSB / SFZP","ple"); |
| 3513 |
< |
leg3->AddEntry(rightofsb,"OFSB / SFZP","ple"); |
| 3514 |
< |
} |
| 3515 |
< |
leg3->AddEntry(bottom,"syst. envelope","l"); |
| 3516 |
< |
leg3->SetX1(0.25);leg3->SetX2(0.6); |
| 3517 |
< |
leg3->SetY1(0.65); |
| 3478 |
> |
tcan->SetLogy(0); |
| 3479 |
> |
rightofzp->Divide(rightsfzp); |
| 3480 |
> |
rightofzp->GetXaxis()->SetRangeUser(0.0,binning[binning.size()-1]); |
| 3481 |
> |
rightofzp->GetYaxis()->SetRangeUser(0.0,2.5); |
| 3482 |
> |
rightofzp->GetYaxis()->SetTitle("#deltaJZB ratio"); |
| 3483 |
> |
rightofzp->Draw(); |
| 3484 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3485 |
> |
rightofsb->Divide(rightsfzp); |
| 3486 |
> |
rightsfsb->Divide(rightsfzp); |
| 3487 |
> |
rightofsb->Draw("same"); |
| 3488 |
> |
rightsfsb->Draw("same"); |
| 3489 |
> |
} |
| 3490 |
> |
|
| 3491 |
> |
TLine *top2 = new TLine(0.0,1.0+linepos,binning[binning.size()-1],1.0+linepos); |
| 3492 |
> |
TLine *center2 = new TLine(0.0,1.0,binning[binning.size()-1],1.0); |
| 3493 |
> |
TLine *bottom2 = new TLine(0.0,1.0-linepos,binning[binning.size()-1],1.0-linepos); |
| 3494 |
> |
|
| 3495 |
> |
top2->SetLineColor(kBlue);top2->SetLineStyle(2); |
| 3496 |
> |
bottom2->SetLineColor(kBlue);bottom2->SetLineStyle(2); |
| 3497 |
> |
center2->SetLineColor(kBlue); |
| 3498 |
> |
|
| 3499 |
> |
top2->Draw("same"); |
| 3500 |
> |
center2->Draw("same"); |
| 3501 |
> |
bottom2->Draw("same"); |
| 3502 |
> |
|
| 3503 |
> |
rightofzp->SetMarkerStyle(21); |
| 3504 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3505 |
> |
rightofsb->SetMarkerStyle(24); |
| 3506 |
> |
rightsfsb->SetMarkerStyle(32); |
| 3507 |
> |
} |
| 3508 |
> |
|
| 3509 |
> |
TLegend *leg3 = make_legend("MC t#bar{t}",0.55,0.75,false); |
| 3510 |
> |
leg3->AddEntry(rightofzp,"OFZP / SFZP","ple"); |
| 3511 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3512 |
> |
leg3->AddEntry(rightsfsb,"SFSB / SFZP","ple"); |
| 3513 |
> |
leg3->AddEntry(rightofsb,"OFSB / SFZP","ple"); |
| 3514 |
> |
} |
| 3515 |
> |
leg3->AddEntry(bottom,"syst. envelope","l"); |
| 3516 |
> |
leg3->SetX1(0.25);leg3->SetX2(0.6); |
| 3517 |
> |
leg3->SetY1(0.65); |
| 3518 |
|
|
| 3519 |
< |
leg3->Draw("same"); |
| 3519 |
> |
leg3->Draw("same"); |
| 3520 |
|
|
| 3521 |
< |
DrawMCPrelim(simulatedlumi); |
| 3522 |
< |
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison_ratio"); |
| 3521 |
> |
DrawMCPrelim(simulatedlumi); |
| 3522 |
> |
CompleteSave(tcan,"Systematics/"+prestring+"/ttbar_deltajzb_comparison_ratio"); |
| 3523 |
> |
|
| 3524 |
> |
} |
| 3525 |
|
|
| 3526 |
|
delete TZem; |
| 3527 |
|
delete nTZem; |
| 3528 |
|
delete TZeemm; |
| 3529 |
|
delete nTZeemm; |
| 3530 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3530 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3531 |
|
delete TSem; |
| 3532 |
|
delete nTSem; |
| 3533 |
|
delete TSeemm; |
| 3536 |
|
|
| 3537 |
|
delete tcan; |
| 3538 |
|
cutWeight=weightbackup; |
| 3539 |
+ |
switch_overunderflow(false); |
| 3540 |
|
} |
| 3541 |
|
|
| 3542 |
|
void ttbar_sidebands_comparison(string mcjzb, vector<float> jzb_binning) { |
| 3569 |
|
nicer_binning.push_back(100); |
| 3570 |
|
nicer_binning.push_back(125); |
| 3571 |
|
nicer_binning.push_back(150); |
| 3572 |
< |
nicer_binning.push_back(175); |
| 3572 |
> |
//nicer_binning.push_back(175); |
| 3573 |
|
nicer_binning.push_back(200); |
| 3574 |
< |
nicer_binning.push_back(230); |
| 3575 |
< |
nicer_binning.push_back(280); |
| 3576 |
< |
nicer_binning.push_back(400); |
| 3574 |
> |
// nicer_binning.push_back(250); |
| 3575 |
> |
// nicer_binning.push_back(300); |
| 3576 |
> |
// nicer_binning.push_back(400); |
| 3577 |
|
|
| 3578 |
|
ttbar_sidebands_comparison(mcjzb,nicer_binning, "ttbar/"); |
| 3579 |
|
} |
| 3580 |
|
|
| 3581 |
|
|
| 3582 |
< |
void zjets_prediction_comparison(string mcjzbWithPU) { |
| 3582 |
> |
void zjets_prediction_comparison(string mcjzbWithPUa, TCut massregioncut, string massregionname) { |
| 3583 |
> |
cout << "****************************************************************************************************************************************************************" << endl; |
| 3584 |
|
TCanvas *zcan = new TCanvas("zcan","zcan"); |
| 3585 |
< |
zcan->SetLogy(1); |
| 3585 |
> |
// zcan->SetLogy(1); |
| 3586 |
|
TCut weightbackup=cutWeight; |
| 3587 |
+ |
|
| 3588 |
+ |
bool UsePURW=true; |
| 3589 |
+ |
|
| 3590 |
+ |
|
| 3591 |
+ |
string mcjzb; |
| 3592 |
+ |
if(UsePURW) { |
| 3593 |
+ |
mcjzb=mcjzbWithPUa; |
| 3594 |
+ |
cout << "Using PURW peak positions" << endl; |
| 3595 |
+ |
} else { |
| 3596 |
+ |
// Do it without PU re-weighting |
| 3597 |
+ |
cutWeight="1.0"; |
| 3598 |
+ |
float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0; |
| 3599 |
+ |
stringstream resultsNoPU; |
| 3600 |
+ |
stringstream noPUdatajzb; |
| 3601 |
+ |
stringstream noPUmcjzb; |
| 3602 |
+ |
|
| 3603 |
+ |
find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,resultsNoPU,false,noPUdatajzb,noPUmcjzb); |
| 3604 |
+ |
dout << "The peak corrected JZB expression for MC without pileup is : " << noPUmcjzb.str() << endl; |
| 3605 |
+ |
|
| 3606 |
+ |
mcjzb = noPUmcjzb.str(); |
| 3607 |
+ |
|
| 3608 |
+ |
} |
| 3609 |
|
|
| 2642 |
– |
/* |
| 2643 |
– |
// Do it without PU re-weighting |
| 2644 |
– |
float MCPeakNoPU=0,MCPeakErrorNoPU=0,DataPeakNoPU=0,DataPeakErrorNoPU=0,MCSigma=0,DataSigma=0; |
| 2645 |
– |
stringstream resultsNoPU; |
| 2646 |
– |
|
| 2647 |
– |
stringstream mcjzbnoPU; |
| 2648 |
– |
find_peaks(MCPeakNoPU,MCPeakErrorNoPU, DataPeakNoPU, DataPeakErrorNoPU,MCSigma,DataSigma,resultsNoPU,true); |
| 2649 |
– |
if(MCPeakNoPU>0) mcjzbnoPU<<"("<<jzbvariablemc<<"-"<<TMath::Abs(MCPeakNoPU)<<")"; |
| 2650 |
– |
else mcjzbnoPU<<"("<<jzbvariablemc<<"+"<<TMath::Abs(MCPeakNoPU)<<")"; |
| 2651 |
– |
|
| 2652 |
– |
string mcjzb = mcjzbnoPU.str(); |
| 2653 |
– |
dout << "The peak corrected JZB expression for MC without pileup is : " << mcjzb << endl; |
| 2654 |
– |
*/ |
| 2655 |
– |
// cutWeight="1.0"; |
| 2656 |
– |
|
| 2657 |
– |
string mcjzb = mcjzbWithPU; // this is with PURW, if you want without it you have to uncomment the part above (and comment out this line) |
| 3610 |
|
|
| 3611 |
|
vector<float> binning; |
| 3612 |
|
binning.push_back(0); |
| 3613 |
|
binning.push_back(10); |
| 3614 |
< |
binning.push_back(30); |
| 3615 |
< |
binning.push_back(50); |
| 3616 |
< |
// binning.push_back(80); |
| 3614 |
> |
binning.push_back(20); |
| 3615 |
> |
binning.push_back(40); |
| 3616 |
> |
binning.push_back(60); |
| 3617 |
> |
// binning.push_back(50); |
| 3618 |
> |
// binning.push_back(60); |
| 3619 |
> |
// binning.push_back(70); |
| 3620 |
> |
// binning.push_back(80); |
| 3621 |
> |
// binning.push_back(90); |
| 3622 |
|
binning.push_back(100); |
| 3623 |
< |
// float sbg_min=0.; |
| 2667 |
< |
// float sbg_max=100.; |
| 2668 |
< |
// int sbg_nbins=5; |
| 3623 |
> |
|
| 3624 |
|
float simulatedlumi = luminosity;//in pb please - adjust to your likings |
| 3625 |
|
|
| 3626 |
|
TCut kPos((mcjzb+">0").c_str()); |
| 3627 |
|
TCut kNeg((mcjzb+"<0").c_str()); |
| 3628 |
|
string var( "abs("+mcjzb+")" ); |
| 3629 |
+ |
|
| 3630 |
+ |
TCut notTau("abs(genMID1)!=15"); |
| 3631 |
+ |
TCut ee_mm_tautau("mll>0"); |
| 3632 |
+ |
|
| 3633 |
|
|
| 3634 |
< |
TCut kcut(cutmass&&cutOSSF&&cutnJets); |
| 3635 |
< |
TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",kcut&&kPos,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3636 |
< |
TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",kcut&&kNeg,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3634 |
> |
TH1F *hJZBpos = systsamples.Draw("hJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&¬Tau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3635 |
> |
TH1F *hJZBneg = systsamples.Draw("hJZBneg",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&cutnJets&¬Tau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3636 |
> |
|
| 3637 |
|
hJZBpos->SetLineColor(kBlack); |
| 3638 |
|
hJZBneg->SetLineColor(kRed); |
| 3639 |
|
|
| 3642 |
|
hJZBneg->Draw("same,hist"); |
| 3643 |
|
hJZBpos->Draw("same,e1"); // So it's on top... |
| 3644 |
|
|
| 3645 |
< |
TLegend *leg = make_legend("MC Z+jets",0.55,0.75,false); |
| 3645 |
> |
TLegend *leg = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.55,0.75,false); |
| 3646 |
|
leg->AddEntry(hJZBpos,"Observed","pe"); |
| 3647 |
|
leg->AddEntry(hJZBneg,"Predicted","l"); |
| 3648 |
|
leg->Draw("same"); |
| 3649 |
|
DrawMCPrelim(simulatedlumi); |
| 3650 |
< |
CompleteSave(zcan,"Systematics/zjets_prediction"); |
| 3650 |
> |
CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction"+any2string(massregionname)); |
| 3651 |
|
|
| 3652 |
|
TH1F* hratio = (TH1F*)hJZBpos->Clone("hratio"); |
| 3653 |
|
hratio->Divide(hJZBneg); |
| 3654 |
|
|
| 3655 |
|
for(int i=1;i<=hJZBpos->GetNbinsX();i++) { |
| 3656 |
< |
cout << "Positive: " << hJZBpos->GetBinContent(i) << " vs Negative : " << hJZBneg->GetBinContent(i) << endl; |
| 3656 |
> |
cout << "Positive: " << hJZBpos->GetBinContent(i) << " vs Negative : " << hJZBneg->GetBinContent(i) << " (ratio : " << hJZBpos->GetBinContent(i) / hJZBneg->GetBinContent(i) << endl; |
| 3657 |
|
} |
| 3658 |
|
|
| 3659 |
< |
zcan->SetLogy(0); |
| 3659 |
> |
// zcan->SetLogy(0); |
| 3660 |
|
hratio->GetYaxis()->SetRangeUser(0,2.5); |
| 3661 |
|
hratio->GetYaxis()->SetTitle("Observed/Predicted"); |
| 3662 |
|
hratio->Draw("e1"); |
| 3674 |
|
center->Draw("same"); |
| 3675 |
|
bottom->Draw("same"); |
| 3676 |
|
|
| 3677 |
< |
TLegend *leg2 = make_legend("MC Z+jets",0.25,0.75,false); |
| 3677 |
> |
TLegend *leg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false); |
| 3678 |
|
leg2->AddEntry(hratio,"obs / pred","pe"); |
| 3679 |
|
leg2->AddEntry(bottom,"syst. envelope","l"); |
| 3680 |
|
leg2->Draw("same"); |
| 3681 |
|
DrawMCPrelim(simulatedlumi); |
| 3682 |
< |
CompleteSave(zcan,"Systematics/zjets_prediction_ratio"); |
| 3682 |
> |
CompleteSave(zcan,"Systematics/ZJets/zjets_eemm_prediction_ratio"+any2string(massregionname)); |
| 3683 |
> |
|
| 3684 |
> |
TCut reducedNJets(cutnJets); |
| 3685 |
> |
|
| 3686 |
> |
TH1F *TAUhJZBpos = systsamples.Draw("TAUhJZBpos",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3687 |
> |
TH1F *LcorrJZBeemm = systsamples.Draw("LcorrJZBeemm",var,binning, "JZB [GeV]", "events",cutmass&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3688 |
> |
TH1F *RcorrJZBem = systsamples.Draw("RcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3689 |
> |
TH1F *LcorrJZBem = systsamples.Draw("LcorrJZBem",var,binning, "JZB [GeV]", "events",cutmass&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3690 |
> |
|
| 3691 |
> |
TH1F *RcorrJZBSBem; |
| 3692 |
> |
TH1F *LcorrJZBSBem; |
| 3693 |
> |
TH1F *RcorrJZBSBeemm; |
| 3694 |
> |
TH1F *LcorrJZBSBeemm; |
| 3695 |
> |
|
| 3696 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3697 |
> |
RcorrJZBSBem = systsamples.Draw("RcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3698 |
> |
LcorrJZBSBem = systsamples.Draw("LcorrJZBSBem",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSOF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3699 |
> |
RcorrJZBSBeemm = systsamples.Draw("RcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kPos&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3700 |
> |
LcorrJZBSBeemm = systsamples.Draw("LcorrJZBSBeemm",var,binning, "JZB [GeV]", "events",sidebandcut&&cutOSSF&&reducedNJets&&ee_mm_tautau&&kNeg&&massregioncut,mc,simulatedlumi,systsamples.FindSample("/DY")); |
| 3701 |
> |
} |
| 3702 |
> |
|
| 3703 |
> |
TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred"); |
| 3704 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3705 |
> |
Bpred->Add(RcorrJZBem,1.0/3.); |
| 3706 |
> |
Bpred->Add(LcorrJZBem,-1.0/3.); |
| 3707 |
> |
Bpred->Add(RcorrJZBSBem,1.0/3.); |
| 3708 |
> |
Bpred->Add(LcorrJZBSBem,-1.0/3.); |
| 3709 |
> |
Bpred->Add(RcorrJZBSBeemm,1.0/3.); |
| 3710 |
> |
Bpred->Add(LcorrJZBSBeemm,-1.0/3.); |
| 3711 |
> |
} else { |
| 3712 |
> |
Bpred->Add(RcorrJZBem,1.0); |
| 3713 |
> |
Bpred->Add(LcorrJZBem,-1.0); |
| 3714 |
> |
} |
| 3715 |
> |
|
| 3716 |
> |
Bpred->SetLineColor(kRed); |
| 3717 |
> |
|
| 3718 |
> |
TAUhJZBpos->SetLineColor(kBlack); |
| 3719 |
> |
Bpred->SetLineColor(kRed); |
| 3720 |
> |
|
| 3721 |
> |
TAUhJZBpos->SetMinimum(1.0); |
| 3722 |
> |
TAUhJZBpos->Draw("e1"); |
| 3723 |
> |
Bpred->Draw("same,hist"); |
| 3724 |
> |
TAUhJZBpos->Draw("same,e1"); |
| 3725 |
> |
|
| 3726 |
> |
TLegend *TAUleg = make_legend("MC Z+jets #rightarrow ee,#mu#mu,#tau#tau",0.55,0.75,false); |
| 3727 |
> |
TAUleg->AddEntry(TAUhJZBpos,"Observed","pe"); |
| 3728 |
> |
TAUleg->AddEntry(Bpred,"Predicted","l"); |
| 3729 |
> |
TAUleg->Draw("same"); |
| 3730 |
> |
DrawMCPrelim(simulatedlumi); |
| 3731 |
> |
CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction__"+any2string(massregionname)); |
| 3732 |
> |
|
| 3733 |
> |
TH1F* TAUhratio = (TH1F*)TAUhJZBpos->Clone("TAUhratio"); |
| 3734 |
> |
TAUhratio->Divide(Bpred); |
| 3735 |
> |
|
| 3736 |
> |
for(int i=1;i<=TAUhJZBpos->GetNbinsX();i++) { |
| 3737 |
> |
cout << "ee/mm/tautau observed: " << TAUhJZBpos->GetBinContent(i) << " vs predicted : " << Bpred->GetBinContent(i) << " (ratio : " << TAUhJZBpos->GetBinContent(i) / Bpred->GetBinContent(i) << endl; |
| 3738 |
> |
} |
| 3739 |
> |
|
| 3740 |
> |
zcan->SetLogy(0); |
| 3741 |
> |
TAUhratio->GetYaxis()->SetRangeUser(0,2.5); |
| 3742 |
> |
TAUhratio->GetYaxis()->SetTitle("Observed/Predicted"); |
| 3743 |
> |
TAUhratio->Draw("e1"); |
| 3744 |
> |
|
| 3745 |
> |
top->Draw("same"); |
| 3746 |
> |
center->Draw("same"); |
| 3747 |
> |
bottom->Draw("same"); |
| 3748 |
> |
|
| 3749 |
> |
TLegend *TAUleg2 = make_legend("MC Z+jets #rightarrow ee,#mu#mu",0.25,0.75,false); |
| 3750 |
> |
TAUleg2->AddEntry(TAUhratio,"obs / pred","pe"); |
| 3751 |
> |
TAUleg2->AddEntry(bottom,"syst. envelope","l"); |
| 3752 |
> |
TAUleg2->Draw("same"); |
| 3753 |
> |
DrawMCPrelim(simulatedlumi); |
| 3754 |
> |
CompleteSave(zcan,"Systematics/ZJets/zjets_eemumutautau_prediction_ratio"+any2string(massregionname)); |
| 3755 |
> |
|
| 3756 |
> |
delete Bpred; |
| 3757 |
> |
delete TAUhJZBpos; |
| 3758 |
> |
delete LcorrJZBeemm; |
| 3759 |
> |
delete RcorrJZBem; |
| 3760 |
> |
delete LcorrJZBem; |
| 3761 |
> |
delete hJZBpos; |
| 3762 |
> |
delete hJZBneg; |
| 3763 |
> |
|
| 3764 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3765 |
> |
delete RcorrJZBSBem; |
| 3766 |
> |
delete LcorrJZBSBem; |
| 3767 |
> |
delete RcorrJZBSBeemm; |
| 3768 |
> |
delete LcorrJZBSBeemm; |
| 3769 |
> |
} |
| 3770 |
> |
|
| 3771 |
|
|
| 3772 |
|
delete zcan; |
| 3773 |
|
cutWeight=weightbackup; |
| 2727 |
– |
|
| 3774 |
|
} |
| 3775 |
|
|
| 3776 |
|
|
| 3777 |
< |
|
| 3777 |
> |
void zjets_prediction_comparison(string mcjzbWithPUa) { |
| 3778 |
> |
zjets_prediction_comparison(mcjzbWithPUa, TCut(""), "nomasscut"); |
| 3779 |
> |
zjets_prediction_comparison(mcjzbWithPUa, TCut("abs(mll-91)<20"), "Zwindow_20"); |
| 3780 |
> |
zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>20&&mll<70"), "LowMassRegion2070"); |
| 3781 |
> |
zjets_prediction_comparison(mcjzbWithPUa, TCut("mll>110"), "HighMassRegion110"); |
| 3782 |
> |
} |
| 3783 |
> |
|
| 3784 |
|
void sideband_assessment(string datajzb, float min=30.0, float max=50.0) { |
| 3785 |
|
tout << endl << endl; |
| 3786 |
|
stringstream bordercut; |
| 3831 |
|
TH1F *RcorrJZBSBeemm; |
| 3832 |
|
TH1F *LcorrJZBSBeemm; |
| 3833 |
|
|
| 3834 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3834 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3835 |
|
RcorrJZBSBem = coll.Draw("RcorrJZBSBem",jzbexpr.c_str(),nbins,low,14000, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity, mysample); |
| 3836 |
|
LcorrJZBSBem = coll.Draw("LcorrJZBSBem",("-("+jzbexpr+")").c_str(),nbins,low,14000, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity, mysample); |
| 3837 |
|
RcorrJZBSBeemm = coll.Draw("RcorrJZBSBeemm",jzbexpr.c_str(),nbins,low,14000, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity, mysample); |
| 3839 |
|
} |
| 3840 |
|
|
| 3841 |
|
TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred"); |
| 3842 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 3842 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 3843 |
|
Bpred->Add(RcorrJZBem,1.0/3.); |
| 3844 |
|
Bpred->Add(LcorrJZBem,-1.0/3.); |
| 3845 |
|
Bpred->Add(RcorrJZBSBem,1.0/3.); |
| 4050 |
|
TH1F *JRcorrJZBSBeemm; |
| 4051 |
|
TH1F *JLcorrJZBSBeemm; |
| 4052 |
|
|
| 4053 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4053 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4054 |
|
RcorrJZBSBem = qcdsamples.Draw("RcorrJZBSBem",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 4055 |
|
LcorrJZBSBem = qcdsamples.Draw("LcorrJZBSBem",("-"+jzb).c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSOF&&cutnJets,is_data, luminosity,use_signal); |
| 4056 |
|
RcorrJZBSBeemm = qcdsamples.Draw("RcorrJZBSBeemm",jzb.c_str(),nbins,low,hi, "JZB [GeV]", "events", sidebandcut&&cutOSSF&&cutnJets,is_data, luminosity,use_signal); |
| 4071 |
|
|
| 4072 |
|
TH1F *Bpred = (TH1F*)LcorrJZBeemm->Clone("Bpred"); |
| 4073 |
|
TH1F *JBpred = (TH1F*)JLcorrJZBeemm->Clone("Bpred"); |
| 4074 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4074 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4075 |
|
Bpred->Add(RcorrJZBem,1.0/3.); |
| 4076 |
|
Bpred->Add(LcorrJZBem,-1.0/3.); |
| 4077 |
|
Bpred->Add(RcorrJZBSBem,1.0/3.); |
| 4141 |
|
dout << "______________________________________________" << endl; |
| 4142 |
|
dout << "How QCD shows up in the different regions: " << endl; |
| 4143 |
|
dout << "OSSF: " << endl; |
| 4144 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4144 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4145 |
|
dout << " Z window: \t" << RcorrJZBeemm->Integral() << " (JZB>0) , " << LcorrJZBeemm->Integral() << " (JZB<0) --> total: " << RcorrJZBeemm->Integral() + LcorrJZBeemm->Integral() << endl; |
| 4146 |
|
dout << " sideband: \t" << RcorrJZBSBeemm->Integral() << " (JZB>0) , " << LcorrJZBSBeemm->Integral() << " (JZB<0) --> total: " << RcorrJZBSBeemm->Integral() + LcorrJZBSBeemm->Integral() << endl; |
| 4147 |
|
} else { |
| 4148 |
|
dout << " " << RcorrJZBeemm->Integral() << " (JZB>0) , " << LcorrJZBeemm->Integral() << " (JZB<0) --> total: " << RcorrJZBeemm->Integral() + LcorrJZBeemm->Integral() << endl; |
| 4149 |
|
} |
| 4150 |
|
dout << "OSOF: " << endl; |
| 4151 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4151 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4152 |
|
dout << " Z window: \t" << RcorrJZBem->Integral() << " (JZB>0) , " << LcorrJZBem->Integral() << " (JZB<0) --> total: " << RcorrJZBem->Integral() + LcorrJZBem->Integral() << endl; |
| 4153 |
|
dout << " sideband: \t" << RcorrJZBSBem->Integral() << " (JZB>0) , " << LcorrJZBSBem->Integral() << " (JZB<0) --> total: " << RcorrJZBSBem->Integral() + LcorrJZBSBem->Integral() << endl; |
| 4154 |
|
} else { |
| 4155 |
|
dout << " Z window: \t" << RcorrJZBem->Integral() << " (JZB>0) , " << LcorrJZBem->Integral() << " (JZB<0) --> total: " << RcorrJZBem->Integral() + LcorrJZBem->Integral() << endl; |
| 4156 |
|
} |
| 4157 |
|
dout << "Therefore: " << endl; |
| 4158 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4158 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4159 |
|
dout << " Prediction increases by : " << LcorrJZBeemm->Integral() << " + (1.0/3)*(" << RcorrJZBSBeemm->Integral() <<"-"<< LcorrJZBSBeemm->Integral() << ") (SFSB) "; |
| 4160 |
|
dout << " + (1.0/3)*(" << RcorrJZBem->Integral() <<"-"<< LcorrJZBem->Integral() << ") (OFZP) "; |
| 4161 |
|
dout << " + (1.0/3)*(" << RcorrJZBSBem->Integral() <<"-"<< LcorrJZBSBem->Integral() << ") (OFSB) "; |
| 4171 |
|
for(int i=0;i<(int)bins.size();i++) { |
| 4172 |
|
dout << " JZB > " << bins[i] << " : " << endl; |
| 4173 |
|
dout << " Observation increases by : " << RcorrJZBeemm->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) << endl; |
| 4174 |
< |
if(PlottingSetup::RestrictToMassPeak) { |
| 4174 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB) { |
| 4175 |
|
dout << " Prediction increases by : " << LcorrJZBeemm->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) + (1.0/3)*(RcorrJZBSBeemm->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) - LcorrJZBSBeemm->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) + RcorrJZBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) - LcorrJZBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) + RcorrJZBSBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) - LcorrJZBSBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX())) << endl; |
| 4176 |
|
} else { |
| 4177 |
|
dout << " Prediction increases by : " << LcorrJZBeemm->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) + RcorrJZBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) - LcorrJZBem->Integral(RcorrJZBeemm->FindBin(bins[i]),RcorrJZBeemm->GetNbinsX()) << endl; |
| 4185 |
|
if(LcorrJZBem) delete LcorrJZBem; |
| 4186 |
|
if(RcorrJZBeemm) delete RcorrJZBeemm; |
| 4187 |
|
if(LcorrJZBeemm) delete LcorrJZBeemm; |
| 4188 |
< |
if(PlottingSetup::RestrictToMassPeak&&RcorrJZBSBem) delete RcorrJZBSBem; |
| 4189 |
< |
if(PlottingSetup::RestrictToMassPeak&&LcorrJZBSBem) delete LcorrJZBSBem; |
| 4190 |
< |
if(PlottingSetup::RestrictToMassPeak&&RcorrJZBSBeemm) delete RcorrJZBSBeemm; |
| 4191 |
< |
if(PlottingSetup::RestrictToMassPeak&&LcorrJZBSBeemm) delete LcorrJZBSBeemm; |
| 4188 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB&&RcorrJZBSBem) delete RcorrJZBSBem; |
| 4189 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB&&LcorrJZBSBem) delete LcorrJZBSBem; |
| 4190 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB&&RcorrJZBSBeemm) delete RcorrJZBSBeemm; |
| 4191 |
> |
if(PlottingSetup::RestrictToMassPeak&&PlottingSetup::UseSidebandsForcJZB&&LcorrJZBSBeemm) delete LcorrJZBSBeemm; |
| 4192 |
|
} |
| 4193 |
|
|
| 4194 |
|
void check_ptsanity() { |
| 4285 |
|
} |
| 4286 |
|
} |
| 4287 |
|
|
| 4288 |
< |
void met_vs_jzb_plots() { |
| 4288 |
> |
void met_vs_jzb_plots(string datajzb, string mcjzb) { |
| 4289 |
|
|
| 4290 |
|
TCanvas *canmetjzb = new TCanvas("canmet","MET vs JZB canvas"); |
| 4291 |
|
canmetjzb->SetRightMargin(0.16); |
| 4294 |
|
findme.push_back("DY"); |
| 4295 |
|
findme.push_back("TTJets"); |
| 4296 |
|
findme.push_back("LM"); |
| 4297 |
< |
|
| 4297 |
> |
/* |
| 4298 |
|
for(int ifind=0;ifind<(int)findme.size();ifind++) { |
| 4299 |
|
vector<int> selsamples = allsamples.FindSample(findme[ifind]); |
| 4300 |
|
TH2F *metvsjzb = new TH2F("metvsjzb","metvsjzb",200,0,100,400,-100,100); |
| 4301 |
|
for(int isel=0;isel<(int)selsamples.size();isel++) { |
| 4302 |
< |
cout << "Producing MET:JZB plot ... working on sample: " << allsamples.collection[selsamples[isel]].filename << endl; |
| 4303 |
< |
allsamples.collection[selsamples[isel]].events->Draw("jzb[1]:met[4]>>+metvsjzb",cutmass&&cutOSSF); |
| 4302 |
> |
dout << "Producing MET:JZB plot ... working on sample: " << allsamples.collection[selsamples[isel]].filename << endl; |
| 4303 |
> |
allsamples.collection[selsamples[isel]].events->Draw("jzb[1]:met[4]>>+metvsjzb",cutmass&&cutOSSF&&cutnJets); |
| 4304 |
|
} |
| 4305 |
|
metvsjzb->Scale(allsamples.collection[selsamples[0]].weight); |
| 4306 |
|
metvsjzb->SetStats(0); |
| 4314 |
|
title->Draw(); |
| 4315 |
|
CompleteSave(canmetjzb,(string)"METvsJZBplots/"+findme[ifind]); |
| 4316 |
|
} |
| 4317 |
+ |
*/ |
| 4318 |
+ |
|
| 4319 |
+ |
dout << "About to produce MET plot for DY split up by JZB" << endl; |
| 4320 |
+ |
|
| 4321 |
+ |
int nbins=14; |
| 4322 |
+ |
float low=0; |
| 4323 |
+ |
float high=140; |
| 4324 |
+ |
|
| 4325 |
+ |
stringstream sLEFT; |
| 4326 |
+ |
sLEFT << "((" << mcjzb << ")<0)"; |
| 4327 |
+ |
TCut LEFT(sLEFT.str().c_str()); |
| 4328 |
+ |
stringstream sRIGHT; |
| 4329 |
+ |
sRIGHT << "((" << mcjzb << ")>0)"; |
| 4330 |
+ |
TCut RIGHT(sRIGHT.str().c_str()); |
| 4331 |
+ |
|
| 4332 |
+ |
TH1F *metleft = allsamples.Draw("metleft","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4333 |
+ |
TH1F *metleftO = allsamples.Draw("metleftO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4334 |
+ |
TH1F *metright = allsamples.Draw("metright","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4335 |
+ |
TH1F *metrightO = allsamples.Draw("metrightO","met[4]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4336 |
+ |
|
| 4337 |
+ |
|
| 4338 |
+ |
TH1F *Bpred = (TH1F*)metleft->Clone("Bpred"); |
| 4339 |
+ |
Bpred->Add(metleftO,-1); |
| 4340 |
+ |
Bpred->Add(metrightO); |
| 4341 |
+ |
TH1F *obs = (TH1F*)metright->Clone("obs"); |
| 4342 |
+ |
|
| 4343 |
+ |
metleft->Add(metleftO,-1); |
| 4344 |
+ |
metright->Add(metrightO,-1); |
| 4345 |
+ |
|
| 4346 |
+ |
metleft->SetLineColor(kRed); |
| 4347 |
+ |
metright->SetLineColor(kBlack); |
| 4348 |
+ |
TPad *metpad = new TPad("metpad","metpad",0,0,1,1); |
| 4349 |
+ |
metpad->cd(); |
| 4350 |
+ |
metpad->SetLogy(1); |
| 4351 |
+ |
metleft->Draw("histo"); |
| 4352 |
+ |
metright->Draw("same"); |
| 4353 |
+ |
TLegend *lg = make_legend(); |
| 4354 |
+ |
lg->SetX1(0.5); |
| 4355 |
+ |
lg->SetY1(0.7); |
| 4356 |
+ |
lg->AddEntry(metright,"JZB>0 (OSOF corrected)","P"); |
| 4357 |
+ |
lg->AddEntry(metleft,"JZB<0 (OSOF corrected)","L"); |
| 4358 |
+ |
lg->SetHeader("DY"); |
| 4359 |
+ |
|
| 4360 |
+ |
lg->Draw(); |
| 4361 |
+ |
save_with_ratio(metright,metleft,metpad->cd(),"METvsJZBplots/ComparingLeftToRightinMETspectrum"); |
| 4362 |
+ |
|
| 4363 |
+ |
TPad *metpad3 = new TPad("metpad3","metpad3",0,0,1,1); |
| 4364 |
+ |
metpad3->cd(); |
| 4365 |
+ |
metpad3->SetLogy(1); |
| 4366 |
+ |
Bpred->SetLineColor(kRed); |
| 4367 |
+ |
Bpred->Draw("histo"); |
| 4368 |
+ |
obs->SetLineColor(kBlack); |
| 4369 |
+ |
obs->Draw("same"); |
| 4370 |
+ |
TLegend *lg2 = make_legend(); |
| 4371 |
+ |
lg2->SetX1(0.5); |
| 4372 |
+ |
lg2->SetY1(0.7); |
| 4373 |
+ |
lg2->AddEntry(obs,"observed","P"); |
| 4374 |
+ |
lg2->AddEntry(Bpred,"predicted","L"); |
| 4375 |
+ |
lg2->SetHeader("DY"); |
| 4376 |
+ |
|
| 4377 |
+ |
lg2->Draw(); |
| 4378 |
+ |
|
| 4379 |
+ |
save_with_ratio(obs,Bpred,metpad3->cd(),"METvsJZBplots/ComparingPredObsinMET"); |
| 4380 |
+ |
|
| 4381 |
+ |
TPad *metpad2 = new TPad("metpad2","metpad2",0,0,1,1); |
| 4382 |
+ |
metpad2->cd(); |
| 4383 |
+ |
metpad2->SetLogy(1); |
| 4384 |
+ |
|
| 4385 |
+ |
TH1F *metlefta = allsamples.Draw("metlefta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4386 |
+ |
TH1F *metleftOa = allsamples.Draw("metleftOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4387 |
+ |
TH1F *metrighta = allsamples.Draw("metrighta","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4388 |
+ |
TH1F *metrightOa = allsamples.Draw("metrightOa","met[2]",nbins,low,high, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity, allsamples.FindSample("DYJets")); |
| 4389 |
+ |
|
| 4390 |
+ |
metlefta->Add(metleftOa,-1); |
| 4391 |
+ |
metrighta->Add(metrightOa,-1); |
| 4392 |
+ |
|
| 4393 |
+ |
metlefta->SetLineColor(kRed); |
| 4394 |
+ |
metpad2->cd(); |
| 4395 |
+ |
metlefta->Draw("histo"); |
| 4396 |
+ |
metrighta->Draw("same"); |
| 4397 |
+ |
lg->Draw(); |
| 4398 |
+ |
save_with_ratio(metrighta,metlefta,metpad2->cd(),"METvsJZBplots/ComparingLeftToRightinMET_type1_spectrum"); |
| 4399 |
+ |
|
| 4400 |
+ |
delete Bpred; |
| 4401 |
+ |
delete obs; |
| 4402 |
+ |
|
| 4403 |
+ |
float newhigh=300; |
| 4404 |
+ |
int newNBins=30; |
| 4405 |
+ |
|
| 4406 |
+ |
TPad *metpad4 = new TPad("metpad4","metpad4",0,0,1,1); |
| 4407 |
+ |
TH1F *Ametleft = allsamples.Draw("Ametleft","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&LEFT,mc, luminosity); |
| 4408 |
+ |
TH1F *AmetleftO = allsamples.Draw("AmetleftO","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&LEFT,mc, luminosity); |
| 4409 |
+ |
TH1F *Ametright = allsamples.Draw("Ametright","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSSF&&cutnJets&&RIGHT,mc, luminosity); |
| 4410 |
+ |
TH1F *AmetrightO = allsamples.Draw("AmetrightO","met[4]",newNBins,low,newhigh, "MET [GeV]", "events", cutmass&&cutOSOF&&cutnJets&&RIGHT,mc, luminosity); |
| 4411 |
+ |
|
| 4412 |
+ |
TH1F *aBpred = (TH1F*)Ametleft->Clone("aBpred"); |
| 4413 |
+ |
aBpred->Add(AmetleftO,-1); |
| 4414 |
+ |
aBpred->Add(AmetrightO); |
| 4415 |
+ |
aBpred->SetLineColor(kRed); |
| 4416 |
+ |
|
| 4417 |
+ |
TH1F *aobs = (TH1F*)Ametright->Clone("aobs"); |
| 4418 |
+ |
metpad4->cd(); |
| 4419 |
+ |
metpad4->SetLogy(1); |
| 4420 |
+ |
aobs->Draw(); |
| 4421 |
+ |
aBpred->Draw("histo,same"); |
| 4422 |
+ |
aobs->Draw("same"); |
| 4423 |
+ |
lg->SetHeader("All MC"); |
| 4424 |
+ |
lg->Draw(); |
| 4425 |
+ |
save_with_ratio(aobs,aBpred,metpad4->cd(),"METvsJZBplots/ComparingPredObsinMET_ALLSAMPLES"); |
| 4426 |
+ |
|
| 4427 |
+ |
|
| 4428 |
+ |
delete lg; |
| 4429 |
+ |
delete canmetjzb; |
| 4430 |
+ |
delete metleft; |
| 4431 |
+ |
delete metleftO; |
| 4432 |
+ |
delete metright; |
| 4433 |
+ |
delete metrightO; |
| 4434 |
|
} |
| 4435 |
|
|
| 4436 |
|
|
