test/python/runPbModel_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("CaloTest")
0004 process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0005
0006 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0007
0008 process.load("SimG4Core.CheckSecondary.PbTarget_cfi")
0009
0010 process.load("SimG4Core.Application.g4SimHits_cfi")
0011
0012 process.MessageLogger = cms.Service("MessageLogger",
0013     destinations = cms.untracked.vstring('cout'),
0014     categories = cms.untracked.vstring('SimG4CoreApplication', 'CheckSecondary'),
0015     cout = cms.untracked.PSet(
0016         threshold = cms.untracked.string('INFO'),
0017         default = cms.untracked.PSet(
0018             limit = cms.untracked.int32(-1)
0019         ),
0020         CheckSecondary = cms.untracked.PSet(
0021             limit = cms.untracked.int32(0)
0022         ),
0023         SimG4CoreApplication = cms.untracked.PSet(
0024             limit = cms.untracked.int32(0)
0025         )
0026     )
0027 )
0028
0029 process.maxEvents = cms.untracked.PSet(
0030     input = cms.untracked.int32(1000000)
0031 )
0032
0033 process.source = cms.Source("EmptySource")
0034
0035 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0036     PGunParameters = cms.PSet(
0037         MinEta = cms.double(0.0),
0038         MaxEta = cms.double(0.0),
0039         MinPhi = cms.double(1.57079632679),
0040         MaxPhi = cms.double(1.57079632679),
0041 #        PartID = cms.vint32(-211),
0042 #        MinE   = cms.double(1.40694),
0043 #        MaxE   = cms.double(1.40694)
0044 #        MinE   = cms.double(5.00195),
0045 #        MaxE   = cms.double(5.00195)
0046         PartID = cms.vint32(2212),
0047 #        MinE   = cms.double(1.3713),
0048 #        MaxE   = cms.double(1.3713)
0049 #        MinE   = cms.double(1.68535),
0050 #        MaxE   = cms.double(1.68535)
0051 #        MinE   = cms.double(2.2092),
0052 #        MaxE   = cms.double(2.2092)
0053 #        MinE   = cms.double(3.1433),
0054 #        MaxE   = cms.double(3.1433)
0055 #        MinE   = cms.double(5.0873),
0056 #        MaxE   = cms.double(5.0873)
0057 #        MinE   = cms.double(6.0729),
0058 #        MaxE   = cms.double(6.0729)
0059 #        MinE   = cms.double(6.5674),
0060 #        MaxE   = cms.double(6.5674)
0061 #        MinE   = cms.double(7.0626),
0062 #        MaxE   = cms.double(7.0626)
0063         MinE   = cms.double(7.5585),
0064         MaxE   = cms.double(7.5585)
0065 #        MinE   = cms.double(8.3032),
0066 #        MaxE   = cms.double(8.3032)
0067 #        MinE   = cms.double(9.0488),
0068 #        MaxE   = cms.double(9.0488)
0069     ),
0070     Verbosity = cms.untracked.int32(0),
0071     AddAntiParticle = cms.bool(False),
0072     firstRun        = cms.untracked.uint32(1)
0073 )
0074
0075 process.Timing = cms.Service("Timing")
0076
0077 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0078     moduleSeeds = cms.PSet(
0079         generator = cms.untracked.uint32(456789),
0080         g4SimHits = cms.untracked.uint32(9876),
0081         VtxSmeared = cms.untracked.uint32(123456789)
0082     ),
0083     sourceSeed = cms.untracked.uint32(135799753)
0084 )
0085
0086 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.g4SimHits)
0087 process.VtxSmeared.SigmaX = 0.00001
0088 process.VtxSmeared.SigmaY = 0.00001
0089 process.VtxSmeared.SigmaZ = 0.00001
0090 process.g4SimHits.UseMagneticField     = False
0091 process.g4SimHits.Physics.type         = 'SimG4Core/Physics/CMSModel'
0092 process.g4SimHits.Physics.Model        = 'Bertini'
0093 process.g4SimHits.Physics.EMPhysics    = False
0094 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0095     CheckSecondary = cms.PSet(
0096         SaveInFile = cms.untracked.string('PbBertini67.5GeV.root'),
0097         Verbosity = cms.untracked.int32(0),
0098         MinimumDeltaE = cms.untracked.double(0.0),
0099         KillAfter = cms.untracked.int32(1)
0100     ),
0101     type = cms.string('CheckSecondary')
0102 ),
0103     cms.PSet(
0104         type = cms.string('KillSecondariesRunAction')
0105     ))
0106