'diel': selectors.dielectron,
'monoel': selectors.monoelectron,
'elmu': selectors.oppflavor,
'eefake': selectors.zee,
'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]))
'gjets': selectors.gammaJets,
'dimu': selectors.dimuon,
'monomu': selectors.monomuon,
'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,
noncollision = [
('halo', selectors.haloMIP(haloNorms[0])),
('haloUp', selectors.haloCSC(haloNorms[0])),
('spikeSieie', selectors.spikeSieie(spikeNorms[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,