Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:32:19

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Eras.Modifier_h2tb_cff import h2tb
0003 
0004 process = cms.Process("PROD",h2tb)
0005 
0006 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0007 
0008 #Geometry
0009 #
0010 process.load("SimG4CMS.HcalTestBeam.TB2007GeometryXML_cfi")
0011 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0012 process.load('Geometry.HcalTestBeamData.hcalDDDSimConstants_cff')
0013 
0014 # Detector simulation (Geant4-based)
0015 #
0016 process.load("SimG4Core.Application.g4SimHits_cfi")
0017 
0018 process.load("IOMC.RandomEngine.IOMC_cff")
0019 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0020 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0021 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0022 
0023 process.load('FWCore.MessageService.MessageLogger_cfi')
0024 if 'MessageLogger' in process.__dict__:
0025     process.MessageLogger.MaterialBudget=dict()
0026 
0027 process.common_beam_direction_parameters = cms.PSet(
0028     MinEta       = cms.double(1.562),
0029     MaxEta       = cms.double(1.562),
0030     MinPhi       = cms.double(0.0436),
0031     MaxPhi       = cms.double(0.0436),
0032     BeamPosition = cms.double(-800.0)
0033 )
0034 
0035 process.source = cms.Source("EmptySource",
0036     firstRun        = cms.untracked.uint32(1),
0037     firstEvent      = cms.untracked.uint32(1)
0038 )
0039 
0040 process.generator = cms.EDProducer("FlatRandomEGunProducer",
0041     PGunParameters = cms.PSet(
0042         process.common_beam_direction_parameters,
0043         PartID = cms.vint32(14),
0044         MinE   = cms.double(10.0),
0045         MaxE   = cms.double(10.0)
0046     ),
0047     AddAntiParticle = cms.bool(False),
0048     Verbosity       = cms.untracked.int32(0)
0049 )
0050 
0051 from IOMC.EventVertexGenerators.VtxSmearedParameters_cfi import *
0052 process.VtxSmeared = cms.EDProducer("BeamProfileVtxGenerator",
0053     process.common_beam_direction_parameters,
0054     VtxSmearedCommon,
0055     BeamMeanX       = cms.double(0.0),
0056     BeamMeanY       = cms.double(0.0),
0057     BeamSigmaX      = cms.double(0.0001),
0058     BeamSigmaY      = cms.double(0.0001),
0059     Psi             = cms.double(999.9),
0060     GaussianProfile = cms.bool(False),
0061     BinX            = cms.int32(50),
0062     BinY            = cms.int32(50),
0063     File            = cms.string('beam.profile'),
0064     UseFile         = cms.bool(False),
0065     TimeOffset      = cms.double(0.)
0066 )
0067 
0068 process.maxEvents = cms.untracked.PSet(
0069     input = cms.untracked.int32(1)
0070 )
0071 
0072 process.TFileService = cms.Service("TFileService",
0073     fileName = cms.string('matbdg_HCAL1.root')
0074 )
0075 
0076 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.g4SimHits)
0077 process.g4SimHits.NonBeamEvent = True
0078 process.g4SimHits.UseMagneticField = False
0079 process.g4SimHits.StackingAction.TrackNeutrino = True
0080 process.g4SimHits.Physics.type = 'SimG4Core/Physics/DummyPhysics'
0081 process.g4SimHits.Physics.DummyEMPhysics = True
0082 process.g4SimHits.Physics.CutsPerRegion = False
0083 process.g4SimHits.ECalSD.TestBeam = True
0084 process.g4SimHits.HCalSD.UseShowerLibrary = False
0085 process.g4SimHits.HCalSD.TestNumberingScheme = True
0086 process.g4SimHits.HCalSD.UseHF = False
0087 process.g4SimHits.HCalSD.ForTBH2 = True
0088 process.g4SimHits.CaloTrkProcessing.TestBeam = True
0089 process.g4SimHits.Watchers = cms.VPSet(cms.PSet(
0090     MaterialBudgetHcal = cms.PSet(
0091         FillHisto    = cms.untracked.bool(False),
0092         PrintSummary = cms.untracked.bool(True),
0093         NBinPhi      = cms.untracked.int32(180),
0094         NBinEta      = cms.untracked.int32(100),
0095         MaxEta       = cms.untracked.double(3.0),
0096         EtaLow       = cms.untracked.double(-3.0),
0097         EtaHigh      = cms.untracked.double(3.0),
0098         RMax         = cms.untracked.double(10.0),
0099         ZMax         = cms.untracked.double(15.0)
0100     ),
0101     type = cms.string('MaterialBudgetHcal')
0102 ))
0103 
0104