import FWCore.ParameterSet.Config as cms

from Configuration.Generator.PythiaUEZ2Settings_cfi import *

generator = cms.EDFilter('Pythia6GeneratorFilter',

	pythiaPylistVerbosity = cms.untracked.int32(0),
	pythiaHepMCVerbosity = cms.untracked.bool(False),
	comEnergy = cms.double(8000.0),
	crossSection = cms.untracked.double(2978915.),     # Given by PYTHIA after running 
	filterEfficiency = cms.untracked.double(0.000295), # Given by PYTHIA after running
	maxEventsToPrint = cms.untracked.int32(0),

	ExternalDecays = cms.PSet(
		EvtGen = cms.untracked.PSet(
			operates_on_particles = cms.vint32(521), 
			use_default_decay = cms.untracked.bool(False),
			decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
			particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
			user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bu_Psi2SK_JpsiPiPiK.dec'),
			list_forced_decays = cms.vstring('MyB+', 'MyB-'),
			),
		parameterSets = cms.vstring('EvtGen')
	),
	PythiaParameters = cms.PSet(
		pythiaUESettingsBlock,
		bbbarSettings = cms.vstring('MSEL = 1'), 
        # This is a vector of ParameterSet names to be read, in this order
	        parameterSets = cms.vstring(
             		'pythiaUESettings',
             		'bbbarSettings')
    )


	)

bfilter = cms.EDFilter(
        "PythiaFilter",
        MaxEta = cms.untracked.double(9999.),
        MinEta = cms.untracked.double(-9999.),
        ParticleID = cms.untracked.int32(521),
	MinPt = cms.untracked.double(7.)
        )

mumugenfilter = cms.EDFilter("MCParticlePairFilter",
	Status = cms.untracked.vint32(1, 1),
	MinPt = cms.untracked.vdouble(0.0, 0.0),
	MaxEta = cms.untracked.vdouble(2.5, 2.5),
	MinEta = cms.untracked.vdouble(-2.5, -2.5),
	ParticleCharge = cms.untracked.int32(-1),
	ParticleID1 = cms.untracked.vint32(13),
	ParticleID2 = cms.untracked.vint32(13)
	)

ProductionFilterSequence = cms.Sequence(generator*bfilter*mumugenfilter)