File indexing completed on 2023-03-17 11:04:55
0001 import FWCore.ParameterSet.Config as cms
0002
0003
0004
0005
0006
0007
0008
0009
0010
0011 process = cms.Process("TEST")
0012 process.load("FWCore.Framework.test.cmsExceptionsFatal_cff")
0013 process.load("SimGeneral.HepPDTESSource.pythiapdt_cfi")
0014
0015
0016 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0017 generator = cms.PSet(
0018 initialSeed = cms.untracked.uint32(123456789),
0019 engineName = cms.untracked.string('HepJamesRandom')
0020 )
0021 )
0022
0023
0024
0025
0026 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0027
0028
0029
0030
0031 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(5))
0032
0033 process.source = cms.Source("EmptySource")
0034
0035 from Configuration.Generator.PythiaUESettings_cfi import *
0036 from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
0037
0038 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0039 pythiaHepMCVerbosity = cms.untracked.bool(False),
0040 maxEventsToPrint = cms.untracked.int32(3),
0041 pythiaPylistVerbosity = cms.untracked.int32(1),
0042
0043
0044
0045 comEnergy = cms.double(10000.0),
0046
0047 ExternalDecays = cms.PSet(
0048
0049 Tauola = cms.untracked.PSet(
0050
0051
0052
0053
0054
0055 UseTauolaPolarization = cms.bool(True),
0056 InputCards = cms.PSet
0057 (
0058 pjak1 = cms.int32(0),
0059 pjak2 = cms.int32(0),
0060 mdtau = cms.int32(202)
0061 )
0062
0063
0064 ),
0065
0066 EvtGen = cms.untracked.PSet(
0067 operates_on_particles = cms.vint32( 0 ),
0068
0069
0070 use_internal_pythia = cms.untracked.bool(False),
0071 use_default_decay = cms.untracked.bool(True),
0072 decay_table = cms.FileInPath('GeneratorInterface/ExternalDecays/data/DECAY_NOLONGLIFE.DEC'),
0073
0074 particle_property_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/evt.pdl'),
0075 user_decay_file = cms.FileInPath('GeneratorInterface/ExternalDecays/data/Onia_mumu.dec'),
0076 list_forced_decays = cms.vstring(),
0077 ),
0078
0079
0080
0081 parameterSets = cms.vstring( 'Tauola',
0082 'EvtGen'
0083 )
0084 ),
0085
0086 PythiaParameters = cms.PSet(
0087
0088 pythiaGeneralSettings = cms.vstring(
0089 'PMAS(25,1)=190.0 !mass of Higgs',
0090 'MSTP(128)=0 !dec.prods out of doc section, point at parents in the main section',
0091 'MSEL=0 !(D=1) to select between full user control (0, then use MSUB) and some preprogrammed alternative: QCD hight pT processes (1, then ISUB=11, 12, 13, 28, 53, 68), QCD low pT processes (2, then ISUB=11, 12, 13, 28, 53, 68, 91, 92, 94, 95)',
0092 'MSTJ(11)=3 !Choice of the fragmentation function',
0093
0094 'MSTP(51)=7 !structure function chosen',
0095
0096
0097
0098
0099 'MSTU(21)=1 !Check on possible errors during program execution',
0100 'MSUB(102)=1 !ggH',
0101 'MSUB(123)=1 !ZZ fusion to H',
0102 'MSUB(124)=1 !WW fusion to H',
0103
0104 'PARP(82)=1.9 !pt cutoff for multiparton interactions',
0105 'PARP(83)=0.5 !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch',
0106 'PARP(84)=0.4 !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch',
0107 'PARP(90)=0.16 !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch',
0108
0109 'CKIN(45)=5. !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch',
0110 'CKIN(46)=150. !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0111 'CKIN(47)=5. !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0112 'CKIN(48)=150. !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch'
0113 ),
0114 pythiaHZZDecays = cms.vstring(
0115 'MDME(174,1)=0 !Z decay into d dbar',
0116 'MDME(175,1)=0 !Z decay into u ubar',
0117 'MDME(176,1)=0 !Z decay into s sbar',
0118 'MDME(177,1)=0 !Z decay into c cbar',
0119
0120 'MDME(178,1)=4 !Z decay into b bbar',
0121
0122 'MDME(179,1)=0 !Z decay into t tbar',
0123 'MDME(182,1)=0 !Z decay into e- e+',
0124 'MDME(183,1)=0 !Z decay into nu_e nu_ebar',
0125 'MDME(184,1)=0 !Z decay into mu- mu+',
0126 'MDME(185,1)=0 !Z decay into nu_mu nu_mubar',
0127
0128 'MDME(186,1)=5 !Z decay into tau- tau+',
0129
0130 'MDME(187,1)=0 !Z decay into nu_tau nu_taubar',
0131 'MDME(210,1)=0 !Higgs decay into dd',
0132 'MDME(211,1)=0 !Higgs decay into uu',
0133 'MDME(212,1)=0 !Higgs decay into ss',
0134 'MDME(213,1)=0 !Higgs decay into cc',
0135 'MDME(214,1)=0 !Higgs decay into bb',
0136 'MDME(215,1)=0 !Higgs decay into tt',
0137 'MDME(216,1)=0 !Higgs decay into',
0138 'MDME(217,1)=0 !Higgs decay into Higgs decay',
0139 'MDME(218,1)=0 !Higgs decay into e nu e',
0140 'MDME(219,1)=0 !Higgs decay into mu nu mu',
0141 'MDME(220,1)=0 !Higgs decay into tau nu tau',
0142 'MDME(221,1)=0 !Higgs decay into Higgs decay',
0143 'MDME(222,1)=0 !Higgs decay into g g',
0144 'MDME(223,1)=0 !Higgs decay into gam gam',
0145 'MDME(224,1)=0 !Higgs decay into gam Z',
0146 'MDME(225,1)=1 !Higgs decay into Z Z',
0147 'MDME(226,1)=0 !Higgs decay into W W'
0148 ),
0149
0150 parameterSets = cms.vstring(
0151
0152 'pythiaGeneralSettings',
0153 'pythiaHZZDecays')
0154 )
0155 )
0156
0157 process.GEN = cms.OutputModule("PoolOutputModule",
0158 fileName = cms.untracked.string('TestHZZ4tau.root')
0159 )
0160
0161 process.p = cms.Path(process.generator)
0162 process.outpath = cms.EndPath(process.GEN)
0163
0164 process.schedule = cms.Schedule(process.p, process.outpath)