ExternalDecays/test/Py8GenFilter_Onia_EvtGen_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
0008 process.source = cms.Source("EmptySource")
0009
0010 process.generator = cms.EDFilter("Pythia8GeneratorFilter",
0011     maxEventsToPrint = cms.untracked.int32(3),
0012     pythiaPylistVerbosity = cms.untracked.int32(1),
0013     filterEfficiency = cms.untracked.double(1.0),
0014     pythiaHepMCVerbosity = cms.untracked.bool(True),
0015     comEnergy = cms.double(7000.),
0016
0017     ExternalDecays = cms.PSet(
0018         EvtGen1 = cms.untracked.PSet(
0019              decay_table = cms.string('GeneratorInterface/EvtGenInterface/data/DECAY_2010.DEC'),
0020              particle_property_file = cms.FileInPath('GeneratorInterface/EvtGenInterface/data/evt.pdl'),
0021              #user_decay_files = cms.vstring('DECAY_2010.DEC'),
0022              user_decay_file = cms.vstring('GeneratorInterface/ExternalDecays/data/Onia_mumu.dec'),
0023              list_forced_decays = cms.vstring('MyUpsilon(2S)'),
0024              operates_on_particles = cms.vint32(0) # 0 (zero) means default list (hardcoded), the list of PDG IDs can be put here
0025              ),
0026         parameterSets = cms.vstring('EvtGen1')
0027     ),
0028
0029     PythiaParameters = cms.PSet(      # MultipleInteractions commands below don't work (obsolete?). Should official settings be used?
0030
0031         py8UESettings = cms.vstring( 'StringZ:usePetersonB = on',      # these 2 together ==
0032                                  'StringZ:usePetersonC = on',      # mstj(11)=3
0033                      'ParticleDecays:limitTau0 = on',  # mstj(22)=2 - decay unstable particles
0034                      'ParticleDecays:tau0Max   = 10.', # parj(71)=10.- for which ctau < 10 mm
0035                      # 'PartonLevel:ISR = on',         # on==D; mstp(61)=1 or 2 ??? Py6 D=2
0036                                  # 'PartonLevel:FSR = on',         # on==D;  mstp(71)=1
0037                      'SigmaProcess:alphaSorder=1',     # D; mstp(2)=1 (D)
0038                      'SigmaProcess:Kfactor=1.',        # D; mstp(33)=0 - no K-factor, i.e. K=1
0039                      # This is the important master switch - I can turn it ON/OFF,
0040                      # and it'll greatly affect the event multiplicity,
0041                      # without messing up too much the principle physics
0042                      'PartonLevel:MPI = off',            # D; mstp(81)=1
0043                      ### Py6 mstp(82)=4 (MI model) does NOT have anything similar in Py8,
0044                      ### because it actually means a new model, and in Py8 old models went away
0045                      ### - only new remain
0046                      #'MultipleInteractions:pT0Ref=1.8387',    # parp(82)=1.8387 - pt cutoff for MI
0047                      #'MultipleInteractions:ecmRef=1960.',     # parp(89)=1960.
0048                      #'MultipleInteractions:coreFraction=0.5', # D; parp(83)=0.5
0049                      #'MultipleInteractions:coreRadius=0.4',   # D; parp(84)=0.4
0050                      ### parp(85) & parp(86) are N/A in Py8
0051                      #'MultipleInteractions:ecmPow=0.16',      # parp(90)=0.16
0052                      'BeamRemnants:primordialKT=on',          # D; mstp(91)=1
0053                      'BeamRemnants:primordialKThard=2.1',     # parp(91)=2.1
0054                      'BeamRemnants:primordialKTremnant=15.',  # parp(93)=15.
0055                      'SpaceShower:pTmin=1.25',                # parp(62)=1.25
0056                      ### parp(64) is N/A in Py8
0057                      'SpaceShower:pTmaxFudge=2.'              # should this be instead of parp(67)=2.5 ???
0058                                                               # In Py6, D==parp(67)=4.
0059                                                               # and in Py8 for this D=1, max=2
0060     ),
0061         py8ProcessSettings = cms.vstring( 'Bottomonium:all = on' # Quarkonia (bottonium/NRQCD) - msel=62
0062 #        py8ProcessSettings = cms.vstring( 'Bottomonium:gg2QQbar[3S1(1)]g = on' # Quarkonia - msel=62, isub=???
0063     ),
0064     parameterSets = cms.vstring('py8UESettings','py8ProcessSettings')
0065     )
0066 )
0067
0068 # The following three lines reduce the clutter of repeated printouts
0069 # of the same exception message.
0070 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0071
0072 process.MessageLogger.cerr.enableStatistics = False
0073
0074
0075 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0076     generator = cms.PSet(
0077         initialSeed = cms.untracked.uint32(123456789),
0078         engineName = cms.untracked.string('HepJamesRandom')
0079     )
0080 )
0081
0082 process.maxEvents = cms.untracked.PSet(
0083     input = cms.untracked.int32(10)
0084 )
0085
0086 process.GEN = cms.OutputModule("PoolOutputModule",
0087     fileName = cms.untracked.string('Py8Onia_EvtGen.root')
0088 )
0089
0090 process.p = cms.Path(process.generator)
0091 process.outpath = cms.EndPath(process.GEN)
0092
0093 process.schedule = cms.Schedule(process.p, process.outpath)
0094