Pythia8Interface/test/Py8_Z2tautau_tauola_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004
0005 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0006 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0007 process.load("PhysicsTools.HepMCCandAlgos.genParticles_cfi")
0008
0009 process.source = cms.Source("EmptySource",
0010      firstLuminosityBlock = cms.untracked.uint32(1),
0011 #     numberEventsInLuminosityBlock = cms.untracked.uint32(200),
0012      firstEvent = cms.untracked.uint32(1),
0013      firstRun = cms.untracked.uint32(1),
0014      numberEventsInRun = cms.untracked.uint32(100)
0015 )
0016
0017 from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
0018
0019 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0020     maxEventsToPrint = cms.untracked.int32(1),
0021     pythiaPylistVerbosity = cms.untracked.int32(1),
0022     filterEfficiency = cms.untracked.double(1.0),
0023     pythiaHepMCVerbosity = cms.untracked.bool(False),
0024     pythiaHepMCVerbosityParticles = cms.untracked.bool(True),
0025     comEnergy = cms.double(8000.),
0026     ExternalDecays = cms.PSet(
0027         Tauola = cms.untracked.PSet(
0028             TauolaPolar,
0029         TauolaDefaultInputCards
0030         ),
0031         parameterSets = cms.vstring('Tauola')
0032     ),
0033
0034     PythiaParameters = cms.PSet(
0035         py8ProcessSettings = cms.vstring( 'WeakSingleBoson:ffbar2gmZ = on'
0036     ),
0037     py8ZDecaySettings = cms.vstring(  '23:onMode = off', # turn OFF all Z decays
0038                       '23:onIfAny = 15'  # turn ON Z->tautau
0039     ),
0040     py8TauDecaySettings = cms.vstring('15:onMode = off', # turn OFF all tau decays
0041     ),
0042         parameterSets = cms.vstring(
0043                                 'py8ProcessSettings',
0044                     'py8ZDecaySettings')
0045     )
0046 )
0047
0048 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0049 #process.MessageLogger = cms.Service("MessageLogger",
0050 #    cout = cms.untracked.PSet(
0051 #        default = cms.untracked.PSet(
0052 #            limit = cms.untracked.int32(100)
0053 #        )
0054 #    ),
0055 #    destinations = cms.untracked.vstring('cout')
0056 #)
0057
0058 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0059     generator = cms.PSet(
0060         initialSeed = cms.untracked.uint32(123456789),
0061         engineName = cms.untracked.string('HepJamesRandom')
0062     )
0063 )
0064
0065 process.maxEvents = cms.untracked.PSet(
0066     input = cms.untracked.int32(100)
0067 )
0068
0069 process.GEN = cms.OutputModule("PoolOutputModule",
0070     fileName = cms.untracked.string('Py8_Z2tautau_tauola.root')
0071 )
0072
0073 process.genParticles = cms.EDProducer("GenParticleProducer",
0074     saveBarCodes = cms.untracked.bool(True),
0075     src = cms.InputTag("generator:unsmeared"),
0076     abortOnUnknownPDGCode = cms.untracked.bool(False)
0077 )
0078
0079 process.printGenParticles = cms.EDAnalyzer("ParticleListDrawer",
0080                                            src = cms.InputTag("genParticles"),
0081                                            maxEventsToPrint = cms.untracked.int32(3) )
0082
0083 process.p = cms.Path(process.generator*process.genParticles*process.printGenParticles)
0084 process.outpath = cms.EndPath(process.GEN)
0085
0086 process.schedule = cms.Schedule(process.p, process.outpath)
0087