import sys
from JMTucker.Tools.BasicAnalyzer_cfg import cms, global_tag, process, file_event_from_argv
from JMTucker.Tools.MiniAOD_cfg import which_global_tag
file_event_from_argv(process)
process.TFileService.fileName = 'trigeff.root'
global_tag(process, which_global_tag(True)) # needed for emulator that checks L1 bits
from JMTucker.MFVNeutralino.TriggerFilter import *
setup_trigger_filter_soup(process, 'pmfvfilt')
setup_trigger_filter_soup(process, 'ponlyhtpt', 'tfonlyhtpt')
setup_trigger_filter_soup(process, 'ponlyht', 'tfonlyht')
setup_trigger_filter_soup(process, 'ponlyht8', 'tfonlyht8')
setup_trigger_filter_soup(process, 'ponlyht9', 'tfonlyht9')
setup_trigger_filter(process, 'pemuht8', need_pat=True)
setup_trigger_filter(process, 'pemuht9', 'emuht900', need_pat=True)
process.emuht900.return_actual = False
process.emuht900.return_ht900 = True
#process.options.wantSummary = True
process.tfonlyhtpt.HLTPaths = [
'HLT_PFHT650_v*',
'HLT_PFHT800_v*',
'HLT_PFHT900_v*',
'HLT_PFHT550_4Jet_v*',
'HLT_PFHT450_SixJet40_v*',
'HLT_PFHT400_SixJet30_v*',
]
files = ['/store/user/tucker/pick1vtx/qcdht1000/pick_%i.root' % i for i in range(50)]
elif args.sample == 'qcdht1500':
files = ['/store/user/tucker/pick1vtx/qcdht1500/161104_153718/0000/pick_%i.root' % i for i in range(1,15) if i != 9]
files += ['/store/user/tucker/pick1vtx/qcdht1500/qcdht1500_%i.root' % i for i in (0,1,2)]
elif args.sample == 'qcdht2000':
files = ['/store/user/tucker/pick1vtx/qcdht2000/pick_%i.root' % i for i in range(23)]
else:
raise ValueError('bad sample %s' % args.sample)
####
process = basic_process('Overlay')
process.source.fileNames = files
process.maxEvents.input = args.max_events
report_every(process, 1000000 if args.batch else 100)
geometry_etc(process, which_global_tag(not args.is_data))
tfileservice(process, args.out_fn)
random_service(process, {'mfvVertices': 12179 + args.which_event,
'mfvOverlayTracks': 12180 + args.which_event})
process.load('CommonTools.ParticleFlow.goodOfflinePrimaryVertices_cfi')
process.load('JMTucker.MFVNeutralino.Vertexer_cfi')
process.mfvVertices.histos = False
process.mfvVertices.verbose = args.debug
process.mfvVertices.use_second_tracks = True
process.mfvVertices.second_track_src = 'mfvOverlayTracks'
if not args.rest_of_event:
process.mfvVertices.disregard_event = True
process.mfvOverlayTracks = cms.EDFilter('MFVOverlayVertexTracks',
#!/usr/bin/env python
import sys
from JMTucker.Tools.BasicAnalyzer_cfg import *
from JMTucker.Tools.MiniAOD_cfg import which_global_tag
process.TFileService.fileName = 'prescales.root'
global_tag(process, which_global_tag(is_mc=False))
process.source.fileNames = ['/store/data/Run2015D/JetHT/AOD/PromptReco-v4/000/260/627/00000/78D8E6A7-6484-E511-89B4-02163E0134F6.root']
from JMTucker.Tools.L1GtUtils_cff import l1GtUtilsTags
add_analyzer(process, 'MFVTriggerPrescales', l1GtUtilsTags)
process.maxEvents.input = -1
if __name__ == '__main__' and hasattr(sys, 'argv') and 'submit' in sys.argv:
import JMTucker.Tools.Samples as Samples
from JMTucker.Tools.CRAB3Submitter import CRABSubmitter
cs = CRABSubmitter('TriggerPrescalesV0',
splitting = 'FileBased',
units_per_job = 20,
total_units = -1,
)
cs.submit_all(Samples.data_samples)
#!/usr/bin/env python
from JMTucker.Tools.MiniAOD_cfg import cms, which_global_tag
from JMTucker.Tools.CMSSWTools import *
is_mc = True
version = 'v9'
batch_name = 'Ntuple' + version.upper()
####
process = basic_process('Ntuple')
report_every(process, 1000000)
registration_warnings(process)
geometry_etc(process, which_global_tag(is_mc))
random_service(process, {'mfvVertices': 1222})
tfileservice(process, 'vertex_histos.root')
output_file(process, 'ntuple.root', [
'drop *',
'keep *_mcStat_*_*',
'keep MFVVertexAuxs_mfvVerticesAux_*_*',
'keep MFVEvent_mfvEvent__*',
])
process.load('CommonTools.ParticleFlow.goodOfflinePrimaryVertices_cfi')
process.load('PhysicsTools.PatAlgos.selectionLayer1.jetSelector_cfi')
process.load('PhysicsTools.PatAlgos.selectionLayer1.muonSelector_cfi')
process.load('PhysicsTools.PatAlgos.selectionLayer1.electronSelector_cfi')
process.load('JMTucker.Tools.PATTupleSelection_cfi')
process.load('JMTucker.MFVNeutralino.Vertexer_cff')
process.load('JMTucker.MFVNeutralino.EventProducer_cfi')
import sys, os
from JMTucker.Tools.BasicAnalyzer_cfg import *
from JMTucker.Tools.MiniAOD_cfg import which_global_tag
simple = False
process.source.fileNames = ['file:ntuple.root']
process.source.noEventSort = cms.untracked.bool(True)
file_event_from_argv(process)
geometry_etc(process, which_global_tag(True))
del process.TFileService
process.MessageLogger.cerr.FwkReport.reportEvery = 1
process.load('JMTucker.MFVNeutralino.VertexSelector_cfi')
process.load('JMTucker.MFVNeutralino.AnalysisCuts_cfi')
printer = cms.EDAnalyzer('MFVPrinter',
print_vertex = cms.bool(False),
print_event = cms.bool(False),
print_vertex_aux = cms.bool(False),
vertex_src = cms.InputTag('mfvVertices'),
event_src = cms.InputTag('mfvEvent'),
vertex_aux_src = cms.InputTag('mfvVerticesAux'),
)
if simple:
process.justPrint = printer.clone(print_event = True, print_vertex_aux = True, vertex_aux_src = 'mfvSelectedVerticesTight')
process.p = cms.Path(process.mfvSelectedVerticesSeq * process.mfvAnalysisCuts * process.justPrint)
else:
process.printRecoVertices = printer.clone(print_vertex = True)
