'name': 'lumi_13TeV_Ghosts',
'type': 'lnN',
'samples': dict((skey, '1.004') for skey in mc if skey not in ['WZ']),
'cuts' : fitcuts
}
#### Theoretical Systematics
# Scale
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
nuisances['QCDscale_VH'] = {
'name' : 'QCDscale_VH',
'samples' : {
'WH_hww' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','WH','125.09','scale','sm'),
'ZH_hww' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ZH','125.09','scale','sm'),
'WH_htt' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','WH','125.09','scale','sm'),
'ZH_htt' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ZH','125.09','scale','sm')
},
'type' : 'lnN',
'cuts' : fitcuts
}
nuisances['QCDscale_ggZH'] = {
'name' : 'QCDscale_ggZH',
'samples' : {
'ggZH_hww': HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ggZH','125.09','scale','sm'),
},
'type' : 'lnN',
'cuts' : fitcuts
'samples' : {
'WWewk' : '1.11',
},
'type' : 'lnN'
}
'''
#### QCD scale uncertainties for Higgs signals other than ggH
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
nuisances['QCDscale_ggH'] = {
'name' : 'QCDscale_ggH',
'samples' : {
'H_htt' : HiggsXS.GetHiggsProdXSNP('YR4prel','13TeV','ggH','125.0','scale','sm'),
},
'type' : 'lnN',
}
nuisances['QCDscale_qqH'] = {
'name' : 'QCDscale_qqH',
'samples' : {
'qqH_hww' : HiggsXS.GetHiggsProdXSNP('YR4prel','13TeV','vbfH','125.0','scale','sm'),
'qqH_htt' : HiggsXS.GetHiggsProdXSNP('YR4prel','13TeV','vbfH','125.0','scale','sm'),
},
'type' : 'lnN',
}
nuisances['QCDscale_VH'] = {
},
}
nuisances['WZRate'] = {
'name': 'WZRate',
'type': 'lnN',
'samples': {
'WZgS_H': '1.16',
},
}
###### pdf uncertainties
## All shifts on signal have to be renormalized at the datacard generation step
valuesggh = HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'ggH', '125.09',
'pdf', 'sm')
valuesggzh = HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'ggZH', '125.09',
'pdf', 'sm')
valuesbbh = HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'bbH', '125.09',
'pdf', 'sm')
nuisances['pdf_Higgs_gg'] = {
'name': 'pdf_Higgs_gg',
'samples': {
'ggH_hww': valuesggh,
'ggH_htt': valuesggh,
'ggZH_hww': valuesggzh,
'bbH_hww': valuesbbh
},
'type': 'lnN',
}
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
#nuisances['QCDscale_qqH'] = {
# 'name': 'QCDscale_qqH',
# 'samples' : {
# 'qqH_htt': HiggsXS.GetHiggsProdXSNP('YR4','13TeV','vbfH','125.09','scale','sm')
# },
# 'type': 'lnN'
#}
nuisances['QCDscale_VH'] = {
'name': 'QCDscale_VH',
'samples': {
'WH_hww':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'WH', '125.09', 'scale',
'sm'),
#'ZH_hww' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ZH','125.09','scale','sm'),
'WH_htt':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'WH', '125.09', 'scale',
'sm'),
'ZH_htt':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'ZH', '125.09', 'scale', 'sm')
},
'type': 'lnN',
}
#nuisances['QCDscale_ggZH'] = {
# 'name' : 'QCDscale_ggZH',
# 'samples' : {
# 'ggZH_hww': HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ggZH','125.09','scale','sm'),
# },
'WH_htt': '1.025',
},
'type': 'lnN',
}
# Theoritical Systematics
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
nuisances['QCDscale_VH'] = {
'name': 'QCDscale_VH',
'samples': {
'WH_hww':
HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'WH', '125.0', 'scale',
'sm'),
'ZH_hww':
HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'ZH', '125.0', 'scale',
'sm'),
'WH_htt':
HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'WH', '125.0', 'scale',
'sm'),
# 'ZH_htt' : HiggsXS.GetHiggsProdXSNP('YR4prel','13TeV','ZH','125.0','scale','sm'),
},
'type': 'lnN',
}
#nuisances['QCDscale_ZH_wh3l'] = {
# 'name' : 'QCDscale_ZH',
# 'samples' : {
# 'ZH_hww' : HiggsXS.GetHiggsProdXSNP('YR4prel','13TeV','ZH','125.0','scale','sm'),
'std_vector_LHE_weight[4]/std_vector_LHE_weight[0]'
],
}
}
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
####################QCD scale uncertainties for Higgs signals other than ggH
nuisances['QCDscale_ggH'] = {
'name': 'QCDscale_ggH',
'samples': {
'ggH_hww':
HiggsXS.GetHiggsProdXSNP(
'YR4prel', '13TeV', 'ggH', '125.0', 'scale',
'sm'), #-> already included in jet bin migration QCD uncertainty?
'ggH_htt':
HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'ggH', '125.0', 'scale',
'sm'),
'H_htt':
HiggsXS.GetHiggsProdXSNP('YR4prel', '13TeV', 'ggH', '125.0', 'scale',
'sm'),
},
'type': 'lnN',
}
nuisances['QCDscale_qqH'] = {
'name': 'QCDscale_qqH',
'samples': {
'qqH_hww':
},
'type': 'lnN',
}
for name in sampleNames:
nuisances['lumi']['samples'].update({name: '1.025'})
# Theoritical Systematics
from LatinoAnalysis.Tools.HiggsXSection import *
HiggsXS = HiggsXSection()
nuisances['QCDscale_VH'] = {
'name': 'QCDscale_VH',
'samples': {
'WH_hww':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'WH', '125.09', 'scale',
'sm'),
'ZH_hww':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'ZH', '125.09', 'scale',
'sm'),
'WH_htt':
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'WH', '125.09', 'scale',
'sm'),
# 'ZH_htt' : HiggsXS.GetHiggsProdXSNP('YR4','13TeV','ZH','125.09','scale','sm'),
},
'type': 'lnN',
}
for name in sampleNames:
if 'WH' in name:
nuisances['QCDscale_VH']['samples'].update({
name:
HiggsXS.GetHiggsProdXSNP('YR4', '13TeV', 'WH', '125.09', 'scale',
