Generator/python/TTbar012Jets_5f_NLO_FXFX_Madgraph_LHE_13TeV_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002 from Configuration.Generator.Pythia8CommonSettings_cfi import *
0003 from Configuration.Generator.Pythia8CUEP8M1Settings_cfi import *
0004 from Configuration.Generator.Pythia8aMCatNLOSettings_cfi import *
0005
0006 externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
0007     args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/2017/13TeV/madgraph/V5_2.6.0/Validation/tt012j_5f_ckm_NLO_FXFX_VMG5_260_slc6_amd64_gcc481_CMSSW_7_1_30_tarball.tar.xz'),
0008     nEvents = cms.untracked.uint32(5000),
0009     numberOfParameters = cms.uint32(1),
0010     outputFile = cms.string('cmsgrid_final.lhe'),
0011     scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh'),
0012     generateConcurrently = cms.untracked.bool(True)
0013 )
0014
0015 #Link to datacards:
0016 #https://github.com/cms-sw/genproductions/tree/c41ab29f3d86c9e53df8b0d76c12cd519adbf013/bin/MadGraph5_aMCatNLO/cards/production/13TeV/tt012j_5f_ckm_NLO_FXFX
0017
0018 generator = cms.EDFilter("Pythia8HadronizerFilter",
0019     maxEventsToPrint = cms.untracked.int32(1),
0020     pythiaPylistVerbosity = cms.untracked.int32(1),
0021     filterEfficiency = cms.untracked.double(1.0),
0022     pythiaHepMCVerbosity = cms.untracked.bool(False),
0023     comEnergy = cms.double(13000.),
0024     ExternalDecays = cms.PSet(
0025         EvtGen130 = cms.untracked.PSet(
0026             decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
0027             particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0028             list_forced_decays = cms.vstring(),
0029             operates_on_particles = cms.vint32()
0030             ),
0031         parameterSets = cms.vstring('EvtGen130')
0032         ),
0033     PythiaParameters = cms.PSet(
0034         pythia8CommonSettingsBlock,
0035         pythia8CUEP8M1SettingsBlock,
0036         pythia8aMCatNLOSettingsBlock,
0037         processParameters = cms.vstring(
0038             'JetMatching:setMad = off',
0039             'JetMatching:scheme = 1',
0040             'JetMatching:merge = on',
0041             'JetMatching:jetAlgorithm = 2',
0042             'JetMatching:etaJetMax = 999.',
0043             'JetMatching:coneRadius = 1.',
0044             'JetMatching:slowJetPower = 1',
0045             'JetMatching:qCut = 40.', #this is the actual merging scale
0046             'JetMatching:doFxFx = on',
0047             'JetMatching:qCutME = 20.',#this must match the ptj cut in the lhe generation step
0048             'JetMatching:nQmatch = 5', #4 corresponds to 4-flavour scheme (no matching of b-quarks), 5 for 5-flavour scheme
0049             'JetMatching:nJetMax = 2', #number of partons in born matrix element for highest multiplicity
0050         ),
0051         parameterSets = cms.vstring('pythia8CommonSettings',
0052                                     'pythia8CUEP8M1Settings',
0053                                     'pythia8aMCatNLOSettings',
0054                                     'processParameters',
0055                                     )
0056     )
0057 )
0058
0059 ProductionFilterSequence = cms.Sequence(externalLHEProducer*generator)