Generator/python/Pythia_H0_pyupda_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.Generator.PythiaUESettings_cfi import *
0004
0005 #
0006 # Example illustrating how to change particle properties/branching ratios
0007 # with PYUPDA cards.
0008 #
0009 # This produces an H0, which decays to a pair of long-lived exotic particles
0010 # which then each decay to a pair of light quarks.
0011 #
0012
0013 generator = cms.EDFilter("Pythia6GeneratorFilter",
0014     pythiaHepMCVerbosity = cms.untracked.bool(True),
0015     maxEventsToPrint = cms.untracked.int32(10),
0016     # Print event
0017     pythiaPylistVerbosity = cms.untracked.int32(3),
0018     # Print decay tables
0019 #    pythiaPylistVerbosity = cms.untracked.int32(12),
0020     filterEfficiency = cms.untracked.double(1.0),
0021     comEnergy = cms.double(10000.0),
0022 #    crossSection = cms.untracked.double(55000000000.),
0023     UseExternalGenerators = cms.untracked.bool(False),
0024 #
0025     PythiaParameters = cms.PSet(
0026         pythiaUESettingsBlock,
0027         pythiaMyParameters = cms.vstring(
0028 # This is needed if you are introducing long-lived exotic particles.
0029           'MSTJ(22)=1    ! Decay ALL unstable particles',
0030
0031           'MSEL=0',
0032 # Request gg -> H0 production
0033           'MSUB(152)=1',
0034           'MWID(35)=2 ! Let me set H0 properties'
0035         ),
0036 #
0037 # This block controls how Pythia interacts with the PYUPDA cards.
0038 #
0039         PYUPDAParameters = cms.vstring(
0040 # Either:
0041 #     1) Specify the location of the PYUPDA table to be read in.
0042          "PYUPDAFILE = 'Configuration/Generator/data/Pythia_H0_pyupda.in'"
0043 #        Optionally ask to call PYUPDA after PYINIT. Don't do this unless you have to.
0044 #        "PYUPDApostPYINIT"
0045 # Or:
0046 #     2) Write current PYUPDA parameters to file (so you can edit them and read back in).
0047 #         "PYUPDAFILE = \'pyupda.out\'"
0048 #         "PYUPDAWRITE"
0049         ),
0050 #
0051         # This is a vector of ParameterSet names to be read, in this order
0052         parameterSets = cms.vstring('pythiaUESettings',
0053                                     'pythiaMyParameters',
0054                                     'PYUPDAParameters')
0055     )
0056 )
0057
0058 # N.B. If your PYUPDA tables introduces new exotic particles, you will need
0059 # to include:
0060 #
0061 # genParticles.abortOnUnknownPDGCode = cms.untracked.bool(False)
0062
0063 # This is no longer necessary, as starting 34X cmsDriver takes care automatically
0064 #
0065 #ProductionFilterSequence = cms.Sequence(generator)