Generator/python/Hydjet_Quenched_MinBias_XeXe_5442GeV_cfi.py

0001 import FWCore.ParameterSet.Config as cms
0002
0003 from Configuration.Generator.Pyquen2015Settings_cff import *
0004 from Configuration.Generator.PythiaUESettings_cfi import *
0005
0006 generator = cms.EDFilter("HydjetGeneratorFilter",
0007                          aBeamTarget = cms.double(129.0), ## beam/target atomic number
0008                          comEnergy = cms.double(5442.0),
0009                          qgpInitialTemperature = cms.double(1.), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
0010                          qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
0011                          hadronFreezoutTemperature = cms.double(0.125),
0012                          doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
0013                          doCollisionalEnLoss = cms.bool(True),
0014                          qgpNumQuarkFlavor = cms.int32(0),  ## number of active quark flavors in qgp; allowed values: 0,1,2,3
0015                          sigmaInelNN = cms.double(70),
0016                          shadowingSwitch = cms.int32(1),
0017                          nMultiplicity = cms.int32(18545),
0018                          fracSoftMultiplicity = cms.double(1.),
0019                          maxLongitudinalRapidity = cms.double(3.75),
0020                          maxTransverseRapidity = cms.double(1.16),
0021                          angularSpectrumSelector = cms.int32(1),
0022                          rotateEventPlane = cms.bool(True),
0023                          allowEmptyEvents = cms.bool(False),
0024                          embeddingMode = cms.int32(0),
0025                          hydjetMode = cms.string('kHydroQJets'),
0026
0027                          PythiaParameters = cms.PSet(
0028                             pythiaUESettingsBlock,
0029                             hydjetPythiaDefault = cms.vstring('MSEL=0   ! user processes',
0030                                                               'CKIN(3)=9.2',# ! ptMin
0031                                                               'MSTP(81)=1'
0032                                                              ),
0033                             myParameters = cms.vstring('MDCY(310,1)=0'),
0034                             pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
0035                                                      'MSUB(12)=1', # q+qbar->q+qbar
0036                                                      'MSUB(13)=1', # q+qbar->g+g
0037                                                      'MSUB(28)=1', # q+g->q+g
0038                                                      'MSUB(53)=1', # g+g->q+qbar
0039                                                      'MSUB(68)=1' # g+g->g+g
0040                                                      ),
0041                             pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
0042                                                               'MSUB(18)=1', # q+qbar->gamma+gamma
0043                                                               'MSUB(29)=1', # q+g->q+gamma
0044                                                               'MSUB(114)=1', # g+g->gamma+gamma
0045                                                               'MSUB(115)=1' # g+g->g+gamma
0046                                                               ),
0047                             pythiaZjets = cms.vstring('MSUB(15)=1',
0048                                                       'MSUB(30)=1'),
0049                             pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
0050                                                                 'MSUB(422) = 1',
0051                                                                 'MSUB(423) = 1',
0052                                                                 'MSUB(424) = 1',
0053                                                                 'MSUB(425) = 1',
0054                                                                 'MSUB(426) = 1',
0055                                                                 'MSUB(427) = 1',
0056                                                                 'MSUB(428) = 1',
0057                                                                 'MSUB(429) = 1',
0058                                                                 'MSUB(430) = 1',
0059                                                                 'MSUB(431) = 1',
0060                                                                 'MSUB(432) = 1',
0061                                                                 'MSUB(433) = 1',
0062                                                                 'MSUB(434) = 1',
0063                                                                 'MSUB(435) = 1',
0064                                                                 'MSUB(436) = 1',
0065                                                                 'MSUB(437) = 1',
0066                                                                 'MSUB(438) = 1',
0067                                                                 'MSUB(439) = 1'
0068                                                                 ),
0069
0070                             pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
0071                                                                  'MSUB(462) = 1',
0072                                                                  'MSUB(463) = 1',
0073                                                                  'MSUB(464) = 1',
0074                                                                  'MSUB(465) = 1',
0075                                                                  'MSUB(466) = 1',
0076                                                                  'MSUB(467) = 1',
0077                                                                  'MSUB(468) = 1',
0078                                                                  'MSUB(469) = 1',
0079                                                                  'MSUB(470) = 1',
0080                                                                  'MSUB(471) = 1',
0081                                                                  'MSUB(472) = 1',
0082                                                                  'MSUB(473) = 1',
0083                                                                  'MSUB(474) = 1',
0084                                                                  'MSUB(475) = 1',
0085                                                                  'MSUB(476) = 1',
0086                                                                  'MSUB(477) = 1',
0087                                                                  'MSUB(478) = 1',
0088                                                                  'MSUB(479) = 1',
0089                                                                  ),
0090                             pythiaWeakBosons = cms.vstring('MSUB(1)=1',
0091                                                            'MSUB(2)=1'),
0092                             pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16',   # Matrix Elements
0093                                                                   'PARP(142)=0.0119',
0094                                                                   'PARP(143)=0.01',
0095                                                                   'PARP(144)=0.01',
0096                                                                   'PARP(145)=0.05',
0097                                                                   'PARP(146)=9.28',
0098                                                                   'PARP(147)=0.15',
0099                                                                   'PARP(148)=0.02',
0100                                                                   'PARP(149)=0.02',
0101                                                                   'PARP(150)=0.085',
0102                                                                   # Meson spin
0103                                                                   'PARJ(13)=0.60',
0104                                                                   'PARJ(14)=0.162',
0105                                                                   'PARJ(15)=0.018',
0106                                                                   'PARJ(16)=0.054',
0107                                                                   # Polarization
0108                                                                   'MSTP(145)=0',
0109                                                                   'MSTP(146)=0',
0110                                                                   'MSTP(147)=0',
0111                                                                   'MSTP(148)=1',
0112                                                                   'MSTP(149)=1',
0113                                                                   # Chi_c branching ratios
0114                                                                   'BRAT(861)=0.202',
0115                                                                   'BRAT(862)=0.798',
0116                                                                   'BRAT(1501)=0.013',
0117                                                                   'BRAT(1502)=0.987',
0118                                                                   'BRAT(1555)=0.356',
0119                                                                   'BRAT(1556)=0.644'
0120                                                                   ),
0121                             # Quarkonia and Weak Bosons added back upon dilepton group's request.
0122                             parameterSets = cms.vstring('pythiaUESettings',
0123                                                         'hydjetPythiaDefault',
0124                                                         'myParameters',
0125                                                         'pythiaJets',
0126                                                         'pythiaPromptPhotons',
0127                                                         'pythiaZjets',
0128                                                         'pythiaBottomoniumNRQCD',
0129                                                         'pythiaCharmoniumNRQCD',
0130                                                         'pythiaQuarkoniaSettings',
0131                                                         'pythiaWeakBosons'
0132                                                         )
0133                          ),
0134                          cFlag = cms.int32(1),
0135                          bMin = cms.double(0),
0136                          bMax = cms.double(30),
0137                          bFixed = cms.double(0)
0138                         )