12 myParList = cms.vdouble(1.006, 1.0, 0.0, 1.82790e+00, 3.66237e+00, 0.965, 1.0)
13 myNameTag = mySample +
"_" + myEnergy +
"_" + mySection
15 process = cms.Process(
"TBSim")
20 process.load(
"FWCore.MessageLogger.MessageLogger_cfi")
23 process.load(
"SimGeneral.HepPDTESSource.pythiapdt_cfi")
25 process.RandomNumberGeneratorService = cms.Service(
"RandomNumberGeneratorService",
26 saveFileName = cms.untracked.string(
''),
28 initialSeed = cms.untracked.uint32(7824367),
29 engineName = cms.untracked.string(
'HepJamesRandom')
32 initialSeed = cms.untracked.uint32(int(mySeed)),
33 engineName = cms.untracked.string(
'HepJamesRandom')
35 VtxSmeared = cms.PSet(
36 initialSeed = cms.untracked.uint32(98765432),
37 engineName = cms.untracked.string(
'HepJamesRandom')
40 initialSeed = cms.untracked.uint32(11),
41 engineName = cms.untracked.string(
'HepJamesRandom')
43 SimEcalTBG4Object = cms.PSet(
44 initialSeed = cms.untracked.uint32(12),
45 engineName = cms.untracked.string(
'HepJamesRandom')
48 initialSeed = cms.untracked.uint32(12345),
49 engineName = cms.untracked.string(
'HepJamesRandom')
51 simEcalUnsuppressedDigis = cms.PSet(
52 initialSeed = cms.untracked.uint32(1234567),
53 engineName = cms.untracked.string(
'HepJamesRandom')
57 process.randomEngineStateProducer = cms.EDProducer(
"RandomEngineStateProducer")
60 process.maxEvents = cms.untracked.PSet(
61 input = cms.untracked.int32(int(myEvent)*10)
67 process.common_beam_direction_parameters = cms.PSet(
68 BeamMeanY = cms.double(0.0),
69 BeamMeanX = cms.double(0.0),
70 MinEta = cms.double(0.221605),
71 MaxEta = cms.double(0.221605),
72 MinPhi = cms.double(0.0467487),
73 MaxPhi = cms.double(0.0467487),
75 BeamPosition = cms.double(-26733.5)
81 process.MuonNumberingInitialization = cms.ESProducer(
"MuonNumberingInitialization")
82 process.XMLIdealGeometryESSource = cms.ESSource(
"XMLIdealGeometryESSource",
83 geomXMLFiles = cms.vstring(
'Geometry/CMSCommonData/data/materials.xml',
84 'Geometry/CMSCommonData/data/rotations.xml',
85 'Geometry/EcalTestBeam/data/ebcon.xml',
86 'Geometry/EcalCommonData/data/ebrot.xml',
87 'Geometry/EcalTestBeam/data/eregalgo.xml',
88 'Geometry/EcalCommonData/data/ebalgo.xml',
89 'Geometry/EcalTestBeam/data/tbrot.xml',
90 'Geometry/EcalTestBeam/data/TBH4.xml',
91 'Geometry/EcalTestBeam/data/TBH4ecalsens.xml',
92 'Geometry/HcalSimData/data/CaloUtil.xml',
95 'SimG4Core/GFlash/TB/gflashTBH4ProdCuts.xml',
96 'Geometry/CMSCommonData/data/FieldParameters.xml'),
97 rootNodeName = cms.string(
'TBH4:OCMS')
101 process.load(
"Geometry.CaloEventSetup.CaloGeometry_cff")
102 process.load(
"Geometry.CaloEventSetup.CaloTopology_cfi")
103 process.load(
"Geometry.EcalMapping.EcalMapping_cfi")
104 process.load(
"Geometry.CaloEventSetup.EcalTrigTowerConstituents_cfi")
105 process.CaloGeometryBuilder.SelectedCalos = [
'EcalBarrel']
111 process.load(
"MagneticField.Engine.uniformMagneticField_cfi")
113 process.source = cms.Source(
"EmptySource")
115 process.generator = cms.EDProducer(
"FlatRandomEGunProducer",
116 PGunParameters = cms.PSet(
117 process.common_beam_direction_parameters,
118 PartID = cms.vint32(11),
119 MinE = cms.double(float(myEnergy)),
120 MaxE = cms.double(float(myEnergy))
122 Verbosity = cms.untracked.int32(0),
124 psethack = cms.string(
'single electron'),
125 AddAntiParticle = cms.bool(
False),
128 process.ProductionFilterSequence = cms.Sequence(process.generator)
136 process.VtxSmeared = cms.EDProducer(
"BeamProfileVtxGenerator",
137 process.common_beam_direction_parameters,
142 BeamSigmaX = cms.double(5.0),
143 BeamSigmaY = cms.double(5.0),
144 Psi = cms.double(999.9),
145 GaussianProfile = cms.bool(
True),
146 BinX = cms.int32(50),
147 BinY = cms.int32(50),
148 File = cms.string(
'beam.profile'),
149 UseFile = cms.bool(
False),
150 TimeOffset = cms.double(0.)
155 process.load(
"SimG4Core.Application.g4SimHits_cfi")
157 process.g4SimHits.UseMagneticField = cms.bool(
False)
158 process.g4SimHits.Physics.DefaultCutValue = 1.
159 process.g4SimHits.NonBeamEvent = cms.bool(
True)
160 process.g4SimHits.Generator.HepMCProductLabel = cms.string(
'generator')
161 process.g4SimHits.Generator.ApplyPCuts = cms.bool(
False)
162 process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(
True)
163 process.g4SimHits.Generator.ApplyPhiCuts = cms.bool(
False)
164 process.g4SimHits.Generator.MaxEtaCut = cms.double(1.5)
165 process.g4SimHits.Generator.MinEtaCut = cms.double(0.0)
166 process.g4SimHits.CaloSD.CorrectTOFBeam = cms.bool(
True)
167 process.g4SimHits.CaloSD.BeamPosition = cms.double(-26733.5)
168 process.g4SimHits.CaloTrkProcessing.TestBeam = cms.bool(
True)
169 process.g4SimHits.StackingAction.MaxTrackTime = cms.double(10000.)
170 process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(10000.)
171 process.g4SimHits.CaloSD.TmaxHit = cms.double(10000.)
172 process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(10000.,10000.,10000.,10000.,10000.)
173 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
174 type = cms.string(
'EcalTBH4Trigger'),
175 verbose = cms.untracked.bool(
False),
177 trigEvents = cms.untracked.int32(int(myEvent))
183 process.SimEcalTBG4Object = cms.EDProducer(
"EcalTBMCInfoProducer",
184 process.common_beam_direction_parameters,
185 CrystalMapFile = cms.FileInPath(
'Geometry/EcalTestBeam/data/BarrelSM1CrystalCenterElectron120GeV.dat'),
186 moduleLabelVtx = cms.untracked.string(
'generator')
191 process.SimEcalTBHodoscope = cms.EDProducer(
"TBHodoActiveVolumeRawInfoProducer")
195 process.SimEcalEventHeader = cms.EDProducer(
"FakeTBEventHeaderProducer",
196 EcalTBInfoLabel = cms.untracked.string(
'SimEcalTBG4Object')
203 process.load(
"SimGeneral.MixingModule.mixNoPU_cfi")
207 process.load(
"CalibCalorimetry.EcalTrivialCondModules.EcalTrivialCondRetrieverTB_cfi")
209 process.load(
"CalibCalorimetry.EcalTrivialCondModules.ESTrivialCondRetriever_cfi")
213 process.load(
"SimCalorimetry.EcalTestBeam.ecaldigi_testbeam_cfi")
216 process.load(
"Configuration.EcalTB.localReco_tbsim_cff")
235 process.ecal_notCont_sim.EBs25notContainment = 1.0
236 process.ecal_notCont_sim.EEs25notContainment = 1.0
237 process.g4SimHits.Physics.type =
'SimG4Core/Physics/GFlash'
238 process.g4SimHits.Physics.GFlash = cms.PSet(
239 bField = cms.double(0.0),
240 GflashEMShowerModel = cms.bool(
True),
241 GflashHadronShowerModel = cms.bool(
True),
242 GflashHistogram = cms.bool(
True),
243 GflashHistogramName = cms.string(
"gflash_histogram_"+myNameTag+
".root"),
244 GflashHadronPhysics = cms.string(
'QGSP_BERT'),
245 GflashHcalOuter = cms.bool(
True),
246 GflashExportToFastSim = cms.bool(
False),
247 watcherOn = cms.bool(
False),
248 Verbosity = cms.untracked.int32(0),
249 tuning_pList = myParList
252 print "physics type : ", process.g4SimHits.Physics.type
255 process.doSimHits = cms.Sequence(process.ProductionFilterSequence*process.VtxSmeared*process.g4SimHits)
256 process.doSimTB = cms.Sequence(process.SimEcalTBG4Object*process.SimEcalTBHodoscope*process.SimEcalEventHeader)
257 process.doEcalDigis = cms.Sequence(process.mix*process.simEcalUnsuppressedDigis)
259 process.p1 = cms.Path(process.doSimHits*process.doSimTB*process.doEcalDigis*process.localReco_tbsim)
265 process.MessageLogger.destinations=cms.untracked.vstring(
'cout'
269 process.MessageLogger.categories=cms.untracked.vstring(
'FwkJob'
280 ,
'BeamProfileVtxGenerator'
286 process.MessageLogger.G4msg = cms.untracked.PSet(
287 noTimeStamps = cms.untracked.bool(
True)
289 ,threshold = cms.untracked.string(
'INFO')
291 ,INFO = cms.untracked.PSet(limit = cms.untracked.int32(0))
292 ,FwkReport = cms.untracked.PSet(limit = cms.untracked.int32(0))
293 ,FwkSummary = cms.untracked.PSet(limit = cms.untracked.int32(0))
294 ,Root_NoDictionary = cms.untracked.PSet(limit = cms.untracked.int32(0))
295 ,FwkJob = cms.untracked.PSet(limit = cms.untracked.int32(0))
296 ,TimeReport = cms.untracked.PSet(limit = cms.untracked.int32(0))
297 ,TimeModule = cms.untracked.PSet(limit = cms.untracked.int32(0))
298 ,TimeEvent = cms.untracked.PSet(limit = cms.untracked.int32(0))
299 ,MemoryCheck = cms.untracked.PSet(limit = cms.untracked.int32(0))
305 ,PhysicsList = cms.untracked.PSet(limit = cms.untracked.int32(-1))
306 ,G4cout = cms.untracked.PSet(limit = cms.untracked.int32(-1))
307 ,G4cerr = cms.untracked.PSet(limit = cms.untracked.int32(-1))
308 ,BeamProfileVtxGenerator = cms.untracked.PSet(limit = cms.untracked.int32(-1))
313 process.options = cms.untracked.PSet(
314 wantSummary = cms.untracked.bool(
True)
319 process.options = cms.untracked.PSet()
320 process.options.SkipEvent = cms.untracked.vstring(
'EventCorruption')