Pythia8Interface/test/biasedTauDecayer.py

0001 #example gen fragment that takes a gridapack which produces
0002 #ttH events with inclusive top decays and undecayed higgs at lhe level
0003 #and selects resonance decays such that events have at least four leptons (electrons, muons, taus)
0004
0005 import FWCore.ParameterSet.Config as cms
0006
0007 process = cms.Process("PROD")
0008
0009 process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
0010
0011 # link to cards:
0012 # https://github.com/cms-sw/genproductions/blob/0d4b4288fa053d9a8aef5c6e123b66bf94c3aee8/bin/Powheg/production/V2/13TeV/Higgs/ttH_inclusive_NNPDF30_13TeV_M125/ttH_inclusive_NNPDF30_13TeV_M125.input
0013
0014 process.source = cms.Source("LHESource",
0015   fileNames = cms.untracked.vstring('file:powheg-Ztautau.lhe')
0016 )
0017
0018 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0019 from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
0020 from Configuration.Generator.Pythia8PowhegEmissionVetoSettings_cfi import *
0021
0022 process.generator = cms.EDFilter("Pythia8HadronizerFilter",
0023                          maxEventsToPrint = cms.untracked.int32(20),
0024                          pythiaPylistVerbosity = cms.untracked.int32(1),
0025                          filterEfficiency = cms.untracked.double(1.0),
0026                          pythiaHepMCVerbosity = cms.untracked.bool(False),
0027                          comEnergy = cms.double(13000.),
0028                          PythiaParameters = cms.PSet(
0029         pythia8CommonSettingsBlock,
0030         pythia8CUEP8M1SettingsBlock,
0031         processParameters = cms.vstring(
0032             'BiasedTauDecayer:filter = on',
0033             'PartonLevel:MPI = off'
0034           ),
0035         parameterSets = cms.vstring('pythia8CommonSettings',
0036                                     'pythia8CUEP8M1Settings',
0037                                     'processParameters'
0038                                     )
0039         )
0040                          )
0041
0042 process.p = cms.Path(process.generator)
0043
0044