################
# Gen fragment #
################


import FWCore.ParameterSet.Config as cms


configurationMetadata = cms.untracked.PSet(
	version = cms.untracked.string('$Revision: 1.8 $'),
	name = cms.untracked.string('$Source: /local/reps/CMSSW/UserCode/Dinardo/B0Analysis/B0KstMuMu/python/PYTHIA6_B0dToJPsiKstMuMuKPi_7TeV_cff.py,v $'),
	annotation = cms.untracked.string('Summer12: Pythia6+EvtGen generation of B0d --> K*0(K pi) J/Psi(Mu+Mu-), 8TeV tune'))


#from Configuration.Generator.PythiaUEZ2Settings_cfi import *
from Configuration.Generator.PythiaUEZ2starSettings_cfi import *


generator = cms.EDFilter("Pythia6GeneratorFilter",
			 pythiaPylistVerbosity = cms.untracked.int32(0),
			 pythiaHepMCVerbosity = cms.untracked.bool(False),
			 #comEnergy = cms.double(7000.0),
			 comEnergy = cms.double(8000.0),
			 crossSection = cms.untracked.double(0.0),     # Given by PYTHIA after running
			 filterEfficiency = cms.untracked.double(1.0), # Given by PYTHIA after running
			 maxEventsToPrint = cms.untracked.int32(0),
			 
			 ExternalDecays = cms.PSet(
	EvtGen = cms.untracked.PSet(
	operates_on_particles = cms.vint32(0),
	use_default_decay = cms.untracked.bool(False),
	decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY.DEC'),
	particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
	user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Bd_JpsiKstar_mumuKpi.dec'),
	list_forced_decays = cms.vstring('MyB0','Myanti-B0')),
	parameterSets = cms.vstring('EvtGen')),
			 
			 PythiaParameters = cms.PSet(
	pythiaUESettingsBlock,
	bbbarSettings = cms.vstring('MSEL = 1'),
	parameterSets = cms.vstring('pythiaUESettings','bbbarSettings')))


###########
# Filters #
###########
# Filter only pp events which produce a B0: efficiency ~6e-3
b0filter = cms.EDFilter("PythiaFilter", ParticleID = cms.untracked.int32(511))

# Filter only pp events which produce a J/Psi OR Psi(2S)
oniafilter = cms.EDFilter("MCSingleParticleFilter",
			  MaxEta = cms.untracked.vdouble(1000.0, 1000.0),
			  Status = cms.untracked.vint32(2, 2),
			  MinEta = cms.untracked.vdouble(-1000.0, -1000.0),
			  MinPt = cms.untracked.vdouble(0.0, 0.0),
			  ParticleID = cms.untracked.vint32(443, 100443))

# Filter on final state muons
mumugenfilter = cms.EDFilter("MCParticlePairFilter",
			     Status = cms.untracked.vint32(1, 1),
			     MinPt = cms.untracked.vdouble(3.3, 3.3),
			     MaxEta = cms.untracked.vdouble(2.3, 2.3),
			     MinEta = cms.untracked.vdouble(-2.3, -2.3),
			     ParticleCharge = cms.untracked.int32(-1),
			     ParticleID1 = cms.untracked.vint32(13),
			     ParticleID2 = cms.untracked.vint32(13))


ProductionFilterSequence = cms.Sequence(generator*b0filter*oniafilter*mumugenfilter)