Line Code
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48
# Copied from https://github.com/cms-sw/genproductions for RelVal June 5, 2014
import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *

generator = cms.EDFilter("Pythia8HadronizerFilter",
                         ExternalDecays = cms.PSet(
    EvtGen = cms.untracked.PSet(
    use_default_decay = cms.untracked.bool(True),
    decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
    particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
    user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Validation.dec'),
    list_forced_decays = cms.vstring(),
    operates_on_particles = cms.vint32(0)
    ),
    parameterSets = cms.vstring('EvtGen')
    ),
                         UseExternalGenerators = cms.untracked.bool(True),
                         maxEventsToPrint = cms.untracked.int32(1),
                         pythiaPylistVerbosity = cms.untracked.int32(1),
                         filterEfficiency = cms.untracked.double(1.0),
                         pythiaHepMCVerbosity = cms.untracked.bool(False),
                         comEnergy = cms.double(13000.),
                         jetMatching = cms.untracked.PSet(
    scheme = cms.string("Madgraph"),
    mode = cms.string("auto"),# soup, or "inclusive" / "exclusive"
    MEMAIN_etaclmax = cms.double(-1),
    MEMAIN_qcut = cms.double(-1),
    MEMAIN_minjets = cms.int32(-1),
    MEMAIN_maxjets = cms.int32(-1),
    MEMAIN_showerkt = cms.double(0), # use 1=yes only for pt-ordered showers !
    MEMAIN_nqmatch = cms.int32(5), #PID of the flavor until which the QCD radiation are kept in the matching procedure;
    # if nqmatch=4, then all showered partons from b's are NOT taken into account
    # Note (JY): I think the default should be 5 (b); anyway, don't try -1 as it'll result in a throw...
    MEMAIN_excres = cms.string(""),
    outTree_flag = cms.int32(0) # 1=yes, write out the tree for future sanity check
    ),
                         PythiaParameters = cms.PSet(
        pythia8CommonSettingsBlock,
        pythia8CUEP8M1SettingsBlock,
        parameterSets = cms.vstring('pythia8CommonSettings',
                                    'pythia8CUEP8M1Settings',
                                    )
        )
                         )

ProductionFilterSequence = cms.Sequence(generator)