Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:07:25

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(23),
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         pythiaDefault = cms.vstring(
0055        'PMAS(5,1)=4.8 ! b quark mass', 
0056            'PMAS(6,1)=172.3 ! t quark mass',
0057            'MSTP(61)=0 ! initial state radiation',
0058        'mstj(41)=1' # per Steve M., instead of mstp(71)... btw, shoult it be 0 or 1 ?
0059            #'MSTP(71)=0 !final state radiation'
0060     ),
0061     pythiaZtoMuons = cms.vstring(
0062        "MDME(174,1)=0",          # !Z decay into d dbar,
0063            "MDME(175,1)=0",          # !Z decay into u ubar,
0064            "MDME(176,1)=0",          # !Z decay into s sbar,
0065            "MDME(177,1)=0",          # !Z decay into c cbar,
0066            "MDME(178,1)=0",          # !Z decay into b bbar,
0067            "MDME(179,1)=0",          # !Z decay into t tbar,
0068            "MDME(182,1)=0",          # !Z decay into e- e+,
0069            "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0070            "MDME(184,1)=1",          # !Z decay into mu- mu+,
0071            "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0072            "MDME(186,1)=0",          # !Z decay into tau- tau+,
0073            "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0074         ),
0075         parameterSets = cms.vstring(
0076             'pythiaDefault',
0077             'pythiaZtoMuons'
0078         )
0079     )
0080 )
0081 
0082 process.FEVT = cms.OutputModule("PoolOutputModule",
0083     process.FEVTSIMEventContent,
0084     fileName = cms.untracked.string('gen_singleTau.root')
0085 )
0086 
0087 process.p = cms.Path(process.generator)
0088 process.outpath = cms.EndPath(process.FEVT)
0089 process.schedule = cms.Schedule(process.p,process.outpath)
0090 
0091