Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2023-03-17 11:04:56

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("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