HiJetAlgos/test/testRegionalJetTracking_cfg.py

0001
0002 import FWCore.ParameterSet.Config as cms
0003 import FWCore.ParameterSet.VarParsing as VarParsing
0004 import os
0005
0006 process = cms.Process("RERECO")
0007
0008 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0009 process.load("Configuration.StandardSequences.Services_cff")
0010 process.load("Configuration.StandardSequences.MagneticField_cff")
0011 process.load("Configuration.StandardSequences.Geometry_cff")
0012
0013 #global tags for conditions data: https://twiki.cern.ch/twiki/bin/view/CMS/SWGuideFrontierConditions
0014 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0015 process.GlobalTag.globaltag = 'MC_3XY_V21::All'
0016
0017 ##################################################################################
0018
0019 # setup 'standard'  options
0020 options = VarParsing.VarParsing ('standard')
0021
0022 # setup any defaults you want
0023 options.output = 'test_out.root'
0024 options.files= '/store/relval/CMSSW_3_5_2/RelValPyquen_DiJet_pt80to120_4TeV/GEN-SIM-RECO/MC_3XY_V21-v1/0001/FA961928-A01F-DF11-85F6-0030487CD7C0.root'
0025 options.maxEvents = 1
0026
0027 # get and parse the command line arguments
0028 options.parseArguments()
0029
0030
0031 ##################################################################################
0032 # Some Services
0033
0034 process.load("FWCore.MessageService.MessageLogger_cfi")
0035 process.MessageLogger.debugModules = ['*']
0036
0037 process.MessageLogger.cerr = cms.untracked.PSet(
0038     threshold = cms.untracked.string('DEBUG'),
0039     DEBUG = cms.untracked.PSet(
0040         limit = cms.untracked.int32(0)
0041     ),
0042     INFO = cms.untracked.PSet(
0043         limit = cms.untracked.int32(0)
0044     ),
0045     HeavyIonVertexing = cms.untracked.PSet(
0046         limit = cms.untracked.int32(-1)
0047     ),
0048     heavyIonHLTVertexing = cms.untracked.PSet(
0049         limit = cms.untracked.int32(-1)
0050     ),
0051     TauRegionalPixelSeedGenerator = cms.untracked.PSet(
0052         limit = cms.untracked.int32(-1)
0053     )
0054 )
0055
0056 process.SimpleMemoryCheck = cms.Service('SimpleMemoryCheck',
0057                                         ignoreTotal=cms.untracked.int32(0),
0058                                         oncePerEventMode = cms.untracked.bool(False)
0059                                         )
0060
0061 process.Timing = cms.Service("Timing")
0062
0063 ##################################################################################
0064 # Input Source
0065 process.source = cms.Source('PoolSource',fileNames = cms.untracked.vstring(options.files))
0066 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(options.maxEvents))
0067
0068 #Reconstruction
0069 process.load("Configuration.StandardSequences.RawToDigi_cff")           # RawToDigi
0070 process.load("Configuration.StandardSequences.ReconstructionHeavyIons_cff") # full heavy ion reconstruction
0071 process.load("RecoHI.HiTracking.secondStep_cff")                            # pair-seeding extension
0072 process.hiNewClusters.Common.maxChi2 = 0 # for testing: remove none of clusters from previous iteration
0073
0074
0075 #### optionally run tracking only seeded by jets
0076 process.goodJets = cms.EDFilter("CaloJetSelector",
0077                                 src = cms.InputTag("iterativeConePu5CaloJets"),
0078                                 cut = cms.string("et > 30.0"),
0079                                 filter = cms.bool(True)
0080                                 )
0081
0082 process.HiTrackingRegionFactoryFromJetsBlock.RegionPSet.JetSrc="goodJets"
0083 process.hiPixel3PrimTracks.RegionFactoryPSet = process.HiTrackingRegionFactoryFromJetsBlock
0084 process.hiNewSeedFromPairs.RegionFactoryPSet = process.HiTrackingRegionFactoryFromJetsBlock
0085
0086 # Output EDM File
0087 process.load("Configuration.EventContent.EventContentHeavyIons_cff")        #load keep/drop output commands
0088 process.output = cms.OutputModule("PoolOutputModule",
0089                                   #process.FEVTDEBUGEventContent,
0090                                   compressionLevel = cms.untracked.int32(2),
0091                                   commitInterval = cms.untracked.uint32(1),
0092                                   fileName = cms.untracked.string(options.output)
0093                                   )
0094
0095 ##################################################################################
0096 # Paths
0097 process.rechits = cms.Sequence(process.siPixelRecHits*process.siStripMatchedRecHits)
0098 process.vtxreco = cms.Sequence(process.offlineBeamSpot * process.trackerlocalreco * process.hiPixelVertices)
0099 process.pxlreco = cms.Sequence(process.vtxreco * process.hiPixel3PrimTracksSequence)
0100 process.trkreco = cms.Sequence(process.offlineBeamSpot * process.trackerlocalreco * process.heavyIonTracking)
0101
0102 process.reco        = cms.Sequence(process.RawToDigi * process.trkreco)
0103 process.rereco      = cms.Sequence(process.rechits * process.heavyIonTracking)
0104 process.rerecojets  = cms.Sequence(process.goodJets * process.rechits * process.heavyIonTracking)
0105
0106 process.path = cms.Path(process.rerecojets   # triplet-seeded regional step
0107                         #* process.secondStep # pair-seeded regional step
0108                         )
0109 process.save = cms.EndPath(process.output)
0110