Configuration/test/RecoToDisplay_batch_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002 import sys
0003
0004 infile_name = sys.argv[1]
0005 outfile_name = '/afs/cern.ch/user/l/lgray/work/public/CMSSW_7_0_0_pre3_singlegconv/src/RecoParticleFlow/Configuration/test/%s/superClusterDump_%i.root'%(sys.argv[4],int(sys.argv[2]))
0006 nevents = int(sys.argv[3])
0007
0008
0009 process = cms.Process("REPROD")
0010 process.load("Configuration.StandardSequences.Reconstruction_cff")
0011 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0012 #process.load("Configuration.StandardSequences.MagneticField_4T_cff")
0013 process.load("Configuration.StandardSequences.GeometryRecoDB_cff")
0014 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0015 from Configuration.AlCa.autoCond import autoCond
0016 process.GlobalTag.globaltag = autoCond['phase1_2022_realistic']
0017
0018 #process.Timing =cms.Service("Timing")
0019 process.maxEvents = cms.untracked.PSet(
0020     input = cms.untracked.int32(nevents)
0021 )
0022
0023 process.source = cms.Source(
0024     "PoolSource",
0025     fileNames = cms.untracked.vstring(
0026     'root://cms-xrd-global.cern.ch/%s'%infile_name
0027     ),
0028     eventsToProcess = cms.untracked.VEventRange(),
0029     #eventsToProcess = cms.untracked.VEventRange('1:1217421-1:1217421'),
0030     #                                             '1:1220344-1:1220344',
0031     #                                             '1:1655912-1:1655912',
0032     #                                             '1:415027-1:415027',
0033     #                                             '1:460640-1:460640',
0034     #                                             '1:2054772-1:2054772'),
0035     secondaryFileNames = cms.untracked.vstring(),
0036     noEventSort = cms.untracked.bool(True),
0037     duplicateCheckMode = cms.untracked.string('noDuplicateCheck')
0038 )
0039
0040 #from RecoParticleFlow.Configuration.reco_QCDForPF_cff import fileNames
0041 #process.source.fileNames = fileNames
0042
0043
0044 # modify reconstruction sequence
0045 #process.hbhereflag = process.hbhereco.clone()
0046 #process.hbhereflag.hbheInput = 'hbhereco'
0047 #process.towerMakerPF.hbheInput = 'hbhereflag'
0048 #process.particleFlowRecHitHCAL.hcalRecHitsHBHE = cms.InputTag("hbhereflag")
0049 process.particleFlowTmp.useHO = False
0050
0051 process.TFileService = cms.Service(
0052     "TFileService",
0053     fileName=cms.string(outfile_name)
0054     )
0055
0056 process.pfElectronSCTree = cms.EDAnalyzer(
0057     "PFSuperClusterTreeMaker",
0058     doGen = cms.untracked.bool(True),
0059     genSrc = cms.InputTag("genParticles"),
0060     primaryVertices = cms.InputTag("offlinePrimaryVertices"),
0061     superClusterSrcEB = cms.InputTag('pfElectronTranslator','pf')
0062     )
0063
0064 process.pfPhotonSCTree = cms.EDAnalyzer(
0065     "PFSuperClusterTreeMaker",
0066     doGen = cms.untracked.bool(True),
0067     genSrc = cms.InputTag("genParticles"),
0068     primaryVertices = cms.InputTag("offlinePrimaryVertices"),
0069     superClusterSrcEB = cms.InputTag('pfPhotonTranslator','pfphot')
0070     )
0071
0072 process.pfMustacheSCTree = cms.EDAnalyzer(
0073     "PFSuperClusterTreeMaker",
0074     doGen = cms.untracked.bool(True),
0075     genSrc = cms.InputTag("genParticles"),
0076     primaryVertices = cms.InputTag("offlinePrimaryVertices"),
0077     superClusterSrcEB = cms.InputTag("particleFlowSuperClusterECAL",
0078                                      "particleFlowSuperClusterECALBarrel"),
0079     superClusterSrcEE = cms.InputTag("particleFlowSuperClusterECAL",
0080                                      "particleFlowSuperClusterECALEndcapWithPreshower")
0081     )
0082
0083 # Local re-reco: Produce tracker rechits, pf rechits and pf clusters
0084 process.localReReco = cms.Sequence(process.siPixelRecHits+
0085                                    process.siStripMatchedRecHits+
0086                                    #process.hbhereflag+
0087                                    process.particleFlowCluster+
0088                                    process.ecalClusters)
0089
0090 # Track re-reco
0091 process.globalReReco =  cms.Sequence(process.offlineBeamSpot+
0092                                      process.recopixelvertexing+
0093                                      process.ckftracks+
0094                                      process.caloTowersRec+
0095                                      process.vertexreco+
0096                                      process.recoJets+
0097                                      process.muonrecoComplete+
0098                                      process.muoncosmicreco+
0099                                      process.egammaGlobalReco+
0100                                      process.pfTrackingGlobalReco+
0101                                      process.egammaHighLevelRecoPrePF+
0102                                      process.muoncosmichighlevelreco+
0103                                      process.metreco)
0104
0105 # Particle Flow re-processing
0106 process.pfReReco = cms.Sequence(process.particleFlowReco+
0107                                 process.egammaHighLevelRecoPostPF+
0108                                 process.muonshighlevelreco+
0109                                 process.particleFlowLinks+
0110                                 process.recoPFJets+
0111                                 process.recoPFMET+
0112                                 process.PFTau+
0113                                 process.pfElectronSCTree+
0114                                 process.pfPhotonSCTree+
0115                                 process.pfMustacheSCTree)
0116
0117
0118 # Gen Info re-processing
0119 process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")
0120 process.load("RecoJets.Configuration.GenJetParticles_cff")
0121 process.load("RecoJets.Configuration.RecoGenJets_cff")
0122 process.load("RecoMET.Configuration.GenMETParticles_cff")
0123 process.load("RecoMET.Configuration.RecoGenMET_cff")
0124 process.load("RecoParticleFlow.PFProducer.particleFlowSimParticle_cff")
0125 process.load("RecoParticleFlow.Configuration.HepMCCopy_cfi")
0126 process.genReReco = cms.Sequence(process.generator+
0127                                  process.genParticles+
0128                                  process.genJetParticles+
0129                                  process.recoGenJets+
0130                                  process.genMETParticles+
0131                                  process.recoGenMET+
0132                                  process.particleFlowSimParticle)
0133
0134 # pf reprocessing
0135 process.p = cms.Path(process.localReReco+
0136                      process.globalReReco+
0137                      process.pfReReco+
0138                      process.genReReco
0139                      #+process.particleFlowCandidateChecker
0140                      )
0141
0142 # And the logger
0143 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0144 #process.MessageLogger.cout = cms.untracked.PSet(
0145 #    threshold = cms.untracked.string('INFO')
0146 #    )
0147 process.options = cms.untracked.PSet(
0148     makeTriggerResults = cms.untracked.bool(True),
0149     #wantSummary = cms.untracked.bool(True),
0150     Rethrow = cms.untracked.vstring('Unknown',
0151         'ProductNotFound',
0152         'DictionaryNotFound',
0153         'InsertFailure',
0154         'Configuration',
0155         'LogicError',
0156         'UnimplementedFeature',
0157         'InvalidReference',
0158         'NullPointerError',
0159         'NoProductSpecified',
0160         'EventTimeout',
0161         'EventCorruption',
0162         'ModuleFailure',
0163         'ScheduleExecutionFailure',
0164         'EventProcessorFailure',
0165         'FileInPathError',
0166         'FatalRootError',
0167         'NotFound')
0168 )
0169
0170 process.MessageLogger.cerr.FwkReport.reportEvery = 1
0171
0172