Geometry/test/runP_Calo_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.source = cms.Source("PoolSource",
0024     fileNames = cms.untracked.vstring('file:single_neutrino_random.root')
0025 )
0026
0027 process.load('FWCore.MessageService.MessageLogger_cfi')
0028 if 'MessageLogger' in process.__dict__:
0029     process.MessageLogger.MaterialBudget=dict()
0030
0031 process.maxEvents = cms.untracked.PSet(
0032     input = cms.untracked.int32(-1)
0033 )
0034
0035 process.TFileService = cms.Service("TFileService",
0036     fileName = cms.string('matbdg_Calo.root')
0037 )
0038
0039 process.p1 = cms.Path(process.g4SimHits)
0040 process.g4SimHits.UseMagneticField = False
0041 process.g4SimHits.StackingAction.TrackNeutrino = True
0042 process.g4SimHits.Physics.type = 'SimG4Core/Physics/DummyPhysics'
0043 process.g4SimHits.Physics.DummyEMPhysics = True
0044 process.g4SimHits.Physics.CutsPerRegion = False
0045 process.g4SimHits.Generator.ApplyEtaCuts = False
0046 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0047     MaterialBudget = cms.PSet(
0048         DetectorTypes = cms.vstring('BeamPipe','Tracker','EM Calorimeter','Hadron Calorimeter','Muon System','Forward Shield'),
0049         Constituents  = cms.vint32(4,2,1,2,1,1),
0050         StackOrder    = cms.vint32(1,2,3,4,5,6),
0051         DetectorNames = cms.vstring('BEAM','BEAM1','BEAM2','BEAM3','Tracker','PLT','ECAL','CALO','VCAL','MUON','OQUA'),
0052         DetectorLevels= cms.vint32(3,3,3,3,3,3,4,3,3,3,3),
0053         EtaBoundaries = cms.vdouble(1.108,2.643,2.780,4.350,4.570,100.0),
0054         RegionTypes   = cms.vint32(0,1,0,1,0,1),
0055         Boundaries    = cms.vdouble(8050.,10860.,1595.,12800.,330.,16006.5),
0056         NBinEta      = cms.untracked.int32(250),
0057         NBinPhi      = cms.untracked.int32(180),
0058         MinEta       = cms.untracked.double(-5.0),
0059         MaxEta       = cms.untracked.double(5.0)
0060     ),
0061     type = cms.string('MaterialBudget')
0062 ))