Pythia8Interface/python/FxFxEWK.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
0004     args = cms.vstring('/afs/cern.ch/user/c/cvicovil/public/TTWJetsToLNu_5f_NLO_FXFX_slc7_amd64_gcc900_CMSSW_12_0_2_tarball.tar.xz'),
0005     nEvents = cms.untracked.uint32(1000),
0006     numberOfParameters = cms.uint32(1),
0007     outputFile = cms.string('cmsgrid_final.lhe'),
0008     scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh')
0009 )
0010
0011 # Link to datacards:
0012 # https://github.com/cms-sw/genproductions/tree/5fb3762c8be13dabb89a8580863856201d56f49c/bin/MadGraph5_aMCatNLO/cards/production/2017/13TeV/TTWJets/TTWJetsToLNu_5f_NLO_FXFX
0013
0014 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0015 from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
0016 from Configuration.Generator.Pythia8aMCatNLOSettings_cfi import *
0017 from Configuration.Generator.PSweightsPythia.PythiaPSweightsSettings_cfi import *
0018
0019 generator = cms.EDFilter("Pythia8HadronizerFilter",
0020     maxEventsToPrint = cms.untracked.int32(1),
0021     pythiaPylistVerbosity = cms.untracked.int32(1),
0022     filterEfficiency = cms.untracked.double(1.0),
0023     pythiaHepMCVerbosity = cms.untracked.bool(False),
0024     comEnergy = cms.double(13000.),
0025     PythiaParameters = cms.PSet(
0026         pythia8CommonSettingsBlock,
0027         pythia8CP5SettingsBlock,
0028         pythia8aMCatNLOSettingsBlock,
0029         pythia8PSweightsSettingsBlock,
0030         processParameters = cms.vstring(
0031             'JetMatching:setMad = off',
0032             'JetMatching:scheme = 1',
0033             'JetMatching:doVeto = off',  # Turn OFF Pythia JetMatching
0034             'JetMatching:merge = on',
0035             'JetMatching:jetAlgorithm = 2',
0036             'JetMatching:etaJetMax = 1000',
0037             'JetMatching:coneRadius = 1.',
0038             'JetMatching:slowJetPower = 1',
0039             'JetMatching:qCut = 42.', #this is the actual merging scale
0040             'JetMatching:doFxFx = on',
0041             'JetMatching:qCutME = 20.',#this must match the ptj cut in the lhe generation step
0042             'JetMatching:nQmatch = 5', #4 corresponds to 4-flavour scheme (no matching of b-quarks), 5 for 5-flavour scheme
0043             'JetMatching:nJetMax = 1', #number of partons in born matrix element for highest multiplicity
0044
0045         ),
0046         parameterSets = cms.vstring('pythia8CommonSettings',
0047                                     'pythia8CP5Settings',
0048                                     'pythia8aMCatNLOSettings',
0049                                     'pythia8PSweightsSettings',
0050                                     'processParameters',
0051                                     )
0052     ),
0053     UserCustomization = cms.VPSet(
0054            cms.PSet(
0055              pluginName = cms.string("JetMatchingEWKFxFx"), # Call the new JetMatching plugin
0056            )
0057         ),
0058 )
0059