pyslha.py

changeset 120
bffefe12df80
parent 118
5ab517422296
child 123
c85e29bc13c4
     1.1 --- a/pyslha.py	Sat Feb 26 02:27:08 2011 +0100
     1.2 +++ b/pyslha.py	Sat Feb 26 23:21:12 2011 +0100
     1.3 @@ -23,7 +23,7 @@
     1.4  """
     1.5  
     1.6  __author__ = "Andy Buckley <andy.buckley@cern.ch"
     1.7 -__version__ = "1.0.7"
     1.8 +__version__ = "1.0.8"
     1.9  
    1.10  
    1.11  def _autotype(var):
    1.12 @@ -252,6 +252,7 @@
    1.13      return rtn
    1.14  
    1.15  
    1.16 +
    1.17  def readISAWIG(isastr, ignorenobr=False):
    1.18      """
    1.19      Read a spectrum definition from a string in the ISAWIG format, returning
    1.20 @@ -266,10 +267,31 @@
    1.21      branching ratio of zero.
    1.22      """
    1.23  
    1.24 -    ## PDG MC ID codes mapped to HERWIG SUSY ID codes, based on
    1.25 +    ## PDG MC ID codes mapped to HERWIG IDHW codes, based on
    1.26      ## http://www.hep.phy.cam.ac.uk/~richardn/HERWIG/ISAWIG/susycodes.html
    1.27 +    ## + the IDPDG array and section 4.13 of the HERWIG manual.
    1.28      HERWIGID2PDGID = {}
    1.29 -    HERWIGID2PDGID[203] =  25 ## HIGGSL0
    1.30 +    HERWIGID2PDGID[7]   = -1
    1.31 +    HERWIGID2PDGID[8]   = -2
    1.32 +    HERWIGID2PDGID[9]   = -3
    1.33 +    HERWIGID2PDGID[10]  = -4
    1.34 +    HERWIGID2PDGID[11]  = -5
    1.35 +    HERWIGID2PDGID[12]  = -6
    1.36 +    HERWIGID2PDGID[13]  =  21
    1.37 +    HERWIGID2PDGID[59]  =  22
    1.38 +    HERWIGID2PDGID[121] =  11
    1.39 +    HERWIGID2PDGID[122] =  12
    1.40 +    HERWIGID2PDGID[123] =  13
    1.41 +    HERWIGID2PDGID[124] =  14
    1.42 +    HERWIGID2PDGID[125] =  15
    1.43 +    HERWIGID2PDGID[126] =  16
    1.44 +    HERWIGID2PDGID[127] = -11
    1.45 +    HERWIGID2PDGID[128] = -12
    1.46 +    HERWIGID2PDGID[129] = -13
    1.47 +    HERWIGID2PDGID[130] = -14
    1.48 +    HERWIGID2PDGID[131] = -15
    1.49 +    HERWIGID2PDGID[132] = -16
    1.50 +    HERWIGID2PDGID[203] =  25 ## HIGGSL0 (== PDG standard in this direction)
    1.51      HERWIGID2PDGID[204] =  35 ## HIGGSH0
    1.52      HERWIGID2PDGID[205] =  36 ## HIGGSA0
    1.53      HERWIGID2PDGID[206] =  37 ## HIGGS+
    1.54 @@ -517,10 +539,31 @@
    1.55      If the ignorenobr parameter is True, do not write decay entries with a
    1.56      branching ratio of zero.
    1.57      """
    1.58 -    ## PDG MC ID codes mapped to HERWIG SUSY ID codes, based on
    1.59 +    ## PDG MC ID codes mapped to HERWIG IDHW codes, based on
    1.60      ## http://www.hep.phy.cam.ac.uk/~richardn/HERWIG/ISAWIG/susycodes.html
    1.61 +    ## + the IDPDG array and section 4.13 of the HERWIG manual.
    1.62      PDGID2HERWIGID = {}
    1.63 -    PDGID2HERWIGID[      25] = 203 ## HIGGSL0 (ADDED)
    1.64 +    PDGID2HERWIGID[      -1] = 7
    1.65 +    PDGID2HERWIGID[      -2] = 8
    1.66 +    PDGID2HERWIGID[      -3] = 9
    1.67 +    PDGID2HERWIGID[      -4] = 10
    1.68 +    PDGID2HERWIGID[      -5] = 11
    1.69 +    PDGID2HERWIGID[      -6] = 12
    1.70 +    PDGID2HERWIGID[      21] = 13
    1.71 +    PDGID2HERWIGID[      22] = 59
    1.72 +    PDGID2HERWIGID[      11] = 121
    1.73 +    PDGID2HERWIGID[      12] = 122
    1.74 +    PDGID2HERWIGID[      13] = 123
    1.75 +    PDGID2HERWIGID[      14] = 124
    1.76 +    PDGID2HERWIGID[      15] = 125
    1.77 +    PDGID2HERWIGID[      16] = 126
    1.78 +    PDGID2HERWIGID[     -11] = 127
    1.79 +    PDGID2HERWIGID[     -12] = 128
    1.80 +    PDGID2HERWIGID[     -13] = 129
    1.81 +    PDGID2HERWIGID[     -14] = 130
    1.82 +    PDGID2HERWIGID[     -15] = 131
    1.83 +    PDGID2HERWIGID[     -16] = 132
    1.84 +    PDGID2HERWIGID[      25] = 203 ## HIGGSL0 (added for PDG standard -> HERWIG IDHW)
    1.85      PDGID2HERWIGID[      26] = 203 ## HIGGSL0
    1.86      PDGID2HERWIGID[      35] = 204 ## HIGGSH0
    1.87      PDGID2HERWIGID[      36] = 205 ## HIGGSA0

mercurial