import FWCore.ParameterSet.Config as cms

externalLHEProducer = cms.EDProducer("ExternalLHEProducer",
    args = cms.vstring('/cvmfs/cms.cern.ch/phys_generator/gridpacks/2017/13TeV/madgraph/V5_2.6.5/HeavyNeutrino_Trilepton_Tau/HeavyNeutrino_trilepton_M-70_V-0p01_tau_LO_slc6_amd64_gcc700_CMSSW_10_2_24_patch1_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')
)
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("Pythia8HadronizerFilter",
    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,
        parameterSets = cms.vstring('pythia8CommonSettings',
                                    'pythia8CP5Settings',
                                    'pythia8PSweightsSettings',
                                    )
    )
)

from PhysicsTools.HepMCCandAlgos.genParticles_cfi import genParticles
from PhysicsTools.JetMCAlgos.TauGenJets_cfi import tauGenJets
from PhysicsTools.JetMCAlgos.TauGenJetsDecayModeSelectorAllHadrons_cfi import tauGenJetsSelectorAllHadrons

genParticles.src = cms.InputTag("generator", "unsmeared")

tauGenJets.GenParticles = cms.InputTag("genParticles")
tauGenJets.includeNeutrinos = cms.bool(False)

genVisTauSelector = cms.EDFilter("CandViewSelector",
    src = cms.InputTag("tauGenJetsSelectorAllHadrons"),
    cut = cms.string("pt > 18")
  )             
        
genVisTauFilter = cms.EDFilter("CandViewCountFilter",
     src = cms.InputTag("genVisTauSelector"),
     minNumber = cms.uint32(1),
  )

lightleptonfilter = cms.EDFilter("MCMultiParticleFilter",
    PtMin = cms.vdouble(15., 15.),
    ParticleID = cms.vint32(11, 13),
    NumRequired = cms.int32(1),
    AcceptMore = cms.bool(True),
    EtaMax = cms.vdouble(3., 3.),
    Status = cms.vint32(1, 1)
)

ProductionFilterSequence = cms.Sequence(generator*genParticles*tauGenJets*tauGenJetsSelectorAllHadrons*genVisTauSelector*genVisTauFilter*lightleptonfilter)


# Link to generator fragment:
# genFragments/Hadronizer/13TeV/Hadronizer_TuneCP5_13TeV_generic_LHE_pythia8_PSweights_tauhFilter_cff.py