test/python/runPbWO4_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.PbWO4Geom_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(20000)
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.010),
0043 #   MaxE   = cms.double(1.010)
0044 #   MinE   = cms.double(2.005),
0045 #   MaxE   = cms.double(2.005)
0046 #   MinE   = cms.double(3.003),
0047 #   MaxE   = cms.double(3.003)
0048 #   MinE   = cms.double(4.002),
0049 #   MaxE   = cms.double(4.002)
0050 #   MinE   = cms.double(5.002),
0051 #   MaxE   = cms.double(5.002)
0052 #   MinE   = cms.double(6.002),
0053 #   MaxE   = cms.double(6.002)
0054 #   MinE   = cms.double(7.001),
0055 #   MaxE   = cms.double(7.001)
0056 #   MinE   = cms.double(8.001),
0057 #   MaxE   = cms.double(8.001)
0058 #   MinE   = cms.double(9.001),
0059 #   MaxE   = cms.double(9.001)
0060 #   MinE   = cms.double(10.001),
0061 #   MaxE   = cms.double(10.001)
0062 #   MinE   = cms.double(15.001),
0063 #   MaxE   = cms.double(15.001)
0064 #   MinE   = cms.double(20.001),
0065 #   MaxE   = cms.double(20.001)
0066 #   MinE   = cms.double(30.001),
0067 #   MaxE   = cms.double(30.001)
0068 #   MinE   = cms.double(50.001),
0069 #   MaxE   = cms.double(50.001)
0070 #   MinE   = cms.double(100.0),
0071 #   MaxE   = cms.double(100.0)
0072 #   MinE   = cms.double(150.0),
0073 #   MaxE   = cms.double(150.0)
0074 #   MinE   = cms.double(200.0),
0075 #   MaxE   = cms.double(200.0)
0076 #   MinE   = cms.double(300.0),
0077 #   MaxE   = cms.double(300.0)
0078         PartID = cms.vint32(2212),
0079     MinE   = cms.double(1.3713),
0080     MaxE   = cms.double(1.3713)
0081 #   MinE   = cms.double(2.2092),
0082 #   MaxE   = cms.double(2.2092)
0083 #   MinE   = cms.double(3.1433),
0084 #   MaxE   = cms.double(3.1433)
0085 #   MinE   = cms.double(4.1086),
0086 #   MaxE   = cms.double(4.1086)
0087 #   MinE   = cms.double(5.0873),
0088 #   MaxE   = cms.double(5.0873)
0089 #   MinE   = cms.double(6.0729),
0090 #   MaxE   = cms.double(6.0729)
0091 #   MinE   = cms.double(7.0626),
0092 #   MaxE   = cms.double(7.0626)
0093 #   MinE   = cms.double(8.0548),
0094 #   MaxE   = cms.double(8.0548)
0095 #   MinE   = cms.double(9.0488),
0096 #   MaxE   = cms.double(9.0488)
0097 #   MinE   = cms.double(10.0439),
0098 #   MaxE   = cms.double(10.0439)
0099 #   MinE   = cms.double(15.0293),
0100 #   MaxE   = cms.double(15.0293)
0101 #   MinE   = cms.double(20.0220),
0102 #   MaxE   = cms.double(20.0220)
0103 #   MinE   = cms.double(30.0147),
0104 #   MaxE   = cms.double(30.0147)
0105 #   MinE   = cms.double(50.0088),
0106 #   MaxE   = cms.double(50.0088)
0107 #   MinE   = cms.double(100.0044),
0108 #   MaxE   = cms.double(100.0044)
0109 #   MinE   = cms.double(150.0029),
0110 #   MaxE   = cms.double(150.0029)
0111 #   MinE   = cms.double(200.0022),
0112 #   MaxE   = cms.double(200.0022)
0113 #   MinE   = cms.double(300.001),
0114 #   MaxE   = cms.double(300.001)
0115     ),
0116     Verbosity = cms.untracked.int32(0),
0117     AddAntiParticle = cms.bool(False),
0118     firstRun        = cms.untracked.uint32(1)
0119 )
0120
0121 process.Timing = cms.Service("Timing")
0122
0123 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0124     moduleSeeds = cms.PSet(
0125         generator = cms.untracked.uint32(456789),
0126         g4SimHits = cms.untracked.uint32(9876),
0127         VtxSmeared = cms.untracked.uint32(123456789)
0128     ),
0129     sourceSeed = cms.untracked.uint32(135799753)
0130 )
0131
0132 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.g4SimHits)
0133 process.VtxSmeared.SigmaX = 0.00001
0134 process.VtxSmeared.SigmaY = 0.00001
0135 process.VtxSmeared.SigmaZ = 0.00001
0136 process.g4SimHits.UseMagneticField     = False
0137 process.g4SimHits.Physics.type         = 'SimG4Core/Physics/QGSP_BERT'
0138 process.g4SimHits.Physics.EMPhysics    = False
0139 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0140     CheckSecondary = cms.PSet(
0141         SaveInFile = cms.untracked.string('PbWO4QGSP_BERT_pro_1.0GeV.root'),
0142         Verbosity = cms.untracked.int32(0),
0143         MinimumDeltaE = cms.untracked.double(0.0),
0144         KillAfter = cms.untracked.int32(1)
0145     ),
0146     type = cms.string('CheckSecondary')
0147 ),
0148     cms.PSet(
0149         type = cms.string('KillSecondariesRunAction')
0150     ))
0151