def L1TReEmulFromRAWLegacyMuon(process): process.load('L1Trigger.Configuration.SimL1Emulator_cff') process.load('L1Trigger.Configuration.CaloTriggerPrimitives_cff') process.simEcalTriggerPrimitiveDigis.Label = 'ecalDigis' process.simHcalTriggerPrimitiveDigis.inputLabel = cms.VInputTag( cms.InputTag('hcalDigis'), cms.InputTag('hcalDigis')) ## - Legacy to upgrade format muon converter process.load('L1Trigger.L1TCommon.muonLegacyInStage2FormatDigis_cfi') process.muonLegacyInStage2FormatDigis.muonSource = cms.InputTag( 'simGmtDigis') ## - DT TP emulator from L1Trigger.DTTrigger.dtTriggerPrimitiveDigis_cfi import dtTriggerPrimitiveDigis process.simDtTriggerPrimitiveDigis = dtTriggerPrimitiveDigis.clone() process.simDtTriggerPrimitiveDigis.digiTag = cms.InputTag('muonDTDigis') ## - TwinMux from L1Trigger.L1TMuonBarrel.simTwinMuxDigis_cfi import simTwinMuxDigis process.simTwinMuxDigisForDttf = simTwinMuxDigis.clone() process.simTwinMuxDigisForDttf.RPC_Source = cms.InputTag('muonRPCDigis') process.simTwinMuxDigisForDttf.DTDigi_Source = cms.InputTag('bmtfDigis') process.simTwinMuxDigisForDttf.DTThetaDigi_Source = cms.InputTag( 'bmtfDigis') ## - CSC TP emulator from L1Trigger.CSCTriggerPrimitives.cscTriggerPrimitiveDigis_cfi import cscTriggerPrimitiveDigis process.simCscTriggerPrimitiveDigis = cscTriggerPrimitiveDigis.clone() process.simCscTriggerPrimitiveDigis.CSCComparatorDigiProducer = cms.InputTag( 'muonCSCDigis', 'MuonCSCComparatorDigi') process.simCscTriggerPrimitiveDigis.CSCWireDigiProducer = cms.InputTag( 'muonCSCDigis', 'MuonCSCWireDigi') # # - CSC Track Finder emulator # from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis process.simCsctfTrackDigis = csctfTrackDigis.clone() process.simCsctfTrackDigis.SectorReceiverInput = cms.untracked.InputTag( 'csctfDigis') process.simCsctfTrackDigis.DTproducer = 'simDtTriggerPrimitiveDigis' from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.simCsctfDigis = csctfDigis.clone() process.simCsctfDigis.CSCTrackProducer = 'simCsctfTrackDigis' ## ## - DT Track Finder emulator ## from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.simDttfDigis = dttfDigis.clone() process.simDttfDigis.DTDigi_Source = 'simTwinMuxDigisForDttf' process.simDttfDigis.CSCStub_Source = 'simCsctfTrackDigis' ## ## - RPC PAC Trigger emulator ## from L1Trigger.RPCTrigger.rpcTriggerDigis_cff import rpcTriggerDigis process.load('L1Trigger.RPCTrigger.RPCConeConfig_cff') process.simRpcTriggerDigis = rpcTriggerDigis.clone() process.simRpcTriggerDigis.label = 'muonRPCDigis' process.simRpcTriggerDigis.RPCTriggerDebug = cms.untracked.int32(1) ## ## - Legacy Global Muon Trigger emulator ## from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.simGmtDigis = gmtDigis.clone() process.simGmtDigis.DTCandidates = cms.InputTag('simDttfDigis', 'DT') process.simGmtDigis.CSCCandidates = cms.InputTag('simCsctfDigis', 'CSC') process.simGmtDigis.RPCbCandidates = cms.InputTag('simRpcTriggerDigis', 'RPCb') process.simGmtDigis.RPCfCandidates = cms.InputTag('simRpcTriggerDigis', 'RPCf') # This is for the upgrade # BMTF process.simBmtfDigis.DTDigi_Source = cms.InputTag('bmtfDigis') process.simBmtfDigis.DTDigi_Theta_Source = cms.InputTag('bmtfDigis') # TwinMux process.simTwinMuxDigis.RPC_Source = cms.InputTag('muonRPCDigis') process.simTwinMuxDigis.DTDigi_Source = cms.InputTag('bmtfDigis') process.simTwinMuxDigis.DTThetaDigi_Source = cms.InputTag('bmtfDigis') # OMTF process.simOmtfDigis.srcRPC = cms.InputTag('muonRPCDigis') process.simOmtfDigis.srcCSC = cms.InputTag('csctfDigis') process.simOmtfDigis.srcDTPh = cms.InputTag('bmtfDigis') process.simOmtfDigis.srcDTTh = cms.InputTag('bmtfDigis') # EMTF process.simEmtfDigis.CSCInput = cms.InputTag('emtfStage2Digis') process.simEmtfDigis.RPCInput = cms.InputTag('muonRPCDigis') # Calo Layer1 process.simCaloStage2Layer1Digis.ecalToken = cms.InputTag( 'ecalDigis:EcalTriggerPrimitives') process.simCaloStage2Layer1Digis.hcalToken = cms.InputTag('hcalDigis:') # - Sequences process.L1TMuonTriggerPrimitives = cms.Sequence( process.simCscTriggerPrimitiveDigis + process.simDtTriggerPrimitiveDigis + process.simTwinMuxDigisForDttf) process.L1TReEmul = cms.Sequence(process.L1TMuonTriggerPrimitives + process.simCsctfTrackDigis + process.simCsctfDigis + process.simDttfDigis + process.simRpcTriggerDigis + process.simGmtDigis + process.muonLegacyInStage2FormatDigis) process.load('L1Trigger.L1TMuon.simMuonQualityAdjusterDigis_cfi') process.L1TReEmul = cms.Sequence(process.L1TReEmul + process.simTwinMuxDigis + process.simBmtfDigis + process.simEmtfDigis + process.simOmtfDigis + process.simGmtCaloSumDigis + process.simMuonQualityAdjusterDigis + process.simGmtStage2Digis) process.L1TReEmul = cms.Sequence(process.L1TReEmul + process.SimL1TechnicalTriggers + process.SimL1TGlobal) process.L1TReEmulPath = cms.Path(process.L1TReEmul) process.schedule.append(process.L1TReEmulPath) print "# L1TReEmul sequence: " print "# {0}".format(process.L1TReEmul) print "# {0}".format(process.schedule) return process
def reEmulation(process, reEmulMuons=True, reEmulCalos=True, patchNtuple=True, runOnPostLS1 = False): print "[L1Menu]: Setting up overall re-emulation" if patchNtuple and hasattr(process,'l1NtupleProducer') and hasattr(process,'l1ExtraTreeProducer') : print "[L1Menu]:\tConfiguring Ntuple to use re-emulated information" ntuple = getattr(process,'l1NtupleProducer') l1ExtraNtuple = getattr(process,'l1ExtraTreeProducer') elif patchNtuple : print "[L1Menu]:\tERROR: FAILED to find ntuple! switching patchNtuple to False" patchNtuple=False process.l1ExtraReEmul = cms.EDProducer( "L1ExtraParticlesProd", muonSource = cms.InputTag("gtDigis"), isolatedEmSource = cms.InputTag("gctDigis","isoEm"), nonIsolatedEmSource = cms.InputTag("gctDigis","nonIsoEm"), forwardJetSource = cms.InputTag("gctDigis","forJets"), centralJetSource = cms.InputTag("gctDigis","cenJets"), tauJetSource = cms.InputTag("gctDigis","tauJets"), etTotalSource = cms.InputTag("gctDigis"), etHadSource = cms.InputTag("gctDigis"), etMissSource = cms.InputTag("gctDigis"), htMissSource = cms.InputTag("gctDigis"), hfRingEtSumsSource = cms.InputTag("gctDigis"), hfRingBitCountsSource = cms.InputTag("gctDigis"), produceMuonParticles = cms.bool(True), produceCaloParticles = cms.bool(True), centralBxOnly = cms.bool(True), ignoreHtMiss = cms.bool(False) ) if reEmulMuons : print "[L1Menu]:\tSetting up muon re-emulation" from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag("csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() if not runOnPostLS1 : from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.csctfReEmulTrackDigis = csctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.readDtDirect = True process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag("csctfDigis") process.csctfReEmulTrackDigis.DtDirectProd = cms.untracked.InputTag("csctfDigis","DT") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") #process.csctfReEmulDigis.SectorProcessor.initializeFromPSet = True process.csctfReEmulSequence = cms.Sequence( process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) else : from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag("dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") process.csctfReEmulTrackDigis.SectorProcessor.PTLUT.PtMethod = cms.untracked.uint32(33) # no triple ganging in ME11a process.csctfReEmulTrackDigis.SectorProcessor.gangedME1a = cms.untracked.bool(False) process.csctfReEmulTrackDigis.SectorProcessor.firmwareSP = cms.uint32(20120730) #core 20120730 process.csctfReEmulTrackDigis.SectorProcessor.initializeFromPSet = cms.bool(True) process.csctfReEmulSequence = cms.Sequence( process.simCscTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) process.load('L1TriggerConfig.GMTConfigProducers.L1MuGMTParameters_cfi') process.L1MuGMTParameters.MergeMethodPtBrl=cms.string("byCombi") process.L1MuGMTParameters.MergeMethodPtFwd=cms.string("byCombi") process.L1MuGMTParameters.VersionSortRankEtaQLUT = cms.uint32(1043) process.L1MuGMTParameters.VersionLUTs = cms.uint32(1) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag("dttfReEmulDigis","DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag("csctfReEmulDigis","CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.l1ExtraReEmul.muonSource = cms.InputTag("gmtReEmulDigis") if patchNtuple : ntuple.gmtSource = cms.InputTag("gmtReEmulDigis") ntuple.dttfSource = cms.InputTag("dttfReEmulDigis") ntuple.csctfTrkSource = cms.InputTag("csctfReEmulTrackDigis") ntuple.csctfStatusSource = cms.InputTag("csctfReEmulTrackDigis") l1ExtraNtuple.muonLabel = cms.untracked.InputTag("l1ExtraReEmul") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis *process.csctfReEmulSequence *process.dttfReEmulDigis *process.gmtReEmulDigis ) if reEmulCalos : print "[L1Menu]:\tSetting up calo re-emulation" print "[L1Menu]:\tWARNING! Just patching GCT sequence, RCT can't be re-emulated after UCT2015 SW patches" # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly # # In MC HCAL need to be re-run as there is no TPG information stored # process.load("SimCalorimetry.HcalSimProducers.hcalUnsuppressedDigis_cfi") # process.load("SimCalorimetry.HcalTrigPrimProducers.hcaltpdigi_cff") # # from L1Trigger.RegionalCaloTrigger.rctDigis_cfi import rctDigis # from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis # # process.hcalReEmulDigis = process.simHcalTriggerPrimitiveDigis.clone() # process.rctReEmulDigis = rctDigis.clone() # process.gctReEmulDigis = gctDigis.clone() # # process.hcalReEmulDigis.inputLabel = cms.VInputTag(cms.InputTag('hcalDigis'), cms.InputTag('hcalDigis')) # #process.HcalTPGCoderULUT.LUTGenerationMode = cms.bool(False) # # process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalDigis:EcalTriggerPrimitives' ) ) # process.rctReEmulDigis.hcalDigis = cms.VInputTag( cms.InputTag( 'hcalReEmulDigis' ) ) # # process.gctReEmulDigis.inputLabel = cms.InputTag("rctReEmulDigis") from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis process.gctReEmulDigis = gctDigis.clone() process.gctReEmulDigis.inputLabel = cms.InputTag("gctDigis") process.l1ExtraReEmul.isolatedEmSource = cms.InputTag("gctReEmulDigis","isoEm") process.l1ExtraReEmul.nonIsolatedEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") process.l1ExtraReEmul.forwardJetSource = cms.InputTag("gctReEmulDigis","forJets") process.l1ExtraReEmul.centralJetSource = cms.InputTag("gctReEmulDigis","cenJets") process.l1ExtraReEmul.tauJetSource = cms.InputTag("gctReEmulDigis","tauJets") process.l1ExtraReEmul.etTotalSource = cms.InputTag("gctDigis") process.l1ExtraReEmul.etHadSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.htMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingEtSumsSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingBitCountsSource = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gctCentralJetsSource = cms.InputTag("gctReEmulDigis","cenJets") ntuple.gctNonIsoEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") ntuple.gctForwardJetsSource = cms.InputTag("gctReEmulDigis","forJets") ntuple.gctIsoEmSource = cms.InputTag("gctReEmulDigis","isoEm") ntuple.gctEnergySumsSource = cms.InputTag("gctReEmulDigis","") ntuple.gctTauJetsSource = cms.InputTag("gctReEmulDigis","tauJets") l1ExtraNtuple.nonIsoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:NonIsolated") l1ExtraNtuple.isoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:Isolated") l1ExtraNtuple.tauJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Tau") l1ExtraNtuple.cenJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Central") l1ExtraNtuple.fwdJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Forward") l1ExtraNtuple.metLabel = cms.untracked.InputTag("l1ExtraReEmul:MET") l1ExtraNtuple.mhtLabel = cms.untracked.InputTag("l1ExtraReEmul:MHT") # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly # ntuple.rctSource = cms.InputTag("rctReEmulDigis") process.reEmulCaloChain = cms.Sequence( # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly #process.hcalReEmulDigis #+ process.rctReEmulDigis process.gctReEmulDigis ) from L1Trigger.GlobalTrigger.gtDigis_cfi import gtDigis process.gtReEmulDigis = gtDigis.clone() if reEmulMuons : process.gtReEmulDigis.GmtInputTag = cms.InputTag("gmtReEmulDigis") if reEmulCalos : process.gtReEmulDigis.GctInputTag = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gtSource = cms.InputTag("gtReEmulDigis") if reEmulMuons and reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulMuons : process.reEmul = cms.Sequence(process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.gtReEmulDigis + process.l1ExtraReEmul) else : process.reEmul = cms.Sequence(process.gtReEmulDigis + process.l1ExtraReEmul)
def L1TReEmulFromRAWLegacyMuon(process): process.load('L1Trigger.Configuration.SimL1Emulator_cff') process.load('L1Trigger.Configuration.CaloTriggerPrimitives_cff') process.simEcalTriggerPrimitiveDigis.Label = 'ecalDigis' process.simHcalTriggerPrimitiveDigis.inputLabel = cms.VInputTag( cms.InputTag('hcalDigis'), cms.InputTag('hcalDigis') ) ## - Legacy to upgrade format muon converter process.load('L1Trigger.L1TCommon.muonLegacyInStage2FormatDigis_cfi') process.muonLegacyInStage2FormatDigis.muonSource = cms.InputTag('simGmtDigis') ## - DT TP emulator from L1Trigger.DTTrigger.dtTriggerPrimitiveDigis_cfi import dtTriggerPrimitiveDigis process.simDtTriggerPrimitiveDigis = dtTriggerPrimitiveDigis.clone() process.simDtTriggerPrimitiveDigis.digiTag = cms.InputTag('muonDTDigis') ## - TwinMux from L1Trigger.L1TTwinMux.simTwinMuxDigis_cfi import simTwinMuxDigis process.simTwinMuxDigisForDttf = simTwinMuxDigis.clone() process.simTwinMuxDigisForDttf.RPC_Source = cms.InputTag('muonRPCDigis') process.simTwinMuxDigisForDttf.DTDigi_Source = cms.InputTag('bmtfDigis') process.simTwinMuxDigisForDttf.DTThetaDigi_Source = cms.InputTag('bmtfDigis') ## - CSC TP emulator from L1Trigger.CSCTriggerPrimitives.cscTriggerPrimitiveDigis_cfi import cscTriggerPrimitiveDigis process.simCscTriggerPrimitiveDigis = cscTriggerPrimitiveDigis.clone() process.simCscTriggerPrimitiveDigis.CSCComparatorDigiProducer = cms.InputTag( 'muonCSCDigis', 'MuonCSCComparatorDigi' ) process.simCscTriggerPrimitiveDigis.CSCWireDigiProducer = cms.InputTag( 'muonCSCDigis', 'MuonCSCWireDigi' ) # # - CSC Track Finder emulator # from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis process.simCsctfTrackDigis = csctfTrackDigis.clone() process.simCsctfTrackDigis.SectorReceiverInput = cms.untracked.InputTag( 'csctfDigis' ) process.simCsctfTrackDigis.DTproducer = 'simDtTriggerPrimitiveDigis' from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.simCsctfDigis = csctfDigis.clone() process.simCsctfDigis.CSCTrackProducer = 'simCsctfTrackDigis' ## ## - DT Track Finder emulator ## from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.simDttfDigis = dttfDigis.clone() process.simDttfDigis.DTDigi_Source = 'simTwinMuxDigisForDttf' process.simDttfDigis.CSCStub_Source = 'simCsctfTrackDigis' ## ## - RPC PAC Trigger emulator ## from L1Trigger.RPCTrigger.rpcTriggerDigis_cff import rpcTriggerDigis process.load('L1Trigger.RPCTrigger.RPCConeConfig_cff') process.simRpcTriggerDigis = rpcTriggerDigis.clone() process.simRpcTriggerDigis.label = 'muonRPCDigis' process.simRpcTriggerDigis.RPCTriggerDebug = cms.untracked.int32(1) ## ## - Legacy Global Muon Trigger emulator ## from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.simGmtDigis = gmtDigis.clone() process.simGmtDigis.DTCandidates = cms.InputTag( 'simDttfDigis', 'DT' ) process.simGmtDigis.CSCCandidates = cms.InputTag( 'simCsctfDigis', 'CSC' ) process.simGmtDigis.RPCbCandidates = cms.InputTag( 'simRpcTriggerDigis', 'RPCb' ) process.simGmtDigis.RPCfCandidates = cms.InputTag( 'simRpcTriggerDigis', 'RPCf' ) # This is for the upgrade # BMTF process.simBmtfDigis.DTDigi_Source = cms.InputTag('bmtfDigis') process.simBmtfDigis.DTDigi_Theta_Source = cms.InputTag('bmtfDigis') # TwinMux process.simTwinMuxDigis.RPC_Source = cms.InputTag('muonRPCDigis') process.simTwinMuxDigis.DTDigi_Source = cms.InputTag('bmtfDigis') process.simTwinMuxDigis.DTThetaDigi_Source = cms.InputTag('bmtfDigis') # OMTF process.simOmtfDigis.srcRPC = cms.InputTag('muonRPCDigis') process.simOmtfDigis.srcCSC = cms.InputTag('csctfDigis') process.simOmtfDigis.srcDTPh = cms.InputTag('bmtfDigis') process.simOmtfDigis.srcDTTh = cms.InputTag('bmtfDigis') # EMTF process.simEmtfDigis.CSCInput = cms.InputTag('emtfStage2Digis') process.simEmtfDigis.RPCInput = cms.InputTag('muonRPCDigis') # Calo Layer1 process.simCaloStage2Layer1Digis.ecalToken = cms.InputTag('ecalDigis:EcalTriggerPrimitives') process.simCaloStage2Layer1Digis.hcalToken = cms.InputTag('hcalDigis:') # - Sequences process.L1TMuonTriggerPrimitives = cms.Sequence(process.simCscTriggerPrimitiveDigis + process.simDtTriggerPrimitiveDigis + process.simTwinMuxDigisForDttf) process.L1TReEmul = cms.Sequence(process.L1TMuonTriggerPrimitives + process.simCsctfTrackDigis + process.simCsctfDigis + process.simDttfDigis + process.simRpcTriggerDigis + process.simGmtDigis + process.muonLegacyInStage2FormatDigis) process.load('L1Trigger.L1TMuon.simMuonQualityAdjusterDigis_cfi') process.L1TReEmul = cms.Sequence( process.L1TReEmul + process.simTwinMuxDigis + process.simBmtfDigis + process.simEmtfDigis + process.simOmtfDigis + process.simGmtCaloSumDigis + process.simMuonQualityAdjusterDigis + process.simGmtStage2Digis) process.L1TReEmul = cms.Sequence( process.L1TReEmul + process.SimL1TechnicalTriggers + process.SimL1TGlobal ) process.L1TReEmulPath = cms.Path(process.L1TReEmul) process.schedule.append(process.L1TReEmulPath) print("# L1TReEmul sequence: ") print("# {0}".format(process.L1TReEmul)) print("# {0}".format(process.schedule)) return process
def customiseL1EmulatorFromRaw(process): # customization fragment to run full emulator chain (TPGs and L1 Muon,Calo and GT emulators) # starting from a RAW file assuming that "RawToDigi_cff" and "SimL1Emulator_cff" # have already been loaded # (GMT digis produced by same module as the GT digis, as GT and GMT have common unpacker) ## process.simRpcTechTrigDigis.RPCDigiLabel = 'muonRPCDigis' ## 2015 L1 Calo Emulator process.load('L1Trigger.L1TCalorimeter.L1TCaloStage1_PPFromRaw_cff') ### 2015 L1 Muon Emulator from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag( "csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() process.load("L1TriggerConfig.RPCTriggerConfig.RPCConeDefinition_cff") process.load("L1Trigger.RPCTrigger.RPCConeConfig_cff") from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator_sim from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator_sim(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag( "dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag( "csctfReEmulTrackDigis") process.csctfReEmulTrackDigis.SectorProcessor.PTLUT.PtMethod = cms.untracked.uint32( 34) # no triple ganging in ME11a process.csctfReEmulTrackDigis.SectorProcessor.gangedME1a = cms.untracked.bool( False) process.csctfReEmulTrackDigis.SectorProcessor.firmwareSP = cms.uint32( 20140515) #core 20120730 process.csctfReEmulTrackDigis.SectorProcessor.initializeFromPSet = cms.bool( True) process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag( "cscReEmulTriggerPrimitiveDigis", "MPCSORTED") process.cscReEmulTriggerPrimitiveDigis = process.simCscTriggerPrimitiveDigis.clone( ) process.cscReEmulTriggerPrimitiveDigis.CSCComparatorDigiProducer = cms.InputTag( "muonCSCDigis", "MuonCSCComparatorDigi") process.cscReEmulTriggerPrimitiveDigis.CSCWireDigiProducer = cms.InputTag( "muonCSCDigis", "MuonCSCWireDigi") process.csctfReEmulSequence = cms.Sequence( process.cscReEmulTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag("dttfReEmulDigis", "DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag( "csctfReEmulDigis", "CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag( "rpcTriggerReEmulDigis", "RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag( "rpcTriggerReEmulDigis", "RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.load('L1TriggerConfig.GMTConfigProducers.L1MuGMTParameters_cfi') process.L1MuGMTParameters.MergeMethodPtBrl = cms.string("byMinPt") process.L1MuGMTParameters.MergeMethodPtFwd = cms.string("byMinPt") process.L1MuGMTParameters.VersionSortRankEtaQLUT = cms.uint32(275) process.L1MuGMTParameters.VersionLUTs = cms.uint32(1) process.es_prefer_gmtConfig = cms.ESPrefer("L1MuGMTParametersProducer", "L1MuGMTParameters") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis * process.csctfReEmulSequence * process.dttfReEmulDigis * process.gmtReEmulDigis) ## remove for 720pre7 and beyond ## use new muon LUTs ## dttfFile='sqlite:dttf_config.db' ## process.GlobalTag.toGet.extend( ## cms.VPSet(cms.PSet(record = cms.string("L1MuDTEtaPatternLutRcd"), ## tag = cms.string("L1MuDTEtaPatternLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTExtLutRcd"), ## tag = cms.string("L1MuDTExtLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTPhiLutRcd"), ## tag = cms.string("L1MuDTPhiLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTPtaLutRcd"), ## tag = cms.string("L1MuDTPtaLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTQualPatternLutRcd"), ## tag = cms.string("L1MuDTQualPatternLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ) ## ) ## ) ## end # GT from L1Trigger.Configuration.SimL1Emulator_cff import simGtDigis simGtDigis.GmtInputTag = 'gmtReEmulDigis' simGtDigis.GctInputTag = 'simCaloStage1LegacyFormatDigis' simGtDigis.TechnicalTriggersInputTags = cms.VInputTag() # run Calo TPGs, L1 GCT, technical triggers, L1 GT SimL1Emulator = cms.Sequence( ## process.reEmulCaloChain + process.L1TCaloStage1_PPFromRaw + process.reEmulMuonChain + process.simGtDigis) # replace the SimL1Emulator in all paths and sequences for iterable in six.itervalues(process.sequences): iterable.replace(process.SimL1Emulator, SimL1Emulator) for iterable in six.itervalues(process.paths): iterable.replace(process.SimL1Emulator, SimL1Emulator) for iterable in six.itervalues(process.endpaths): iterable.replace(process.SimL1Emulator, SimL1Emulator) process.SimL1Emulator = SimL1Emulator return process
def reEmulation(process, reEmulMuons=True, reEmulCalos=True, patchNtuple=True, runOnPostLS1 = False, useStage1Layer2=False): print "[L1Menu]: Setting up overall re-emulation" if patchNtuple and hasattr(process,'l1NtupleProducer') and hasattr(process,'l1ExtraTreeProducer') : print "[L1Menu]:\tConfiguring Ntuple to use re-emulated information" ntuple = getattr(process,'l1NtupleProducer') l1ExtraNtuple = getattr(process,'l1ExtraTreeProducer') elif patchNtuple : print "[L1Menu]:\tERROR: FAILED to find ntuple! switching patchNtuple to False" patchNtuple=False process.l1ExtraReEmul = cms.EDProducer( "L1ExtraParticlesProd", muonSource = cms.InputTag("gtDigis"), isolatedEmSource = cms.InputTag("gctDigis","isoEm"), nonIsolatedEmSource = cms.InputTag("gctDigis","nonIsoEm"), forwardJetSource = cms.InputTag("gctDigis","forJets"), centralJetSource = cms.InputTag("gctDigis","cenJets"), tauJetSource = cms.InputTag("gctDigis","tauJets"), isoTauJetSource = cms.InputTag("gctDigis","isoTauJets"), etTotalSource = cms.InputTag("gctDigis"), etHadSource = cms.InputTag("gctDigis"), etMissSource = cms.InputTag("gctDigis"), htMissSource = cms.InputTag("gctDigis"), hfRingEtSumsSource = cms.InputTag("gctDigis"), hfRingBitCountsSource = cms.InputTag("gctDigis"), produceMuonParticles = cms.bool(True), produceCaloParticles = cms.bool(True), centralBxOnly = cms.bool(True), ignoreHtMiss = cms.bool(False) ) if reEmulMuons : print "[L1Menu]:\tSetting up muon re-emulation" from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag("csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() if not runOnPostLS1 : from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.csctfReEmulTrackDigis = csctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.readDtDirect = True process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag("csctfDigis") process.csctfReEmulTrackDigis.DtDirectProd = cms.untracked.InputTag("csctfDigis","DT") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") #process.csctfReEmulDigis.SectorProcessor.initializeFromPSet = True process.csctfReEmulSequence = cms.Sequence( process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) else : from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator_sim from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator_sim(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag("dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") process.csctfReEmulTrackDigis.SectorProcessor.PTLUT.PtMethod = cms.untracked.uint32(34) # no triple ganging in ME11a process.csctfReEmulTrackDigis.SectorProcessor.gangedME1a = cms.untracked.bool(False) process.csctfReEmulTrackDigis.SectorProcessor.firmwareSP = cms.uint32(20140515) #core 20120730 process.csctfReEmulTrackDigis.SectorProcessor.initializeFromPSet = cms.bool(True) process.csctfReEmulSequence = cms.Sequence( process.simCscTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) process.load('L1TriggerConfig.GMTConfigProducers.L1MuGMTParameters_cfi') process.L1MuGMTParameters.MergeMethodPtBrl=cms.string("byCombi") process.L1MuGMTParameters.MergeMethodPtFwd=cms.string("byCombi") process.L1MuGMTParameters.VersionSortRankEtaQLUT = cms.uint32(1043) process.L1MuGMTParameters.VersionLUTs = cms.uint32(1) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag("dttfReEmulDigis","DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag("csctfReEmulDigis","CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.l1ExtraReEmul.muonSource = cms.InputTag("gmtReEmulDigis") if patchNtuple : ntuple.gmtSource = cms.InputTag("gmtReEmulDigis") ntuple.dttfSource = cms.InputTag("dttfReEmulDigis") ntuple.csctfTrkSource = cms.InputTag("csctfReEmulTrackDigis") ntuple.csctfStatusSource = cms.InputTag("csctfReEmulTrackDigis") l1ExtraNtuple.muonLabel = cms.untracked.InputTag("l1ExtraReEmul") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis *process.csctfReEmulSequence *process.dttfReEmulDigis *process.gmtReEmulDigis ) if reEmulCalos : print "[L1Menu]:\tSetting up calo re-emulation" ## redo ECAL TP's ## #from SimCalorimetry.EcalTrigPrimProducers.ecalTriggerPrimitiveDigis_cff import simEcalTriggerPrimitiveDigis #process.ecalReEmulDigis = simEcalTriggerPrimitiveDigis.clone() #process.ecalReEmulDigis.InstanceEB = cms.string('ebDigis') #process.ecalReEmulDigis.InstanceEE = cms.string('eeDigis') #process.ecalReEmulDigis.Label = cms.string('ecalDigis') ### from L1Trigger.Configuration.SimL1Emulator_cff import simRctDigis process.rctReEmulDigis = process.simRctDigis.clone() process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalDigis:EcalTriggerPrimitives' ) ) #process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalReEmulDigis' ) ) process.rctReEmulDigis.hcalDigis = cms.VInputTag( cms.InputTag( 'hcalDigis' ) ) from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis process.gctReEmulDigis = gctDigis.clone() process.gctReEmulDigis.inputLabel = cms.InputTag("gctDigis") process.l1ExtraReEmul.isolatedEmSource = cms.InputTag("gctReEmulDigis","isoEm") process.l1ExtraReEmul.nonIsolatedEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") process.l1ExtraReEmul.forwardJetSource = cms.InputTag("gctReEmulDigis","forJets") process.l1ExtraReEmul.centralJetSource = cms.InputTag("gctReEmulDigis","cenJets") process.l1ExtraReEmul.tauJetSource = cms.InputTag("gctReEmulDigis","tauJets") process.l1ExtraReEmul.isoTauJetSource = cms.InputTag("gctReEmulDigis","isoTauJets") process.l1ExtraReEmul.etTotalSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etHadSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.htMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingEtSumsSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingBitCountsSource = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gctCentralJetsSource = cms.InputTag("gctReEmulDigis","cenJets") ntuple.gctNonIsoEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") ntuple.gctForwardJetsSource = cms.InputTag("gctReEmulDigis","forJets") ntuple.gctIsoEmSource = cms.InputTag("gctReEmulDigis","isoEm") ntuple.gctEnergySumsSource = cms.InputTag("gctReEmulDigis","") ntuple.gctTauJetsSource = cms.InputTag("gctReEmulDigis","tauJets") if useStage1Layer2: ntuple.gctIsoTauJetsSource = cms.InputTag("gctReEmulDigis","isoTauJets") else: ntuple.gctIsoTauJetsSource = cms.InputTag("none","isoTauJets") l1ExtraNtuple.nonIsoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:NonIsolated") l1ExtraNtuple.isoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:Isolated") l1ExtraNtuple.tauJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Tau") l1ExtraNtuple.isoTauJetLabel = cms.untracked.InputTag("l1ExtraReEmul:IsoTau") l1ExtraNtuple.cenJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Central") l1ExtraNtuple.fwdJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Forward") l1ExtraNtuple.metLabel = cms.untracked.InputTag("l1ExtraReEmul:MET") l1ExtraNtuple.mhtLabel = cms.untracked.InputTag("l1ExtraReEmul:MHT") # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly # ntuple.rctSource = cms.InputTag("rctReEmulDigis") process.reEmulCaloChain = cms.Sequence( # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly #process.hcalReEmulDigis #process.ecalReEmulDigis process.rctReEmulDigis +process.gctReEmulDigis ) from L1Trigger.GlobalTrigger.gtDigis_cfi import gtDigis process.gtReEmulDigis = gtDigis.clone() if reEmulMuons : process.gtReEmulDigis.GmtInputTag = cms.InputTag("gmtReEmulDigis") if reEmulCalos : process.gctReEmulDigis.inputLabel = cms.InputTag("rctReEmulDigis") process.gtReEmulDigis.GctInputTag = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gtSource = cms.InputTag("gtReEmulDigis") if reEmulMuons and reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulMuons : process.reEmul = cms.Sequence(process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.gtReEmulDigis + process.l1ExtraReEmul) else : process.reEmul = cms.Sequence(process.gtReEmulDigis + process.l1ExtraReEmul)
def customiseL1EmulatorFromRaw(process): # customization fragment to run full emulator chain (TPGs and L1 Muon,Calo and GT emulators) # starting from a RAW file assuming that "RawToDigi_cff" and "SimL1Emulator_cff" # have already been loaded # (GMT digis produced by same module as the GT digis, as GT and GMT have common unpacker) ## process.simRpcTechTrigDigis.RPCDigiLabel = 'muonRPCDigis' ## 2015 L1 Calo Emulator process.load('L1Trigger.L1TCalorimeter.L1TCaloStage1_PPFromRaw_cff') ### 2015 L1 Muon Emulator from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag("csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() process.load("L1TriggerConfig.RPCTriggerConfig.RPCConeDefinition_cff") process.load("L1Trigger.RPCTrigger.RPCConeConfig_cff") from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator_sim from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator_sim(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag("dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") process.csctfReEmulTrackDigis.SectorProcessor.PTLUT.PtMethod = cms.untracked.uint32(34) # no triple ganging in ME11a process.csctfReEmulTrackDigis.SectorProcessor.gangedME1a = cms.untracked.bool(False) process.csctfReEmulTrackDigis.SectorProcessor.firmwareSP = cms.uint32(20140515) #core 20120730 process.csctfReEmulTrackDigis.SectorProcessor.initializeFromPSet = cms.bool(True) process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag("cscReEmulTriggerPrimitiveDigis","MPCSORTED") process.cscReEmulTriggerPrimitiveDigis = process.simCscTriggerPrimitiveDigis.clone() process.cscReEmulTriggerPrimitiveDigis.CSCComparatorDigiProducer = cms.InputTag("muonCSCDigis","MuonCSCComparatorDigi") process.cscReEmulTriggerPrimitiveDigis.CSCWireDigiProducer = cms.InputTag("muonCSCDigis","MuonCSCWireDigi") process.csctfReEmulSequence = cms.Sequence( process.cscReEmulTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag("dttfReEmulDigis","DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag("csctfReEmulDigis","CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.load('L1TriggerConfig.GMTConfigProducers.L1MuGMTParameters_cfi') process.L1MuGMTParameters.MergeMethodPtBrl=cms.string("byMinPt") process.L1MuGMTParameters.MergeMethodPtFwd=cms.string("byMinPt") process.L1MuGMTParameters.VersionSortRankEtaQLUT = cms.uint32(275) process.L1MuGMTParameters.VersionLUTs = cms.uint32(1) process.es_prefer_gmtConfig = cms.ESPrefer("L1MuGMTParametersProducer","L1MuGMTParameters") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis *process.csctfReEmulSequence *process.dttfReEmulDigis *process.gmtReEmulDigis ) ## remove for 720pre7 and beyond ## use new muon LUTs ## dttfFile='sqlite:dttf_config.db' ## process.GlobalTag.toGet.extend( ## cms.VPSet(cms.PSet(record = cms.string("L1MuDTEtaPatternLutRcd"), ## tag = cms.string("L1MuDTEtaPatternLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTExtLutRcd"), ## tag = cms.string("L1MuDTExtLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTPhiLutRcd"), ## tag = cms.string("L1MuDTPhiLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTPtaLutRcd"), ## tag = cms.string("L1MuDTPtaLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ), ## cms.PSet(record = cms.string("L1MuDTQualPatternLutRcd"), ## tag = cms.string("L1MuDTQualPatternLut_CRAFT09_hlt"), ## connect = cms.untracked.string(dttfFile) ## ) ## ) ## ) ## end # GT from L1Trigger.Configuration.SimL1Emulator_cff import simGtDigis simGtDigis.GmtInputTag = 'gmtReEmulDigis' simGtDigis.GctInputTag = 'simCaloStage1LegacyFormatDigis' simGtDigis.TechnicalTriggersInputTags = cms.VInputTag( ) # run Calo TPGs, L1 GCT, technical triggers, L1 GT SimL1Emulator = cms.Sequence( ## process.reEmulCaloChain + process.L1TCaloStage1_PPFromRaw + process.reEmulMuonChain + process.simGtDigis ) # replace the SimL1Emulator in all paths and sequences for iterable in process.sequences.itervalues(): iterable.replace( process.SimL1Emulator, SimL1Emulator) for iterable in process.paths.itervalues(): iterable.replace( process.SimL1Emulator, SimL1Emulator) for iterable in process.endpaths.itervalues(): iterable.replace( process.SimL1Emulator, SimL1Emulator) process.SimL1Emulator = SimL1Emulator return process
def reEmulation(process, reEmulMuons=True, reEmulCalos=True, patchNtuple=True, runOnPostLS1 = True, useStage1Layer2=True, reEmulRCT = True): print "[L1Menu]: Setting up overall re-emulation" if patchNtuple and hasattr(process,'l1NtupleProducer') and hasattr(process,'l1ExtraTreeProducer') : print "[L1Menu]:\tConfiguring Ntuple to use re-emulated information" ntuple = getattr(process,'l1NtupleProducer') l1ExtraNtuple = getattr(process,'l1ExtraTreeProducer') elif patchNtuple : print "[L1Menu]:\tERROR: FAILED to find ntuple! switching patchNtuple to False" patchNtuple=False process.l1ExtraReEmul = cms.EDProducer( "L1ExtraParticlesProd", muonSource = cms.InputTag("gtDigis"), isolatedEmSource = cms.InputTag("gctDigis","isoEm"), nonIsolatedEmSource = cms.InputTag("gctDigis","nonIsoEm"), forwardJetSource = cms.InputTag("gctDigis","forJets"), centralJetSource = cms.InputTag("gctDigis","cenJets"), tauJetSource = cms.InputTag("gctDigis","tauJets"), isoTauJetSource = cms.InputTag("gctDigis","isoTauJets"), etTotalSource = cms.InputTag("gctDigis"), etHadSource = cms.InputTag("gctDigis"), etMissSource = cms.InputTag("gctDigis"), htMissSource = cms.InputTag("gctDigis"), hfRingEtSumsSource = cms.InputTag("gctDigis"), hfRingBitCountsSource = cms.InputTag("gctDigis"), produceMuonParticles = cms.bool(True), produceCaloParticles = cms.bool(True), centralBxOnly = cms.bool(True), ignoreHtMiss = cms.bool(False) ) if reEmulMuons : print "[L1Menu]:\tSetting up muon re-emulation" from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag("csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() if not runOnPostLS1 : from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.csctfReEmulTrackDigis = csctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.readDtDirect = True process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag("csctfDigis") process.csctfReEmulTrackDigis.DtDirectProd = cms.untracked.InputTag("csctfDigis","DT") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") process.csctfReEmulSequence = cms.Sequence( process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) else : from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator_sim from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator_sim(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag("dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag("csctfReEmulTrackDigis") process.csctfReEmulTrackDigis.SectorProcessor.PTLUT.PtMethod = cms.untracked.uint32(34) # no triple ganging in ME11a process.csctfReEmulTrackDigis.SectorProcessor.gangedME1a = cms.untracked.bool(False) process.csctfReEmulTrackDigis.SectorProcessor.firmwareSP = cms.uint32(20140515) #core 20120730 process.csctfReEmulTrackDigis.SectorProcessor.initializeFromPSet = cms.bool(True) process.csctfReEmulSequence = cms.Sequence( process.simCscTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis ) process.load('L1TriggerConfig.GMTConfigProducers.L1MuGMTParameters_cfi') process.L1MuGMTParameters.MergeMethodPtBrl=cms.string("byCombi") process.L1MuGMTParameters.MergeMethodPtFwd=cms.string("byCombi") process.L1MuGMTParameters.VersionSortRankEtaQLUT = cms.uint32(1043) process.L1MuGMTParameters.VersionLUTs = cms.uint32(1) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag("dttfReEmulDigis","DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag("csctfReEmulDigis","CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag("rpcTriggerReEmulDigis","RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.l1ExtraReEmul.muonSource = cms.InputTag("gmtReEmulDigis") if patchNtuple : ntuple.gmtSource = cms.InputTag("gmtReEmulDigis") ntuple.dttfSource = cms.InputTag("dttfReEmulDigis") ntuple.csctfTrkSource = cms.InputTag("csctfReEmulTrackDigis") ntuple.csctfStatusSource = cms.InputTag("csctfReEmulTrackDigis") l1ExtraNtuple.muonLabel = cms.untracked.InputTag("l1ExtraReEmul") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis *process.csctfReEmulSequence *process.dttfReEmulDigis *process.gmtReEmulDigis ) if reEmulCalos : print "[L1Menu]:\tSetting up calo re-emulation" ## redo ECAL TP's ## #from SimCalorimetry.EcalTrigPrimProducers.ecalTriggerPrimitiveDigis_cff import simEcalTriggerPrimitiveDigis #process.ecalReEmulDigis = simEcalTriggerPrimitiveDigis.clone() #process.ecalReEmulDigis.InstanceEB = cms.string('ebDigis') #process.ecalReEmulDigis.InstanceEE = cms.string('eeDigis') #process.ecalReEmulDigis.Label = cms.string('ecalDigis') if reEmulRCT : from L1Trigger.Configuration.SimL1Emulator_cff import simRctDigis process.rctReEmulDigis = process.simRctDigis.clone() process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalDigis:EcalTriggerPrimitives' ) ) ##process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalReEmulDigis' ) ) process.rctReEmulDigis.hcalDigis = cms.VInputTag( cms.InputTag( 'hcalDigis' ) ) from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis process.gctReEmulDigis = gctDigis.clone() process.gctReEmulDigis.inputLabel = cms.InputTag("gctDigis") process.l1ExtraReEmul.isolatedEmSource = cms.InputTag("gctReEmulDigis","isoEm") process.l1ExtraReEmul.nonIsolatedEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") process.l1ExtraReEmul.forwardJetSource = cms.InputTag("gctReEmulDigis","forJets") process.l1ExtraReEmul.centralJetSource = cms.InputTag("gctReEmulDigis","cenJets") process.l1ExtraReEmul.tauJetSource = cms.InputTag("gctReEmulDigis","tauJets") process.l1ExtraReEmul.isoTauJetSource = cms.InputTag("gctReEmulDigis","isoTauJets") process.l1ExtraReEmul.etTotalSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etHadSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.htMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingEtSumsSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingBitCountsSource = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gctCentralJetsSource = cms.InputTag("gctReEmulDigis","cenJets") ntuple.gctNonIsoEmSource = cms.InputTag("gctReEmulDigis","nonIsoEm") ntuple.gctForwardJetsSource = cms.InputTag("gctReEmulDigis","forJets") ntuple.gctIsoEmSource = cms.InputTag("gctReEmulDigis","isoEm") ntuple.gctEnergySumsSource = cms.InputTag("gctReEmulDigis","") ntuple.gctTauJetsSource = cms.InputTag("gctReEmulDigis","tauJets") if useStage1Layer2: ntuple.gctIsoTauJetsSource = cms.InputTag("gctReEmulDigis","isoTauJets") else: ntuple.gctIsoTauJetsSource = cms.InputTag("none","isoTauJets") l1ExtraNtuple.nonIsoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:NonIsolated") l1ExtraNtuple.isoEmLabel = cms.untracked.InputTag("l1ExtraReEmul:Isolated") l1ExtraNtuple.tauJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Tau") l1ExtraNtuple.isoTauJetLabel = cms.untracked.InputTag("l1ExtraReEmul:IsoTau") l1ExtraNtuple.cenJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Central") l1ExtraNtuple.fwdJetLabel = cms.untracked.InputTag("l1ExtraReEmul:Forward") l1ExtraNtuple.metLabel = cms.untracked.InputTag("l1ExtraReEmul:MET") l1ExtraNtuple.mhtLabel = cms.untracked.InputTag("l1ExtraReEmul:MHT") if reEmulRCT : process.gctReEmulDigis.inputLabel = cms.InputTag("rctReEmulDigis") process.reEmulCaloChain = cms.Sequence( #process.hcalReEmulDigis #process.ecalReEmulDigis process.rctReEmulDigis *process.gctReEmulDigis ) else : process.reEmulCaloChain = cms.Sequence( process.gctReEmulDigis ) from L1Trigger.GlobalTrigger.gtDigis_cfi import gtDigis process.gtReEmulDigis = gtDigis.clone() if reEmulMuons : process.gtReEmulDigis.GmtInputTag = cms.InputTag("gmtReEmulDigis") if reEmulCalos : process.gtReEmulDigis.GctInputTag = cms.InputTag("gctReEmulDigis") if patchNtuple : ntuple.gtSource = cms.InputTag("gtReEmulDigis") if reEmulMuons and reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulMuons : process.reEmul = cms.Sequence(process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulCalos : process.reEmul = cms.Sequence(process.reEmulCaloChain + process.gtReEmulDigis + process.l1ExtraReEmul) else : process.reEmul = cms.Sequence(process.gtReEmulDigis + process.l1ExtraReEmul)
def reEmulation(process, reEmulMuons=True, reEmulCalos=True, patchNtuple=True, runOnPostLS1=False): print "[L1Menu]: Setting up overall re-emulation" if patchNtuple and hasattr(process, 'l1NtupleProducer') and hasattr( process, 'l1ExtraTreeProducer'): print "[L1Menu]:\tConfiguring Ntuple to use re-emulated information" ntuple = getattr(process, 'l1NtupleProducer') l1ExtraNtuple = getattr(process, 'l1ExtraTreeProducer') elif patchNtuple: print "[L1Menu]:\tERROR: FAILED to find ntuple! switching patchNtuple to False" patchNtuple = False process.l1ExtraReEmul = cms.EDProducer( "L1ExtraParticlesProd", muonSource=cms.InputTag("gtDigis"), isolatedEmSource=cms.InputTag("gctDigis", "isoEm"), nonIsolatedEmSource=cms.InputTag("gctDigis", "nonIsoEm"), forwardJetSource=cms.InputTag("gctDigis", "forJets"), centralJetSource=cms.InputTag("gctDigis", "cenJets"), tauJetSource=cms.InputTag("gctDigis", "tauJets"), etTotalSource=cms.InputTag("gctDigis"), etHadSource=cms.InputTag("gctDigis"), etMissSource=cms.InputTag("gctDigis"), htMissSource=cms.InputTag("gctDigis"), hfRingEtSumsSource=cms.InputTag("gctDigis"), hfRingBitCountsSource=cms.InputTag("gctDigis"), produceMuonParticles=cms.bool(True), produceCaloParticles=cms.bool(True), centralBxOnly=cms.bool(True), ignoreHtMiss=cms.bool(False)) if reEmulMuons: print "[L1Menu]:\tSetting up muon re-emulation" from L1Trigger.DTTrackFinder.dttfDigis_cfi import dttfDigis process.dttfReEmulDigis = dttfDigis.clone() process.dttfReEmulDigis.DTDigi_Source = cms.InputTag("dttfDigis") process.dttfReEmulDigis.CSCStub_Source = cms.InputTag( "csctfReEmulTrackDigis") from L1Trigger.RPCTrigger.rpcTriggerDigis_cfi import rpcTriggerDigis process.rpcTriggerReEmulDigis = rpcTriggerDigis.clone() if not runOnPostLS1: from L1Trigger.CSCTrackFinder.csctfTrackDigis_cfi import csctfTrackDigis from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis process.csctfReEmulTrackDigis = csctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.readDtDirect = True process.csctfReEmulTrackDigis.SectorReceiverInput = cms.untracked.InputTag( "csctfDigis") process.csctfReEmulTrackDigis.DtDirectProd = cms.untracked.InputTag( "csctfDigis", "DT") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag( "csctfReEmulTrackDigis") #process.csctfReEmulDigis.SectorProcessor.initializeFromPSet = True process.csctfReEmulSequence = cms.Sequence( process.csctfReEmulTrackDigis * process.csctfReEmulDigis) else: from SLHCUpgradeSimulations.Configuration.muonCustoms import customise_csc_L1Emulator from L1Trigger.CSCTrackFinder.csctfDigis_cfi import csctfDigis customise_csc_L1Emulator(process) process.csctfReEmulTrackDigis = process.simCsctfTrackDigis.clone() process.csctfReEmulDigis = csctfDigis.clone() process.csctfReEmulTrackDigis.DTproducer = cms.untracked.InputTag( "dttfDigis") process.csctfReEmulDigis.CSCTrackProducer = cms.untracked.InputTag( "csctfReEmulTrackDigis") process.csctfReEmulSequence = cms.Sequence( process.simCscTriggerPrimitiveDigis * process.csctfReEmulTrackDigis * process.csctfReEmulDigis) from L1Trigger.GlobalMuonTrigger.gmtDigis_cfi import gmtDigis process.gmtReEmulDigis = gmtDigis.clone() process.gmtReEmulDigis.DTCandidates = cms.InputTag( "dttfReEmulDigis", "DT") process.gmtReEmulDigis.CSCCandidates = cms.InputTag( "csctfReEmulDigis", "CSC") process.gmtReEmulDigis.RPCbCandidates = cms.InputTag( "rpcTriggerReEmulDigis", "RPCb") process.gmtReEmulDigis.RPCfCandidates = cms.InputTag( "rpcTriggerReEmulDigis", "RPCf") process.gmtReEmulDigis.MipIsoData = cms.InputTag("none") process.l1ExtraReEmul.muonSource = cms.InputTag("gmtReEmulDigis") if patchNtuple: ntuple.gmtSource = cms.InputTag("gmtReEmulDigis") ntuple.dttfSource = cms.InputTag("dttfReEmulDigis") ntuple.csctfTrkSource = cms.InputTag("csctfReEmulTrackDigis") ntuple.csctfStatusSource = cms.InputTag("csctfReEmulTrackDigis") l1ExtraNtuple.muonLabel = cms.untracked.InputTag("l1ExtraReEmul") process.reEmulMuonChain = cms.Sequence( process.rpcTriggerReEmulDigis * process.csctfReEmulSequence * process.dttfReEmulDigis * process.gmtReEmulDigis) if reEmulCalos: print "[L1Menu]:\tSetting up calo re-emulation" print "[L1Menu]:\tWARNING! Just patching GCT sequence, RCT can't be re-emulated after UCT2015 SW patches" # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly # # In MC HCAL need to be re-run as there is no TPG information stored # process.load("SimCalorimetry.HcalSimProducers.hcalUnsuppressedDigis_cfi") # process.load("SimCalorimetry.HcalTrigPrimProducers.hcaltpdigi_cff") # # from L1Trigger.RegionalCaloTrigger.rctDigis_cfi import rctDigis # from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis # # process.hcalReEmulDigis = process.simHcalTriggerPrimitiveDigis.clone() # process.rctReEmulDigis = rctDigis.clone() # process.gctReEmulDigis = gctDigis.clone() # # process.hcalReEmulDigis.inputLabel = cms.VInputTag(cms.InputTag('hcalDigis'), cms.InputTag('hcalDigis')) # #process.HcalTPGCoderULUT.LUTGenerationMode = cms.bool(False) # # process.rctReEmulDigis.ecalDigis = cms.VInputTag( cms.InputTag( 'ecalDigis:EcalTriggerPrimitives' ) ) # process.rctReEmulDigis.hcalDigis = cms.VInputTag( cms.InputTag( 'hcalReEmulDigis' ) ) # # process.gctReEmulDigis.inputLabel = cms.InputTag("rctReEmulDigis") from L1Trigger.GlobalCaloTrigger.gctDigis_cfi import gctDigis process.gctReEmulDigis = gctDigis.clone() process.gctReEmulDigis.inputLabel = cms.InputTag("gctDigis") process.l1ExtraReEmul.isolatedEmSource = cms.InputTag( "gctReEmulDigis", "isoEm") process.l1ExtraReEmul.nonIsolatedEmSource = cms.InputTag( "gctReEmulDigis", "nonIsoEm") process.l1ExtraReEmul.forwardJetSource = cms.InputTag( "gctReEmulDigis", "forJets") process.l1ExtraReEmul.centralJetSource = cms.InputTag( "gctReEmulDigis", "cenJets") process.l1ExtraReEmul.tauJetSource = cms.InputTag( "gctReEmulDigis", "tauJets") process.l1ExtraReEmul.etTotalSource = cms.InputTag("gctDigis") process.l1ExtraReEmul.etHadSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.etMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.htMissSource = cms.InputTag("gctReEmulDigis") process.l1ExtraReEmul.hfRingEtSumsSource = cms.InputTag( "gctReEmulDigis") process.l1ExtraReEmul.hfRingBitCountsSource = cms.InputTag( "gctReEmulDigis") if patchNtuple: ntuple.gctCentralJetsSource = cms.InputTag("gctReEmulDigis", "cenJets") ntuple.gctNonIsoEmSource = cms.InputTag("gctReEmulDigis", "nonIsoEm") ntuple.gctForwardJetsSource = cms.InputTag("gctReEmulDigis", "forJets") ntuple.gctIsoEmSource = cms.InputTag("gctReEmulDigis", "isoEm") ntuple.gctEnergySumsSource = cms.InputTag("gctReEmulDigis", "") ntuple.gctTauJetsSource = cms.InputTag("gctReEmulDigis", "tauJets") l1ExtraNtuple.nonIsoEmLabel = cms.untracked.InputTag( "l1ExtraReEmul:NonIsolated") l1ExtraNtuple.isoEmLabel = cms.untracked.InputTag( "l1ExtraReEmul:Isolated") l1ExtraNtuple.tauJetLabel = cms.untracked.InputTag( "l1ExtraReEmul:Tau") l1ExtraNtuple.cenJetLabel = cms.untracked.InputTag( "l1ExtraReEmul:Central") l1ExtraNtuple.fwdJetLabel = cms.untracked.InputTag( "l1ExtraReEmul:Forward") l1ExtraNtuple.metLabel = cms.untracked.InputTag( "l1ExtraReEmul:MET") l1ExtraNtuple.mhtLabel = cms.untracked.InputTag( "l1ExtraReEmul:MHT") # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly # ntuple.rctSource = cms.InputTag("rctReEmulDigis") process.reEmulCaloChain = cms.Sequence( # Need to have RCT emulator configurable and not UCT 2015 patches # in order to run 2012 RCT emulator correctly #process.hcalReEmulDigis #+ process.rctReEmulDigis process.gctReEmulDigis) from L1Trigger.GlobalTrigger.gtDigis_cfi import gtDigis process.gtReEmulDigis = gtDigis.clone() if reEmulMuons: process.gtReEmulDigis.GmtInputTag = cms.InputTag("gmtReEmulDigis") if reEmulCalos: process.gtReEmulDigis.GctInputTag = cms.InputTag("gctReEmulDigis") if patchNtuple: ntuple.gtSource = cms.InputTag("gtReEmulDigis") if reEmulMuons and reEmulCalos: process.reEmul = cms.Sequence(process.reEmulCaloChain + process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulMuons: process.reEmul = cms.Sequence(process.reEmulMuonChain + process.gtReEmulDigis + process.l1ExtraReEmul) elif reEmulCalos: process.reEmul = cms.Sequence(process.reEmulCaloChain + process.gtReEmulDigis + process.l1ExtraReEmul) else: process.reEmul = cms.Sequence(process.gtReEmulDigis + process.l1ExtraReEmul)