def event_filter(process, mode, settings, output_commands, **kwargs):
if mode:
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
setup_event_filter(process,
input_is_miniaod=settings.is_miniaod,
mode=mode,
**kwargs)
def event_filter(process, mode, settings, output_commands, **kwargs):
if mode:
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
# for temp including events that don't pass MET trigger and MET lower than 150
#setup_event_filter(process, input_is_miniaod=settings.is_miniaod, mode=mode, trigger_filter = False, event_filter = False, event_filter_require_vertex = False, **kwargs)
setup_event_filter(process,
input_is_miniaod=settings.is_miniaod,
mode=mode,
event_filter_require_vertex=True,
**kwargs)
process.tt = cms.EDAnalyzer(
'TrackingTreer',
process.jmtNtupleFillerMiniAOD,
track_cut_level=cms.int32(
0
), # -1 = all, 0 = pt & pix & strip, 1 = 0 + min_r, 2 = 1 + nsigmadxybs
)
process.tt.track_ref_getter.tracks_maps_srcs = []
process.p = cms.Path(process.tt)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
if use_btag_triggers:
setup_event_filter(process,
input_is_miniaod=True,
mode='bjets OR displaced dijet veto HT novtx',
event_filter_jes_mult=0)
elif use_lepton_triggers:
setup_event_filter(process,
input_is_miniaod=True,
mode='leptons',
event_filter_jes_mult=0)
elif use_displ_lepton_triggers:
setup_event_filter(process,
input_is_miniaod=True,
mode='displaced leptons',
event_filter_jes_mult=0)
else:
setup_event_filter(process,
input_is_miniaod=True,
mode='jets only novtx',
process.tt = cms.EDAnalyzer(
'TrackingTreer',
process.jmtNtupleFillerMiniAOD,
track_cut_level=cms.int32(
1
), # -1 = all, 0 = pt & pix & strip, 1 = 0 + min_r, 2 = 1 + nsigmadxybs
)
process.tt.track_ref_getter.tracks_maps_srcs = []
process.p = cms.Path(process.tt)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
if use_btag_triggers:
setup_event_filter(process,
input_is_miniaod=True,
mode='bjets OR displaced dijet veto HT novtx',
event_filter_jes_mult=0)
elif use_MET_triggers:
setup_event_filter(process,
input_is_miniaod=True,
mode='met only',
event_filter_jes_mult=0,
event_filter_require_vertex=False)
#setup_event_filter(process, input_is_miniaod=True, mode='jets only novtx', event_filter_jes_mult=0)
else:
setup_event_filter(process,
input_is_miniaod=True,
mode='jets only novtx',
event_filter_jes_mult=0)
ReferencedTagsTaskAdder(process)('p')
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import *
settings = CMSSWSettings()
settings.is_mc = True
settings.is_miniaod = True
settings.cross = '' # 2017to2018' # 2017to2017p8'
sample_files(process, 'qcdht1000_2017', 'miniaod')
geometry_etc(process, settings)
tfileservice(process, 'filtercheck.root')
cmssw_from_argv(process)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
sef = lambda *a, **kwa: setup_event_filter(
process, *a, input_is_miniaod=True, **kwa)
sef('pTrigger', mode='trigger jets only')
sef('pTriggerBjets', mode='trigger bjets only', name_ex='bjets')
sef('pTriggerDispDijet',
mode='trigger displaced dijet only',
name_ex='displaced_dijet')
sef('pTriggerOR',
mode='trigger HT OR bjets OR displaced dijet',
name_ex='HT_OR_bjets_OR_displaced_dijet')
sef('pJets', mode='jets only novtx', name_ex='NoVtx'
) # be sure to generate a different name for each subsequent use
sef('pNoJESUp',
mode='jets only novtx',
name_ex='NoJESUp',
event_filter_jes_mult=0)
#sef('pFull', mode = 'jets only', name_ex = 'Full') # uncomment to get efficiency of ntuple-level vertex filter
def event_filter(process, mode, settings, output_commands, **kwargs):
if mode:
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
setup_event_filter(process, input_is_miniaod=settings.is_miniaod, mode=mode, **kwargs)
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import *
settings = CMSSWSettings()
settings.is_mc = True
settings.is_miniaod = True
settings.cross = '' # 2017to2018' # 2017to2017p8'
sample_files(process, 'qcdht1000_2017', 'miniaod')
geometry_etc(process, settings)
tfileservice(process, 'filtercheck.root')
cmssw_from_argv(process)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
sef = lambda *a,**kwa: setup_event_filter(process, *a, input_is_miniaod=True, **kwa)
sef('pTrigger', mode = 'trigger jets only')
sef('pJets', mode = 'jets only novtx', name_ex = 'NoVtx') # be sure to generate a different name for each subsequent use
sef('pNoJESUp', mode = 'jets only novtx', name_ex = 'NoJESUp', event_filter_jes_mult = 0)
#sef('pFull', mode = 'jets only', name_ex = 'Full') # uncomment to get efficiency of ntuple-level vertex filter
if len(process.mfvTriggerFilter.HLTPaths) > 1:
for x in process.mfvTriggerFilter.HLTPaths:
filt_name = ''.join(x.split('_')[1:-1])
sef('p%s' % filt_name, name_ex=x)
getattr(process, filt_name).HLTPaths = [x]
import JMTucker.Tools.SimpleTriggerEfficiency_cfi as SimpleTriggerEfficiency
SimpleTriggerEfficiency.setup_endpath(process)
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
process.TrackerMapper = cms.EDAnalyzer('TrackerMapper',
track_src = cms.InputTag('jmtUnpackedCandidateTracks'),
heavy_flavor_src = cms.InputTag(''),
beamspot_src = cms.InputTag('offlineBeamSpot'),
primary_vertex_src = cms.InputTag('goodOfflinePrimaryVertices'),
weight_src = cms.InputTag('jmtWeightMiniAOD'),
use_duplicateMerge = cms.int32(-1),
old_stlayers_cut = cms.bool(False),
)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
if use_btag_triggers :
event_filter = setup_event_filter(process,
path_name = '',
trigger_filter = 'bjets OR displaced dijet veto HT',
event_filter = 'bjets OR displaced dijet veto HT',
event_filter_jes_mult = 0,
event_filter_require_vertex = False,
input_is_miniaod = True)
else :
event_filter = setup_event_filter(process,
path_name = '',
trigger_filter = 'jets only',
event_filter = 'jets only',
event_filter_jes_mult = 0,
event_filter_require_vertex = False,
input_is_miniaod = True)
common = cms.Sequence(event_filter * process.goodOfflinePrimaryVertices * process.jmtUnpackedCandidateTracks * process.jmtWeightMiniAOD)
process.load('JMTucker.Tools.WeightProducer_cfi')
process.selectedPatJets.src = 'updatedJetsMiniAOD'
process.selectedPatJets.cut = process.jtupleParams.jetCut
process.tt = cms.EDAnalyzer('TrackingTreer',
process.jmtNtupleFillerMiniAOD,
track_sel = cms.bool(True),
)
process.tt.track_ref_getter.tracks_maps_srcs = []
process.p = cms.Path(process.tt)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
setup_event_filter(process, input_is_miniaod=True, mode='jets only novtx', event_filter_jes_mult=0)
ReferencedTagsTaskAdder(process)('p')
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
from JMTucker.Tools.MetaSubmitter import *
import JMTucker.Tools.Samples as Samples
from JMTucker.Tools.Year import year
if year == 2017:
samples = Samples.qcd_samples_2017 + Samples.data_samples_2017
elif year == 2018:
samples = Samples.qcd_samples_2018 + Samples.data_samples_2018
set_splitting(samples, 'miniaod', 'default', json_path('ana_2017p8_1pc.json'), 16)
'JMTBTagEfficiency',
weight_src=cms.InputTag('jmtWeightMiniAOD'),
jets_src=cms.InputTag('selectedPatJets'),
jet_pt_min=cms.double(20),
old=cms.bool(False),
)
process.p = cms.Path(process.JMTBTagEfficiency)
if year == 2017:
process.JMTBTagEfficiencyOld = process.JMTBTagEfficiency.clone(old=True)
process.p *= process.JMTBTagEfficiencyOld
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
setup_event_filter(process,
input_is_miniaod=True,
mode='jets only novtx',
event_filter_jes_mult=0)
ReferencedTagsTaskAdder(process)('p')
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
from JMTucker.Tools.MetaSubmitter import *
import JMTucker.Tools.Samples as Samples
if year == 2017:
samples = Samples.qcd_samples_2017 + Samples.ttbar_samples_2017
elif year == 2018:
samples = Samples.qcd_samples_2018 + Samples.ttbar_samples_2018
set_splitting(samples, 'miniaod', 'default', default_files_per=16)
process.TrackerMapper = cms.EDAnalyzer('TrackerMapper',
track_src = cms.InputTag('jmtUnpackedCandidateTracks'),
heavy_flavor_src = cms.InputTag(''),
beamspot_src = cms.InputTag('offlineBeamSpot'),
primary_vertex_src = cms.InputTag('goodOfflinePrimaryVertices'),
weight_src = cms.InputTag('jmtWeightMiniAOD'),
use_duplicateMerge = cms.int32(-1),
old_stlayers_cut = cms.bool(False),
)
from JMTucker.MFVNeutralino.EventFilter import setup_event_filter
jetsOnly = setup_event_filter(process,
path_name = '',
trigger_filter = 'jets only',
event_filter = 'jets only',
event_filter_jes_mult = 0,
event_filter_require_vertex = False,
input_is_miniaod = True)
common = cms.Sequence(jetsOnly * process.goodOfflinePrimaryVertices * process.jmtUnpackedCandidateTracks * process.jmtWeightMiniAOD)
if False:
process.load('JMTucker.Tools.RescaledTracks_cfi')
process.jmtRescaledTracks.tracks_src = 'jmtUnpackedCandidateTracks'
common *= process.jmtRescaledTracks
process.TrackerMapper.track_src = 'jmtRescaledTracks'
process.p = cms.Path(common * process.TrackerMapper)
for name, filt in ('LightFlavor', process.lightFlavor), ('HeavyFlavor', process.heavyFlavor): #, ('BFlavor', process.bFlavor), ('DisplacedGenPV', process.displacedGenPV):
