Mon, 15 Jul 2013 16:50:38 +0200
Add credit to ChangeLog
1 #! /usr/bin/env python
3 import pyslha
5 blocks, decays = pyslha.readSLHA("""\
6 BLOCK DCINFO # Decay Program information
7 1 SDECAY # decay calculator
8 2 1.1a # version number
9 #
10 BLOCK SPINFO # Spectrum calculator information
11 1 SOFTSUSY # spectrum calculator
12 2 2.0.5 # version number
13 #
14 BLOCK MODSEL # Model selection
15 1 1 sugra
17 BLOCK NMIX # Neutralino Mixing Matrix
18 1 1 9.86364430E-01 # N_11
19 1 2 -5.31103553E-02 # N_12
20 1 4 -5.31186117E-02 # N_14
21 1 3 1.46433995E-01 # N_13
22 2 1 9.93505358E-02 # N_21
23 2 2 9.44949299E-01 # N_22
24 2 3 -2.69846720E-01 # N_23
25 2 4 1.56150698E-01 # N_24
27 BLOCK MASS # Mass Spectrum
28 # PDG code mass particle
29 5 4.88991651E+00 # b-quark pole mass calculated from mb(mb)_Msbar
30 6 1.75000000E+02 # mt pole mass (not read by ME)
31 24 7.98290131E+01 # W+
32 25 1.10899057E+02 # h
33 35 3.99960116E+02 # H
34 36 3.99583917E+02 # A
35 37 4.07879012E+02 # H+
36 1000001 5.68441109E+02 # ~d_L
37 """)
40 for block in blocks.values():
41 print block
43 print pyslha.writeSLHA(blocks, decays)