import FWCore.ParameterSet.Config as cms

generator = cms.EDFilter("Herwig6HadronizerFilter",
	comEnergy                = cms.double(8000.0),
	crossSection             = cms.untracked.double(211.1),
	doMPInteraction          = cms.bool(True),
	emulatePythiaStatusCodes = cms.untracked.bool(True),
	filterEfficiency         = cms.untracked.double(1.0),
	herwigHepMCVerbosity     = cms.untracked.bool(False),
	herwigVerbosity          = cms.untracked.int32(0),
	lhapdfSetPath            = cms.untracked.string(''),
	maxEventsToPrint         = cms.untracked.int32(0),
	printCards               = cms.untracked.bool(False),
	useJimmy                 = cms.bool(True),

	HerwigParameters = cms.PSet(
		herwigUEsettings = cms.vstring(
			'JMUEO     = 1       ! multiparton interaction model',
			'PTJIM     = 4.189   ! 2.8x(sqrt(s)/1.8TeV)^0.27 @ 8 TeV',
			'JMRAD(73) = 1.8     ! inverse proton radius squared',
			'PRSOF     = 0.0     ! prob. of a soft underlying event',
			'MAXER     = 1000000 ! max error'
		),
                herwigMcatnlo = cms.vstring(
			'PTMIN      = 0.5    ! minimum pt in hadronic jet'
		),
		parameterSets = cms.vstring('herwigUEsettings',
                                            'herwigMcatnlo')
	)
)

dileptonFilter = cms.EDFilter("MCParticlePairFilter",
    ParticleID1 = cms.untracked.vint32(11,13,15,-11,-13,-15),
    ParticleID2 = cms.untracked.vint32(11,13,15,-11,-13,-15),
    MinPt = cms.untracked.vdouble(2*[15]),
    Status = cms.untracked.vint32(2*[3])
)

ProductionFilterSequence = cms.Sequence(generator * dileptonFilter )