Configuration/test/RecoToDisplay_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 import sys
0003
0004 process = cms.Process("REPROD")
0005 process.load("Configuration.StandardSequences.Reconstruction_cff")
0006 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0007 #process.load("Configuration.StandardSequences.MagneticField_4T_cff")
0008 process.load("Configuration.StandardSequences.GeometryDB_cff")
0009 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0010 from Configuration.AlCa.autoCond import autoCond
0011 process.GlobalTag.globaltag = autoCond['phase1_2022_realistic']
0012
0013 #process.Timing =cms.Service("Timing")
0014 process.maxEvents = cms.untracked.PSet(
0015     input = cms.untracked.int32(-1)
0016 )
0017
0018 process.source = cms.Source(
0019     "PoolSource",
0020     fileNames = cms.untracked.vstring(
0021     '/store/relval/CMSSW_7_1_0_pre5/RelValTTbar_13/GEN-SIM-RECO/PU25ns_POSTLS171_V1-v2/00000/1A198FA1-B3BC-E311-963F-02163E00CD6B.root'
0022     ),
0023     eventsToProcess = cms.untracked.VEventRange(),
0024     #eventsToProcess = cms.untracked.VEventRange('1:1217421-1:1217421'),
0025     #                                             '1:1220344-1:1220344',
0026     #                                             '1:1655912-1:1655912',
0027     #                                             '1:415027-1:415027',
0028     #                                             '1:460640-1:460640',
0029     #                                             '1:2054772-1:2054772'),
0030     secondaryFileNames = cms.untracked.vstring(),
0031     noEventSort = cms.untracked.bool(True),
0032     duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
0033 )
0034
0035 #from RecoParticleFlow.Configuration.reco_QCDForPF_cff import fileNames
0036 #process.source.fileNames = fileNames
0037
0038 process.dump = cms.EDAnalyzer("EventContentAnalyzer")
0039
0040
0041 process.load("RecoParticleFlow.Configuration.ReDisplay_EventContent_cff")
0042 process.display = cms.OutputModule("PoolOutputModule",
0043     process.DisplayEventContent,
0044     #outputCommands = cms.untracked.vstring('keep *'),
0045     fileName = cms.untracked.string('display.root')
0046 )
0047
0048 process.load("Configuration.EventContent.EventContent_cff")
0049 process.reco = cms.OutputModule("PoolOutputModule",
0050     process.RECOSIMEventContent,
0051     fileName = cms.untracked.string('/afs/cern.ch/user/l/lgray/work/public/phogun_35GeV.root')
0052 )
0053
0054 # modify reconstruction sequence
0055 #process.hbhereflag = process.hbhereco.clone()
0056 #process.hbhereflag.hbheInput = 'hbhereco'
0057 #process.towerMakerPF.hbheInput = 'hbhereflag'
0058 #process.particleFlowRecHitHCAL.hcalRecHitsHBHE = cms.InputTag("hbhereflag")
0059 process.particleFlowTmp.useHO = False
0060
0061 # Local re-reco: Produce tracker rechits, pf rechits and pf clusters
0062 process.localReReco = cms.Sequence(process.siPixelRecHits+
0063                                    process.siStripMatchedRecHits+
0064                                    #process.hbhereflag+
0065                                    process.particleFlowCluster+
0066                                    process.ecalClusters)
0067
0068 # Track re-reco
0069 process.globalReReco =  cms.Sequence(process.offlineBeamSpot+
0070                                      process.MeasurementTrackerEvent+
0071                                      process.recopixelvertexing+
0072                                      process.ckftracks+
0073                                      process.caloTowersRec+
0074                                      process.vertexreco+
0075                                      process.recoJets+
0076                                      process.muonrecoComplete+
0077                                      process.muoncosmicreco+
0078                                      process.egammaGlobalReco+
0079                                      process.pfTrackingGlobalReco+
0080                                      process.egammaHighLevelRecoPrePF+
0081                                      process.muoncosmichighlevelreco+
0082                                      process.metreco)
0083
0084 # Particle Flow re-processing
0085 process.pfReReco = cms.Sequence(process.particleFlowReco+
0086                                 process.egammaHighLevelRecoPostPF+
0087                                 process.muonshighlevelreco+
0088                                 process.particleFlowLinks+
0089                                 process.recoPFJets+
0090                                 process.recoPFMET+
0091                                 process.PFTau)
0092 #pfeg switch
0093 def switch_on_pfeg(the_process):
0094     the_process.particleFlowTmp.useEGammaFilters = cms.bool(True)
0095     the_process.particleFlow.GsfElectrons = cms.InputTag('gedGsfElectrons')
0096     the_process.particleFlow.Photons = cms.InputTag('gedPhotons')
0097     the_process.particleFlowReco.remove(the_process.pfElectronTranslatorSequence)
0098     the_process.particleFlowReco.remove(the_process.pfPhotonTranslatorSequence)
0099     the_process.egammaHighLevelRecoPostPF.remove(the_process.gsfElectronMergingSequence)
0100     the_process.interestingEleIsoDetIdEB.emObjectLabel = \
0101                                               cms.InputTag('gedGsfElectrons')
0102     the_process.interestingEleIsoDetIdEE.emObjectLabel = \
0103                                               cms.InputTag('gedGsfElectrons')
0104     the_process.interestingGamIsoDetIdEB.emObjectLabel = \
0105                                               cms.InputTag('gedPhotons')
0106     the_process.interestingGamIsoDetIdEE.emObjectLabel = \
0107                                               cms.InputTag('gedPhotons')
0108     the_process.PhotonIDProd.photonProducer = cms.string('gedPhotons')
0109     the_process.eidRobustLoose.src = cms.InputTag('gedGsfElectrons')
0110     the_process.eidRobustTight.src = cms.InputTag('gedGsfElectrons')
0111     the_process.eidRobustHighEnergy.src = cms.InputTag('gedGsfElectrons')
0112     the_process.eidLoose.src = cms.InputTag('gedGsfElectrons')
0113     the_process.eidTight.src = cms.InputTag('gedGsfElectrons')
0114 switch_on_pfeg(process)
0115
0116 # Gen Info re-processing
0117 process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")
0118 process.load("RecoJets.Configuration.GenJetParticles_cff")
0119 process.load("RecoJets.Configuration.RecoGenJets_cff")
0120 process.load("RecoMET.Configuration.GenMETParticles_cff")
0121 process.load("RecoMET.Configuration.RecoGenMET_cff")
0122 process.load("RecoParticleFlow.PFProducer.particleFlowSimParticle_cff")
0123 process.load("RecoParticleFlow.Configuration.HepMCCopy_cfi")
0124 process.genReReco = cms.Sequence(process.generator+
0125                                  process.genParticles+
0126                                  process.genJetParticles+
0127                                  process.recoGenJets+
0128                                  process.genMETParticles+
0129                                  process.recoGenMET+
0130                                  process.particleFlowSimParticle)
0131
0132 #process.load("RecoParticleFlow.PFProducer.particleFlowCandidateChecker_cfi")
0133 #process.particleFlowCandidateChecker.pfCandidatesReco = cms.InputTag("particleFlow","","REPROD")
0134 #process.particleFlowCandidateChecker.pfCandidatesReReco = cms.InputTag("particleFlow","","REPROD2")
0135 #process.particleFlowCandidateChecker.pfJetsReco = cms.InputTag("ak5PFJets","","REPROD")
0136 #process.particleFlowCandidateChecker.pfJetsReReco = cms.InputTag("ak5PFJets","","REPROD2")
0137 # The complete reprocessing
0138 process.p = cms.Path(process.localReReco+
0139                      process.globalReReco+
0140                      process.pfReReco+
0141                      process.genReReco
0142                      #+process.particleFlowEGammaCandidateChecker
0143                      )
0144
0145 # And the output.
0146 process.outpath = cms.EndPath(
0147     process.reco +
0148     process.display
0149 )
0150
0151 # And the monitoring
0152 process.SimpleMemoryCheck = cms.Service("SimpleMemoryCheck",
0153                                         ignoreTotal=cms.untracked.int32(1),
0154                                         jobReportOutputOnly = cms.untracked.bool(True)
0155                                         )
0156 process.Timing = cms.Service("Timing",
0157                              summaryOnly = cms.untracked.bool(True)
0158                              )
0159
0160 # And the logger
0161 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0162 process.options = cms.untracked.PSet(
0163     makeTriggerResults = cms.untracked.bool(True),
0164     wantSummary = cms.untracked.bool(True),
0165     Rethrow = cms.untracked.vstring('Unknown',
0166         'ProductNotFound',
0167         'DictionaryNotFound',
0168         'InsertFailure',
0169         'Configuration',
0170         'LogicError',
0171         'UnimplementedFeature',
0172         'InvalidReference',
0173         'NullPointerError',
0174         'NoProductSpecified',
0175         'EventTimeout',
0176         'EventCorruption',
0177         'ModuleFailure',
0178         'ScheduleExecutionFailure',
0179         'EventProcessorFailure',
0180         'FileInPathError',
0181         'FatalRootError',
0182         'NotFound')
0183 )
0184
0185 process.MessageLogger.cerr.FwkReport.reportEvery = 1
0186 #process.MessageLogger.cout = cms.untracked.PSet(
0187 #    threshold = cms.untracked.string('INFO')
0188 #    )
0189
0190