4 if hasattr(process.VtxSmeared,
"X0"):
5 VertexX = process.VtxSmeared.X0
6 VertexY = process.VtxSmeared.Y0
7 VertexZ = process.VtxSmeared.Z0
9 if hasattr(process.VtxSmeared,
"MeanX"):
10 VertexX = process.VtxSmeared.MeanX
11 VertexY = process.VtxSmeared.MeanY
12 VertexZ = process.VtxSmeared.MeanZ
14 process.load(
"SimG4Core.Application.g4SimHits_cfi")
15 process.g4SimHits.Generator.HepMCProductLabel =
'LHCTransport' 16 process.g4SimHits.Generator.MinEtaCut = -13.0
17 process.g4SimHits.Generator.MaxEtaCut = 13.0
19 process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(2000.0)
20 process.g4SimHits.StackingAction.MaxTrackTime = cms.double(2000.0)
22 process.load(
'SimTransport.HectorProducer.HectorTransportCTPPS_cfi')
23 process.LHCTransport.CTPPSHector.VtxMeanX = VertexX
24 process.LHCTransport.CTPPSHector.VtxMeanY = VertexY
25 process.LHCTransport.CTPPSHector.VtxMeanZ = VertexZ
27 process.transport_step = cms.Path(process.LHCTransport)
28 process.psim = cms.Sequence(cms.SequencePlaceholder(
"randomEngineStateProducer")*process.LHCTransport*process.g4SimHits)
31
return(e1-e2)*(e1-e2)+dp *dp