File indexing completed on 2024-04-06 12:29:56
0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("PROD")
0004
0005 process.load("SimG4CMS.Calo.pythiapdt_cfi")
0006 process.load('Configuration.StandardSequences.Generator_cff')
0007 process.load("IOMC.EventVertexGenerators.VtxSmearedGauss_cfi")
0008 process.load("Geometry.CMSCommonData.cmsSimIdealGeometryXML_cfi")
0009 process.load("Geometry.TrackerNumberingBuilder.trackerNumberingGeometry_cfi")
0010 process.load("Geometry.EcalCommonData.ecalSimulationParameters_cff")
0011 process.load("Geometry.HcalCommonData.hcalDDDSimConstants_cff")
0012 process.load("Configuration.StandardSequences.MagneticField_cff")
0013
0014 process.load("Configuration.EventContent.EventContent_cff")
0015
0016 process.load("Configuration.StandardSequences.FrontierConditions_GlobalTag_cff")
0017 from Configuration.AlCa.autoCond import autoCond
0018 process.GlobalTag.globaltag = autoCond['run1_mc']
0019
0020 process.load("SimG4Core.Application.g4SimHits_cfi")
0021
0022 process.MessageLogger = cms.Service("MessageLogger",
0023 cerr = cms.untracked.PSet(
0024 enable = cms.untracked.bool(False)
0025 ),
0026 cout = cms.untracked.PSet(
0027 DEBUG = cms.untracked.PSet(
0028 limit = cms.untracked.int32(0)
0029 ),
0030 G4cerr = cms.untracked.PSet(
0031 limit = cms.untracked.int32(-1)
0032 ),
0033 G4cout = cms.untracked.PSet(
0034 limit = cms.untracked.int32(-1)
0035 ),
0036 INFO = cms.untracked.PSet(
0037 limit = cms.untracked.int32(-1)
0038 ),
0039 SimG4CoreWatcher = cms.untracked.PSet(
0040 limit = cms.untracked.int32(-1)
0041 ),
0042 enable = cms.untracked.bool(True)
0043 ),
0044 debugModules = cms.untracked.vstring('*')
0045 )
0046
0047 process.load("IOMC.RandomEngine.IOMC_cff")
0048 process.RandomNumberGeneratorService.generator.initialSeed = 456789
0049 process.RandomNumberGeneratorService.g4SimHits.initialSeed = 9876
0050 process.RandomNumberGeneratorService.VtxSmeared.initialSeed = 123456789
0051
0052 process.maxEvents = cms.untracked.PSet(
0053 input = cms.untracked.int32(500)
0054 )
0055
0056 process.source = cms.Source("LHESource",
0057 fileNames = cms.untracked.vstring('file:../data/monopole.lhe')
0058 )
0059
0060 process.generator = cms.EDFilter("Pythia6HadronizerFilter",
0061 pythiaPylistVerbosity = cms.untracked.int32(1),
0062 pythiaHepMCVerbosity = cms.untracked.bool(True),
0063 comEnergy = cms.double(900.0),
0064 maxEventsToPrint = cms.untracked.int32(0),
0065 PythiaParameters = cms.PSet(
0066 pythiaUESettings = cms.vstring('MSTJ(11)=3 ! Choice of the fragmentation function',
0067 'MSTJ(22)=2 ! Decay those unstable particles',
0068 'PARJ(71)=10. ! for which ctau 10 mm',
0069 'MSTP(2)=1 ! which order running alphaS',
0070 'MSTP(33)=0 ! no K factors in hard cross sections',
0071 'MSTP(51)=10042 ! structure function chosen (external PDF CTEQ6L1)',
0072 'MSTP(52)=2 ! work with LHAPDF',
0073 'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default',
0074 'MSTP(82)=4 ! Defines the multi-parton model',
0075 'MSTU(21)=1 ! Check on possible errors during program execution',
0076 'PARP(82)=1.8387 ! pt cutoff for multiparton interactions',
0077 'PARP(89)=1960. ! sqrts for which PARP82 is set',
0078 'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter',
0079 'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter',
0080 'PARP(90)=0.16 ! Multiple interactions: rescaling power',
0081 'PARP(67)=2.5 ! amount of initial-state radiation',
0082 'PARP(85)=1.0 ! gluon prod. mechanism in MI',
0083 'PARP(86)=1.0 ! gluon prod. mechanism in MI',
0084 'PARP(62)=1.25 ! ',
0085 'PARP(64)=0.2 ! ',
0086 'MSTP(91)=1 !',
0087 'PARP(91)=2.1 ! kt distribution',
0088 'PARP(93)=15.0 ! '),
0089 processParameters = cms.vstring('MSEL=0 ! User defined processes',
0090 'PMAS(5,1)=4.8 ! b quark mass',
0091 'PMAS(6,1)=172.5 ! t quark mass',
0092 'MSTJ(1)=1 ! Fragmentation/hadronization on or off',
0093 'MSTP(61)=1 ! Parton showering on or off'),
0094 parameterSets = cms.vstring('pythiaUESettings',
0095 'processParameters')
0096 )
0097 )
0098
0099 process.o1 = cms.OutputModule("PoolOutputModule",
0100 process.FEVTSIMEventContent,
0101 fileName = cms.untracked.string('simevent.root')
0102 )
0103
0104 process.Timing = cms.Service("Timing")
0105
0106 process.Tracer = cms.Service("Tracer")
0107
0108 process.p1 = cms.Path(process.generator*process.VtxSmeared*process.g4SimHits)
0109 process.outpath = cms.EndPath(process.o1)
0110 process.g4SimHits.Physics.type = 'SimG4Core/Physics/QGSP_BERT_EML'
0111 process.g4SimHits.Physics.MonopoleCharge = 1
0112 process.g4SimHits.Physics.MonopoleTransport = False
0113 process.g4SimHits.Physics.Verbosity = 2
0114 process.g4SimHits.G4Commands = ['/tracking/verbose 1']
0115 process.common_maximum_timex = cms.PSet(
0116 MaxTrackTime = cms.double(1000.0),
0117 MaxTimeNames = cms.vstring(),
0118 MaxTrackTimes = cms.vdouble(),
0119 DeadRegions = cms.vstring(),
0120 CriticalEnergyForVacuum = cms.double(2.0),
0121 CriticalDensity = cms.double(1e-15)
0122 )
0123 process.g4SimHits.StackingAction = cms.PSet(
0124 process.common_heavy_suppression,
0125 process.common_maximum_timex,
0126 TrackNeutrino = cms.bool(False),
0127 KillDeltaRay = cms.bool(False),
0128 KillHeavy = cms.bool(False),
0129 KillGamma = cms.bool(False),
0130 GammaThreshold= cms.double(0.0001),
0131 SaveFirstLevelSecondary = cms.untracked.bool(True),
0132 SavePrimaryDecayProductsAndConversionsInTracker = cms.untracked.bool(True),
0133 SavePrimaryDecayProductsAndConversionsInCalo = cms.untracked.bool(True),
0134 SavePrimaryDecayProductsAndConversionsInMuon = cms.untracked.bool(True),
0135 RusRoGammaEnergyLimit = cms.double(5.0),
0136 RusRoEcalGamma = cms.double(0.3),
0137 RusRoHcalGamma = cms.double(0.3),
0138 RusRoMuonIronGamma = cms.double(0.3),
0139 RusRoPreShowerGamma = cms.double(0.3),
0140 RusRoCastorGamma = cms.double(0.3),
0141 RusRoWorldGamma = cms.double(0.3),
0142 RusRoNeutronEnergyLimit= cms.double(10.0),
0143 RusRoEcalNeutron = cms.double(0.1),
0144 RusRoHcalNeutron = cms.double(0.1),
0145 RusRoMuonIronNeutron = cms.double(0.1),
0146 RusRoPreShowerNeutron = cms.double(0.1),
0147 RusRoCastorNeutron = cms.double(0.1),
0148 RusRoWorldNeutron = cms.double(0.1),
0149 RusRoProtonEnergyLimit = cms.double(0.0),
0150 RusRoEcalProton = cms.double(1.0),
0151 RusRoHcalProton = cms.double(1.0),
0152 RusRoMuonIronProton = cms.double(1.0),
0153 RusRoPreShowerProton = cms.double(1.0),
0154 RusRoCastorProton = cms.double(1.0),
0155 RusRoWorldProton = cms.double(1.0)
0156 )
0157 process.g4SimHits.SteppingAction = cms.PSet(
0158 process.common_maximum_timex,
0159 EkinNames = cms.vstring(),
0160 EkinThresholds = cms.vdouble(),
0161 EkinParticles = cms.vstring(),
0162 Verbosity = cms.untracked.int32(2)
0163 )
0164 process.g4SimHits.Watchers = cms.VPSet(
0165 cms.PSet(
0166 CheckForHighEtPhotons = cms.untracked.bool(False),
0167 G4Verbose = cms.untracked.bool(True),
0168 EventMin = cms.untracked.int32(0),
0169 EventMax = cms.untracked.int32(1),
0170 EventStep = cms.untracked.int32(1),
0171 TrackMin = cms.untracked.int32(0),
0172 TrackMax = cms.untracked.int32(99999999),
0173 TrackStep = cms.untracked.int32(1),
0174 PDGids = cms.untracked.vint32(4110000,-4110000),
0175 VerboseLevel = cms.untracked.int32(2),
0176 DEBUG = cms.untracked.bool(False),
0177 type = cms.string('TrackingVerboseAction')
0178 ),
0179 cms.PSet(
0180 ChangeFromFirstStep = cms.untracked.bool(True),
0181 type = cms.string('MonopoleSteppingAction')
0182 )
0183 )
0184
