b_discriminator_tag=cms.double(0.97),
njets=cms.uint32(njets),
nbjets=cms.uint32(nbjets),
tau=cms.double(1),
sig_theta=cms.double(0.2),
sig_phi=cms.double(0.2),
)
process.mfvMovedTree = cms.EDAnalyzer(
"MFVMovedTracksTreer",
event_src=cms.InputTag("mfvEvent"),
weight_src=cms.InputTag("mfvWeight"),
mover_src=cms.string("mfvMovedTracks"),
vertices_src=cms.InputTag("mfvVerticesAux"),
max_dist2move=cms.double(0.02),
apply_presel=cms.bool(True),
njets_req=cms.uint32(njets),
nbjets_req=cms.uint32(nbjets),
for_mctruth=cms.bool(False),
)
process.p *= (
process.mfvMovedTracks * process.mfvVertexSequence * process.mfvEvent * process.mfvWeight * process.mfvMovedTree
)
process.options.wantSummary = True
process.maxEvents.input = 100
file_event_from_argv(process)
if __name__ == "__main__" and hasattr(sys, "argv") and "submit" in sys.argv:
from JMTucker.Tools.CRAB3Submitter import CRABSubmitter
def customize_before_unscheduled(process):
process.load('JMTucker.MFVNeutralino.Vertexer_cff')
process.load('JMTucker.MFVNeutralino.EventProducer_cfi')
process.p = cms.Path(process.mfvVertexSequence * process.mfvEvent)
tfileservice(process, 'vertex_histos.root')
random_service(process, {'mfvVertices': 1222})
process.out.outputCommands = output_commands = [
'drop *',
'keep *_mcStat_*_*',
'keep MFVVertexAuxs_mfvVerticesAux_*_*',
'keep MFVEvent_mfvEvent__*',
]
if prepare_vis:
process.mfvSelectedVerticesTight.produce_vertices = True
process.mfvSelectedVerticesTight.produce_tracks = True
process.mfvSelectedVerticesTightLargeErr.produce_vertices = True
process.mfvSelectedVerticesTightLargeErr.produce_tracks = True
process.load('JMTucker.MFVNeutralino.VertexRefitter_cfi')
process.mfvVertexRefitsDrop0 = process.mfvVertexRefits.clone(n_tracks_to_drop = 0)
process.mfvVertexRefitsDrop2 = process.mfvVertexRefits.clone(n_tracks_to_drop = 2)
process.mfvVertexRefitsLargeErr = process.mfvVertexRefits.clone(vertex_src = 'mfvSelectedVerticesTightLargeErr')
process.mfvVertexRefitsLargeErrDrop2 = process.mfvVertexRefits.clone(vertex_src = 'mfvSelectedVerticesTightLargeErr', n_tracks_to_drop = 2)
process.p *= process.mfvVertexRefits * process.mfvVertexRefitsDrop2 * process.mfvVertexRefitsDrop0 * process.mfvVertexRefitsLargeErr * process.mfvVertexRefitsLargeErrDrop2
output_commands += [
'keep *_mfvVertices_*_*',
'keep *_mfvSelectedVerticesTight_*_*',
'keep *_mfvVertexRefits_*_*',
'keep *_mfvVertexRefitsDrop2_*_*',
'keep *_mfvVertexRefitsDrop0_*_*',
'keep *_mfvVertexRefitsLargeErr_*_*',
'keep *_mfvVertexRefitsLargeErrDrop2_*_*',
]
if is_mc:
process.mfvGenParticles = cms.EDProducer('MFVGenParticles',
gen_particles_src = cms.InputTag('genParticles'),
print_info = cms.bool(True),
)
process.p *= process.mfvGenParticles
output_commands += ['keep *_mfvGenParticles_*_*']
if keep_all:
process.mfvEvent.skip_event_filter = ''
def dedrop(l):
return [x for x in l if not x.strip().startswith('drop')]
if is_mc:
process.out.outputCommands += \
process.AODSIMEventContent.outputCommands + \
dedrop(process.MINIAODSIMEventContent.outputCommands) + \
dedrop(output_commands)
else:
process.out.outputCommands += \
process.AODEventContent.outputCommands + \
dedrop(process.MINIAODEventContent.outputCommands) + \
dedrop(output_commands)
