Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 14:29:09

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),  ## (MeV)
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), ## (MeV)
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), ## (MeV)
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