File indexing completed on 2023-03-17 11:04:56
0001
0002
0003
0004
0005 import FWCore.ParameterSet.Config as cms
0006
0007 process = cms.Process("Gen")
0008
0009
0010
0011
0012
0013
0014
0015
0016
0017
0018 process.load("FWCore.MessageService.MessageLogger_cfi")
0019
0020 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0021
0022
0023 process.load("Configuration.EventContent.EventContent_cff")
0024
0025 process.maxEvents = cms.untracked.PSet(
0026 input = cms.untracked.int32(5)
0027 )
0028
0029 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0030 generator = cms.PSet(
0031 initialSeed = cms.untracked.uint32(123456789),
0032 engineName = cms.untracked.string('HepJamesRandom')
0033 )
0034 )
0035
0036
0037 process.source = cms.Source("EmptySource")
0038
0039 process.generator = cms.EDProducer("Pythia6PtGun",
0040 maxEventsToPrint = cms.untracked.int32(5),
0041 pythiaPylistVerbosity = cms.untracked.int32(1),
0042 pythiaHepMCVerbosity = cms.untracked.bool(True),
0043 PGunParameters = cms.PSet(
0044 ParticleID = cms.vint32(-15),
0045 AddAntiParticle = cms.bool(False),
0046 MinPhi = cms.double(-3.14159265359),
0047 MaxPhi = cms.double(3.14159265359),
0048 MinPt = cms.double(50.0),
0049 MaxPt = cms.double(50.0001),
0050 MinEta = cms.double(-2.4),
0051 MaxEta = cms.double(2.4)
0052 ),
0053 PythiaParameters = cms.PSet(
0054 pythiaTauJets = cms.vstring(
0055 'MDME(89,1)=0 ! no tau->electron',
0056 'MDME(90,1)=0 ! no tau->muon'
0057 ),
0058 pythiaUESettings = cms.vstring(
0059 'MSTJ(11)=3 ! Choice of the fragmentation function',
0060 'MSTJ(22)=2 ! Decay those unstable particles',
0061 'PARJ(71)=10 . ! for which ctau 10 mm',
0062 'MSTP(2)=1 ! which order running alphaS',
0063 'MSTP(33)=0 ! no K factors in hard cross sections',
0064 'MSTP(51)=7 ! structure function chosen',
0065 'MSTP(81)=1 ! multiple parton interactions 1 is Pythia default',
0066 'MSTP(82)=4 ! Defines the multi-parton model',
0067 'MSTU(21)=1 ! Check on possible errors during program execution',
0068 'PARP(82)=1.9409 ! pt cutoff for multiparton interactions',
0069 'PARP(89)=1960. ! sqrts for which PARP82 is set',
0070 'PARP(83)=0.5 ! Multiple interactions: matter distrbn parameter',
0071 'PARP(84)=0.4 ! Multiple interactions: matter distribution parameter',
0072 'PARP(90)=0.16 ! Multiple interactions: rescaling power',
0073 'PARP(67)=2.5 ! amount of initial-state radiation',
0074 'PARP(85)=1.0 ! gluon prod. mechanism in MI',
0075 'PARP(86)=1.0 ! gluon prod. mechanism in MI',
0076 'PARP(62)=1.25 ! ',
0077 'PARP(64)=0.2 ! ',
0078 'MSTP(91)=1 !',
0079 'PARP(91)=2.1 ! kt distribution',
0080 'PARP(93)=15.0 ! '
0081 ),
0082 parameterSets = cms.vstring(
0083 'pythiaUESettings',
0084 'pythiaTauJets'
0085 )
0086 )
0087 )
0088
0089 process.FEVT = cms.OutputModule("PoolOutputModule",
0090 process.FEVTSIMEventContent,
0091 fileName = cms.untracked.string('gen_singleTau.root')
0092 )
0093
0094 process.p = cms.Path(process.generator)
0095 process.outpath = cms.EndPath(process.FEVT)
0096 process.schedule = cms.Schedule(process.p,process.outpath)
0097
0098