import FWCore.ParameterSet.Config as cms

externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/UL/13TeV/madgraph/V5_2.6.5/gammaJets/GJets_Pta-40_slc7_amd64_gcc700_CMSSW_10_6_0_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')
)


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

generator = cms.EDFilter("Pythia8HadronizerFilter",
                         HepMCFilter = cms.PSet(
                              filterName = cms.string('PythiaHepMCFilterGammaGamma'),
                              filterParameters = cms.PSet(
                              AcceptPrompts = cms.bool(True),
                              PtSeedThr = cms.double(5.0),
                              NTkConeSum = cms.int32(4),
                              InvMassMin = cms.double(80.0),
                              InvMassMax = cms.double(14000.0),
                              EtaElThr = cms.double(2.8),
                              EtaSeedThr = cms.double(2.8),
                              EtaMaxCandidate = cms.double(3.0),
                              EnergyCut = cms.double(1.0),
                              dPhiSeedMax = cms.double(0.2),
                              dEtaSeedMax = cms.double(0.12),
                              EtaTkThr = cms.double(2.2),
                              dRNarrowCone = cms.double(0.02),
                              PtElThr = cms.double(2.0),
                              NTkConeMax = cms.int32(2),
                              EtaGammaThr = cms.double(2.8),
                              PromptPtThreshold = cms.double(15.0),
                              PtMinCandidate1 = cms.double(15.0),
                              PtMinCandidate2 = cms.double(15.0),
                              PtGammaThr = cms.double(0.0),
                              PtTkThr = cms.double(1.6),
                              dRSeedMax = cms.double(0.0),
                              dRTkMax = cms.double(0.2),
                        )
    ),
	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.),
    nAttempts = cms.uint32(20), # extra
    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 = 19.', #this is the actual merging scale
            'JetMatching:nQmatch = 5', #5 for 5-flavour scheme (matching of b-quarks)
            'JetMatching:nJetMax = 4', #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',
                                    )
    )
)

# Link to generator fragment:
#https://raw.githubusercontent.com/cms-sw/genproductions/a10c4661b08ff764b90a8b1947689c52c5c992fb/python/ThirteenTeV/Hadronizer/Hadronizer_TuneCP5_13TeV_MLM_5f_max4j_qCut19_LHE_pythia8_cff.py