'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',