Pythia6Interface/test/Py6GenFilter_stophadrons_cfg.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 process = cms.Process("TEST")
0004 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0005 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0006 #process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0007
0008
0009 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0010     generator = cms.PSet(
0011         initialSeed = cms.untracked.uint32(123456789),
0012         engineName = cms.untracked.string('HepJamesRandom')
0013     )
0014 )
0015
0016
0017 # The following three lines reduce the clutter of repeated printouts
0018 # of the same exception message.
0019 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0020
0021
0022
0023
0024 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(1))
0025
0026 process.source = cms.Source("EmptySource")
0027
0028 from Configuration.Generator.PythiaUESettings_cfi import *
0029
0030 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0031     pythiaHepMCVerbosity = cms.untracked.bool(True),
0032     maxEventsToPrint = cms.untracked.int32(2),
0033     pythiaPylistVerbosity = cms.untracked.int32(1),
0034     comEnergy = cms.double(10000.0),
0035
0036     stopHadrons = cms.bool(True),
0037     gluinoHadrons = cms.bool(False),
0038
0039     PythiaParameters = cms.PSet(
0040
0041         pythiaUESettingsBlock,
0042     #
0043         # settings for stop-hadrons have been taken from original example Py6 add-on
0044     # http://projects.hepforge.org/pythia6/examples/main78.f
0045     #
0046     pythiaStopHadrons = cms.vstring(
0047            'MSEL=0           !  User defined processes',
0048            'IMSS(1)=1        !  brute force',
0049            'MSUB(261)=1      !  subprocess',
0050            'MSUB(264)=1      !  subprocess',
0051            'IMSS(3)=1',
0052            'RMSS(3)=1000.',
0053            'RMSS(1)=1000.',
0054            'RMSS(2)=1000.',
0055            'RMSS(4)=5000.',
0056            'MDCY(302,1)=0    ! set stop stable',
0057        'MWID(302)=0',
0058            'IMSS(5)=1',
0059            'RMSS(12)=200.',
0060        'RMSS(10)=300.'
0061            ### 'MSTJ(14)=-1', # this is actually hardcoded in Py6Had class
0062            ### 'MSTP(111)=0'  # this is a MANDATORY, thus it's hardcoded in Py6Had class
0063     ),
0064         # This is a vector of ParameterSet names to be read, in this order
0065         parameterSets = cms.vstring(
0066         'pythiaUESettings',
0067             'pythiaStopHadrons')
0068     )
0069 )
0070
0071 process.GEN = cms.OutputModule("PoolOutputModule",
0072     fileName = cms.untracked.string('stopHadrons.root')
0073 )
0074
0075 process.p = cms.Path(process.generator)
0076 process.outpath = cms.EndPath(process.GEN)
0077
0078 process.schedule = cms.Schedule(process.p, process.outpath)