| 1 |
peller |
1.1 |
[Analysis]
|
| 2 |
|
|
tag = 8TeV
|
| 3 |
|
|
|
| 4 |
|
|
# If this Flag is used then 50% of MC samples will go in training and the rest will be scaled up by a factor of 2.
|
| 5 |
|
|
# if you train on 50% but then set it to 'Fakse' while writing DC, then you get your biased training samples in the limit!!!
|
| 6 |
|
|
# be sure to have this right !!!
|
| 7 |
|
|
# >>>>>>>>>>>>>>>>>>>
|
| 8 |
|
|
# >>>>>>>>>>>>>>>>>>>
|
| 9 |
|
|
TrainFlag = True
|
| 10 |
|
|
#Mjj-> False
|
| 11 |
|
|
# <<<<<<<<<<<<<<<<<<<
|
| 12 |
|
|
# <<<<<<<<<<<<<<<<<<<
|
| 13 |
|
|
|
| 14 |
|
|
[VHbbNameSpace]
|
| 15 |
|
|
library=../interface/VHbbNameSpace_h.so
|
| 16 |
|
|
|
| 17 |
|
|
###########################################################################################
|
| 18 |
|
|
# REGRESSION
|
| 19 |
|
|
[Regression]
|
| 20 |
|
|
#regWeight = ../data/MVA_BDT_REG_Jun16.weights.xml
|
| 21 |
|
|
#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'}
|
| 22 |
|
|
#regVars = ["Jet_pt","Jet_eta","Jet_e","Jet_ptRaw", "Jet_chf","Jet_nconstituents", "Jet_vtxPt", "Jet_vtx3dL", "Jet_vtx3deL", "Jet_MET_dPhi", "METet"]
|
| 23 |
|
|
regWeight = ../data/MVA_BDT_REG_Matt.weights.xml
|
| 24 |
|
|
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"}
|
| 25 |
|
|
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)"]
|
| 26 |
|
|
useMET = True
|
| 27 |
|
|
usePtRaw = True
|
| 28 |
|
|
useRho25 = False
|
| 29 |
|
|
useMt = False
|
| 30 |
|
|
|
| 31 |
|
|
[BTagReshaping]
|
| 32 |
|
|
library=<!Directories|vhbbpath!>/interface/BTagReshaping_h.so
|
| 33 |
|
|
|
| 34 |
|
|
[VHbbNameSpace]
|
| 35 |
|
|
library=<!Directories|vhbbpath!>/interface/VHbbNameSpace_h.so
|
| 36 |
|
|
|
| 37 |
|
|
[AngularLike]
|
| 38 |
|
|
backgrounds = ['Zbb','Zudsg','ttbar','ZZ']
|
| 39 |
|
|
yields = {'Zbb':584.8,'Zudsg':881.61,'ttbar':135.7,'ZZ':44.04,'Sig':5.13}
|
| 40 |
|
|
|
| 41 |
|
|
##############################################################################################
|
| 42 |
|
|
# SYSTEMATICS:
|
| 43 |
|
|
|
| 44 |
|
|
[systematics]
|
| 45 |
|
|
systematics: Nominal JER_up JER_down JES_up JES_down beff_up beff_down bmis_up bmis_down beff1_up beff1_down
|
| 46 |
|
|
|
| 47 |
|
|
|
| 48 |
|
|
[treeVarsICHEP]
|
| 49 |
|
|
Nominal: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 50 |
|
|
JER_up: H_JER.mass_up H_JER.pt_up HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt_JER_up[0] hJet_pt_JER_up[1] deltaPullAngle deltaPullAngle2
|
| 51 |
|
|
JER_down: H_JER.mass_down H_JER.pt_down HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt_JER_down[0] hJet_pt_JER_down[1] deltaPullAngle deltaPullAngle2
|
| 52 |
|
|
JES_up: H_JES.mass_up H_JES.pt_up HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt_JES_up[0] hJet_pt_JES_up[1] deltaPullAngle deltaPullAngle2
|
| 53 |
|
|
JES_down: H_JES.mass_down H_JES.pt_down HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt_JES_down[0] hJet_pt_JES_down[1] deltaPullAngle deltaPullAngle2
|
| 54 |
|
|
beff_up: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csvUp[0],hJet_csvUp[1]) min(hJet_csvUp[0],hJet_csvUp[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 55 |
|
|
beff_down: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csvDown[0],hJet_csvDown[1]) min(hJet_csvDown[0],hJet_csvDown[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 56 |
|
|
bmis_up: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csvFUp[0],hJet_csvFUp[1]) min(hJet_csvFUp[0],hJet_csvFUp[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 57 |
|
|
bmis_down: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csvFDown[0],hJet_csvFDown[1]) min(hJet_csvFDown[0],hJet_csvFDown[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 58 |
|
|
beff1_up: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csv_upBC4p[0],hJet_csv_upBC4p[1]) min(hJet_csv_upBC4p[0],hJet_csv_upBC4p[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 59 |
|
|
beff1_down: H.mass H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csv_downBC4p[0],hJet_csv_downBC4p[1]) min(hJet_csv_downBC4p[0],hJet_csv_downBC4p[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 60 |
|
|
|
| 61 |
|
|
[treeVarsNoMass]
|
| 62 |
|
|
Nominal: H.pt HVdPhi V.mass V.pt H.dEta max(hJet_csv[0],hJet_csv[1]) min(hJet_csv[0],hJet_csv[1]) Sum$(aJet_pt>20.&abs(aJet_eta)<2.4) H.dR hJet_pt[0] hJet_pt[1] deltaPullAngle deltaPullAngle2
|
| 63 |
|
|
|
| 64 |
|
|
|
| 65 |
|
|
###########################################################################################
|
| 66 |
|
|
# other things:
|
| 67 |
|
|
|
| 68 |
|
|
[Weights]
|
| 69 |
|
|
#5fb:
|
| 70 |
|
|
#weightF: (PUweightAB*weightTrig2012AB*lheWeight)
|
| 71 |
|
|
#weightF_sys_UP: (PUweightP*weightTrig2012AB*lheWeight)
|
| 72 |
|
|
#weightF_sys_DOWN: (PUweightM*weightTrig2012AB*lheWeight)
|
| 73 |
|
|
|
| 74 |
|
|
weightF: (PUweight*weightTrig2012*lheWeight)
|
| 75 |
|
|
weightF_sys_UP: (PUweightP*weightTrig2012*lheWeight)
|
| 76 |
|
|
weightF_sys_DOWN: (PUweightM*weightTrig2012*lheWeight)
|
| 77 |
|
|
|
| 78 |
|
|
#weightF: (PUweight*lheWeight)
|
