| 14 |
|
|
| 15 |
|
#how and which groups in the plot |
| 16 |
|
setup: ZH,VV,DYb,DYc,DYlight,TT,ST |
| 17 |
< |
setupLog: DYb,DYc,DYlight,TT,ST,VV,ZH |
| 17 |
> |
setupLog: ZH,DYb,DYc,DYlight,TT,ST,VV |
| 18 |
> |
#setupLog: DYb,DYc,DYlight,TT,ST,VV,ZH |
| 19 |
|
|
| 20 |
|
#colors |
| 21 |
|
colorDict: {'DYlight':401,'DYc':41,'DYb':5,'TT':596,'ST':602,'VV':17,'ZH':2} |
| 22 |
< |
typLegendDict: {'DYlight':'Z+udsg','DYc':'Z + c#bar{c}','DYb':'Z + b#bar{b}','TT':'t#bar{t}','ST':'Single top','VV':'VV','ZH':'ZH(125 GeV)'} |
| 22 |
> |
typLegendDict: {'DYlight':'Z+udsg','DYc':'Z + c#bar{c}','DYb':'Z + b#bar{b}','TT':'t#bar{t}','ST':'Single top','VV':'VV','ZH':'ZH','Overlay':'ZH (125 GeV)'} |
| 23 |
|
|
| 24 |
|
addOverFlow = True |
| 25 |
|
|
| 26 |
|
#To be processed |
| 26 |
– |
#List: Light_RTight,bb_RTight,Top_RTight,Light_RMed,bb_RMed,Top_RMed |
| 27 |
– |
#List: Light_RTight,Light_RMed |
| 27 |
|
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm |
| 29 |
– |
|
| 28 |
|
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,HighPt_Zee,LowPt_Zee,HighPtLooseBTag_Zee |
| 31 |
– |
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm |
| 32 |
– |
|
| 33 |
– |
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm |
| 34 |
– |
|
| 29 |
|
#List: HighPt_Zmm,LowPt_Zmm,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_Zee,LowPt_Zee,HighPtLooseBTag_Zee,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 30 |
< |
List: HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 30 |
> |
#List: HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 31 |
> |
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 32 |
> |
#List: HighPt_Zmm,LowPt_Zmm,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_Zee,LowPt_Zee,HighPtLooseBTag_Zee,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 33 |
|
|
| 34 |
< |
#List: HighPt_Zmm,LowPt_Zmm,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm |
| 35 |
< |
#List: HighPt_Zee,LowPt_Zee,HighPtLooseBTag_Zee,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 40 |
< |
#List: RTight,RMed |
| 34 |
> |
#ALL PLOTS |
| 35 |
> |
#List: HighPt_Zmm,LowPt_Zmm,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_Zee,LowPt_Zee,HighPtLooseBTag_Zee,HighPt_MJJ_Zee,LowPt_MJJ_Zee,ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm,ZLight_HighPt_beforeReg,ZLight_HighPt,ZLight_LowPt_beforeReg,ZLight_LowPt |
| 36 |
|
|
| 37 |
< |
#,Light |
| 37 |
> |
#List: HighPt_Zmm |
| 38 |
> |
List: HighPtLooseBTag_Zmm |
| 39 |
> |
# Regression |
| 40 |
> |
#List: ZLight,ZLight_beforeReg,Zbb,Zbb_beforeReg,TTbar,TTbar_beforeReg |
| 41 |
|
|
| 42 |
< |
#List: ZLight_HighPt_Zee,ZLight_LowPt_Zee,Zbb_Zee,TTbar_Zee,ZLight_HighPt_Zmm,ZLight_LowPt_Zmm,Zbb_Zmm,TTbar_Zmm,HighPt_Zee,HighPt_Zmm,LowPt_Zee,LowPt_Zmm,HighPtLooseBTag_Zee,HighPtLooseBTag_Zmm,HighPt_MJJ_Zmm,LowPt_MJJ_Zmm,HighPt_MJJ_Zee,LowPt_MJJ_Zee |
| 42 |
> |
#List: ZLight_HighPt_beforeReg,ZLight_HighPt,ZLight_LowPt_beforeReg,ZLight_LowPt |
| 43 |
|
|
| 44 |
|
#var definition |
| 47 |
– |
#BDT_vars: Hmass,Hpt,Vpt,HVdPhi,Vmass,HdR,csv0,csv1,AK7,Sum |
| 45 |
|
BDT_vars: Hmass,Hpt,Vpt,HVdPhi,Vmass,HdR,HdEta,minCSV,maxCSV,hJet0pt,hJet1pt,PullAngle,PullAngle2,rHVpt,Sum |
| 46 |
|
BDTs: BDT8_RTight,BDT8_RMed,BDT8_RTightLooseBTag |
| 47 |
|
|
| 197 |
|
log: False |
| 198 |
|
blind: False |
| 199 |
|
|
| 200 |
+ |
[Plot:ZLight_LowPt] |
| 201 |
+ |
Datas: Zmm Zee |
| 202 |
+ |
Normalize: True |
| 203 |
+ |
vars: Hmass |
| 204 |
+ |
Signal: ZH125 |
| 205 |
+ |
log: False |
| 206 |
+ |
blind: False |
| 207 |
+ |
|
| 208 |
+ |
[Plot:ZLight_LowPt_beforeReg] |
| 209 |
+ |
Datas: Zmm Zee |
| 210 |
+ |
Normalize: True |
| 211 |
+ |
vars: Hmass_NoReg |
| 212 |
+ |
Signal: ZH125 |
| 213 |
+ |
log: False |
| 214 |
+ |
blind: False |
| 215 |
+ |
|
| 216 |
+ |
[Plot:ZLight_HighPt] |
| 217 |
+ |
Datas: Zmm Zee |
| 218 |
+ |
Normalize: True |
| 219 |
+ |
vars: Hmass |
| 220 |
+ |
Signal: ZH125 |
| 221 |
+ |
log: False |
| 222 |
+ |
blind: False |
| 223 |
+ |
|
| 224 |
+ |
[Plot:ZLight_HighPt_beforeReg] |
| 225 |
+ |
Datas: Zmm Zee |
| 226 |
+ |
Normalize: True |
| 227 |
+ |
vars: Hmass_NoReg |
| 228 |
+ |
Signal: ZH125 |
| 229 |
+ |
log: False |
| 230 |
+ |
blind: False |
| 231 |
+ |
|
| 232 |
+ |
[Plot:Zbb_beforeReg] |
| 233 |
+ |
Datas: Zmm Zee |
| 234 |
+ |
Normalize: True |
| 235 |
+ |
vars: Hmass_NoReg |
| 236 |
+ |
Signal: ZH125 |
| 237 |
+ |
log: False |
| 238 |
+ |
blind: False |
| 239 |
+ |
|
| 240 |
+ |
[Plot:Zbb] |
| 241 |
+ |
Datas: Zmm Zee |
| 242 |
+ |
Normalize: True |
| 243 |
+ |
vars: Hmass |
| 244 |
+ |
Signal: ZH125 |
| 245 |
+ |
log: False |
| 246 |
+ |
blind: False |
| 247 |
+ |
|
| 248 |
+ |
[Plot:TTbar_beforeReg] |
| 249 |
+ |
Datas: Zmm Zee |
| 250 |
+ |
Normalize: True |
| 251 |
+ |
vars: Hmass_NoReg |
| 252 |
+ |
Signal: ZH125 |
| 253 |
+ |
log: False |
| 254 |
+ |
blind: False |
| 255 |
+ |
|
| 256 |
+ |
[Plot:TTbar] |
| 257 |
+ |
Datas: Zmm Zee |
| 258 |
+ |
Normalize: True |
| 259 |
+ |
vars: Hmass |
| 260 |
+ |
Signal: ZH125 |
| 261 |
+ |
log: False |
| 262 |
+ |
blind: False |
| 263 |
+ |
|
| 264 |
+ |
|
| 265 |
+ |
|
| 266 |
|
|
| 267 |
|
|
| 268 |
|
########################################################################################### |
| 269 |
|
# REGRESSION |
| 270 |
|
[Regression] |
| 271 |
< |
regWeight = ../data/MVA_BDT_REG_Jun16.weights.xml |
| 272 |
< |
regDict = {"Jet_pt": "hJet_pt[0]", "Jet_eta": "hJet_eta[0]", "Jet_e": "hJet_e[0]", "Jet_ptRaw": "Jet_ptRaw", "Jet_chf": "hJet_chf[0]","Jet_nconstituents": "hJet_nconstituents[0]", "Jet_vtxPt": "hJet_vtxPt[0]", "Jet_vtx3dL": "hJet_vtx3dL[0]", "Jet_vtx3deL": "hJet_vtx3deL[0]", "Jet_MET_dPhi": "Jet_MET_dPhi", 'METet': 'METet'} |
| 273 |
< |
regVars = ["Jet_pt","Jet_eta","Jet_e","Jet_ptRaw", "Jet_chf","Jet_nconstituents", "Jet_vtxPt", "Jet_vtx3dL", "Jet_vtx3deL", "Jet_MET_dPhi", "METet"] |
| 271 |
> |
#regWeight = ../data/MVA_BDT_REG_Jun16.weights.xml |
| 272 |
> |
#regDict = {"Jet_pt": "hJet_pt[0]", "Jet_eta": "hJet_eta[0]", "Jet_e": "hJet_e[0]", "Jet_ptRaw": "Jet_ptRaw", "Jet_chf": "hJet_chf[0]","Jet_nconstituents": "hJet_nconstituents[0]", "Jet_vtxPt": "hJet_vtxPt[0]", "Jet_vtx3dL": "hJet_vtx3dL[0]", "Jet_vtx3deL": "hJet_vtx3deL[0]", "Jet_MET_dPhi": "Jet_MET_dPhi", 'METet': 'METet'} |
| 273 |
> |
#regVars = ["Jet_pt","Jet_eta","Jet_e","Jet_ptRaw", "Jet_chf","Jet_nconstituents", "Jet_vtxPt", "Jet_vtx3dL", "Jet_vtx3deL", "Jet_MET_dPhi", "METet"] |
| 274 |
> |
regWeight = ../data/MVA_BDT_REG_Matt.weights.xml |
| 275 |
> |
regDict = {"hJet_pt":"hJet_pt[0]","evalJERBias(hJet_ptRaw,hJet_genPt,hJet_eta)":"Jet_ptRaw","rho25":"rho25","hJet_JECUnc":"hJet_JECUnc[0]","evalEt(hJet_pt,hJet_eta,hJet_phi,hJet_e)":"Jet_et","evalMt(hJet_pt,hJet_eta,hJet_phi,hJet_e)":"Jet_mt","hJet_ptLeadTrack":"hJet_ptLeadTrack[0]","max(0,(hJet_SoftLeptptRel*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))":"max(0,(hJet_SoftLeptptRel[0]*(hJet_SoftLeptIdlooseMu[0]==1||hJet_SoftLeptId95[0]==1)))","max(0,(hJet_SoftLeptdR*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))":"max(0,(hJet_SoftLeptdR[0]*(hJet_SoftLeptIdlooseMu[0]==1||hJet_SoftLeptId95[0]==1)))","max(0,(hJet_SoftLeptPt*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))":"max(0,(hJet_SoftLeptPt[0]*(hJet_SoftLeptIdlooseMu[0]==1||hJet_SoftLeptId95[0]==1)))","hJet_cef":"hJet_cef[0]","hJet_nconstituents":"hJet_nconstituents[0]","max(0,hJet_vtxPt)":"max(0,hJet_vtxPt[0])","max(0,hJet_vtxMass)":"max(0,hJet_vtxMass[0])","max(0,hJet_vtx3dL)":"max(0,hJet_vtx3dL[0])","max(0,hJet_vtx3deL)":"max(0,hJet_vtx3deL[0])","METnoPU.et":"METet","deltaPhi(METnoPU.phi,hJet_phi)":"Jet_MET_dPhi"} |
| 276 |
> |
regVars = ["hJet_pt","evalJERBias(hJet_ptRaw,hJet_genPt,hJet_eta)","rho25","hJet_JECUnc","evalEt(hJet_pt,hJet_eta,hJet_phi,hJet_e)","evalMt(hJet_pt,hJet_eta,hJet_phi,hJet_e)","hJet_ptLeadTrack","max(0,(hJet_SoftLeptptRel*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))","max(0,(hJet_SoftLeptdR*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))","max(0,(hJet_SoftLeptPt*(hJet_SoftLeptIdlooseMu==1||hJet_SoftLeptId95==1)))","hJet_cef","hJet_nconstituents","max(0,hJet_vtxPt)","max(0,hJet_vtxMass)","max(0,hJet_vtx3dL)","max(0,hJet_vtx3deL)","METnoPU.et","deltaPhi(METnoPU.phi,hJet_phi)"] |
| 277 |
|
useMET = True |
| 278 |
|
usePtRaw = True |
| 279 |
|
useRho25 = False |
| 286 |
|
# LIMIT |
| 287 |
|
[LimitGeneral] |
| 288 |
|
varName = CMS_vhbb_BDT_Zll |
| 289 |
< |
range = BDT output,10,-1,1 |
| 289 |
> |
range = BDT output,20,-1,1 |
| 290 |
|
|
| 291 |
|
systematicsnaming7TeV={'JER':'CMS_res_j','JES':'CMS_scale_j','beff':'CMS_eff_b','bmis':'CMS_fake_b_7TeV'} |
| 292 |
|
systematicsnaming8TeV={'JER':'CMS_res_j','JES':'CMS_scale_j','beff':'CMS_eff_b','bmis':'CMS_fake_b_8TeV'} |
| 846 |
|
ntBtag = !(hJet_csv[0] > <!ZHbb|CSVT!> | hJet_csv[1] > <!ZHbb|CSVT!>) |
| 847 |
|
twoCSV0 = hJet_csv[0] > 0. & hJet_csv[1] > 0. |
| 848 |
|
preselection = hJet_pt[0] > 20. & hJet_pt[1] > 20. & abs(hJet_eta[0]) < 2.4 & abs(hJet_eta[1]) < 2.4 & hbhe & hJet_puJetIdL[0] > 0. & hJet_puJetIdL[1] > 0. |
| 849 |
+ |
preselection_beforeReg = (hJet_pt[0]/hJet_regWeight[0]) > 20. & (hJet_pt[1]*hJet_regWeight[1]) > 20. & abs(hJet_eta[0]) < 2.4 & abs(hJet_eta[1]) < 2.4 & hbhe & hJet_puJetIdL[0] > 0. & hJet_puJetIdL[1] > 0. |
| 850 |
|
zWindow = V.mass > 75. & V.mass < 105 |
| 851 |
|
zVeto = (V.mass > 105 | V.mass < 75.) |
| 852 |
|
noAddJet = Sum$(aJet_pt > 20 & abs(aJet_eta) < 2.4) == 0 |
| 853 |
|
max1AddJet = Sum$(aJet_pt > 20 & abs(aJet_eta) < 2.4) < 2. |
| 854 |
+ |
vetoHMass_beforeReg = (HNoReg.mass < 90. | HNoReg.mass > 145.) |
| 855 |
|
vetoHMass = (H.mass < 90. | H.mass > 145.) |
| 856 |
|
pullAngle = deltaPullAngle < 10. & deltaPullAngle2 < 10. |
| 857 |
|
rTight = V.pt > 100. |
| 859 |
|
rMed = V.pt > 50. & V.pt < 100. |
| 860 |
|
upperHMass = H.mass < 250. |
| 861 |
|
looseHMass = H.mass > 80. & H.mass < 150. |
| 862 |
+ |
upperHMass_beforeReg = HNoReg.mass < 250. |
| 863 |
+ |
looseHMass_beforeReg = HNoReg.mass > 80. & HNoReg.mass < 150. |
| 864 |
|
dPhiVH = abs( HVdPhi ) > 2.9 |
| 865 |
|
looseBTag = !(hJet_csv[0] > <!ZHbb|CSVL!> & hJet_csv[1] > <!ZHbb|CSVL!>) & (hJet_csv[0] > <!ZHbb|CSVL!> | hJet_csv[1] > <!ZHbb|CSVL!>) |
| 866 |
|
|
| 889 |
|
TTbar = <!ZHbb|preselection!> & <!ZHbb|upperHMass!> & <!ZHbb|tlBtag!> & <!ZHbb|zVeto!> & H.pt > 100. |
| 890 |
|
Zbb = <!ZHbb|preselection!> & <!ZHbb|vetoHMass!> & <!ZHbb|upperHMass!> & <!ZHbb|zWindow!> & <!ZHbb|tlBtag!> & <!ZHbb|max1AddJet!> & <!ZHbb|dPhiVH!> |
| 891 |
|
|
| 892 |
+ |
ZLight_beforeReg: <!ZHbb|preselection_beforeReg!> & <!ZHbb|noAddJet!> & <!ZHbb|zWindow!> & <!ZHbb|ntBtag!> & <!ZHbb|upperHMass_beforeReg!> & <!ZHbb|twoCSV0!> & <!ZHbb|dPhiVH!> & <!ZHbb|maxBtag!> > <!ZHbb|CSVL!> |
| 893 |
+ |
TTbar_beforeReg = <!ZHbb|preselection_beforeReg!> & <!ZHbb|upperHMass_beforeReg!> & <!ZHbb|tlBtag!> & <!ZHbb|zVeto!> & HNoReg.pt > 100. |
| 894 |
+ |
Zbb_beforeReg = <!ZHbb|preselection_beforeReg!> & <!ZHbb|vetoHMass_beforeReg!> & <!ZHbb|upperHMass_beforeReg!> & <!ZHbb|zWindow!> & <!ZHbb|tlBtag!> & <!ZHbb|max1AddJet!> & <!ZHbb|dPhiVH!> |
| 895 |
+ |
|
| 896 |
|
ZLight_HighPt = <!Cuts|ZLight!> & <!ZHbb|rTight!> |
| 897 |
+ |
ZLight_HighPt_beforeReg = <!Cuts|ZLight_beforeReg!> & <!ZHbb|rTight!> |
| 898 |
|
TTbar_HighPt = <!Cuts|TTbar!> & <!ZHbb|rTight!> |
| 899 |
|
Zbb_HighPt = <!Cuts|Zbb!> & <!ZHbb|rTight!> |
| 900 |
|
|
| 901 |
|
ZLight_LowPt = <!Cuts|ZLight!> & <!ZHbb|rMed!> |
| 902 |
+ |
ZLight_LowPt_beforeReg = <!Cuts|ZLight_beforeReg!> & <!ZHbb|rMed!> |
| 903 |
|
TTbar_LowPt = <!Cuts|TTbar!> & <!ZHbb|rMed!> |
| 904 |
|
Zbb_LowPt = <!Cuts|Zbb!> & <!ZHbb|rMed!> |
| 905 |
|
|
| 945 |
|
# other things: |
| 946 |
|
|
| 947 |
|
[Weights] |
| 948 |
< |
weightF: (PUweight*weightTrig2012A*lheWeight) |
| 948 |
> |
weightF: (PUweight*weightTrig2012*lheWeight) |
| 949 |
|
|
| 950 |
|
# LocalWords: Mjj |
