Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2021-02-14 13:29:55

0001 import FWCore.ParameterSet.Config as cms
0002 
0003 #from Configuration.GenProduction.PythiaUESettings_cfi import "
0004 from Configuration.Generator.PythiaUESettings_cfi import *
0005 
0006 process = cms.Process("TEST")
0007 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0008 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0009 #process.load("SimGeneral.HepPDTESSource.pdt_cfi")
0010 
0011 process.load("Configuration.StandardSequences.Services_cff")
0012 
0013 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0014     generator = cms.PSet(
0015         initialSeed = cms.untracked.uint32(123456789),
0016         engineName = cms.untracked.string('HepJamesRandom')
0017     )
0018 )
0019 
0020 process.randomEngineStateProducer = cms.EDProducer("RandomEngineStateProducer")
0021 
0022 # The following three lines reduce the clutter of repeated printouts
0023 # of the same exception message.
0024 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0025 
0026 process.MessageLogger.cerr.enableStatistics = False
0027 
0028 
0029 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(50))
0030 
0031 process.source = cms.Source("LHESource",
0032     fileNames = cms.untracked.vstring('file:edde.lhe')
0033 )
0034 
0035 process.generator = cms.EDFilter("Pythia6HadronizerFilter",
0036     pythiaHepMCVerbosity = cms.untracked.bool(True),
0037     maxEventsToPrint = cms.untracked.int32(0),
0038     pythiaPylistVerbosity = cms.untracked.int32(1),
0039     comEnergy = cms.double(10000.0),
0040     PythiaParameters = cms.PSet(
0041         pythiaUESettingsBlock,
0042         processParameters = cms.vstring('MSEL=0         ! User defined processes', 
0043 #=          'MSTJ(1)=1       ! Fragmentation/hadronization on or off',
0044             'MSTP(61)=1      ! Parton showering on or off',
0045 #           'MSTU(53) = 0            ! no smearing energy',
0046 #           'MSTU(54) = 3            ! massless jets',
0047             'MSTP(71) =1             ! Final-state QCD and QED radiation',
0048             'MSTP(81) =1             ! multiple interaction',
0049             'MSTP(111)=1             ! fragmentation and decay',
0050             'MSTP(122)=0             ! switch off X section print out',
0051 #...Higgs decay definition...
0052             'MDME(210,1) =0           ! h0 -> d dbar',
0053             'MDME(211,1) =0           ! h0 -> u ubar',
0054                     'MDME(212,1) =0           ! h0 -> s sbar',
0055             'MDME(213,1) =0           ! h0 -> c cbar',
0056             'MDME(214,1) =1           ! h0 -> b bbar',
0057             'MDME(215,1) =0           ! h0 -> t tbar',
0058             'MDME(216,1) =-1          ! h0 -> bprime bbar',
0059             'MDME(217,1) =-1          ! h0 -> tprime tbar',
0060             'MDME(218,1) =0           ! h0 -> e+e-',
0061             'MDME(219,1) =0           ! h0 -> mu+mu-',
0062             'MDME(220,1) =0           ! h0 -> tau+tau-',
0063             'MDME(221,1) =-1          ! h0 -> tauprime+ tauprime-',
0064             'MDME(222,1) =0           ! h0 ->  gg',
0065             'MDME(223,1) =0           ! h0-> gamma gamma',
0066             'MDME(224,1) =0           ! h0 -> gamma Z0',
0067             'MDME(225,1) =0           ! h0 -> Z0 Z0',
0068             'MDME(226,1) =0           ! h0 -> W+W-'
0069     ),
0070     # This is a vector of ParameterSet names to be read, in this order
0071     parameterSets = cms.vstring('pythiaUESettings', 'processParameters')
0072     )
0073 )
0074 
0075 process.GEN = cms.OutputModule("PoolOutputModule",
0076     fileName = cms.untracked.string('EDDE.root')
0077 )
0078 
0079 process.p = cms.Path(process.generator)
0080 process.p1 = cms.Path(process.randomEngineStateProducer)
0081 process.outpath = cms.EndPath(process.GEN)
0082 
0083 process.schedule = cms.Schedule(process.p, process.p1, process.outpath)