HGCalTestBeam/test/HGCalTBGenSim_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process('SIM')
0004
0005 # import of standard configurations
0006 process.load('Configuration.StandardSequences.Services_cff')
0007 process.load('SimGeneral.HepPDTESSource.pythiapdt_cfi')
0008 process.load('Configuration.EventContent.EventContent_cff')
0009 process.load('SimGeneral.MixingModule.mixNoPU_cfi')
0010 process.load('SimG4CMS.HGCalTestBeam.HGCalTB160Module16XML_cfi')
0011 process.load('Geometry.HGCalCommonData.hgcalNumberingInitialization_cfi')
0012 process.load('Geometry.HGCalCommonData.hgcalParametersInitialization_cfi')
0013 process.load('Configuration.StandardSequences.MagneticField_0T_cff')
0014 process.load('Configuration.StandardSequences.Generator_cff')
0015 process.load('IOMC.EventVertexGenerators.VtxSmearedGauss_cfi')
0016 process.load('GeneratorInterface.Core.genFilterSummary_cff')
0017 process.load('Configuration.StandardSequences.SimIdeal_cff')
0018 process.load('Configuration.StandardSequences.EndOfProcess_cff')
0019 process.load('Configuration.StandardSequences.FrontierConditions_GlobalTag_cff')
0020 process.load('SimG4CMS.HGCalTestBeam.HGCalTBAnalyzer_cfi')
0021
0022 process.maxEvents = cms.untracked.PSet(
0023     input = cms.untracked.int32(1000)
0024 )
0025
0026 process.MessageLogger = cms.Service("MessageLogger",
0027     cout = cms.untracked.PSet(
0028         HGCSim = cms.untracked.PSet(
0029             limit = cms.untracked.int32(-1)
0030         ),
0031         default = cms.untracked.PSet(
0032             limit = cms.untracked.int32(0)
0033         ),
0034         enable = cms.untracked.bool(True)
0035     )
0036 )
0037
0038 # Input source
0039 process.source = cms.Source("EmptySource")
0040
0041 process.options = cms.untracked.PSet(
0042
0043 )
0044
0045 # Production Info
0046 process.configurationMetadata = cms.untracked.PSet(
0047     annotation = cms.untracked.string('SingleElectronE4_cfi'),
0048     name = cms.untracked.string('Applications'),
0049     version = cms.untracked.string('$Revision: 1.19 $')
0050 )
0051
0052 # Output definition
0053
0054 process.RAWSIMoutput = cms.OutputModule("PoolOutputModule",
0055     SelectEvents = cms.untracked.PSet(
0056         SelectEvents = cms.vstring('generation_step')
0057     ),
0058     dataset = cms.untracked.PSet(
0059         dataTier = cms.untracked.string('GEN-SIM'),
0060         filterName = cms.untracked.string('')
0061     ),
0062     eventAutoFlushCompressedSize = cms.untracked.int32(5242880),
0063     fileName = cms.untracked.string('file:gensim.root'),
0064     outputCommands = process.RAWSIMEventContent.outputCommands,
0065     splitLevel = cms.untracked.int32(0)
0066 )
0067
0068 # Additional output definition
0069 process.TFileService = cms.Service("TFileService",
0070                                    fileName = cms.string('TBGenSim.root')
0071                                    )
0072
0073 # Other statements
0074 process.genstepfilter.triggerConditions=cms.vstring("generation_step")
0075 from Configuration.AlCa.GlobalTag import GlobalTag
0076 process.GlobalTag = GlobalTag(process.GlobalTag, 'auto:run2_mc', '')
0077
0078 process.generator = cms.EDProducer("GaussRandomPThetaGunProducer",
0079     AddAntiParticle = cms.bool(False),
0080     PGunParameters = cms.PSet(
0081         MeanP    = cms.double(4.0),
0082         SigmaP   = cms.double(0.108),
0083         MinTheta = cms.double(0.0),
0084         MaxTheta = cms.double(0.0),
0085         MinPhi = cms.double(-3.14159265359),
0086         MaxPhi = cms.double(3.14159265359),
0087         PartID = cms.vint32(11)
0088     ),
0089     Verbosity = cms.untracked.int32(0),
0090     firstRun = cms.untracked.uint32(1),
0091     psethack = cms.string('single electron E 4')
0092 )
0093 process.VtxSmeared.MeanX  = 0.0
0094 process.VtxSmeared.SigmaX = 0.55
0095 process.VtxSmeared.MeanY  = 0.0
0096 process.VtxSmeared.SigmaY = 0.65
0097 process.VtxSmeared.MeanZ  = -3500.0
0098 process.VtxSmeared.SigmaZ = 0
0099 process.HGCalTBAnalyzer.doDigis = False
0100 process.HGCalTBAnalyzer.doRecHits = False
0101 process.g4SimHits.OnlySDs = ['HGCSensitiveDetector',
0102                              'HGCalTB1601SensitiveDetector']
0103 process.g4SimHits.HGCSD.Detectors = 1
0104 process.g4SimHits.HGCSD.UseDetector = 1
0105
0106 # Path and EndPath definitions
0107 process.generation_step = cms.Path(process.pgen)
0108 process.simulation_step = cms.Path(process.psim)
0109 process.genfiltersummary_step = cms.EndPath(process.genFilterSummary)
0110 process.analysis_step = cms.Path(process.HGCalTBAnalyzer)
0111 process.endjob_step = cms.EndPath(process.endOfProcess)
0112 process.RAWSIMoutput_step = cms.EndPath(process.RAWSIMoutput)
0113
0114 # Schedule definition
0115 process.schedule = cms.Schedule(process.generation_step,process.genfiltersummary_step,process.simulation_step,process.analysis_step,process.endjob_step,process.RAWSIMoutput_step)
0116 # filter all path with the production filter sequence
0117 for path in process.paths:
0118     getattr(process,path)._seq = process.generator * getattr(process,path)._seq
0119
0120