|
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
|
import FWCore.ParameterSet.Config as cms
process = cms.Process("TEST")
process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
process.source = cms.Source("EmptySource")
process.generator = cms.EDFilter("Pythia8GeneratorFilter",
maxEventsToPrint = cms.untracked.int32(0),
pythiaPylistVerbosity = cms.untracked.int32(1),
filterEfficiency = cms.untracked.double(1.0),
pythiaHepMCVerbosity = cms.untracked.bool(True),
comEnergy = cms.double(7000.),
ExternalDecays = cms.PSet(
Tauola = cms.untracked.PSet(
UseTauolaPolarization = cms.bool(True),
InputCards = cms.PSet
(
pjak1 = cms.int32(1), # electron mode
pjak2 = cms.int32(2), # muon mode
mdtau = cms.int32(0) # does NOT matter at this point...
)
),
parameterSets = cms.vstring('Tauola')
),
PythiaParameters = cms.PSet(
py8GeneralSettings = cms.vstring( 'Main:timesAllowErrors = 10000',
'ParticleDecays:limitTau0 = on',
'ParticleDecays:tauMax = 10'
),
py8UESettings = cms.vstring( 'Tune:pp 2',
'Tune:ee 3'
),
py8SpecialSettings = cms.vstring('HadronLevel:Hadronize = off', # mstp(111)=1
#' ParticleDecays:limitTau0 = on', # mstj(22)=2 - decay unstable particles
# 'ParticleDecays:tau0Max = 10.', # parj(71)=10.- for which ctau < 10 mm
'PartonLevel:ISR = off',
'PartonLevel:FSR = off',
'PartonLevel:MPI = off' # mstp(81)=0
# for pp intractions, SpaceShower things are not important,
# but TimeShower settings are relevat
'TimeShower:QEDshowerByL = off',
'TimeShower:QEDshowerByQ = off',
'SpaceShower:QEDshowerByL = off',
'SpaceShower:QEDshowerByQ = off',
'PhaseSpace:pTHatMin = 5.',
'PhaseSpace:pTHatMax = 150.',
'StringZ:usePetersonB = on', # these 2 together ==
'StringZ:usePetersonC = on' # mstj(11)=3
),
py8ProcessSettings = cms.vstring( # 'HiggsSM:gg2H = on',
# 'HiggsSM:ff2Hff(t:ZZ) = on',
# 'HiggsSM:ff2Hff(t:WW) = on'
# 'WeakDoubleBoson:ffbar2ZW = on'
'WeakSingleBoson:ffbar2gmZ = on' # msel=11
),
py8ZDecaySettings = cms.vstring( # '24:onMode = off', '24:onIfAny = 15',
# '25:onMode = off', '25:onIfAny = 23',
'23:onMode = off', # turn OFF all Z decays
# '23:onIfAny = 13', # turn ON Z->mumu
'23:onIfAny = 15' # turn ON Z->tautau
),
parameterSets = cms.vstring('py8GeneralSettings',
'py8UESettings',
'py8ProcessSettings',
'py8ZDecaySettings'
)
)
)
# The following three lines reduce the clutter of repeated printouts
# of the same exception message.
process.load("FWCore.MessageLogger.MessageLogger_cfi")
process.MessageLogger.cerr.enableStatistics = False
process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
generator = cms.PSet(
initialSeed = cms.untracked.uint32(123456789),
engineName = cms.untracked.string('HepJamesRandom')
)
)
process.maxEvents = cms.untracked.PSet(
input = cms.untracked.int32(20)
)
process.GEN = cms.OutputModule("PoolOutputModule",
fileName = cms.untracked.string('Py8_Z2tautau.root')
)
process.p = cms.Path(process.generator)
process.outpath = cms.EndPath(process.GEN)
process.schedule = cms.Schedule(process.p, process.outpath)
#process.schedule = cms.Schedule(process.p, process.p1)
|