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)