00001 import FWCore.ParameterSet.Config as cms
00002
00003 process = cms.Process("Sim")
00004 process.load("FWCore.MessageLogger.MessageLogger_cfi")
00005
00006 process.load("SimG4Core.GFlash.TB.h4TB2006GeometryXML_cfi")
00007
00008 process.load("Configuration.EcalTB.simulation_tbsim_cff")
00009
00010 process.load("Configuration.EcalTB.digitization_tbsim_cff")
00011
00012 process.load("Configuration.EcalTB.localReco_tbsim_cff")
00013
00014 process.maxEvents = cms.untracked.PSet(
00015 input = cms.untracked.int32(10)
00016 )
00017 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
00018 moduleSeeds = cms.PSet(
00019 g4SimHits = cms.untracked.uint32(9876),
00020 SimEcalTBG4Object = cms.untracked.uint32(5432),
00021 ecalUnsuppressedDigis = cms.untracked.uint32(54321),
00022 VtxSmeared = cms.untracked.uint32(12345)
00023 ),
00024 sourceSeed = cms.untracked.uint32(98765)
00025 )
00026
00027 process.common_beam_direction_parameters = cms.PSet(
00028 BeamMeanY = cms.untracked.double(0.0),
00029 BeamMeanX = cms.untracked.double(0.0),
00030 MaxEta = cms.untracked.double(0.221525),
00031 MaxPhi = cms.untracked.double(0.0467617),
00032 MinEta = cms.untracked.double(0.221525),
00033 BeamPosition = cms.untracked.double(-26733.5),
00034 MinPhi = cms.untracked.double(0.0467617)
00035 )
00036 process.source = cms.Source("FlatRandomEGunSource",
00037 PGunParameters = cms.untracked.PSet(
00038 process.common_beam_direction_parameters,
00039 MaxE = cms.untracked.double(20.0),
00040 MinE = cms.untracked.double(20.0),
00041 PartID = cms.untracked.vint32(11)
00042 ),
00043 Verbosity = cms.untracked.int32(0)
00044 )
00045
00046 process.VtxSmeared = cms.EDProducer("BeamProfileVtxGenerator",
00047 process.common_beam_direction_parameters,
00048 BeamSigmaX = cms.untracked.double(2.4),
00049 BeamSigmaY = cms.untracked.double(2.4),
00050 GaussianProfile = cms.untracked.bool(False)
00051 )
00052
00053 process.TreeProducerCalibSimul = cms.EDAnalyzer("TreeProducerCalibSimul",
00054 rootfile = cms.untracked.string('treeTB_gf.root'),
00055 eventHeaderCollection = cms.string(''),
00056 eventHeaderProducer = cms.string('SimEcalEventHeader'),
00057 txtfile = cms.untracked.string('treeTB_gf.txt'),
00058 EBRecHitCollection = cms.string('EcalRecHitsEB'),
00059 tdcRecInfoCollection = cms.string('EcalTBTDCRecInfo'),
00060 xtalInBeam = cms.untracked.int32(248),
00061 hodoRecInfoProducer = cms.string('ecalTBSimHodoscopeReconstructor'),
00062 hodoRecInfoCollection = cms.string('EcalTBHodoscopeRecInfo'),
00063 RecHitProducer = cms.string('ecalTBSimRecHit'),
00064 tdcRecInfoProducer = cms.string('ecalTBSimTDCReconstructor')
00065 )
00066
00067 process.o1 = cms.OutputModule("PoolOutputModule",
00068 outputCommands = cms.untracked.vstring('keep *',
00069 'drop PSimHits_g4SimHits_*_Sim',
00070 'keep PCaloHits_g4SimHits_EcalHitsEB_Sim',
00071 'keep PCaloHits_g4SimHits_CaloHitsTk_Sim',
00072 'keep PCaloHits_g4SimHits_EcalTBH4BeamHits_Sim'),
00073 fileName = cms.untracked.string('recoTB_gf.root')
00074 )
00075
00076 process.doSimHits = cms.Sequence(process.VtxSmeared*process.g4SimHits)
00077 process.doSimTB = cms.Sequence(process.SimEcalTBG4Object*process.SimEcalTBHodoscope*process.SimEcalEventHeader)
00078 process.doEcalDigis = cms.Sequence(process.mix*process.ecalUnsuppressedDigis)
00079 process.p1 = cms.Path(process.doSimHits*process.doSimTB*process.doEcalDigis*process.localReco_tbsim*process.TreeProducerCalibSimul)
00080 process.outpath = cms.EndPath(process.o1)
00081 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
00082 type = cms.string('EcalTBH4Trigger'),
00083 verbose = cms.untracked.bool(False),
00084 trigEvents = cms.untracked.int32(100000)
00085 ))
00086 process.g4SimHits.Physics.type = 'SimG4Core/Physics/GFlash'
00087 process.g4SimHits.Physics.GFlash = cms.PSet(
00088 bField = cms.double(0.0),
00089 GflashEMShowerModel = cms.bool(True),
00090 GflashHadronShowerModel = cms.bool(True),
00091 GflashHistogram = cms.bool(True),
00092 GflashHistogramName = cms.string('gflash_histogram_h4.root'),
00093 GflashHadronPhysics = cms.string('QGSP_BERT'),
00094 tuning_pList = cms.vdouble()
00095 )
00096 process.ecal_notCont_sim.EBs25notContainment = 1.0
00097 process.ecal_notCont_sim.EEs25notContainment = 1.0
00098
00099
1.7.3