CMS 3D CMS Logo

Functions
TB2006Analysis_cfi Namespace Reference

Functions

def testbeam2006 (process)
 

Function Documentation

def TB2006Analysis_cfi.testbeam2006 (   process)

Definition at line 5 of file TB2006Analysis_cfi.py.

References reco.return().

5 def testbeam2006(process):
6 
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')
17 
18  process.RandomNumberGeneratorService.generator.initialSeed = 456789
19  process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
20  process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
21 
22  process.source = cms.Source("EmptySource",
23  firstRun = cms.untracked.uint32(1),
24  firstEvent = cms.untracked.uint32(1)
25  )
26 
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)
36  )
37 
38  process.generator = cms.EDProducer("FlatRandomEGunProducer",
39  PGunParameters = cms.PSet(
40  process.common_beam_parameters
41  ),
42  Verbosity = cms.untracked.int32(0),
43  AddAntiParticle = cms.bool(False)
44  )
45 
46  process.VtxSmeared = cms.EDProducer("BeamProfileVtxGenerator",
47  process.common_beam_parameters,
48  VtxSmearedCommon,
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),
55  BinX = cms.int32(50),
56  BinY = cms.int32(50),
57  File = cms.string('beam.profile'),
58  UseFile = cms.bool(False),
59  TimeOffset = cms.double(0.)
60  )
61 
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)
76  )
77  )
78 
79  process.p1 = cms.Path(process.generator*process.VtxSmeared*process.generatorSmeared*process.g4SimHits*process.testbeam)
80 
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()
85 
86  process.g4SimHits.NonBeamEvent = True
87  process.g4SimHits.UseMagneticField = False
88 
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)
92 
93  process.g4SimHits.Generator.ApplyEtaCuts = cms.bool(False)
94 
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()
99 
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()
104 
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)
107 
108  process.g4SimHits.HCalSD.UseShowerLibrary = False
109  process.g4SimHits.HCalSD.UseHF = False
110  process.g4SimHits.HCalSD.ForTBH2 = False
111 
112  return(process)
113 
def testbeam2006(process)
return(e1-e2)*(e1-e2)+dp *dp