File indexing completed on 2024-04-06 12:03:36
0001
0002
0003 import FWCore.ParameterSet.Config as cms
0004
0005 from Configuration.Generator.PythiaUESettings_cfi import *
0006
0007 collisionParameters4TeV = cms.PSet(aBeamTarget = cms.double(208.0),
0008 comEnergy = cms.double(4000.0)
0009 )
0010
0011 collisionParameters2760GeV = cms.PSet(aBeamTarget = cms.double(208.0),
0012 comEnergy = cms.double(2760.0)
0013 )
0014
0015 collisionParameters5362GeV = cms.PSet(aBeamTarget = cms.double(208.0),
0016 comEnergy = cms.double(5362.0)
0017 )
0018
0019 collisionParameters = collisionParameters2760GeV.clone()
0020
0021 qgpParameters = cms.PSet(qgpInitialTemperature = cms.double(1.0),
0022 qgpProperTimeFormation = cms.double(0.1),
0023 hadronFreezoutTemperature = cms.double(0.14),
0024 doRadiativeEnLoss = cms.bool(True),
0025 doCollisionalEnLoss = cms.bool(False),
0026 qgpNumQuarkFlavor = cms.int32(0),
0027 )
0028
0029 pyquenParameters = cms.PSet(doIsospin = cms.bool(True),
0030 angularSpectrumSelector = cms.int32(0),
0031 embeddingMode = cms.int32(0),
0032 )
0033
0034 hydjetParameters = cms.PSet(sigmaInelNN = cms.double(58),
0035 shadowingSwitch = cms.int32(0),
0036 nMultiplicity = cms.int32(21500),
0037 fracSoftMultiplicity = cms.double(1.),
0038 maxLongitudinalRapidity = cms.double(4.5),
0039 maxTransverseRapidity = cms.double(1.),
0040 rotateEventPlane = cms.bool(True),
0041 allowEmptyEvents = cms.bool(False),
0042 embeddingMode = cms.int32(0)
0043 )
0044
0045 pyquenPythiaDefaultBlock = cms.PSet(
0046 pythiaUESettingsBlock,
0047 hydjetPythiaDefault = cms.vstring('MSEL=0 ! user processes',
0048 'CKIN(3)=6.',
0049 'MSTP(81)=0'
0050 ),
0051 ppDefault = cms.vstring('MSEL=1 ! QCD hight pT processes (only jets)',
0052 'CKIN(3)=6.',
0053 'MSTP(81)=0'
0054 ),
0055 pythiaHirootDefault = cms.vstring('MSEL=0',
0056 'MSTU(21)=1',
0057 'PARU(14)=1.',
0058 'MSTP(81)=0',
0059 'PMAS(5,1)=4.8',
0060 'PMAS(6,1)=175.0',
0061 'CKIN(3)=7.',
0062 'MSTJ(22)=2',
0063 'PARJ(71)=10.',
0064 'PARP(67)=1.',
0065 'PARP(82)=1.9',
0066 'PARP(85)=0.33',
0067 'PARP(86)=0.66',
0068 'PARP(89)=1000.',
0069 'PARP(91)=1.0',
0070 'MSTJ(11)=3',
0071 'MSTJ(22)=2'
0072 ),
0073 ppJets = cms.vstring('MSEL=1 ! QCD hight pT processes'),
0074 customProcesses = cms.vstring('MSEL=0 ! User processes'),
0075 pythiaJets = cms.vstring('MSUB(11)=1',
0076 'MSUB(12)=1',
0077 'MSUB(13)=1',
0078 'MSUB(28)=1',
0079 'MSUB(53)=1',
0080 'MSUB(68)=1'
0081 ),
0082 pythiaPromptPhotons = cms.vstring('MSUB(14)=1',
0083 'MSUB(18)=1',
0084 'MSUB(29)=1',
0085 'MSUB(114)=1',
0086 'MSUB(115)=1'
0087 ),
0088
0089 pythiaWeakBosons = cms.vstring('MSUB(1)=1',
0090 'MSUB(2)=1'),
0091
0092 pythiaZjets = cms.vstring('MSUB(15)=1',
0093 'MSUB(30)=1'),
0094
0095 pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
0096 'MSUB(422) = 1',
0097 'MSUB(423) = 1',
0098 'MSUB(424) = 1',
0099 'MSUB(425) = 1',
0100 'MSUB(426) = 1',
0101 'MSUB(427) = 1',
0102 'MSUB(428) = 1',
0103 'MSUB(429) = 1',
0104 'MSUB(430) = 1',
0105 'MSUB(431) = 1',
0106 'MSUB(432) = 1',
0107 'MSUB(433) = 1',
0108 'MSUB(434) = 1',
0109 'MSUB(435) = 1',
0110 'MSUB(436) = 1',
0111 'MSUB(437) = 1',
0112 'MSUB(438) = 1',
0113 'MSUB(439) = 1'
0114 ),
0115
0116 pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
0117 'MSUB(462) = 1',
0118 'MSUB(463) = 1',
0119 'MSUB(464) = 1',
0120 'MSUB(465) = 1',
0121 'MSUB(466) = 1',
0122 'MSUB(467) = 1',
0123 'MSUB(468) = 1',
0124 'MSUB(469) = 1',
0125 'MSUB(470) = 1',
0126 'MSUB(471) = 1',
0127 'MSUB(472) = 1',
0128 'MSUB(473) = 1',
0129 'MSUB(474) = 1',
0130 'MSUB(475) = 1',
0131 'MSUB(476) = 1',
0132 'MSUB(477) = 1',
0133 'MSUB(478) = 1',
0134 'MSUB(479) = 1',
0135 ),
0136
0137 pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16',
0138 'PARP(142)=0.0119',
0139 'PARP(143)=0.01',
0140 'PARP(144)=0.01',
0141 'PARP(145)=0.05',
0142 'PARP(146)=9.28',
0143 'PARP(147)=0.15',
0144 'PARP(148)=0.02',
0145 'PARP(149)=0.02',
0146 'PARP(150)=0.085',
0147
0148 'PARJ(13)=0.60',
0149 'PARJ(14)=0.162',
0150 'PARJ(15)=0.018',
0151 'PARJ(16)=0.054',
0152
0153 'MSTP(145)=0',
0154 'MSTP(146)=0',
0155 'MSTP(147)=0',
0156 'MSTP(148)=1',
0157 'MSTP(149)=1',
0158
0159 'BRAT(861)=0.202',
0160 'BRAT(862)=0.798',
0161 'BRAT(1501)=0.013',
0162 'BRAT(1502)=0.987',
0163 'BRAT(1555)=0.356',
0164 'BRAT(1556)=0.644'
0165 ),
0166
0167 pythiaZtoMuons = cms.vstring("MDME(174,1)=0",
0168 "MDME(175,1)=0",
0169 "MDME(176,1)=0",
0170 "MDME(177,1)=0",
0171 "MDME(178,1)=0",
0172 "MDME(179,1)=0",
0173 "MDME(182,1)=0",
0174 "MDME(183,1)=0",
0175 "MDME(184,1)=1",
0176 "MDME(185,1)=0",
0177 "MDME(186,1)=0",
0178 "MDME(187,1)=0"
0179 ),
0180
0181
0182 pythiaZtoElectrons = cms.vstring("MDME(174,1)=0",
0183 "MDME(175,1)=0",
0184 "MDME(176,1)=0",
0185 "MDME(177,1)=0",
0186 "MDME(178,1)=0",
0187 "MDME(179,1)=0",
0188 "MDME(182,1)=1",
0189 "MDME(183,1)=0",
0190 "MDME(184,1)=0",
0191 "MDME(185,1)=0",
0192 "MDME(186,1)=0",
0193 "MDME(187,1)=0"
0194 ),
0195
0196 pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0",
0197 "MDME(175,1)=0",
0198 "MDME(176,1)=0",
0199 "MDME(177,1)=0",
0200 "MDME(178,1)=0",
0201 "MDME(179,1)=0",
0202 "MDME(182,1)=1",
0203 "MDME(183,1)=0",
0204 "MDME(184,1)=1",
0205 "MDME(185,1)=0",
0206 "MDME(186,1)=0",
0207 "MDME(187,1)=0"
0208 ),
0209
0210 pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ',
0211 'BRAT(1035) = 1 ',
0212 'BRAT(1036) = 0 ',
0213 'BRAT(1037) = 0 ',
0214 'BRAT(1038) = 0 ',
0215 'BRAT(1039) = 0 ',
0216 'BRAT(1040) = 0 ',
0217 'BRAT(1041) = 0 ',
0218 'BRAT(1042) = 0 ',
0219 'MDME(1034,1) = 0 ',
0220 'MDME(1035,1) = 1 ',
0221 'MDME(1036,1) = 0 ',
0222 'MDME(1037,1) = 0 ',
0223 'MDME(1038,1) = 0 ',
0224 'MDME(1039,1) = 0 ',
0225 'MDME(1040,1) = 0 ',
0226 'MDME(1041,1) = 0 ',
0227 'MDME(1042,1) = 0 ',
0228 ),
0229
0230 pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ',
0231 'BRAT(859) = 1 ',
0232 'BRAT(860) = 0 ',
0233 'MDME(858,1) = 0 ',
0234 'MDME(859,1) = 1 ',
0235 'MDME(860,1) = 0 ',
0236 ),
0237
0238 pythiaMuonCandidates = cms.vstring(
0239 'CKIN(3)=20',
0240 'MSTJ(22)=2',
0241 'PARJ(71)=40.'
0242 ),
0243 myParameters = cms.vstring(
0244 )
0245
0246 )
0247
0248
0249 impactParameters = cms.PSet(cFlag = cms.int32(1),
0250 bFixed = cms.double(0),
0251 bMin = cms.double(0),
0252 bMax = cms.double(30)
0253 )
0254