PyquenInterface/python/pyquenPythiaDefault_cff.py

0001 #Default Pythia Paramters for Hydjet & Pyquen
0002
0003 import FWCore.ParameterSet.Config as cms
0004
0005 from Configuration.Generator.PythiaUESettings_cfi import *
0006 pyquenPythiaDefaultBlock = cms.PSet(
0007     pythiaUESettingsBlock,
0008     pyquenMain = cms.vstring(
0009                             'MSTP(122)=0', # ! no printout of Pythia initialization information hereinafter
0010                             'MSTP(128)=2'  # ! instructs pythia not to put multiple copies of resonances in the event record
0011                             ),
0012
0013     ppDefault = cms.vstring('MSEL=1   ! QCD hight pT processes',
0014                             'CKIN(3)=7.',# ! ptMin
0015                             'MSTP(81)=0'
0016                             ),
0017     pythiaHirootDefault = cms.vstring('MSEL=0', # ! Only user defined processes,
0018                                       'MSTU(21)=1', # ! to avoid stopping run',
0019                                       'PARU(14)=1.', # ! tolerance parameter to adjust fragmentation',
0020                                       'MSTP(81)=0', # ! pp multiple scattering off',
0021                                       'PMAS(5,1)=4.8', # ! b quark mass',
0022                                       'PMAS(6,1)=175.0', # ! t quark mass'
0023                                       'CKIN(3)=7.',# ! ptMin
0024                                       'MSTJ(22)=2',
0025                                       'PARJ(71)=10.', # Decays only if life time < 10mm
0026                                       'PARP(67)=1.',
0027                                       'PARP(82)=1.9',
0028                                       'PARP(85)=0.33',
0029                                       'PARP(86)=0.66',
0030                                       'PARP(89)=1000.',
0031                                       'PARP(91)=1.0',
0032                                       'MSTJ(11)=3',
0033                                       'MSTJ(22)=2'
0034                                       ),
0035
0036     ppJets = cms.vstring('MSEL=1   ! QCD hight pT processes',
0037              'MSTP(122)=0' # ! no printout of Pythia initialization information hereinafter
0038             ),
0039     pythiaJets = cms.vstring('MSUB(11)=1', # q+q->q+q
0040                              'MSUB(12)=1', # q+qbar->q+qbar
0041                              'MSUB(13)=1', # q+qbar->g+g
0042                              'MSUB(28)=1', # q+g->q+g
0043                              'MSUB(53)=1', # g+g->q+qbar
0044                              'MSUB(68)=1' # g+g->g+g
0045                              ),
0046     pythiaPromptPhotons = cms.vstring('MSUB(14)=1', # q+qbar->g+gamma
0047                                       'MSUB(18)=1', # q+qbar->gamma+gamma
0048                                       'MSUB(29)=1', # q+g->q+gamma
0049                                       'MSUB(114)=1', # g+g->gamma+gamma
0050                                       'MSUB(115)=1' # g+g->g+gamma
0051                                       ),
0052
0053     pythiaWeakBosons = cms.vstring('MSUB(1)=1',
0054                                    'MSUB(2)=1'),
0055
0056     pythiaZjets = cms.vstring('MSUB(15)=1',
0057                               'MSUB(30)=1'),
0058
0059     pythiaCharmoniumNRQCD = cms.vstring('MSUB(421) = 1',
0060                                         'MSUB(422) = 1',
0061                                         'MSUB(423) = 1',
0062                                         'MSUB(424) = 1',
0063                                         'MSUB(425) = 1',
0064                                         'MSUB(426) = 1',
0065                                         'MSUB(427) = 1',
0066                                         'MSUB(428) = 1',
0067                                         'MSUB(429) = 1',
0068                                         'MSUB(430) = 1',
0069                                         'MSUB(431) = 1',
0070                                         'MSUB(432) = 1',
0071                                         'MSUB(433) = 1',
0072                                         'MSUB(434) = 1',
0073                                         'MSUB(435) = 1',
0074                                         'MSUB(436) = 1',
0075                                         'MSUB(437) = 1',
0076                                         'MSUB(438) = 1',
0077                                         'MSUB(439) = 1'
0078                                         ),
0079
0080     pythiaBottomoniumNRQCD = cms.vstring('MSUB(461) = 1',
0081                                          'MSUB(462) = 1',
0082                                          'MSUB(463) = 1',
0083                                          'MSUB(464) = 1',
0084                                          'MSUB(465) = 1',
0085                                          'MSUB(466) = 1',
0086                                          'MSUB(467) = 1',
0087                                          'MSUB(468) = 1',
0088                                          'MSUB(469) = 1',
0089                                          'MSUB(470) = 1',
0090                                          'MSUB(471) = 1',
0091                                          'MSUB(472) = 1',
0092                                          'MSUB(473) = 1',
0093                                          'MSUB(474) = 1',
0094                                          'MSUB(475) = 1',
0095                                          'MSUB(476) = 1',
0096                                          'MSUB(477) = 1',
0097                                          'MSUB(478) = 1',
0098                                          'MSUB(479) = 1',
0099                                          ),
0100     pythiaQuarkoniaSettings = cms.vstring('PARP(141)=1.16',   # Matrix Elements
0101                                           'PARP(142)=0.0119',
0102                                           'PARP(143)=0.01',
0103                                           'PARP(144)=0.01',
0104                                           'PARP(145)=0.05',
0105                                           'PARP(146)=9.28',
0106                                           'PARP(147)=0.15',
0107                                           'PARP(148)=0.02',
0108                                           'PARP(149)=0.02',
0109                                           'PARP(150)=0.085',
0110                                           # Meson spin
0111                                           'PARJ(13)=0.60',
0112                                           'PARJ(14)=0.162',
0113                                           'PARJ(15)=0.018',
0114                                           'PARJ(16)=0.054',
0115                                           # Polarization
0116                                           'MSTP(145)=0',
0117                                           'MSTP(146)=0',
0118                                           'MSTP(147)=0',
0119                                           'MSTP(148)=1',
0120                                           'MSTP(149)=1',
0121                                           # Chi_c branching ratios
0122                                           'BRAT(861)=0.202',
0123                                           'BRAT(862)=0.798',
0124                                           'BRAT(1501)=0.013',
0125                                           'BRAT(1502)=0.987',
0126                                           'BRAT(1555)=0.356',
0127                                           'BRAT(1556)=0.644'
0128                                           ),
0129     pythiaZtoMuons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0130                                  "MDME(175,1)=0",          # !Z decay into u ubar,
0131                                  "MDME(176,1)=0",          # !Z decay into s sbar,
0132                                  "MDME(177,1)=0",          # !Z decay into c cbar,
0133                                  "MDME(178,1)=0",          # !Z decay into b bbar,
0134                                  "MDME(179,1)=0",          # !Z decay into t tbar,
0135                                  "MDME(182,1)=0",          # !Z decay into e- e+,
0136                                  "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0137                                  "MDME(184,1)=1",          # !Z decay into mu- mu+,
0138                                  "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0139                                  "MDME(186,1)=0",          # !Z decay into tau- tau+,
0140                                  "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0141                                  ),
0142
0143
0144     pythiaZtoElectrons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0145                                      "MDME(175,1)=0",          # !Z decay into u ubar,
0146                                      "MDME(176,1)=0",          # !Z decay into s sbar,
0147                                      "MDME(177,1)=0",          # !Z decay into c cbar,
0148                                      "MDME(178,1)=0",          # !Z decay into b bbar,
0149                                      "MDME(179,1)=0",          # !Z decay into t tbar,
0150                                      "MDME(182,1)=1",          # !Z decay into e- e+,
0151                                      "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0152                                      "MDME(184,1)=0",          # !Z decay into mu- mu+,
0153                                      "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0154                                      "MDME(186,1)=0",          # !Z decay into tau- tau+,
0155                                      "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0156                                      ),
0157
0158     pythiaZtoMuonsAndElectrons = cms.vstring("MDME(174,1)=0",          # !Z decay into d dbar,
0159                                              "MDME(175,1)=0",          # !Z decay into u ubar,
0160                                              "MDME(176,1)=0",          # !Z decay into s sbar,
0161                                              "MDME(177,1)=0",          # !Z decay into c cbar,
0162                                              "MDME(178,1)=0",          # !Z decay into b bbar,
0163                                              "MDME(179,1)=0",          # !Z decay into t tbar,
0164                                              "MDME(182,1)=1",          # !Z decay into e- e+,
0165                                              "MDME(183,1)=0",          # !Z decay into nu_e nu_ebar,
0166                                              "MDME(184,1)=1",          # !Z decay into mu- mu+,
0167                                              "MDME(185,1)=0",          # !Z decay into nu_mu nu_mubar,
0168                                              "MDME(186,1)=0",          # !Z decay into tau- tau+,
0169                                              "MDME(187,1)=0"           # !Z decay into nu_tau nu_taubar
0170                                              ),
0171
0172     pythiaUpsilonToMuons = cms.vstring('BRAT(1034) = 0 ',  # switch off',
0173                                        'BRAT(1035) = 1 ',  # switch on',
0174                                        'BRAT(1036) = 0 ',  # switch off',
0175                                        'BRAT(1037) = 0 ',  # switch off',
0176                                        'BRAT(1038) = 0 ',  # switch off',
0177                                        'BRAT(1039) = 0 ',  # switch off',
0178                                        'BRAT(1040) = 0 ',  # switch off',
0179                                        'BRAT(1041) = 0 ',  # switch off',
0180                                        'BRAT(1042) = 0 ',  # switch off',
0181                                        'MDME(1034,1) = 0 ',  # switch off',
0182                                        'MDME(1035,1) = 1 ',  # switch on',
0183                                        'MDME(1036,1) = 0 ',  # switch off',
0184                                        'MDME(1037,1) = 0 ',  # switch off',
0185                                        'MDME(1038,1) = 0 ',  # switch off',
0186                                        'MDME(1039,1) = 0 ',  # switch off',
0187                                        'MDME(1040,1) = 0 ',  # switch off',
0188                                        'MDME(1041,1) = 0 ',  # switch off',
0189                                        'MDME(1042,1) = 0 ',  # switch off'
0190                                        ),
0191
0192     pythiaJpsiToMuons = cms.vstring('BRAT(858) = 0 ',  # switch off',
0193                                     'BRAT(859) = 1 ',  # switch on',
0194                                     'BRAT(860) = 0 ',  # switch off',
0195                                     'MDME(858,1) = 0 ',  # switch off',
0196                                     'MDME(859,1) = 1 ',  # switch on',
0197                                     'MDME(860,1) = 0 ',  # switch off'
0198                                     ),
0199
0200     pythiaMuonCandidates = cms.vstring(
0201     'CKIN(3)=20',
0202     'MSTJ(22)=2',
0203     'PARJ(71)=40.'
0204     ),
0205     myParameters = cms.vstring(
0206     )
0207
0208     )