Exemplo n.º 1
0
    'wenu': selectors.wenuall
}

data_sph = ['monoph', 'efake', 'hfake', 'hfakeUp', 'hfakeDown', 'purity', 'purityUp', 'purityDown', 'gjets']
data_smu = ['dimu', 'monomu', 'elmu']
data_sel = ['diel', 'monoel', 'eefake']
mc_cand = ['monoph']
mc_qcd = ['hfake', 'hfakeUp', 'hfakeDown', 'purity', 'purityUp', 'purityDown', 'gjets'] 
mc_sig = ['monoph', 'signalRaw']
mc_lep = ['monomu', 'monoel']
mc_dilep = ['dimu', 'diel', 'elmu']
mc_vgcand = [(region, selectors.kfactor(defaults[region])) for region in mc_cand]
mc_vglep = [(region, selectors.kfactor(defaults[region])) for region in mc_lep]
mc_vgdilep = [(region, selectors.kfactor(defaults[region])) for region in mc_dilep]
mc_gj = [('monoph', selectors.kfactor(selectors.gjSmeared)), ('purity', selectors.kfactor(selectors.purity))]
mc_wlnu = [(region, selectors.wlnu(defaults[region])) for region in mc_cand] + ['wenu']

sphLumi = allsamples['sph-d3'].lumi + allsamples['sph-d4'].lumi
haloNorms = [ 5.9 * allsamples[sph].lumi / sphLumi for sph in ['sph-d3', 'sph-d4'] ]

selectors = {
    # Data
    'sph-d3': data_sph + [('halo', selectors.haloCSC(haloNorms[0]))
                          ,('haloUp', selectors.haloMIP(haloNorms[0]))
                          ,('haloDown', selectors.haloSieie(haloNorms[0]))
                          ],
    'sph-d4': data_sph + [('halo', selectors.haloCSC(haloNorms[1]))
                          ,('haloUp', selectors.haloMIP(haloNorms[1]))
                          ,('haloDown', selectors.haloSieie(haloNorms[1]))
                          ],
    'smu-d3': data_smu,
Exemplo n.º 2
0
    'diel': selectors.dielectron,
    'monoel': selectors.monoelectron,
    'elmu': selectors.oppflavor,
    'eefake': selectors.zee
}

data_sph = ['monoph', 'efake', 'hfake', 'hfakeUp', 'hfakeDown', 'gjets']
data_smu = ['dimu', 'monomu', 'elmu']
data_sel = ['diel', 'monoel', 'eefake']
mc_cand = ['monoph']
mc_lep = ['monomu', 'monoel']
mc_dilep = ['dimu', 'diel', 'elmu']
mc_vgcand = [(region, selectors.kfactor(defaults[region])) for region in ['monoph']]
mc_vglep = [(region, selectors.kfactor(defaults[region])) for region in ['monomu', 'monoel']]
mc_vgdilep = [(region, selectors.kfactor(defaults[region])) for region in ['dimu', 'diel', 'elmu']]
mc_wlnu = [(region, selectors.wlnu(defaults[region])) for region in ['monoph']]

selectors = {
    # Data
    'sph-d3': data_sph,
    'sph-d4': data_sph,
    'smu-d3': data_smu,
    'smu-d4': data_smu,
    'sel-d3': data_sel,
    'sel-d4': data_sel,
    # MC for signal region
    'znng-130': mc_vgcand,
    'wnlg-130': mc_vgcand + mc_vglep,
    'g-40': mc_cand,
    'g-100': mc_cand,
    'g-200': mc_cand,
Exemplo n.º 3
0
data_15 = []
data_sph = ['monoph', 'efake', 'hfake', 'hfakeUp', 'hfakeDown', 'purity', 'purityUp', 'purityDown', 'lowmt', 'lowmtEfake', 'lowmtHfake', 'gjets', 'dimu', 'diel', 'monomu', 'monoel'] + noncollision
data_smu = ['dimu', 'monomu', 'monomuHfake', 'elmu', 'zmmJets'] # are SinglePhoton triggers in this PD? (do the samples know about them, obviously they are not used to define it)
data_sel = ['diel', 'monoel', 'monoelHfake', 'eefake', 'zeeJets'] # are SinglePhoton triggers in this PD? (do the samples know about them, obviously they are not used to define it)
mc_cand = ['monoph', 'purity']
mc_qcd = ['hfake', 'hfakeUp', 'hfakeDown', 'purity', 'purityUp', 'purityDown', 'gjets'] 
mc_sig = ['monoph', 'purity', 'signalRaw']
mc_lep = ['monomu', 'monoel']
mc_dilep = ['dimu', 'diel', 'elmu', 'zmmJets', 'zeeJets']
mc_vgcand = [(region, selectors.kfactor(defaults[region])) for region in mc_cand]
mc_vglep = [(region, selectors.kfactor(defaults[region])) for region in mc_lep]
mc_vgdilep = [(region, selectors.kfactor(defaults[region])) for region in mc_dilep]
#mc_gj = [('raw', selectors.kfactor(defaults['monoph'])), ('monoph', selectors.kfactor(selectors.gjSmeared)), ('purity', selectors.kfactor(selectors.purity))]
#mc_gj = [('raw', selectors.kfactor(defaults['monoph'])), ('monoph', selectors.kfactor(defaults['monoph'])), ('purity', selectors.kfactor(selectors.purity))]
mc_gj = [(region, selectors.kfactor(defaults[region])) for region in  mc_qcd + mc_cand]
mc_wlnu = [(region, selectors.wlnu(defaults[region])) for region in mc_cand] + ['wenu', 'zmmJets', 'zeeJets']
mc_lowmt = ['lowmt']
mc_vglowmt = [(region, selectors.kfactor(defaults[region])) for region in mc_lowmt]

selectors = {
    # Data 2016
    'sph-16b-r': data_sph,
    'sph-16c-r': data_sph,
    'sph-16d-r': data_sph,
    'sph-16e-r': data_sph,
    'sph-16f-r': data_sph,
    'sph-16g-r': data_sph,
    'sph-16h': data_sph,
    # MC for signal region
    'znng-130': mc_vgcand + mc_vglowmt,
    'wnlg-130': mc_vgcand + mc_vglep + mc_vglowmt,