Back to home page

Project CMSSW displayed by LXR

 
 

    


File indexing completed on 2024-08-19 22:31:09

0001 #Default Pythia Paramters for Hydjet & Pyquen
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), ## beam/target atomic number
0008                                comEnergy = cms.double(4000.0)
0009                                )
0010 
0011 collisionParameters2760GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
0012                                    comEnergy = cms.double(2760.0)
0013                                    )
0014 
0015 collisionParameters5020GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
0016                                    comEnergy = cms.double(5020.0)
0017                                    )
0018 
0019 collisionParameters5362GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
0020                                    comEnergy = cms.double(5362.0)
0021                                    )
0022 
0023 collisionParameters5519GeV = cms.PSet(aBeamTarget = cms.double(208.0), ## beam/target atomic number
0024                                    comEnergy = cms.double(5519.0)
0025                                    )
0026 
0027 qgpParameters = cms.PSet(qgpInitialTemperature = cms.double(1.1), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
0028                          qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
0029                          hadronFreezoutTemperature = cms.double(0.125),
0030                          doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
0031                          doCollisionalEnLoss = cms.bool(True),
0032                          qgpNumQuarkFlavor = cms.int32(0),  ## number of active quark flavors in qgp; allowed values: 0,1,2,3
0033                          )
0034 
0035 qgpParameters2023 = cms.PSet(qgpInitialTemperature = cms.double(1.0), ## initial temperature of QGP; allowed range [0.2,2.0]GeV;
0036                              qgpProperTimeFormation = cms.double(0.1), ## proper time of QGP formation; allowed range [0.01,10.0]fm/c;
0037                              hadronFreezoutTemperature = cms.double(0.125),
0038                              doRadiativeEnLoss = cms.bool(True), ## if true, perform partonic radiative en loss
0039                              doCollisionalEnLoss = cms.bool(True),
0040                              qgpNumQuarkFlavor = cms.int32(0),  ## number of active quark flavors in qgp; allowed values: 0,1,2,3
0041                          )
0042 
0043 pyquenParameters  = cms.PSet(doIsospin = cms.bool(True),
0044                              angularSpectrumSelector = cms.int32(0), ## angular emitted gluon spectrum
0045                              embeddingMode = cms.int32(0),
0046                              )
0047 
0048 hydjetParameters = cms.PSet(sigmaInelNN = cms.double(70),
0049                             shadowingSwitch = cms.int32(1),
0050                             nMultiplicity = cms.int32(25000),
0051                             fracSoftMultiplicity = cms.double(1.),
0052                             maxLongitudinalRapidity = cms.double(4.5),
0053                             maxTransverseRapidity = cms.double(1.25),
0054                             rotateEventPlane = cms.bool(True),
0055                             allowEmptyEvents = cms.bool(False),
0056                             angularSpectrumSelector = cms.int32(1), ## angular emitted gluon spectrum
0057                             embeddingMode = cms.int32(0)
0058                             )
0059 
0060 hydjetParameters2023 = cms.PSet(sigmaInelNN = cms.double(70),
0061                                 shadowingSwitch = cms.int32(1),
0062                                 nMultiplicity = cms.int32(9000),
0063                                 fracSoftMultiplicity = cms.double(1.),
0064                                 maxLongitudinalRapidity = cms.double(2.3),
0065                                 maxTransverseRapidity = cms.double(1.15),
0066                                 rotateEventPlane = cms.bool(True),
0067                                 allowEmptyEvents = cms.bool(False),
0068                                 angularSpectrumSelector = cms.int32(0), ## angular emitted gluon spectrum
0069                                 embeddingMode = cms.int32(0)
0070                             )
0071 
0072 pyquenPythiaDefaultBlock = cms.PSet(
0073     pythiaUESettingsBlock,
0074     hydjetPythiaDefault = cms.vstring('MSEL=0   ! user processes',
0075                                       'CKIN(3)=10.',# ! ptMin
0076                                       'MSTP(81)=1'
0077                                       ),
0078     ppDefault = cms.vstring('MSEL=1   ! QCD hight pT processes (only jets)',
0079                             'CKIN(3)=10.',# ! ptMin
0080                             'MSTP(81)=1'
0081                             ),
0082     pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
0083                                 'MSTU(21)=1', # ! to avoid stopping run',
0084                                 'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
0085                                 'MSTP(81)=0', # ! pp multiple scattering off',
0086                                 'PMAS(5,1)=4.8', # ! b quark mass',
0087                                 'PMAS(6,1)=175.0', # ! t quark mass'
0088                                 'CKIN(3)=7.',# ! ptMin
0089                                 'MSTJ(22)=2',
0090                                 'PARJ(71)=10.', # Decays only if life time < 10mm
0091                                 'PARP(67)=1.',
0092                                 'PARP(82)=1.9',
0093                                 'PARP(85)=0.33',
0094                                 'PARP(86)=0.66',
0095                                 'PARP(89)=1000.',
0096                                 'PARP(91)=1.0',
0097                                 'MSTJ(11)=3',
0098                                 'MSTJ(22)=2'
0099                                 ),
0100     ppJets = cms.vstring('MSEL=1   ! QCD hight pT processes'),
0101     customProcesses = cms.vstring('MSEL=0   ! User processes'),
0102     pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
0103                              'MSUB(12)=1', # q+qbar->q+qbar
0104                              'MSUB(13)=1', # q+qbar->g+g
0105                              'MSUB(28)=1', # q+g->q+g
0106                              'MSUB(53)=1', # g+g->q+qbar
0107                              'MSUB(68)=1' # g+g->g+g
0108                              ),
0109     pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
0110                                       'MSUB(18)=1', # q+qbar->gamma+gamma
0111                                       'MSUB(29)=1', # q+g->q+gamma
0112                                       'MSUB(114)=1', # g+g->gamma+gamma
0113                                       'MSUB(115)=1' # g+g->g+gamma
0114                                       ),
0115 
0116     pythiaWeakBosons = cms.vstring('MSUB(1)=1',
0117                                    'MSUB(2)=1'),
0118 
0119     pythiaZjets = cms.vstring('MSUB(15)=1',
0120                               'MSUB(30)=1'),
0121 
0122     pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
0123                                         'MSUB(422) = 1',
0124                                         'MSUB(423) = 1',
0125                                         'MSUB(424) = 1',
0126                                         'MSUB(425) = 1',
0127                                         'MSUB(426) = 1',
0128                                         'MSUB(427) = 1',
0129                                         'MSUB(428) = 1',
0130                                         'MSUB(429) = 1',
0131                                         'MSUB(430) = 1',
0132                                         'MSUB(431) = 1',
0133                                         'MSUB(432) = 1',
0134                                         'MSUB(433) = 1',
0135                                         'MSUB(434) = 1',
0136                                         'MSUB(435) = 1',
0137                                         'MSUB(436) = 1',
0138                                         'MSUB(437) = 1',
0139                                         'MSUB(438) = 1',
0140                                         'MSUB(439) = 1'
0141                                         ),
0142 
0143     pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
0144                                          'MSUB(462) = 1',
0145                                          'MSUB(463) = 1',
0146                                          'MSUB(464) = 1',
0147                                          'MSUB(465) = 1',
0148                                          'MSUB(466) = 1',
0149                                          'MSUB(467) = 1',
0150                                          'MSUB(468) = 1',
0151                                          'MSUB(469) = 1',
0152                                          'MSUB(470) = 1',
0153                                          'MSUB(471) = 1',
0154                                          'MSUB(472) = 1',
0155                                          'MSUB(473) = 1',
0156                                          'MSUB(474) = 1',
0157                                          'MSUB(475) = 1',
0158                                          'MSUB(476) = 1',
0159                                          'MSUB(477) = 1',
0160                                          'MSUB(478) = 1',
0161                                          'MSUB(479) = 1',
0162                                          ),
0163 
0164     pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16',   # Matrix Elements
0165                                           'PARP(142)=0.0119',
0166                                           'PARP(143)=0.01',
0167                                           'PARP(144)=0.01',
0168                                           'PARP(145)=0.05',
0169                                           'PARP(146)=9.28',
0170                                           'PARP(147)=0.15',
0171                                           'PARP(148)=0.02',
0172                                           'PARP(149)=0.02',
0173                                           'PARP(150)=0.085',
0174                                           # Meson spin
0175                                           'PARJ(13)=0.60',
0176                                           'PARJ(14)=0.162',
0177                                           'PARJ(15)=0.018',
0178                                           'PARJ(16)=0.054',
0179                                           # Polarization
0180                                           'MSTP(145)=0',
0181                                           'MSTP(146)=0',
0182                                           'MSTP(147)=0',
0183                                           'MSTP(148)=1',
0184                                           'MSTP(149)=1',
0185                                           # Chi_c branching ratios
0186                                           'BRAT(861)=0.202',
0187                                           'BRAT(862)=0.798',
0188                                           'BRAT(1501)=0.013',
0189                                           'BRAT(1502)=0.987',
0190                                           'BRAT(1555)=0.356',
0191                                           'BRAT(1556)=0.644'
0192                                           ),
0193 
0194     pythiaZtoMuons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0195                                  "MDME(175,1)=0",          # !Z decay into u ubar,
0196                                  "MDME(176,1)=0",          # !Z decay into s sbar,
0197                                  "MDME(177,1)=0",          # !Z decay into c cbar,
0198                                  "MDME(178,1)=0",          # !Z decay into b bbar,
0199                                  "MDME(179,1)=0",          # !Z decay into t tbar,
0200                                  "MDME(182,1)=0",          # !Z decay into e- e+,
0201                                  "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0202                                  "MDME(184,1)=1",          # !Z decay into mu- mu+,
0203                                  "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0204                                  "MDME(186,1)=0",          # !Z decay into tau- tau+,
0205                                  "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0206                                  ),
0207 
0208 
0209     pythiaZtoElectrons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0210                                      "MDME(175,1)=0",          # !Z decay into u ubar,
0211                                      "MDME(176,1)=0",          # !Z decay into s sbar,
0212                                      "MDME(177,1)=0",          # !Z decay into c cbar,
0213                                      "MDME(178,1)=0",          # !Z decay into b bbar,
0214                                      "MDME(179,1)=0",          # !Z decay into t tbar,
0215                                      "MDME(182,1)=1",          # !Z decay into e- e+,
0216                                      "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0217                                      "MDME(184,1)=0",          # !Z decay into mu- mu+,
0218                                      "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0219                                      "MDME(186,1)=0",          # !Z decay into tau- tau+,
0220                                      "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0221                                      ),
0222 
0223     pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0224                                              "MDME(175,1)=0",          # !Z decay into u ubar,
0225                                              "MDME(176,1)=0",          # !Z decay into s sbar,
0226                                              "MDME(177,1)=0",          # !Z decay into c cbar,
0227                                              "MDME(178,1)=0",          # !Z decay into b bbar,
0228                                              "MDME(179,1)=0",          # !Z decay into t tbar,
0229                                              "MDME(182,1)=1",          # !Z decay into e- e+,
0230                                              "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0231                                              "MDME(184,1)=1",          # !Z decay into mu- mu+,
0232                                              "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0233                                              "MDME(186,1)=0",          # !Z decay into tau- tau+,
0234                                              "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0235                                              ),
0236 
0237     pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ',  # switch off',
0238                                        'BRAT(1035) = 1 ',  # switch on',
0239                                        'BRAT(1036) = 0 ',  # switch off',
0240                                        'BRAT(1037) = 0 ',  # switch off',
0241                                        'BRAT(1038) = 0 ',  # switch off',
0242                                        'BRAT(1039) = 0 ',  # switch off',
0243                                        'BRAT(1040) = 0 ',  # switch off',
0244                                        'BRAT(1041) = 0 ',  # switch off',
0245                                        'BRAT(1042) = 0 ',  # switch off',
0246                                        'MDME(1034,1) = 0 ',  # switch off',
0247                                        'MDME(1035,1) = 1 ',  # switch on',
0248                                        'MDME(1036,1) = 0 ',  # switch off',
0249                                        'MDME(1037,1) = 0 ',  # switch off',
0250                                        'MDME(1038,1) = 0 ',  # switch off',
0251                                        'MDME(1039,1) = 0 ',  # switch off',
0252                                        'MDME(1040,1) = 0 ',  # switch off',
0253                                        'MDME(1041,1) = 0 ',  # switch off',
0254                                        'MDME(1042,1) = 0 ',  # switch off'
0255                                        ),
0256 
0257    pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ',  # switch off',
0258                                    'BRAT(859) = 1 ',  # switch on',
0259                                    'BRAT(860) = 0 ',  # switch off',
0260                                    'MDME(858,1) = 0 ',  # switch off',
0261                                    'MDME(859,1) = 1 ',  # switch on',
0262                                    'MDME(860,1) = 0 ',  # switch off'
0263                                    ),
0264 
0265     pythiaMuonCandidates = cms.vstring(
0266     'CKIN(3)=20',
0267     'MSTJ(22)=2',
0268     'PARJ(71)=40.'
0269     ),
0270     myParameters = cms.vstring('MDCY(310,1)=0')
0271 
0272 )