File indexing completed on 2024-04-06 12:13:55
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
0007
0008 process.RandomNumberGeneratorService = cms.Service("RandomNumberGeneratorService",
0009 generator = cms.PSet(
0010 initialSeed = cms.untracked.uint32(123456789),
0011 engineName = cms.untracked.string('HepJamesRandom')
0012 )
0013 )
0014
0015
0016
0017
0018 process.load("FWCore.MessageLogger.MessageLogger_cfi")
0019
0020
0021
0022
0023 process.maxEvents = cms.untracked.PSet(input = cms.untracked.int32(20))
0024
0025
0026 process.source = cms.Source("EmptySource",
0027 firstLuminosityBlock = cms.untracked.uint32(1),
0028 numberEventsInLuminosityBlock = cms.untracked.uint32(5),
0029 firstEvent = cms.untracked.uint32(1),
0030 firstRun = cms.untracked.uint32(1),
0031 numberEventsInRun = cms.untracked.uint32(10)
0032 )
0033
0034 from Configuration.Generator.PythiaUESettings_cfi import *
0035 from GeneratorInterface.ExternalDecays.TauolaSettings_cff import *
0036
0037 process.generator = cms.EDFilter("Pythia6GeneratorFilter",
0038 pythiaHepMCVerbosity = cms.untracked.bool(True),
0039 maxEventsToPrint = cms.untracked.int32(3),
0040 pythiaPylistVerbosity = cms.untracked.int32(1),
0041
0042
0043 displayPythiaCards = cms.untracked.bool(True),
0044 comEnergy = cms.double(10000.0),
0045 ExternalDecays = cms.PSet(
0046 Tauola = cms.untracked.PSet(
0047
0048
0049
0050
0051
0052 UseTauolaPolarization = cms.bool(True),
0053 InputCards = cms.PSet
0054 (
0055 pjak1 = cms.int32(0),
0056 pjak2 = cms.int32(0),
0057 mdtau = cms.int32(240)
0058 )
0059
0060
0061 ),
0062 parameterSets = cms.vstring('Tauola')
0063 ),
0064 PythiaParameters = cms.PSet(
0065
0066 pythiaHZZ4tau = cms.vstring('PMAS(25,1)=190.0 !mass of Higgs',
0067 '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)',
0068 'MSTJ(11)=3 !Choice of the fragmentation function',
0069 'MSTJ(41)=3 !Switch off Pythia QED bremsshtrahlung off leptons (only)',
0070 'MSTP(51)=7 !structure function chosen',
0071 'MSTP(61)=0 ! no initial-state showers',
0072 'MSTP(71)=0 ! no final-state showers',
0073 'MSTP(81)=0 ! no multiple interactions',
0074 'MSTP(111)=0 ! no hadronization',
0075 'MSTU(21)=1 !Check on possible errors during program execution',
0076 'MSUB(102)=1 !ggH',
0077 'MSUB(123)=1 !ZZ fusion to H',
0078 'MSUB(124)=1 !WW fusion to H',
0079 'PARP(82)=1.9 !pt cutoff for multiparton interactions',
0080 'PARP(83)=0.5 !Multiple interactions: matter distrbn parameter Registered by Chris.Seez@cern.ch',
0081 'PARP(84)=0.4 !Multiple interactions: matter distribution parameter Registered by Chris.Seez@cern.ch',
0082 'PARP(90)=0.16 !Multiple interactions: rescaling power Registered by Chris.Seez@cern.ch',
0083 'CKIN(45)=5. !high mass cut on m2 in 2 to 2 process Registered by Chris.Seez@cern.ch',
0084 'CKIN(46)=150. !high mass cut on secondary resonance m1 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0085 'CKIN(47)=5. !low mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0086 'CKIN(48)=150. !high mass cut on secondary resonance m2 in 2->1->2 process Registered by Alexandre.Nikitenko@cern.ch',
0087 'MDME(174,1)=0 !Z decay into d dbar',
0088 'MDME(175,1)=0 !Z decay into u ubar',
0089 'MDME(176,1)=0 !Z decay into s sbar',
0090 'MDME(177,1)=0 !Z decay into c cbar',
0091 'MDME(178,1)=0 !Z decay into b bbar',
0092 'MDME(179,1)=0 !Z decay into t tbar',
0093 'MDME(182,1)=0 !Z decay into e- e+',
0094 'MDME(183,1)=0 !Z decay into nu_e nu_ebar',
0095 'MDME(184,1)=0 !Z decay into mu- mu+',
0096 'MDME(185,1)=0 !Z decay into nu_mu nu_mubar',
0097 'MDME(186,1)=1 !Z decay into tau- tau+',
0098 'MDME(187,1)=0 !Z decay into nu_tau nu_taubar',
0099 'MDME(210,1)=0 !Higgs decay into dd',
0100 'MDME(211,1)=0 !Higgs decay into uu',
0101 'MDME(212,1)=0 !Higgs decay into ss',
0102 'MDME(213,1)=0 !Higgs decay into cc',
0103 'MDME(214,1)=0 !Higgs decay into bb',
0104 'MDME(215,1)=0 !Higgs decay into tt',
0105 'MDME(216,1)=0 !Higgs decay into',
0106 'MDME(217,1)=0 !Higgs decay into Higgs decay',
0107 'MDME(218,1)=0 !Higgs decay into e nu e',
0108 'MDME(219,1)=0 !Higgs decay into mu nu mu',
0109 'MDME(220,1)=0 !Higgs decay into tau nu tau',
0110 'MDME(221,1)=0 !Higgs decay into Higgs decay',
0111 'MDME(222,1)=0 !Higgs decay into g g',
0112 'MDME(223,1)=0 !Higgs decay into gam gam',
0113 'MDME(224,1)=0 !Higgs decay into gam Z',
0114 'MDME(225,1)=1 !Higgs decay into Z Z',
0115 'MDME(226,1)=0 !Higgs decay into W W',
0116 'MSTP(128)=0 !dec.prods out of doc section, point at parents in the main section'),
0117
0118
0119 parameterSets = cms.vstring(
0120
0121 'pythiaHZZ4tau')
0122 )
0123 )
0124
0125 process.GEN = cms.OutputModule("PoolOutputModule",
0126 fileName = cms.untracked.string('Py6_HZZ4tau.root')
0127 )
0128
0129
0130
0131
0132
0133
0134
0135 process.p = cms.Path(process.generator)
0136
0137 process.outpath = cms.EndPath(process.GEN)
0138
0139
0140 process.schedule = cms.Schedule(process.p, process.outpath)