| 1 | (* Automatically generated from HDecay output *
|
|---|
| 2 | * (c.f. "HDECAY: A Program for Higgs boson decays in the standard model and *
|
|---|
| 3 | * its supersymmetric extension.", A. Djouadi, J. Kalinowski, M. Spira, *
|
|---|
| 4 | * Comput.Phys.Commun.108:56-74,1998) *)
|
|---|
| 5 |
|
|---|
| 6 | BeginPackage ["Hdecay`"];
|
|---|
| 7 | Begin ["Hdecay`p`"];
|
|---|
| 8 |
|
|---|
| 9 | data = {{0, 0}, {1., 2.251*^-7}, {6.0201, 0.00002011}, {11.0402, 0.0001738}, {16.0603, 0.0004549},
|
|---|
| 10 | {21.0804, 0.0006405}, {26.1005, 0.0007918}, {31.1206, 0.0009277}, {36.1407, 0.001055},
|
|---|
| 11 | {41.1608, 0.001178}, {46.1809, 0.001298}, {51.201, 0.001416}, {56.2211, 0.001533},
|
|---|
| 12 | {61.2412, 0.001649}, {66.2613, 0.001766}, {71.2814, 0.001883}, {76.3015, 0.002},
|
|---|
| 13 | {81.3216, 0.002119}, {86.3417, 0.00224}, {91.3618, 0.002364}, {96.3819, 0.002495},
|
|---|
| 14 | {101.402, 0.002642}, {106.422, 0.002816}, {111.442, 0.003041}, {116.462, 0.003347},
|
|---|
| 15 | {121.482, 0.003782}, {126.503, 0.00442}, {131.523, 0.005373}, {136.543, 0.00683},
|
|---|
| 16 | {141.563, 0.009133}, {146.583, 0.013}, {151.603, 0.02027}, {156.623, 0.03847},
|
|---|
| 17 | {161.643, 0.1316}, {166.663, 0.2903}, {171.683, 0.4168}, {176.704, 0.5386},
|
|---|
| 18 | {181.724, 0.6883}, {186.744, 0.9035}, {191.764, 1.105}, {196.784, 1.3}, {201.804, 1.499},
|
|---|
| 19 | {206.824, 1.705}, {211.844, 1.922}, {216.864, 2.151}, {221.884, 2.394}, {226.905, 2.651},
|
|---|
| 20 | {231.925, 2.924}, {236.945, 3.213}, {241.965, 3.519}, {246.985, 3.843}, {252.005, 4.184},
|
|---|
| 21 | {257.025, 4.545}, {262.045, 4.925}, {267.065, 5.324}, {272.085, 5.742}, {277.106, 6.182},
|
|---|
| 22 | {282.126, 6.642}, {287.146, 7.123}, {292.166, 7.625}, {297.186, 8.15}, {302.206, 8.696},
|
|---|
| 23 | {307.226, 9.265}, {312.246, 9.856}, {317.266, 10.47}, {322.286, 11.1}, {327.307, 11.76},
|
|---|
| 24 | {332.327, 12.44}, {337.347, 13.13}, {342.367, 13.83}, {347.387, 14.61}, {352.407, 15.7},
|
|---|
| 25 | {357.427, 16.91}, {362.447, 18.18}, {367.467, 19.51}, {372.487, 20.89}, {377.508, 22.31},
|
|---|
| 26 | {382.528, 23.77}, {387.548, 25.27}, {392.568, 26.81}, {397.588, 28.38}, {402.608, 29.99},
|
|---|
| 27 | {407.628, 31.63}, {412.648, 33.3}, {417.668, 35.01}, {422.688, 36.76}, {427.709, 38.53},
|
|---|
| 28 | {432.729, 40.34}, {437.749, 42.19}, {442.769, 44.06}, {447.789, 45.97}, {452.809, 47.92},
|
|---|
| 29 | {457.829, 49.9}, {462.849, 51.91}, {467.869, 53.96}, {472.889, 56.05}, {477.91, 58.17},
|
|---|
| 30 | {482.93, 60.33}, {487.95, 62.52}, {492.97, 64.76},
|
|---|
| 31 | {497.99, 67.03}};
|
|---|
| 32 |
|
|---|
| 33 | interpolation[x_, l_] := Plus @@ ((#[[1]] * x^#[[2]])& /@ l);
|
|---|
| 34 |
|
|---|
| 35 | params = {{-0.26293901351164234, 1}, {0.001368182199849752, 2}, {-2.3048363155706936*^-6, 3},
|
|---|
| 36 | {3.1498414039397794*^-9, 4}, {-1.977234012939224*^-12, 5},
|
|---|
| 37 | {6.665973758092558*^-16, 6}};
|
|---|
| 38 |
|
|---|
| 39 | Hdecay`Hdecay[mh_] := If[mh < 500, Max[Interpolation[data, InterpolationOrder->1][mh], 0], interpolation[mh, params]]/;(mh >= 0);
|
|---|
| 40 | Hdecay`Hdecay::usage = "usage: Hdecay [mh]";
|
|---|
| 41 |
|
|---|
| 42 | End[];
|
|---|
| 43 | Protect[Hdecay];
|
|---|
| 44 | EndPackage[];
|
|---|
