import FWCore.ParameterSet.Config as cms

#Made with AutoClone script: /afs/cern.ch/user/x/xulyu/work/public/mcrequest/btWprimeToBottomTop.py
externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
 args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/2017/13TeV/madgraph/V5_2.6.5/WPrime_tbOnly_FL/WPrime_tbOnly_LL_M800_slc7_amd64_gcc700_CMSSW_10_6_19_tarball.tar.xz'),
 nEvents = cms.untracked.uint32(5000),
 numberOfParameters = cms.uint32(1),
 outputFile = cms.string('cmsgrid_final.lhe'),
 scriptName = cms.FileInPath('GeneratorInterface/LHEInterface/data/run_generic_tarball_cvmfs.sh'),
    generateConcurrently = cms.untracked.bool(True)
)

#Link to datacards:
#https://github.com/cms-sw/genproductions/tree/fbe8787313a801b1cecd6e071a083a2710658d2a/bin/MadGraph5_aMCatNLO/cards/production/13TeV/WPrime_tbOnly
#PR: https://github.com/cms-sw/genproductions/pull/3200

import FWCore.ParameterSet.Config as cms

from Configuration.Generator.Pythia8CommonSettings_cfi import *
from Configuration.Generator.MCTunes2017.PythiaCP5Settings_cfi import *
from Configuration.Generator.PSweightsPythia.PythiaPSweightsSettings_cfi import *

generator = cms.EDFilter('Pythia8ConcurrentHadronizerFilter',
                          	
                            
    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.),

    PythiaParameters = cms.PSet(
        pythia8CommonSettingsBlock,
        pythia8CP5SettingsBlock,
        pythia8PSweightsSettingsBlock,
        processParameters = cms.vstring(
            'JetMatching:setMad = off',
            'JetMatching:scheme = 1',
            'JetMatching:merge = on',
           'JetMatching:jetAlgorithm = 2',
           'JetMatching:etaJetMax = 5.',
           'JetMatching:coneRadius = 1.',
           'JetMatching:slowJetPower = 1',
           'JetMatching:qCut = 80.', #this is the actual merging scale
           'JetMatching:nQmatch = 5', #4 corresponds to 4-flavour scheme (no matching of b-quarks), 5 for 5-flavour scheme
           'JetMatching:nJetMax = 1', #number of partons in born matrix element for highest multiplicity
           'JetMatching:doShowerKt = off', #off for MLM matching, turn on for shower-kT matching
        ),

        
        parameterSets = cms.vstring('pythia8CommonSettings',
                                    'pythia8CP5Settings',
                                    'pythia8PSweightsSettings',
                                    'processParameters',
                                    ),
               )
             )