test/python/testHGCalSimWatcherV14_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.Eras.Era_Phase2C11I13M9_cff import Phase2C11I13M9
0004 process = cms.Process('testHGCalRecoLocal',Phase2C11I13M9)
0005
0006 # import of standard configurations
0007 process.load('Configuration.StandardSequences.Services_cff')
0008 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0009 process.load('FWCore.MessageService.MessageLogger_cfi')
0010 process.load('Configuration.EventContent.EventContent_cff')
0011 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0012 process.load('Configuration.Geometry.GeometryExtended2026D77Reco_cff')
0013 process.load('Configuration.Geometry.GeometryExtended2026D77_cff')
0014 process.load('Configuration.StandardSequences.MagneticField_cff')
0015 process.load('Configuration.StandardSequences.Generator_cff')
0016 process.load('IOMC.EventVertexGenerators.VtxSmearedRealistic50ns13TeVCollision_cfi')
0017 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0018 process.load('Configuration.StandardSequences.SimIdeal_cff')
0019 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0020 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0021 process.load('Configuration.StandardSequences.Digi_cff')
0022 process.load('Configuration.StandardSequences.SimL1Emulator_cff')
0023 process.load('Configuration.StandardSequences.L1TrackTrigger_cff')
0024 process.load('Configuration.StandardSequences.DigiToRaw_cff')
0025 process.load('HLTrigger.Configuration.HLT_Fake2_cff')
0026 process.load('Configuration.StandardSequences.RawToDigi_cff')
0027 process.load('Configuration.StandardSequences.L1Reco_cff')
0028 process.load('Configuration.StandardSequences.Reconstruction_cff')
0029 process.load('Configuration.StandardSequences.RecoSim_cff')
0030 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0031
0032 process.maxEvents = cms.untracked.PSet(
0033     input = cms.untracked.int32(1000)
0034 )
0035
0036 process.MessageLogger.cerr.FwkReport.reportEvery = 5
0037 if hasattr(process,'MessageLogger'):
0038     process.MessageLogger.ValidHGCal=dict()
0039     process.MessageLogger.HGCalGeom=dict()
0040
0041 # Input source
0042 process.source = cms.Source("EmptySource")
0043
0044 process.options = cms.untracked.PSet(
0045     wantSummary = cms.untracked.bool(True),
0046     numberOfConcurrentRuns = cms.untracked.uint32(1),
0047     numberOfStreams = cms.untracked.uint32(0),
0048     numberOfThreads = cms.untracked.uint32(1),
0049     printDependencies = cms.untracked.bool(False),
0050     sizeOfStackForThreadsInKB = cms.optional.untracked.uint32,
0051 )
0052
0053 # Production Info
0054 process.configurationMetadata = cms.untracked.PSet(
0055     version = cms.untracked.string(''),
0056     annotation = cms.untracked.string(''),
0057     name = cms.untracked.string('Applications')
0058 )
0059
0060
0061 # Output definition
0062 process.output = cms.OutputModule("PoolOutputModule",
0063     splitLevel = cms.untracked.int32(0),
0064     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0065     outputCommands = cms.untracked.vstring(
0066         'keep *_*hbhe*_*_*',
0067     'keep *_g4SimHits_*_*',
0068 #       'keep *_mix_*_*',
0069     'keep *_*HGC*_*_*',
0070         ),
0071     fileName = cms.untracked.string('file:testHGCalSimWatcherV14.root'),
0072     dataset = cms.untracked.PSet(
0073         filterName = cms.untracked.string(''),
0074         dataTier = cms.untracked.string('GEN-SIM-DIGI-RAW-RECO')
0075     ),
0076     SelectEvents = cms.untracked.PSet(
0077         SelectEvents = cms.vstring('generation_step')
0078     )
0079 )
0080
0081 # Additional output definition
0082 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0083     SimG4HGCalValidation = cms.PSet(
0084         Names = cms.vstring(
0085         'HGCalEESensitive',
0086         'HGCalHESensitive',
0087                 'HGCalHESiliconSensitive',
0088         'HGCalHEScintillatorSensitive',
0089         ),
0090     Types          = cms.vint32(0,0,0,0),
0091     DetTypes       = cms.vint32(0,1,1,2),
0092     LabelLayerInfo = cms.string("HGCalInfoLayer"),
0093     Verbosity      = cms.untracked.int32(0),
0094     ),
0095     type = cms.string('SimG4HGCalValidation')
0096 ))
0097
0098 # Other statements
0099 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0100 from Configuration.AlCa.GlobalTag import GlobalTag
0101 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:phase2_realistic_T21', '')
0102
0103 process.generator = cms.EDProducer("FlatRandomPtGunProducer",
0104     PGunParameters = cms.PSet(
0105         MaxPt = cms.double(20.0),
0106         MinPt = cms.double(20.0),
0107         #PartID = cms.vint32(11), #--->electron
0108         PartID = cms.vint32(13), #--->muon
0109         #PartID = cms.vint32(211), #--->pion
0110         MaxEta = cms.double(3.0),
0111         MaxPhi = cms.double(3.14159265359),
0112         MinEta = cms.double(1.2),
0113         MinPhi = cms.double(-3.14159265359)
0114     ),
0115     Verbosity = cms.untracked.int32(0),
0116     psethack = cms.string('single muon pt 35'),
0117     AddAntiParticle = cms.bool(False),
0118     firstRun = cms.untracked.uint32(1)
0119 )
0120
0121
0122 #Modified to produce hgceedigis
0123 process.mix.digitizers = cms.PSet(process.theDigitizersValid)
0124
0125 process.ProductionFilterSequence = cms.Sequence(process.generator)
0126
0127 # Path and EndPath definitions
0128 process.generation_step = cms.Path(process.pgen)
0129 process.simulation_step = cms.Path(process.psim)
0130 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0131 process.digitisation_step = cms.Path(process.pdigi_valid)
0132 process.L1simulation_step = cms.Path(process.SimL1Emulator)
0133 process.L1TrackTrigger_step = cms.Path(process.L1TrackTrigger)
0134 process.digi2raw_step = cms.Path(process.DigiToRaw)
0135 process.raw2digi_step = cms.Path(process.RawToDigi)
0136 process.L1Reco_step = cms.Path(process.L1Reco)
0137 process.reconstruction_step = cms.Path(process.localreco)
0138 process.recosim_step = cms.Path(process.recosim)
0139 process.out_step = cms.EndPath(process.output)
0140
0141 # Schedule definition
0142 process.schedule = cms.Schedule(process.generation_step,
0143                 process.simulation_step,
0144                 process.digitisation_step,
0145                                 process.L1simulation_step,
0146                                 process.L1TrackTrigger_step,
0147                                 process.digi2raw_step,
0148                                 process.raw2digi_step,
0149                                 process.L1Reco_step,
0150                                 process.reconstruction_step,
0151                                 process.recosim_step,
0152                                 process.out_step
0153                 )
0154
0155 # filter all path with the production filter sequence
0156 for path in process.paths:
0157         if getattr(process,path)._seq is not None: getattr(process,path)._seq = process.ProductionFilterSequence * getattr(process,path)._seq