Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-04-06 12:13:55

0001 # The following comments couldn't be translated into the new config version:
0002 
0003 # Single tau(+decays) ptgun
0004 
0005 import FWCore.ParameterSet.Config as cms
0006 
0007 process = cms.Process("Gen")
0008 # this example configuration offers some minimum
0009 # annotation, to help users get through; please
0010 # don't hesitate to read through the comments
0011 # use MessageLogger to redirect/suppress multiple
0012 # service messages coming from the system
0013 #
0014 # in this config below, we use the replace option to make
0015 # the logger let out messages of severity ERROR (INFO level
0016 # will be suppressed), and we want to limit the number to 10
0017 #
0018 process.load("FWCore.MessageService.MessageLogger_cfi")
0019 
0020 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0021 
0022 # Event output
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("Pythia6PartonPtGun",
0040     maxEventsToPrint = cms.untracked.int32(5),
0041     pythiaPylistVerbosity = cms.untracked.int32(1),
0042     pythiaHepMCVerbosity = cms.untracked.bool(True),    
0043     PGunParameters = cms.PSet(
0044         PartonID = cms.int32(21),
0045         MinPhi = cms.double(-3.14159265359),
0046         MaxPhi = cms.double(3.14159265359),
0047         MinPt = cms.double(50.0),
0048         MaxPt = cms.double(50.0001),
0049         MinEta = cms.double(-1.0),
0050         MaxEta = cms.double(1.0)
0051     ),
0052     PythiaParameters = cms.PSet(
0053         pythiaJets = cms.vstring(),
0054         pythiaUESettings = cms.vstring(
0055             'MSTJ(11)=3     ! Choice of the fragmentation function',
0056             'MSTJ(22)=2     ! Decay those unstable particles',
0057             'PARJ(71)=10 .  ! for which ctau  10 mm',
0058             'MSTP(2)=1      ! which order running alphaS',
0059             'MSTP(33)=0     ! no K factors in hard cross sections',
0060             'MSTP(51)=7     ! structure function chosen',
0061             'MSTP(81)=1     ! multiple parton interactions 1 is Pythia default',
0062             'MSTP(82)=4     ! Defines the multi-parton model',
0063             'MSTU(21)=1     ! Check on possible errors during program execution',
0064             'PARP(82)=1.9409   ! pt cutoff for multiparton interactions',
0065             'PARP(89)=1960. ! sqrts for which PARP82 is set',
0066             'PARP(83)=0.5   ! Multiple interactions: matter distrbn parameter',
0067             'PARP(84)=0.4   ! Multiple interactions: matter distribution parameter',
0068             'PARP(90)=0.16  ! Multiple interactions: rescaling power',
0069             'PARP(67)=2.5    ! amount of initial-state radiation',
0070             'PARP(85)=1.0  ! gluon prod. mechanism in MI',
0071             'PARP(86)=1.0  ! gluon prod. mechanism in MI',
0072             'PARP(62)=1.25   ! ',
0073             'PARP(64)=0.2    ! ',
0074             'MSTP(91)=1     !',
0075             'PARP(91)=2.1   ! kt distribution',
0076             'PARP(93)=15.0  ! '
0077         ),
0078         parameterSets = cms.vstring(
0079             'pythiaUESettings',
0080             'pythiaJets'
0081         )
0082     )
0083 )
0084 
0085 process.FEVT = cms.OutputModule("PoolOutputModule",
0086     process.FEVTSIMEventContent,
0087     fileName = cms.untracked.string('gen_TopGun.root')
0088 )
0089 
0090 process.p = cms.Path(process.generator)
0091 process.outpath = cms.EndPath(process.FEVT)
0092 process.schedule = cms.Schedule(process.p,process.outpath)
0093 
0094