Geometry/test/runP_Forward_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PROD")
0004
0005 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0006
0007 #Geometry
0008 #
0009 process.load("Configuration.Geometry.GeometryExtendedReco_cff")
0010
0011 #Magnetic Field
0012 #
0013 process.load("Configuration.StandardSequences.MagneticField_38T_cff")
0014
0015 # Detector simulation (Geant4-based)
0016 #
0017 process.load("SimG4Core.Application.g4SimHits_cfi")
0018
0019 process.load("IOMC.RandomEngine.IOMC_cff")
0020 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0021 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0022
0023 process.load('FWCore.MessageService.MessageLogger_cfi')
0024 if 'MessageLogger' in process.__dict__:
0025     process.MessageLogger.MaterialBudget=dict()
0026
0027 process.source = cms.Source("EmptySource",
0028     firstRun        = cms.untracked.uint32(1),
0029     firstEvent      = cms.untracked.uint32(1)
0030 )
0031
0032 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0033     PGunParameters = cms.PSet(
0034         PartID = cms.vint32(14),
0035         MinEta = cms.double(-7.0),
0036         MaxEta = cms.double(-4.5),
0037         MinPhi = cms.double(-3.14159265359),
0038         MaxPhi = cms.double(3.14159265359),
0039         MinE   = cms.double(10.0),
0040         MaxE   = cms.double(10.0)
0041     ),
0042     AddAntiParticle = cms.bool(False),
0043     Verbosity       = cms.untracked.int32(0)
0044 )
0045
0046 process.maxEvents = cms.untracked.PSet(
0047     input = cms.untracked.int32(100000)
0048 )
0049
0050 process.TFileService = cms.Service("TFileService",
0051     fileName = cms.string('matbdg_Fwd.root')
0052 )
0053
0054 process.p1 = cms.Path(process.generator*process.g4SimHits)
0055 process.g4SimHits.UseMagneticField = False
0056 process.g4SimHits.StackingAction.TrackNeutrino = True
0057 process.g4SimHits.Physics.type = 'SimG4Core/Physics/DummyPhysics'
0058 process.g4SimHits.Physics.DummyEMPhysics = True
0059 process.g4SimHits.Physics.CutsPerRegion = False
0060 process.g4SimHits.Generator.ApplyEtaCuts = False
0061
0062 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0063     MaterialBudgetForward = cms.PSet(
0064         DetectorTypes = cms.vstring('BeamPipe','Tracker','EM Calorimeter','Hadron Calorimeter','Forward Calorimeter','TOTEM','CASTOR','Forward Shield','Muon System'),
0065         Constituents  = cms.vint32(4,2,1,1,1,2,2,1,1),
0066         StackOrder    = cms.vint32(1,2,3,4,6,8,9,7,5),
0067         DetectorNames = cms.vstring('BEAM','BEAM1','BEAM2','BEAM3','Tracker','PLT','ECAL','CALO','VCAL','TotemT1','TotemT2','CastorF','CastorB','OQUA','MUON'),
0068         DetectorLevels= cms.vint32(3,3,3,3,3,3,4,3,3,3,3,3,3,3,3),
0069         EtaBoundaries = cms.vdouble(1.108,2.643,2.780,4.350,4.570,100.0),
0070         RegionTypes   = cms.vint32(0,1,0,1,0,1),
0071         Boundaries    = cms.vdouble(8050.,10860.,1595.,12800.,330.,16006.5),
0072         NBinEta      = cms.untracked.int32(250),
0073         NBinPhi      = cms.untracked.int32(180),
0074         MinEta       = cms.untracked.double(-7.0),
0075         MaxEta       = cms.untracked.double(-4.5)
0076     ),
0077     type = cms.string('MaterialBudgetForward')
0078 ))
0079
0080