import FWCore.ParameterSet.Config as cms

process = cms.Process("PROD")

process.load("Configuration.StandardSequences.SimulationRandomNumberGeneratorSeeds_cff")
process.printTree = cms.EDAnalyzer("ParticleListDrawer",
                                     maxEventsToPrint = cms.untracked.int32(5),
                                     printVertex = cms.untracked.bool(False),
                                     src = cms.InputTag("genParticles")
                                   )


#process.source = cms.Source("EmptySource")
process.source = cms.Source("LHESource",
    fileNames = cms.untracked.vstring('file:fourtop500_v4.lhe')
)

process.generator = cms.EDFilter("Pythia8HadronizerFilter",
    maxEventsToPrint = cms.untracked.int32(5),
    pythiaPylistVerbosity = cms.untracked.int32(1),
    filterEfficiency = cms.untracked.double(1.0),
    pythiaHepMCVerbosity = cms.untracked.bool(False),
    comEnergy = cms.double(7000.),
    LHEInputFileName = cms.untracked.string('fourtop500_v4.lhe'),
#    PythiaParameters = cms.PSet(
#        pythia8_example03 = cms.vstring(''),
#        parameterSets = cms.vstring('pythia8_example03')    
#    )
    PythiaParameters = cms.PSet(
	 pythia8_black = cms.vstring( ## see details on http://home.thep.lu.se/~torbjorn/php8135/ExtraDimensionalProcesses.php?filepath=files/
#    'Top:gg2ttbar = on',
#    'Top:qqbar2ttbar = on',
    'Tune:pp 5',
    'PDF:pSet = 5'    
#    'PartonLevel:MI = on'
#    'PartonLevel:ISR = on',
#    'PartonLevel:FSR = on'
         ),
 	 parameterSets = cms.vstring('pythia8_black')
    )
)

process.MessageLogger = cms.Service("MessageLogger",
    cout = cms.untracked.PSet(
        default = cms.untracked.PSet(
            limit = cms.untracked.int32(0)
        )
    ),
    destinations = cms.untracked.vstring('cout')
)

process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
    moduleSeeds = cms.PSet(
        generator = cms.untracked.uint32(123456),
        g4SimHits = cms.untracked.uint32(123456788),
        VtxSmeared = cms.untracked.uint32(123456789)
    ),
)

process.maxEvents = cms.untracked.PSet(
    input = cms.untracked.int32(10)
)

process.GEN = cms.OutputModule("PoolOutputModule",
    fileName = cms.untracked.string('fourtop500_v4.root')
)

process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")

process.p = cms.Path(process.generator*process.genParticles)

process.outpath = cms.EndPath(process.GEN)

process.schedule = cms.Schedule(process.p, process.outpath)