from Configuration.Generator.HerwigppDefaults_cfi import *
from Configuration.Generator.HerwigppUE_CUETHS1_cfi import *
from Configuration.Generator.HerwigppPDF_CTEQ6_LO_cfi import *
from Configuration.Generator.HerwigppEnergy_13TeV_cfi import *

generator = cms.EDFilter("ThePEGGeneratorFilter",
	herwigDefaultsBlock,
	herwigppUESettingsBlock,
	herwigppPDFSettingsBlock,
	herwigppEnergySettingsBlock,

    crossSection = cms.untracked.double(1),
    filterEfficiency = cms.untracked.double(1),
    configFiles = cms.vstring(),
    
    parameterSets = cms.vstring(
      'hwpp_cmsDefaults',
      'hwpp_ue_CUETHS1',
      'hwpp_pdf_CTEQ6L1_CUETHS1',
      'hwpp_cm_13TeV',
      'productionParameters',
    ),
    productionParameters = cms.vstring(                                                                    
      'cd /Herwig/MatrixElements/',
      'insert SimpleQCD:MatrixElements[0] MEQCD2to2',
      'cd /',
      'set /Herwig/Cuts/JetKtCut:MinKT 10*GeV',
      'set /Herwig/Cuts/JetKtCut:MaxKT 35*GeV',
      'set /Herwig/Cuts/QCDCuts:MHatMin 0.0*GeV',
      'set /Herwig/UnderlyingEvent/MPIHandler:IdenticalToUE 0',
    ),
)

from PhysicsTools.HepMCCandAlgos.genParticles_cfi import genParticles
from RecoJets.Configuration.GenJetParticles_cff import genParticlesForJets
from RecoJets.JetProducers.ak4GenJets_cfi import ak4GenJets
 
 
backwardJetSelector = cms.EDFilter("CandViewSelector",
      src = cms.InputTag("ak4GenJets"),
      cut = cms.string('pt > 10 & eta < -2.8')
)
 
backwardJetFilter = cms.EDFilter("CandViewCountFilter",
    src = cms.InputTag("backwardJetSelector"),
    minNumber = cms.uint32(1)
)
 
forwardJetSelector = cms.EDFilter("CandViewSelector",
      src = cms.InputTag("ak4GenJets"),
      cut = cms.string('pt > 10 & eta > 2.8')
)
 
forwardJetFilter = cms.EDFilter("CandViewCountFilter",
    src = cms.InputTag("forwardJetSelector"),
    minNumber = cms.uint32(1)
)
 
ProductionFilterSequence = cms.Sequence(generator+genParticles+genParticlesForJets+ak4GenJets+
                                        backwardJetSelector+backwardJetFilter+
                                        forwardJetSelector+forwardJetFilter)