7 process.load(
'FWCore.MessageService.MessageLogger_cfi')
8 process.load(
'Configuration.StandardSequences.Services_cff')
9 process.load(
'SimGeneral.HepPDTESSource.pdt_cfi')
10 process.load(
'Configuration.EventContent.EventContent_cff')
11 process.load(
'Geometry.HcalCommonData.hcalParameters_cfi')
12 process.load(
'Geometry.HcalCommonData.hcalDDDSimConstants_cfi')
13 process.load(
'IOMC.EventVertexGenerators.VtxSmearedFlat_cfi')
14 process.load(
'GeneratorInterface.Core.generatorSmeared_cfi')
15 process.load(
'SimG4Core.Application.g4SimHits_cfi')
16 process.load(
'IOMC.RandomEngine.IOMC_cff')
18 process.RandomNumberGeneratorService.generator.initialSeed = 456789
19 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
20 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
22 process.source = cms.Source(
"EmptySource",
23 firstRun = cms.untracked.uint32(1),
24 firstEvent = cms.untracked.uint32(1)
27 process.common_beam_parameters = cms.PSet(
28 MinE = cms.double(50.0),
29 MaxE = cms.double(50.0),
30 PartID = cms.vint32(-211),
31 MinEta = cms.double(0.2175),
32 MaxEta = cms.double(0.2175),
33 MinPhi = cms.double(-0.1309),
34 MaxPhi = cms.double(-0.1309),
35 BeamPosition = cms.double(-800.0)
38 process.generator = cms.EDProducer(
"FlatRandomEGunProducer",
39 PGunParameters = cms.PSet(
40 process.common_beam_parameters
42 Verbosity = cms.untracked.int32(0),
43 AddAntiParticle = cms.bool(
False)
46 process.VtxSmeared = cms.EDProducer(
"BeamProfileVtxGenerator",
47 process.common_beam_parameters,
49 BeamMeanX = cms.double(0.0),
50 BeamMeanY = cms.double(0.0),
51 BeamSigmaX = cms.double(0.0001),
52 BeamSigmaY = cms.double(0.0001),
53 Psi = cms.double(999.9),
54 GaussianProfile = cms.bool(
False),
57 File = cms.string(
'beam.profile'),
58 UseFile = cms.bool(
False),
59 TimeOffset = cms.double(0.)
62 process.testbeam = cms.EDAnalyzer(
"HcalTB06Analysis",
63 process.common_beam_parameters,
64 ECAL = cms.bool(
True),
65 TestBeamAnalysis = cms.PSet(
66 Verbose = cms.untracked.bool(
False),
67 ETtotMax = cms.untracked.double(400.),
68 EHCalMax = cms.untracked.double(400.),
69 beamEnergy = cms.untracked.double(50.),
70 TimeLimit = cms.double(180.),
71 EcalWidth = cms.double(0.362),
72 HcalWidth = cms.double(0.640),
73 EcalFactor = cms.double(1.),
74 HcalFactor = cms.double(100.),
75 MakeTree = cms.untracked.bool(
False)
79 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.testbeam)
81 process.common_maximum_time.MaxTrackTime = cms.double(1000.0)
82 process.common_maximum_time.MaxTimeNames = cms.vstring()
83 process.common_maximum_time.MaxTrackTimes = cms.vstring()
84 process.common_maximum_time.DeadRegions = cms.vstring()
86 process.g4SimHits.NonBeamEvent =
True 87 process.g4SimHits.UseMagneticField =
False 89 process.g4SimHits.Physics.type =
'SimG4Core/Physics/QGSP_FTFP_BERT_EML' 90 process.g4SimHits.Physics.Region =
'HcalRegion' 91 process.g4SimHits.Physics.MaxTrackTime = cms.double(1000.0)
93 process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(
False)
95 process.g4SimHits.StackingAction.MaxTrackTime = cms.double(1000.0)
96 process.g4SimHits.StackingAction.MaxTimeNames = cms.vstring()
97 process.g4SimHits.StackingAction.MaxTrackTimes = cms.vdouble()
98 process.g4SimHits.StackingAction.DeadRegions = cms.vstring()
100 process.g4SimHits.SteppingAction.MaxTrackTime = cms.double(1000.0)
101 process.g4SimHits.SteppingAction.MaxTimeNames = cms.vstring()
102 process.g4SimHits.SteppingAction.MaxTrackTimes = cms.vdouble()
103 process.g4SimHits.SteppingAction.DeadRegions = cms.vstring()
105 process.g4SimHits.CaloSD.EminHits = cms.vdouble(0.0,0.0,0.0,0.0)
106 process.g4SimHits.CaloSD.TmaxHits = cms.vdouble(1000.0,1000.0,1000.0,1000.0)
108 process.g4SimHits.HCalSD.UseShowerLibrary =
False 109 process.g4SimHits.HCalSD.UseHF =
False 110 process.g4SimHits.HCalSD.ForTBH2 =
False 113 def testbeam2006(process)
return(e1-e2)*(e1-e2)+dp *dp