8 process.load(
'FWCore.MessageService.MessageLogger_cfi')
9 process.load(
'Configuration.StandardSequences.Services_cff')
10 process.load(
'SimGeneral.HepPDTESSource.pdt_cfi')
11 process.load(
'Configuration.EventContent.EventContent_cff')
12 process.load(
'Geometry.HcalCommonData.hcalParameters_cfi')
13 process.load(
'Geometry.HcalCommonData.hcalDDDSimConstants_cfi')
14 process.load(
'IOMC.EventVertexGenerators.VtxSmearedFlat_cfi')
15 process.load(
'GeneratorInterface.Core.generatorSmeared_cfi')
16 process.load(
'SimG4Core.Application.g4SimHits_cfi')
17 process.load(
'IOMC.RandomEngine.IOMC_cff')
19 process.RandomNumberGeneratorService.generator.initialSeed = 456789
20 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
21 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
23 process.source = cms.Source(
"EmptySource",
24 firstRun = cms.untracked.uint32(1),
25 firstEvent = cms.untracked.uint32(1)
28 process.common_beam_parameters = cms.PSet(
29 MinE = cms.double(50.0),
30 MaxE = cms.double(50.0),
31 PartID = cms.vint32(-211),
32 MinEta = cms.double(0.2175),
33 MaxEta = cms.double(0.2175),
34 MinPhi = cms.double(-0.1309),
35 MaxPhi = cms.double(-0.1309),
36 BeamPosition = cms.double(-800.0)
39 process.generator = cms.EDProducer(
"FlatRandomEGunProducer",
40 PGunParameters = cms.PSet(
41 process.common_beam_parameters
43 Verbosity = cms.untracked.int32(0),
44 AddAntiParticle = cms.bool(
False)
47 process.VtxSmeared = cms.EDProducer(
"BeamProfileVtxGenerator",
48 process.common_beam_parameters,
50 BeamMeanX = cms.double(0.0),
51 BeamMeanY = cms.double(0.0),
52 BeamSigmaX = cms.double(0.0001),
53 BeamSigmaY = cms.double(0.0001),
54 Psi = cms.double(999.9),
55 GaussianProfile = cms.bool(
False),
58 File = cms.string(
'beam.profile'),
59 UseFile = cms.bool(
False),
60 TimeOffset = cms.double(0.)
63 process.testbeam = cms.EDAnalyzer(
"HcalTB06Analysis",
64 process.common_beam_parameters,
65 ECAL = cms.bool(
True),
66 TestBeamAnalysis = cms.PSet(
67 Verbose = cms.untracked.bool(
False),
68 ETtotMax = cms.untracked.double(400.),
69 EHCalMax = cms.untracked.double(400.),
70 EcalWidth = cms.double(0.1),
71 EcalFactor = cms.double(1.),
72 HcalFactor = cms.double(100.)
76 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.testbeam)
78 process.common_maximum_time.MaxTrackTime = cms.double(1000.0)
79 process.common_maximum_time.MaxTimeNames = cms.vstring()
80 process.common_maximum_time.MaxTrackTimes = cms.vstring()
81 process.common_maximum_time.DeadRegions = cms.vstring()
83 process.g4SimHits.NonBeamEvent =
True
84 process.g4SimHits.UseMagneticField =
False
86 process.g4SimHits.Physics.type =
'SimG4Core/Physics/QGSP_FTFP_BERT_EML'
87 process.g4SimHits.Physics.Region =
'HcalRegion'
88 process.g4SimHits.Physics.MaxTrackTime = cms.double(1000.0)
90 process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(
False)
92 process.g4SimHits.StackingAction.MaxTrackTime = cms.double(1000.0)
93 process.g4SimHits.StackingAction.MaxTimeNames = cms.vstring()
94 process.g4SimHits.StackingAction.MaxTrackTimes = cms.vdouble()
95 process.g4SimHits.StackingAction.DeadRegions = cms.vstring()
96 process.g4SimHits.StackingAction.RusRoGammaEnergyLimit = cms.double(0.0)
97 process.g4SimHits.StackingAction.RusRoNeutronEnergyLimit = cms.double(0.0)
99 process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(1000.0)
100 process.g4SimHits.SteppingAction.MaxTimeNames = cms.vstring()
101 process.g4SimHits.SteppingAction.MaxTrackTimes = cms.vdouble()
102 process.g4SimHits.SteppingAction.DeadRegions = cms.vstring()
104 process.g4SimHits.CaloSD.EminHits = cms.vdouble(0.0,0.0,0.0,0.0)
105 process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(1000.0,1000.0,1000.0,1000.0)
107 process.g4SimHits.HCalSD.UseShowerLibrary =
False
108 process.g4SimHits.HCalSD.UseHF =
False
109 process.g4SimHits.HCalSD.ForTBH2 =
False
return(e1-e2)*(e1-e2)+dp *dp