SoftSUSY is hosted by Hepforge, IPPP Durham
SOFTSUSY  4.1
twoBodyDecays.h
1 
11 #ifndef TWOBODYDECAYS_H
12 #define TWOBODYDECAYS_H
13 
14 #include "nmssmsoftsusy.h"
15 #include "decays.h"
16 #include "softsusy.h"
17 #include "physpars.h"
18 #include "lowe.h"
19 #include "softpars.h"
20 #include "softsusy.h"
21 #include "flavoursoft.h"
22 #include "susy.h"
23 #include "particle.h"
24 #include <iostream>
25 #include <cstring>
26 #include <stdlib.h>
27 #include <cmath>
28 #include <stdlib.h>
29 #include <vector>
30 #include <iomanip>
31 #include <complex>
32 
33 using namespace std;
34 const double fpi = 0.13041 / root2, mpiplus = 0.13957018, mpi0 = 0.1349766;
35 
37 double charginoToNeutralino1pion(const MssmSoftsusy * m);
38 
39 double gluinoamplitudedecay (double m1, double m2, double m3,
40  double alphastrong);
41 double gluinoamplitudedecaymix (double m1, double m2, double m3,
42  double alphastrong, double squarkmix,
43  double theta);
44 double squarkamplitudedecaygluino (double m1, double m2, double m3,
45  double alphastrong);
46 double squarkamplitudedecaygluinomix (double m1, double m2, double m3,
47  double alphastrong, double squarkmix,
48  double theta);
49 double squarkamplitudedecaycharginoW1 (double m1, double m2, double m3,
50  double g, double gamma);
51 double squarkamplitudedecaycharginoW2 (double m1, double m2, double m3,
52  double g, double gamma);
54 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
55  double theta, double beta, double mWboson, double runmt, double runmb,
56  double torb);
58 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
59  double theta, double beta, double mWboson, double runmt, double runmb,
60  double torb);
62 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
63  double theta, double beta, double mWboson, double mup, double mdown,
64  double torb);
66 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
67  double theta, double beta, double mWboson, double mup, double mdown,
68  double torb);
70 (double m1, double m2, double m3, double g, double gprime,
71  DoubleMatrix & mixNeut, int neutralino, int uord );
73 (double m1, double m2, double m3, double g, double gprime,
74  DoubleMatrix & mixNeut, int neutralino, int uord );
76 (double m1, double m2, double m3, double mWboson, double theta, double beta,
77  DoubleMatrix & mixNeut, double g, double gp, double runmt, int squark,
78  int oneortwo, int neutralino);
80 (double m1, double m2, double m3, double g, double thetat, double thetab,
81  int m1torb, int m1oneortwo, int m3torb, int m3oneortwo);
83 (double m1, double m2, double m3, double g, double mWboson, double beta,
84  double thetat, double thetab, double greekmu, double At, double Ab,
85  double mt, double mb, int t1or2, int b1or2);
86  double squark32amplitudedecayneutralHiggssquark3 (double m1, double m2, double m3, double g, double gp, double mWboson, double beta, double alpha, double thetat, double thetab, double greekmu, double At, double Ab, double mt, double mb, int torb, char phi);
87  double squark32amplitudedecaysquark3Zboson (double m1, double m2, double m3, double g, double gp, double theta);
88  double sleptonamplitudedecayleptonneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, char LorR, int neutralino);
89  double sneutrinoamplitudedecayneutrinoneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
90  double sleptonamplitudedecaychargino (double m1, double m2, double m3, double g, double theta, int chargino);
91  double stauamplitudedecaytauneutralino (double m1, double m2, double m3, double g, double gp, double mWboson, DoubleMatrix & mixNeut, double theta, double beta, int oneortwo, int neutralino);
92  double stausneutrinoamplitudedecaytauneutrinoneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
93  double stauamplitudedecaynutauchargino (double m1, double m2, double m3, double g, double mWboson, double theta, double thetaL, double beta, double mtau, int oneortwo, int chargino);
94  double stausneutrinoamplitudedecaytauchargino (double m1, double m2, double m3, double g, double mWboson, double beta, double thetaL, double thetaR, int chargino);
95  double stauamplitudedecaysnustauHminus (double m1, double m2, double m3, double g, double mWboson, double beta, double thetatau, double mtau, double mu, double Atau, int oneortwo);
96  double stauamplitudedecaysnustauWboson (double m1, double m2, double m3, double g, double thetatau, int oneortwo);
97  double stau2amplitudedecaystau1Zboson (double m1, double m2, double m3, double g, double gp, double thetatau);
98  double stau2amplitudedecaystau1phi (double m1, double m2, double m3, double g, double gp, double thetatau, double beta, double alpha, double mWboson, double mtau, double greekmu, double Atau, char phi);
99  double charginoamplitudedecayquarksquarkL (double m1, double m2, double m3, double g, double theta, int chargino);
100  double charginoamplitudedecayquarksquarkmix (double m1, double m2, double m3, double g, double theta, double thetaL, double thetaR, double beta, double runmt, double runmb, double mWboson, int chargino, int upordowntypesquark, int oneortwo);
101  double charginoamplitudedecayleptonsleptonL (double m1, double m2, double m3, double g, double thetaLorR, int chargino);
102  double charginoamplitudedecaysnutautau (double m1, double m2, double m3, double g, double thetaL, double thetaR, double beta, double mWboson, int chargino);
103  double charginoamplitudedecaystaunutau (double m1, double m2, double m3, double g, double thetaL, double thetaR, double thetatau, double beta, double mWboson, double mtau, int oneortwo, int chargino);
104  double charginoamplitudedecayWbosonneutralino (double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int chargino, int neutralino);
105  double charginoamplitudedecayHminusneutralino (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, DoubleMatrix & mixNeut, int chargino, int neutralino);
106  double chargino2amplitudedecaychargino1Zboson (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR);
107  double chargino2amplitudedecaychargino1neutHiggs (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, double alpha, char phi);
108  double neutralinoamplitudedecayquarksquarkLorR (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int uordtype , char LorR, int neutralino);
109  double neutralinoamplitudedecayleptonsleptonLorR (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, char LorR, int neutralino);
110  double neutralinoamplitudedecayneutrinosneutrinoL (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
111  double neutralinoamplitudedecaysquark3quarkmix (double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix & mixNeut, double g, double gp, double runmq, int squark , int oneortwo, int neutralino);
112  double neutralinoamplitudedecaystautau (double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix & mixNeut, double g, double gp, int oneortwo, int neutralino);
113  double neutralinoamplitudedecaycharginoWboson (double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
114  double neutralinoamplitudedecaycharginoHplus (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
115  double neutralinoamplitudedecayneutralinoZboson (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int ineutralino, int fneutralino);
116  double neutralinoamplitudedecayneutralinoneutHiggs (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, double mixingangle, int ineutralino, int fneutralino, char phi);
117 
118  double higgslorHamplitudedecayquarkantiquark (double m1, double m2, double g, double alpha, double beta, double mWboson, int uord, char lorH, DoubleMatrix & CPEMix, bool nmssmIsIt, bool QCD, double alphas);
119  double higgsAamplitudedecayquarkantiquark (double m1, double m2, double g, double beta, double mWboson, int uord, bool QCD, double alphas);
120  double higgsAamplitudedecayquarkantiquarkNMSSM (double m1, double m2, double beta, DoubleMatrix & CPOMix, int uord, int higgs, bool QCD, double alphas);
121  double higgsphiamplitudedecayneutralinoneutralino (double m1, double m2, double m3, double g, double tanthetaW, double mixingangle, DoubleMatrix & mixNeut, int ineutralino, int fneutralino, char phi);
122  double higgsphiamplitudedecaysamechargino (double m1, double m2, double g, double thetaL, double thetaR, double alpha, double beta, int chargino, char phi);
123  double higgsphiamplitudedecaydifchargino (double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double beta, char phi);
124  double higgshamplitudedecayAA (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
125  double higgsHamplitudedecayhh (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
126  double higgsHamplitudedecayAA (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
127  double higgsHamplitudedecayHplusHminus (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
128  double higgshamplitudedecayhiggsAZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
129  double higgsHamplitudedecayhiggsAZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
130  double higgsAamplitudedecayhiggshZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
131  double higgsAamplitudedecayhiggsHZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
132  double higgshamplitudedecay2squarksamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mupq, double mdownq, int sq);
133  double higgshamplitudedecay2squarkdiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mupq, double mdownq, double greekmu, double Aup, double Adown, int sq);
134  double higgsHamplitudedecay2squarksamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mupq, double mdownq, int sq);
135  double higgsHamplitudedecay2squarkdiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mupq, double mdownq, double greekmu, double Aup, double Adown, int sq);
136  double higgshamplitudedecay2sleptonsamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl);
137  double higgshamplitudedecay2sleptondiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl);
138  double higgsHamplitudedecay2sleptonsamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl);
139  double higgsHamplitudedecay2sleptondiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl);
140  double higgshamplitudedecaystop1stop1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
141  double higgshamplitudedecaystop2stop2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
142  double higgshamplitudedecaystop1stop2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
143  double higgshamplitudedecaysbottom1sbottom1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
144  double higgshamplitudedecaysbottom2sbottom2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
145  double higgshamplitudedecaysbottom1sbottom2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
146  double higgsHamplitudedecaystop1stop1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
147  double higgsHamplitudedecaystop2stop2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
148  double higgsHamplitudedecaystop1stop2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
149  double higgsHamplitudedecaysbottom1sbottom1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
150  double higgsHamplitudedecaysbottom2sbottom2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
151  double higgsHamplitudedecaysbottom1sbottom2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta);
152  double higgshamplitudedecaystau1stau1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
153  double higgshamplitudedecaystau2stau2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
154  double higgshamplitudedecaystau1stau2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
155  double higgshtestamplitudedecaystau1stau1 (double m1, double m2, double thetatau, double g, double gp, double mWboson, double alpha, double beta, double mtau, double greekmu, double Atau);
156  double higgsHamplitudedecaystau1stau1 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
157  double higgsHamplitudedecaystau2stau2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
158  double higgsHamplitudedecaystau1stau2 (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta);
159  double higgsAamplitudedecaysfermions (double m1, double m2, double m3, double g, double mWboson, double mf, double greekmu, double Asf, double beta, char uord);
160  double higgsHplusamplitudedecayquarkantiquark (double m1, double m2, double m3, double g, double mWboson, double beta, DoubleMatrix & VCKM, int quark, int antiquark);
161  double higgsHplusamplitudedecayneutralinochargino (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
162  double higgsHplusamplitudedecayneutralinocharginosusyhitway (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
163  double higgsHplusamplitudedecayWbosonhiggsh (double m1, double m2, double m3, double g, double alpha, double beta);
164  DoubleVector higgsHplusamplitudedecaysquarksquark (double m1, double m2, double m3, double g, double beta, double mWboson, double mup, double mdown, double greekmu, double Aup, double Adown);
165  DoubleVector higgsHplusamplitudedecaysquarksquarkmix (double m1, double m2, double m3, double g, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double thetatop, double thetabottom);
166  double higgsesamplitudedecaygammagammatotal(double m1, double g, double gprime, double alphaEmrun, double mWboson, double polemw, double alpha, double beta, double mtop, double mbottom, double mcharm, double mtau, double mHpm, double mstop1, double mstop2, double msbottom1, double msbottom2, double mstau1, double mstau2, double mchargino1, double mchargino2, double thetaL, double thetaR, double thetat, double thetab, double thetatau, double greekmu, double Atop, double Abottom, double Atau, char higgstype);
167  double higgsesamplitudedecaygluongluontotal(double m1, double g, double gs, double gprime, double mWboson, double alpha, double beta, double mtop, double mbottom, double mcharm, double mstop1, double mstop2, double msbottom1, double msbottom2, double thetat, double thetab, double greekmu, double Atop, double Abottom, double mstrange, double mscharmL, double mscharmR, double msstrangeL, double msstrangeR, double Acharm, double Astrange, double mup, double mdown, double msupL, double msupR, double msdownL, double msdownR, double Aup, double Adown, char higgstype, bool QCD);
169  DoubleVector higgshamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix & CPEMix, bool nmssmIsIt);
170  DoubleVector higgsHamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix & CPEMix, bool nmssmIsIt);
171  DoubleVector higgsH3amplitudedecayVVNMSSM(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix & CPEMix, bool nmssmIsIt);
172  double higgsAamplitudedecayHpmWboson(double m1, double mWboson, double mHpm, double g, double thetaA, int pseudoscalar, bool nmssmIsIt);
173  double stop2amplitudedecaystop1CPevenhiggsNMSSM (double mst2, double mst1, double mh, double mt , double thetat, DoubleMatrix & CPEMix, double beta, double mWboson, double g, double gp, double At, double mueff, double lam, int higgs);
174  double stop2amplitudedecaystop1CPoddhiggsNMSSM (double mst2, double mst1, double ma, double mt , double thetat, DoubleMatrix & CPOMix, double beta, double mWboson, double g, double At, double mueff, double lam, int higgs);
175  double sbottom2amplitudedecaysbottom1CPevenhiggsNMSSM (double msb2, double msb1, double mh, double mb , double thetab, DoubleMatrix & CPEMix, double beta, double mWboson, double g, double gp, double Ab, double mueff, double lam, int higgs);
176  double sbottom2amplitudedecaysbottom1CPoddhiggsNMSSM (double msb2, double msb1, double ma, double mb , double thetab, DoubleMatrix & CPOMix, double beta, double mWboson, double g , double Ab, double mueff, double lam, int higgs);
177  double stau2amplitudedecaystau1CPevenhiggsNMSSM (double mstau2, double mstau1, double mh, double mtau, double thetatau, DoubleMatrix & CPEMix, double beta, double mWboson, double g, double gp, double Atau, double mueff, double lam, int higgs);
178  double stau2amplitudedecaystau1CPoddhiggsNMSSM (double mstau2, double mstau1, double ma, double mtau, double thetatau, DoubleMatrix & CPOMix, double beta, double mWboson, double g , double Atau, double mueff, double lam, int higgs);
179  double chargino2amplitudedecaychargino1CPevenhiggsNMSSM (double mchar2, double mchar1, double mh, double g, double lam, double thetaL, double thetaR, DoubleMatrix & CPEMix, int higgs);
180  double chargino2amplitudedecaychargino1CPoddhiggsNMSSM (double mchar2, double mchar1, double mA, double g, double lam, double thetaL, double thetaR, DoubleMatrix & CPOMix, int higgs);
181  double neutralinoamplitudedecaycharginoWNMSSM (double mneut, double mchar, double mWboson, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
182  double neutralinoamplitudedecayneutralinoZNMSSM (double mneuti, double mneutj, double mZboson, double g, double gp, DoubleMatrix & mixNeut, int neutralinoi, int neutralinoj);
183  double neutralinoamplitudecharginoHpmNMSSM (double mneut, double mchar, double mHp, double g, double gp, double beta, double thetaL, double thetaR, double lam, DoubleMatrix & mixNeut, int neutralino, int chargino);
184  double neutralinoamplitudedecayneutralinoCPevenhiggsNMSSM (double mneuti, double mneutj, double mhiggs, double g, double gp, double lam, double kappa, DoubleMatrix & mixNeut, DoubleMatrix & CPEMix, int neutralinoi, int neutralinoj, int higgs);
185  double neutralinoamplitudedecayneutralinoCPoddhiggsNMSSM (double mneuti, double mneutj, double ma, double g, double gp, double lam, double kappa, DoubleMatrix & mixNeut, DoubleMatrix & CPOMix, int neuti, int neutj, int higgsa);
186  double neutralinoamplitudedecaysfermionfermionfirst2genNMSSM (double mneut, double msf, double mf, double g, double gp, DoubleMatrix & mixNeut, int neut, char type, char LorR);
187  double neutralinoamplitudestoptopNMSSM (double mneut, double mst, double mt, double g, double gp, double thetat, double beta, double mWboson, DoubleMatrix & mixNeut, double runmt, int neut, int stop);
188  double neutralinoamplitudesbottombottomNMSSM (double mneut, double msb, double mb, double g, double gp, double thetab, double beta, double mWboson, DoubleMatrix & mixNeut, double runmb, int neut, int sbottom);
189  double neutralinoamplitudestautauNMSSM (double mneut, double mstau, double mtau, double g, double gp, double thetatau, double beta, double mWboson, DoubleMatrix & mixNeut, double runmtau, int neut, int stau);
190  double neutralinoamplitudestauneutrinotauneutrinoNMSSM (double mneut, double mstaunu, double mtaunu, double g, double gp, DoubleMatrix & mixNeut, int neut);
191 
192  double squarkamplitudedecayquarkneutralinoNMSSM (double m1, double mq, double mneut, double g, double gp, DoubleMatrix & mixNeut, char uord, char LorR, int neut);
193  double sleptonamplitudedecayleptonneutralinoNMSSM (double m1, double ml, double mneut, double g, double gp, DoubleMatrix & mixNeut, char uord, char LorR, int neut);
194  double stopamplitudedecaytopneutralinoNMSSM (double m1, double mt, double mneut, double g, double gp, double thetat, DoubleMatrix & mixNeut, double runmt, double mWboson, double beta, int stop, int neut);
195  double sbottomamplitudedecaybottomneutralinoNMSSM (double m1, double mb, double mneut, double g, double gp, double thetab, DoubleMatrix & mixNeut, double runmb, double mWboson, double beta, int sbottom, int neut);
196  double stauamplitudedecaytauneutralinoNMSSM (double m1, double mtau, double mneut, double g, double gp, double thetatau, DoubleMatrix & mixNeut, double runmtau, double mWboson, double beta, int stau, int neut);
197  double charginoiamplitudedecayneutralinojHpmNMSSM (double mchar, double mneut, double mHpm, double g, double gp, double thetaL, double thetaR, double beta, DoubleMatrix & mixNeut, double lam, int chargino, int neut);
198  double charginoiamplitudedecayneutralinojWNMSSM (double mchar, double mneut, double mWboson, double g, double gp, double thetaL, double thetaR, DoubleMatrix & mixNeut, int chargino, int neut);
199  double HpmamplitudecharginojneutralinoiNMSSM (double mHp, double mchar, double mneut, double g, double gp, double beta, double thetaL, double thetaR, double lam, DoubleMatrix & mixNeut, int neutralino, int chargino);
200  double snutauamplitudedecaynutauneutralinoNMSSM (double m1, double mneut, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
201 
202  double higgsesamplitudedecayZbosonphotontotal(double m1, double mZboson, double g, double gprime, double alphaEmrun, double polemw, double runmw, double alpha, double beta, double mtop, double mbottom, double mcharm, double mstrange, double mstop1, double mstop2, double msbottom1, double msbottom2, double mHplus, double thetat, double thetab, double greekmu, double Atop, double Abottom, char higgstype);
203 double gluinoamplitudedecaydgausscharginoqqpbarfirsttwogen (double mgluino, double mchargino, double mquark, double mquarkp, double msqL, double msqpL, double g, double thetaL, double thetaR, double alphas, int charg, bool onetothree);
205 double gluinoamplitudedecay1to3neutfirsttwogen (double m1, double m2, double m3, double m4, double m5, double g, double gp, DoubleMatrix & mixNeut, double alphas, char uord, int neut, int Nsteps, int adaptive, bool onetothree, double approx);
206  double gluinoamplitudedecay1to3charfirsttwogen (double m1, double m2, double m3, double m4, double m5, double m6, double g, double thetaL, double thetaR, double alphas, int charg, int Nsteps, int adaptive, bool onetothree, double approx);
207  double gluinoamplitudedecay1to3neutttbar (double m1, double m2, double m3, double m4, double m5, double mw, double g, double gp, double thetat, double beta,double alphas, DoubleMatrix & mixNeut, double runmq, int neutralino, int Nsteps, int adaptive, bool onetothree, double approx);
208  double gluinoamplitudedecay1to3neutbbbar (double m1, double m2, double m3, double m4, double m5, double mw, double g, double gp, double thetab, double beta,double alphas, DoubleMatrix & mixNeut, double runmq, int neutralino, int Nsteps, int adaptive, bool onetothree, double approx);
209  double gluinoamplitudedecaychartbbar (double m1, double m2, double m3, double m4, double m5, double m6, double m7, double m8, double alphas, double thetat, double thetab, double mw, double g, double gp, double gammaL, double gammaR, double beta, double runmt, double runmb, int chargino, int Nsteps, int adaptive, bool onetothree, double approx);
210  double neutralinoamplitudedecayneutffbar (double m1, double m2, double mf, double msf1, double msf2, double mz, double mh, double mH, double mA, double runmf, double mw, double thetaf, double beta, double alpha, double g, double gp, DoubleMatrix & mixNeut, int neutralinoj, int neutralinoi, char qorl, char uord, int Nsteps, int adaptive, bool onetothree, double approx);
211  double gluinoamplitudedecaygravitino (double m1, double mgrav, double MPlreduced, int gravonoff, int gluNLSP);
212  double squarkamplitudedecaygravitino(double m1, double mgrav, double mquark, double MPlreduced, int gravonoff, int squNLSP);
213  double neutralinoamplitudedecayphotongravitino(double m1, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double g, double gp, int neutralino, int gravonoff, int neutNLSP);
214  double neutralinoamplitudedecayZgravitino(double m1, double mZ, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double g, double gp, double beta, int neutralino, int gravonoff, int neutNLSP);
215  double neutralinoamplitudedecayphigravitino(double m1, double mphi, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double alpha, double beta, int neutralino, int gravonoff, char phi, int neutNLSP);
216 
218  double higgsAamplitudedecaysfermionsNMSSM (double m1, double m2, double m3, double g, double mWboson, double mf, double Asf, double beta, double lam, double mueff, DoubleMatrix & CPOMix, char uord, int pseudoscalar);
219  double higgsAamplitudedecaysamecharginoNMSSM (double m1, double m2, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix & CPOMix, int chargino, int pseudoscalar);
220  double higgsAamplitudedecaydifcharginoNMSSM (double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix & CPOMix, int pseudoscalar);
221  double higgsAamplitudedecayneutralinoneutralinoNMSSM(double m1, double m2, double m3, double g, double tanthetaW, double lam, double kappa, DoubleMatrix & CPOMix, DoubleMatrix & mixNeut, int ineutralino, int fneutralino, int pseudoscalar);
222  double higgsAamplitudedecayhiggshorHZbosonNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double thetaA, DoubleMatrix & CPEMix, int pseudoscalar, int higgs);
223  double higgslHamplitudedecayquarkantiquarkNMSSM (double MSSMamplitude, DoubleMatrix & CPEMix, double alpha, int higgs, char uord);
224  double higgshamplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
225  double higgshamplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
226  double higgsHamplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
227  double higgsHamplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
228  double higgsH3amplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
229  double higgsH3amplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
230  double higgshamplitudedecay2squarkdiffhandNMSSM (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson,double mq, double Aq, double mueff, double lam, DoubleMatrix & CPEMix, int sq, int higgs);
231  double higgsphiamplitudedecaysamecharginoNMSSM (double m1, double m2, double g, double thetaL, double thetaR, double lam, DoubleMatrix & CPEMix, int chargino, int higgs);
232  double higgsphiamplitudedecaydiffcharginoNMSSM (double m1, double m2, double m3, double g, double thetaL, double thetaR, double lam, DoubleMatrix & CPEMix, int higgs);
233  double higgsAamplitudedecaygammagammaNMSSM (double m1, double g, double gprime, double alpha, double mWboson, DoubleMatrix & CPOMix, double beta, double mtop, double mbottom, double mcharm, double mtau, double mch1, double mch2, double thetaL, double thetaR, double lam, int higgs);
234  double higgsAamplitudedecaygluongluonNMSSM (double m1, double g, double gs, double alphas, double mWboson, DoubleMatrix & CPOMix, double beta, double mtop, double mbottom, double mcharm, double lam, int higgs, bool QCD);
235  double higgsAamplitudedecayZgammaNMSSM (double m1, double g, double gp, double alpha, double mWboson, double mZboson, DoubleMatrix & CPOMix, double beta, double mtop, double mbottom, double mcharm, double mch1, double mch2, double thetaL, double thetaR, double lam, int higgs);
236  double higgsCPevenamplitudedecaygammagammaNMSSM(double m1, double mtop, double mbottom, double mcharm, double mtau, double mWboson, double mHpm, double mchar1, double mchar2, double mscharmL, double mscharmR, double mstop1, double mstop2, double msstrangeL, double msstrangeR, double msbottom1, double msbottom2, double msmuonL, double msmuonR, double mstau1, double mstau2, DoubleMatrix & CPEMix, double beta, double g, double gp, double alpha, double thetat, double thetab, double thetatau, double thetaL, double thetaR, double At, double Ab, double Atau, double mu, double mueff, double lam, double kappa, double Alambda, int higgs);
237  double higgsCPevenamplitudedecaygluongluonNMSSM(double m1, double mtop, double mbottom, double mcharm, double mWboson, double mscharmL, double mscharmR, double mstop1, double mstop2, double msstrangeL, double msstrangeR, double msbottom1, double msbottom2, double msupL, double msupR, double msdownL, double msdownR, double runmt, double runmb, DoubleMatrix & CPEMix, double beta, double g, double gp, double gs, double alphas, double thetat, double thetab, double thetaL, double thetaR, double At, double Ab, double mu, double mueff, double lam, double kappa, double Alambda, int higgs, bool QCD);
238  double higgshamplitudedecayZgammaNMSSM (double m1, double g, double gp, double alpha, double mWboson, double mZboson, double mHpm, DoubleMatrix & CPEMix, double beta, double mtop, double mbottom, double mcharm, double mch1, double mch2, double thetaL, double thetaR, double lam, double kappa, double Alambda, double greekmu, double mueff, int higgs);
239  double higgshamplitudedecayneutineutjNMSSM (double m1, double mneuti, double mneutj, double g, double gp, DoubleMatrix & CPEMix, DoubleMatrix & mixNeut, double lam, double kappa, int neuti, int neutj, int higgs);
240  double higgsCPevenamplitudedecayAANMSSM(double m1, double mA1, double mA2, double mWboson, double runmt, double runmb, double g, double gp, double beta, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, double lam, double kappa, double Alambda, double Akappa, double mueff, int higgs, int pseudoscalar1, int pseudoscalar2);
241  double higgsCPevenamplitudedecaypseudoscalarZNMSSM (double m1, double mA, double mZboson, double g, double gp, double beta, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, int higgs, int pseudoscalar);
242  double higgsCPevenamplitudedecayHpHmNMSSM (double m1, double mHpm, double mWboson, double g, double gp, double mtop, double mbottom, double beta, double lam, double mueff, double kappa, double Alambda, DoubleMatrix & CPEMix, int higgs);
243  double higgsCPevenamplitudedecayhhorhHorHH(double m1, double mh1, double mh2, double g, double gp, double runmw, double beta, double lam, double Alambda, double kappa, double Akappa, double mueff, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, int higgs1, int higgs2);
244  double higgsCPevenamplitudedecayhhorhHorHHNMSSM(double m1, double mh1, double mh2, double g, double gp, double mWboson, double mtop, double mbottom, double beta, double lam, double Alambda, double kappa, double Akappa, double mueff, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, int higgs1, int higgs2, int starthiggs);
245  double higgsA2amplitudedecayA1CPevenNMSSM(double m1, double mA1, double mh, double mWboson, double runmt, double runmb, double g, double gp, double beta, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, double lam, double kappa, double Alambda, double Akappa, double mueff, int higgs);
246  double higgsCPevenamplitudedecayWHpmNMSSM (double m1, double mWboson, double mHpm, double beta, double g, DoubleMatrix & CPEMix, int higgs);
247  double higgsCPevenamplitudedecaystopistopiNMSSM (double m1, double mstopi, double thetat, double runmt, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double At, double mueff, double lam, int stop, int higgs);
248  double higgsCPevenamplitudedecaystopistopjNMSSM (double m1, double mstopi, double mstopj, double thetat, double runmt, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double At, double mueff, double lam, int higgs);
249  double higgsCPevenamplitudedecaysbottomisbottomiNMSSM (double m1, double msbottomi, double thetab, double runmb, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double Ab, double mueff, double lam, int sbottom, int higgs);
250  double higgsCPevenamplitudedecaysbottomisbottomjNMSSM (double m1, double msbottomi, double msbottomj, double thetab, double runmb, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double Ab, double mueff, double lam, int higgs);
251  double higgsCPevenamplitudedecaystauistauiNMSSM (double m1, double mstaui, double thetatau, double runmtau, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double Atau, double mueff, double lam, int stau, int higgs);
252  double higgsCPevenamplitudedecaystauistaujNMSSM (double m1, double mstaui, double mstauj, double thetatau, double runmtau, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, double Atau, double mueff, double lam, int higgs);
253  double higgsCPevenamplitudedecaysnusnuNMSSM (double m1, double msnu, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, int higgs);
254 
255  double higgsesamplitudedecaygammagammatotal(double m1, double g, double gprime, double alphaEmrun, double mWboson, double polemw, double alpha, double beta, double mtop, double mbottom, double mcharm, double mtau, double mHpm, double mstop1, double mstop2, double msbottom1, double msbottom2, double mstau1, double mstau2, double mchargino1, double mchargino2, double thetaL, double thetaR, double thetat, double thetab, double thetatau, double greekmu, double Atop, double Abottom, double Atau, char higgstype);
257 void slhaDecays(ostream & fout, vector<Particle> & decayTable, bool outputPartialWidths);
258 DoubleVector hggQCDcorrections(double amplitudeW, double alphas, int Nf, char higgs, double prefactor, double SMtotr, double SMtoti, double sqtotr, double sqtoti);
259  DoubleVector higgsmatrixelementgammagammaviatops (double m1, double mtop, double alpha, double beta, char higgstype);
260  DoubleVector higgsmatrixelementgammagammaviastops (double m1, double mstop1, double mstop2, double mtop, double mbottom, double mWboson, double thetat, double g, double gprime, double alpha, double beta, double greekmu, double Atop, double Abottom, char higgstype);
261  DoubleVector higgsmatrixelementgammagammaviabottoms (double m1, double mbottom, double alpha, double beta, char higgstype);
262  DoubleVector higgsmatrixelementgammagammaviasbottoms (double m1, double msbottom1, double msbottom2, double mbottom, double mtop, double mWboson, double thetab, double g, double gprime, double alpha, double beta, double Atop, double Abottom, double greekmu, char higgstype);
263  DoubleVector higgsmatrixelementgammagammaviastaus (double m1, double mstau1, double mstau2, double mtau, double mWboson, double thetatau, double g, double gprime, double alpha, double beta, double greekmu, double Atau, char higgstype);
264  DoubleVector higgsmatrixelementgammagammaviaWbosons (double m1, double mWboson, double alpha, double beta, double g, double gprime, char higgstype);
265  DoubleVector higgsmatrixelementgammagammaviaHpms (double m1, double mHpm, double mWboson, double alpha, double beta, double g, double gprime, char higgstype);
266  DoubleVector higgsmatrixelementgammagammaviachargino1s (double m1, double mchargino1, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype);
267  DoubleVector higgsmatrixelementgammagammaviachargino2s (double m1, double mchargino2, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype);
268  DoubleVector higgsmatrixelementgammagammaviacharms (double m1, double mcharm, double alpha, double beta, char higgstype);
269  DoubleVector higgsmatrixelementgammagammaviataus (double m1, double mtau, double alpha, double beta, char higgstype);
270 #endif
double higgsHamplitudedecaystau1stau2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for H->stau1 antistau2 or H->stau2 antistau1 as these widths have the same v...
Definition: twoBodyDecays.cpp:3635
DoubleVector higgsmatrixelementgammagammaviacharms(double m1, double mcharm, double alpha, double beta, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4222
double higgsesamplitudedecaygammagammatotal(double m1, double g, double gprime, double alphaEmrun, double mWboson, double polemw, double alpha, double beta, double mtop, double mbottom, double mcharm, double mtau, double mHpm, double mstop1, double mstop2, double msbottom1, double msbottom2, double mstau1, double mstau2, double mchargino1, double mchargino2, double thetaL, double thetaR, double thetat, double thetab, double thetatau, double greekmu, double Atop, double Abottom, double Atau, char higgstype)
function that calculates the partial width for h->gamma gamma at 1-loop (forbidden at tree-level) ...
Definition: twoBodyDecays.cpp:3908
double higgsAamplitudedecayZgammaNMSSM(double m1, double g, double gp, double alpha, double mWboson, double mZboson, DoubleMatrix &CPOMix, double beta, double mtop, double mbottom, double mcharm, double mch1, double mch2, double thetaL, double thetaR, double lam, int higgs)
Definition: threeBodyDecays.cpp:4275
DoubleVector higgsHplusamplitudedecaysquarksquarkmix(double m1, double m2, double m3, double g, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double thetatop, double thetabottom)
calculates partial width for Hplus->squarki antisquarkj where i,j are each 1/2 so mixing included...
Definition: twoBodyDecays.cpp:3866
double squarkRamplitudedecayneutralino(double m1, double m2, double m3, double g, double gprime, DoubleMatrix &mixNeut, int neutralino, int uord)
neutralino takes values of 1, 2, 3, or 4 and denotes which neutralino mass eigenstate we decay into...
Definition: twoBodyDecays.cpp:232
DoubleVector higgsHplusamplitudedecaysquarksquark(double m1, double m2, double m3, double g, double beta, double mWboson, double mup, double mdown, double greekmu, double Aup, double Adown)
Definition: twoBodyDecays.cpp:3834
Two loop corrections to the bottom mass. This file has been generated at Fri 7 Apr 2017 21:00:53 with...
double higgshamplitudedecay2sleptonsamehand(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl)
calculates partial width for h->slepton slepton with no mixing and sleptons of same handedness...
Definition: twoBodyDecays.cpp:2823
double neutralinoamplitudedecayphotongravitino(double m1, double mgrav, double MPlreduced, DoubleMatrix &mixNeut, double g, double gp, int neutralino, int gravonoff, int neutNLSP)
Function that calculates the decays to gravitino and photon of a neutralino if int gravonoff == 1...
Definition: twoBodyDecays.cpp:5187
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35
double neutralinoamplitudedecayneutrinosneutrinoL(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
m2 must be neutrino mass (i.e. 0)
Definition: twoBodyDecays.cpp:1364
double higgshamplitudedecay2sleptonsamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sl)
Definition: twoBodyDecays.cpp:2476
double charginoToNeutralino1pion(const MssmSoftsusy *m)
First, do hadronic decays.
Definition: twoBodyDecays.cpp:14
double charginoamplitudedecayquarksquarkL(double m1, double m2, double m3, double g, double theta, int chargino)
quark mass is m2
Definition: twoBodyDecays.cpp:897
void slhaDecays(ostream &fout, vector< Particle > &decayTable, bool outputPartialWidths)
outputs a decay table in SLHA format
Definition: twoBodyDecays.cpp:5243
double higgsAamplitudedecaysfermionsNMSSM(double m1, double m2, double m3, double g, double mWboson, double mf, double Asf, double beta, double lam, double mueff, DoubleMatrix &CPOMix, char uord, int pseudoscalar)
mf here should be runmf and mWboson should be runmw
Definition: twoBodyDecays.cpp:3706
DoubleVector higgsHamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix &CPEMix, bool nmssmIsIt)
Function that calculates the Heavy Higgs decays to two vector bosons, assuming both are on-shell...
Definition: twoBodyDecays.cpp:4818
double snutauamplitudedecaynutauneutralinoNMSSM(double m1, double mneut, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: threeBodyDecays.cpp:5894
double higgsAamplitudedecaysamecharginoNMSSM(double m1, double m2, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix &CPOMix, int chargino, int pseudoscalar)
Definition: twoBodyDecays.cpp:2020
double higgshamplitudedecaysbottom1sbottom2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->sbottom1 antisbottom2 or sbottom2 antisbottom1 as they are the same ...
Definition: twoBodyDecays.cpp:3164
double higgsHamplitudedecaystop2stop2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->stop2 antistop2
Definition: twoBodyDecays.cpp:3242
DoubleVector higgsmatrixelementgammagammaviachargino1s(double m1, double mchargino1, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4473
double higgsesamplitudedecayZbosonphotontotal(double m1, double mZboson, double g, double gprime, double alphaEmrun, double polemw, double runmw, double alpha, double beta, double mtop, double mbottom, double mcharm, double mstrange, double mstop1, double mstop2, double msbottom1, double msbottom2, double mHplus, double thetat, double thetab, double greekmu, double Atop, double Abottom, char higgstype)
function that calculates the partial width for h->Z gamma at 1-loop (forbidden at tree-level) ...
Definition: twoBodyDecays.cpp:4952
double chargino2amplitudedecaychargino1neutHiggs(double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, double alpha, char phi)
m2 must be neutral Higgs mass
Definition: twoBodyDecays.cpp:1225
double higgsHamplitudedecaystop1stop2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->stop1 antistop2 or stop2 antistop1 as they are the same ...
Definition: twoBodyDecays.cpp:3281
double higgsCPevenamplitudedecaygluongluonNMSSM(double m1, double mtop, double mbottom, double mcharm, double mWboson, double mscharmL, double mscharmR, double mstop1, double mstop2, double msstrangeL, double msstrangeR, double msbottom1, double msbottom2, double msupL, double msupR, double msdownL, double msdownR, double runmt, double runmb, DoubleMatrix &CPEMix, double beta, double g, double gp, double gs, double alphas, double thetat, double thetab, double thetaL, double thetaR, double At, double Ab, double mu, double mueff, double lam, double kappa, double Alambda, int higgs, bool QCD)
Definition: threeBodyDecays.cpp:4405
double higgsAamplitudedecayhiggsHZboson(double m1, double m2, double m3, double g, double gp, double alpha, double beta)
calculates partial width for A->HZ, m2 must be Zboson mass
Definition: twoBodyDecays.cpp:2324
double higgsHamplitudedecay2squarkdiffhand(double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mupq, double mdownq, double greekmu, double Aup, double Adown, int sq)
calculates partial width for H->squark squark with no mixing and squarks of different handedness...
Definition: twoBodyDecays.cpp:2786
DoubleVector higgsH3amplitudedecayVVNMSSM(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix &CPEMix, bool nmssmIsIt)
Function that calculates the Heavy Higgs decays to two vector bosons, assuming both are on-shell...
Definition: twoBodyDecays.cpp:4885
double squark1amplitudedecaycharginoW1mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double runmt, double runmb, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:145
double higgsHplusamplitudedecayneutralinochargino(double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix &mixNeut, int neutralino, int chargino)
Calculates the partial width for decays of the charged H+ higgs to a neutralino and a chragino Wtilda...
Definition: twoBodyDecays.cpp:3767
double neutralinoamplitudedecayphigravitino(double m1, double mphi, double mgrav, double MPlreduced, DoubleMatrix &mixNeut, double alpha, double beta, int neutralino, int gravonoff, char phi, int neutNLSP)
Function that calculates the decays to gravitino and Z boson of a neutralino if int gravonoff == 1...
Definition: twoBodyDecays.cpp:5220
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...
double squark3amplitudedecaychargedHiggssquark3(double m1, double m2, double m3, double g, double mWboson, double beta, double thetat, double thetab, double greekmu, double At, double Ab, double mt, double mb, int t1or2, int b1or2)
Definition: twoBodyDecays.cpp:392
double neutralinoamplitudedecayleptonsleptonLorR(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, char LorR, int neutralino)
m2 must be lepton mass
Definition: twoBodyDecays.cpp:1323
double higgsAamplitudedecayhiggshorHZbosonNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double thetaA, DoubleMatrix &CPEMix, int pseudoscalar, int higgs)
m2 must be Z mass
Definition: twoBodyDecays.cpp:2349
double charginoamplitudedecayquarksquarkmix(double m1, double m2, double m3, double g, double theta, double thetaL, double thetaR, double beta, double runmt, double runmb, double mWboson, int chargino, int upordowntypesquark, int oneortwo)
quark mass is m2
Definition: twoBodyDecays.cpp:931
double higgsHamplitudedecaysbottom1sbottom1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->sbottom1 antisbottom1
Definition: twoBodyDecays.cpp:3320
double squark1amplitudedecaycharginoW2mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double runmt, double runmb, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:162
DoubleVector higgsmatrixelementgammagammaviachargino2s(double m1, double mchargino2, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4525
double higgsHamplitudedecaystau2stau2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for H->stau2 antistau2
Definition: twoBodyDecays.cpp:3595
double higgsHplusamplitudedecayWbosonhiggsh(double m1, double m2, double m3, double g, double alpha, double beta)
Calculates the partial width for a charged Higgs boson H+ to decay to a Wboson and a neutral light sc...
Definition: twoBodyDecays.cpp:3808
double higgsAamplitudedecayquarkantiquark(double m1, double m2, double g, double beta, double mWboson, int uord, bool QCD, double alphas)
uord indicates if it&#39;s an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1759
double higgshamplitudedecay2squarksamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sq)
takes MSSM amplitude and applies relevant NMSSM prefactor to get NMSSM amplitude
Definition: twoBodyDecays.cpp:2436
double higgshamplitudedecayhiggsAZboson(double m1, double m2, double m3, double g, double gp, double alpha, double beta)
calculates partial width for h->AZ, m2 must be Zboson mass
Definition: twoBodyDecays.cpp:2246
double charginoamplitudedecayWbosonneutralino(double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix &mixNeut, int chargino, int neutralino)
m2 must be Wboson mass
Definition: twoBodyDecays.cpp:1117
double higgsHamplitudedecaystau1stau1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for H->stau1 antistau1
Definition: twoBodyDecays.cpp:3555
double higgsH3amplitudedecay2squarksamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sq)
Definition: twoBodyDecays.cpp:2745
double stauamplitudedecaynutauchargino(double m1, double m2, double m3, double g, double mWboson, double theta, double thetaL, double beta, double mtau, int oneortwo, int chargino)
Definition: twoBodyDecays.cpp:682
MssmSusy contains all SUSY couplings and tan beta, as well as their beta functions.
double gluinoamplitudedecaygravitino(double m1, double mgrav, double MPlreduced, int gravonoff, int gluNLSP)
Decays to gravitinos.
Definition: twoBodyDecays.cpp:5161
double squark3amplitudedecayneutralino(double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix &mixNeut, double g, double gp, double runmq, int squark, int oneortwo, int neutralino)
Definition: twoBodyDecays.cpp:254
double squark2amplitudedecaycharginoW2mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double runmt, double runmb, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:196
double higgshamplitudedecayneutineutjNMSSM(double m1, double mneuti, double mneutj, double g, double gp, DoubleMatrix &CPEMix, DoubleMatrix &mixNeut, double lam, double kappa, int neuti, int neutj, int higgs)
Definition: threeBodyDecays.cpp:4582
double neutralinoamplitudedecaystautau(double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix &mixNeut, double g, double gp, int oneortwo, int neutralino)
m2 must be tau (i.e. lepton) mass
Definition: twoBodyDecays.cpp:1465
DoubleVector higgsmatrixelementgammagammaviataus(double m1, double mtau, double alpha, double beta, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4263
double higgslorHamplitudedecayquarkantiquark(double m1, double m2, double g, double alpha, double beta, double mWboson, int uord, char lorH, DoubleMatrix &CPEMix, bool NMSSMmodel, bool QCD, double alphas)
uord indicates if it&#39;s an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1659
double squarkLamplitudedecayneutralino(double m1, double m2, double m3, double g, double gprime, DoubleMatrix &mixNeut, int neutralino, int uord)
neutralino takes values of 1, 2, 3, or 4 and denotes which neutralino mass eigenstate we decay into...
Definition: twoBodyDecays.cpp:215
double higgsphiamplitudedecaysamechargino(double m1, double m2, double g, double thetaL, double thetaR, double alpha, double beta, int chargino, char phi)
phi tells it whether a "h", "H" or "A" is decaying
Definition: twoBodyDecays.cpp:1931
double higgsAamplitudedecaygluongluonNMSSM(double m1, double g, double gs, double alphas, double mWboson, DoubleMatrix &CPOMix, double beta, double mtop, double mbottom, double mcharm, double lam, int higgs, bool QCD)
Definition: threeBodyDecays.cpp:4222
double higgsAamplitudedecaysfermions(double m1, double m2, double m3, double g, double mWboson, double mf, double greekmu, double Asf, double beta, char uord)
calculates partial width for A->sfermion1 sfermion2, these sfermions must be of the same type...
Definition: twoBodyDecays.cpp:3675
double neutralinoamplitudedecayquarksquarkLorR(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int uordtype, char LorR, int neutralino)
m2 must be quark mass
Definition: twoBodyDecays.cpp:1272
double higgsHamplitudedecaystop1stop1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->stop1 antistop1
Definition: twoBodyDecays.cpp:3203
double higgshamplitudedecaystop1stop2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->stop1 antistop2 or stop2 antistop1 as they are the same ...
Definition: twoBodyDecays.cpp:3047
double squark2amplitudedecaycharginoW1mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double runmt, double runmb, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:180
DoubleVector higgsmatrixelementgammagammaviastaus(double m1, double mstau1, double mstau2, double mtau, double mWboson, double thetatau, double g, double gprime, double alpha, double beta, double greekmu, double Atau, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4303
double higgsH3amplitudedecay2sleptonsamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sl)
Definition: twoBodyDecays.cpp:2552
double higgshamplitudedecayZgammaNMSSM(double m1, double g, double gp, double alpha, double mWboson, double mZboson, double mHpm, DoubleMatrix &CPEMix, double beta, double mtop, double mbottom, double mcharm, double mch1, double mch2, double thetaL, double thetaR, double lam, double kappa, double Alambda, double greekmu, double mueff, int higgs)
Definition: threeBodyDecays.cpp:4509
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
double charginoamplitudedecayleptonsleptonL(double m1, double m2, double m3, double g, double thetaLorR, int chargino)
lepton mass is m2, use thetaR for decays to sneutrinos and leptons, use thetaL for decays to sleptons...
Definition: twoBodyDecays.cpp:986
double higgshamplitudedecayAA(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for h->AA
Definition: twoBodyDecays.cpp:2146
DoubleVector hggQCDcorrections(double amplitudeW, double alphas, int Nf, char higgs, double prefactor, double SMtotr, double SMtoti, double sqtotr, double sqtoti)
Function for QCD corrections to h->gg.
Definition: twoBodyDecays.cpp:4661
DoubleVector higgsmatrixelementgammagammaviaWbosons(double m1, double mWboson, double alpha, double beta, double g, double gprime, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4394
double stau2amplitudedecaystau1Zboson(double m1, double m2, double m3, double g, double gp, double thetatau)
m3 must be mz here
Definition: twoBodyDecays.cpp:834
double higgshamplitudedecaystau1stau1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for h->stau1 antistau1
Definition: twoBodyDecays.cpp:3437
double squarkamplitudedecaygravitino(double m1, double mgrav, double mquark, double MPlreduced, int gravonoff, int squNLSP)
Function that calculates the decays to gravitinos of squark if int gravonoff == 1, otherwise decays to gravitinos are off.
Definition: twoBodyDecays.cpp:5174
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:211
double higgshamplitudedecaystop2stop2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->stop2 antistop2
Definition: twoBodyDecays.cpp:3008
double higgsCPevenamplitudedecaygammagammaNMSSM(double m1, double mtop, double mbottom, double mcharm, double mtau, double mWboson, double mHpm, double mchar1, double mchar2, double mscharmL, double mscharmR, double mstop1, double mstop2, double msstrangeL, double msstrangeR, double msbottom1, double msbottom2, double msmuonL, double msmuonR, double mstau1, double mstau2, DoubleMatrix &CPEMix, double beta, double g, double gp, double alpha, double thetat, double thetab, double thetatau, double thetaL, double thetaR, double At, double Ab, double Atau, double mu, double mueff, double lam, double kappa, double Alambda, int higgs)
Definition: threeBodyDecays.cpp:4322
double squark32amplitudedecayneutralHiggssquark3(double m1, double m2, double m3, double g, double gp, double mWboson, double beta, double alpha, double thetat, double thetab, double greekmu, double At, double Ab, double mt, double mb, int torb, char phi)
Definition: twoBodyDecays.cpp:450
Code calculates decay modes and prints out an SLHA format file with them in. For R-parity conserving ...
DoubleVector higgsmatrixelementgammagammaviaHpms(double m1, double mHpm, double mWboson, double alpha, double beta, double g, double gprime, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4431
DoubleVector higgsmatrixelementgammagammaviastops(double m1, double mstop1, double mstop2, double mtop, double mbottom, double mWboson, double thetat, double g, double gprime, double alpha, double beta, double greekmu, double Atop, double Abottom, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:3997
double stauamplitudedecaysnustauWboson(double m1, double m2, double m3, double g, double thetatau, int oneortwo)
m3 must be mw here
Definition: twoBodyDecays.cpp:802
double sleptonamplitudedecayleptonneutralinoNMSSM(double m1, double ml, double mneut, double g, double gp, DoubleMatrix &mixNeut, char uord, char LorR, int neut)
Definition: threeBodyDecays.cpp:5597
NmssmSoftsusy object contains all the NMSSM functionality: SUSY breaking and preserving parameters...
double higgsesamplitudedecaygluongluontotal(double m1, double g, double gs, double gprime, double mWboson, double alpha, double beta, double mtop, double mbottom, double mcharm, double mstop1, double mstop2, double msbottom1, double msbottom2, double thetat, double thetab, double greekmu, double Atop, double Abottom, double mstrange, double mscharmL, double mscharmR, double msstrangeL, double msstrangeR, double Acharm, double Astrange, double mup, double mdown, double msupL, double msupR, double msdownL, double msdownR, double Aup, double Adown, char higgstype, bool QCD)
function that calculates the partial width for h->gluon gluon at 1-loop (forbidden at tree-level) ...
Definition: twoBodyDecays.cpp:4582
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...
double charginoamplitudedecaystaunutau(double m1, double m2, double m3, double g, double thetaL, double thetaR, double thetatau, double beta, double mWboson, double mtau, int oneortwo, int chargino)
m2 must be nutau mass (i.e. 0)
Definition: twoBodyDecays.cpp:1057
double squark3amplitudedecaysquark3Wboson(double m1, double m2, double m3, double g, double thetat, double thetab, int m1torb, int m1oneortwo, int m3torb, int m3oneortwo)
m1torb tells the function if the initial squark is a stop or sbottom, m1oneortwo tells it whether it&#39;...
Definition: twoBodyDecays.cpp:318
double neutralinoamplitudedecayZgravitino(double m1, double mZ, double mgrav, double MPlreduced, DoubleMatrix &mixNeut, double g, double gp, double beta, int neutralino, int gravonoff, int neutNLSP)
Function that calculates the decays to gravitino and Z boson of a neutralino if int gravonoff == 1...
Definition: twoBodyDecays.cpp:5203
double higgsAamplitudedecayhiggshZboson(double m1, double m2, double m3, double g, double gp, double alpha, double beta)
calculates partial width for A->hZ, m2 must be Zboson mass
Definition: twoBodyDecays.cpp:2297
double charginoamplitudedecayHminusneutralino(double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, DoubleMatrix &mixNeut, int chargino, int neutralino)
m2 must be Hminus mass
Definition: twoBodyDecays.cpp:1154
double neutralinoamplitudedecaysquark3quarkmix(double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix &mixNeut, double g, double gp, double runmq, int squark, int oneortwo, int neutralino)
m2 must be quark mass
Definition: twoBodyDecays.cpp:1393
double higgshamplitudedecay2squarkdiffhand(double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mupq, double mdownq, double greekmu, double Aup, double Adown, int sq)
calculates partial width for h->squark squark with no mixing and squarks of different handedness...
Definition: twoBodyDecays.cpp:2589
Soft SUSY breaking parameters.
double m1
decay global variable declarations
Definition: decays.cpp:14
double higgsHamplitudedecayhh(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for h->hh
Definition: twoBodyDecays.cpp:2171
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:77
DoubleVector higgsmatrixelementgammagammaviatops(double m1, double mtop, double alpha, double beta, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:3956
double higgshamplitudedecaystop1stop1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->stop1 antistop1
Definition: twoBodyDecays.cpp:2968
double higgsHamplitudedecayhiggsAZboson(double m1, double m2, double m3, double g, double gp, double alpha, double beta)
calculates partial width for H->AZ, m2 must be Zboson mass
Definition: twoBodyDecays.cpp:2271
double higgshamplitudedecaysbottom2sbottom2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->sbottom2 antisbottom2
Definition: twoBodyDecays.cpp:3125
double higgsHamplitudedecayHplusHminus(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for H->H+H-
Definition: twoBodyDecays.cpp:2221
double higgshamplitudedecay2sleptondiffhand(double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl)
calculates partial width for h->slepton slepton with no mixing and sleptons of different handedness...
Definition: twoBodyDecays.cpp:2862
double stausneutrinoamplitudedecaytauchargino(double m1, double m2, double m3, double g, double mWboson, double beta, double thetaL, double thetaR, int chargino)
Definition: twoBodyDecays.cpp:731
double higgshamplitudedecay2squarkdiffhandNMSSM(double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mq, double Aq, double mueff, double lam, DoubleMatrix &CPEMix, int sq, int higgs)
Definition: twoBodyDecays.cpp:2627
double higgsAamplitudedecaygammagammaNMSSM(double m1, double g, double gprime, double alpha, double mWboson, DoubleMatrix &CPOMix, double beta, double mtop, double mbottom, double mcharm, double mtau, double mch1, double mch2, double thetaL, double thetaR, double lam, int higgs)
Definition: threeBodyDecays.cpp:4180
double higgsHamplitudedecay2sleptonsamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sl)
Definition: twoBodyDecays.cpp:2514
double higgsHplusamplitudedecayquarkantiquark(double m1, double m2, double m3, double g, double mWboson, double beta, DoubleMatrix &VCKM, int quark, int antiquark)
calculates partial width for H+ ->quark1 antiquark2, this is the same as H- -> antiquark1 quark2...
Definition: twoBodyDecays.cpp:3740
double higgsHamplitudedecay2sleptonsamehand(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl)
calculates partial width for H->slepton slepton with no mixing and sleptons of same handedness...
Definition: twoBodyDecays.cpp:2895
double charginoamplitudedecaysnutautau(double m1, double m2, double m3, double g, double thetaL, double thetaR, double beta, double mWboson, int chargino)
m2 must be tau mass
Definition: twoBodyDecays.cpp:1020
double sneutrinoamplitudedecayneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: twoBodyDecays.cpp:572
double higgshamplitudedecay2squarksamehand(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mupq, double mdownq, int sq)
calculates partial width for h->squark squark with no mixing and squarks of same handedness, therefore for first two generations, int sq tells the function which squarks it&#39;s decaying into - uL, dL, uR, dR for sq = 1,2,3,4 respectively
Definition: twoBodyDecays.cpp:2394
double sleptonamplitudedecaychargino(double m1, double m2, double m3, double g, double theta, int chargino)
for both sleptonL decays to charginos + neutrinos and for sneutrino decays to lepton + charginos - ju...
Definition: twoBodyDecays.cpp:595
Definitions of particles and container for decay widths.
double stauamplitudedecaytauneutralino(double m1, double m2, double m3, double g, double gp, double mWboson, DoubleMatrix &mixNeut, double theta, double beta, int oneortwo, int neutralino)
Definition: twoBodyDecays.cpp:626
double neutralinoamplitudedecaycharginoWboson(double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix &mixNeut, int neutralino, int chargino)
m2 must be W boson mass (expect to be pole mass 80.4...), neutralino is i in T&B whilst chargino is j...
Definition: twoBodyDecays.cpp:1509
DoubleVector higgsmatrixelementgammagammaviabottoms(double m1, double mbottom, double alpha, double beta, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4090
double higgsphiamplitudedecaydifchargino(double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double beta, char phi)
phi tells it whether a "h", "H" or "A" is decaying
Definition: twoBodyDecays.cpp:2076
double higgshamplitudedecaysbottom1sbottom1(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for h->sbottom1 antisbottom1
Definition: twoBodyDecays.cpp:3086
DoubleVector higgshamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix &CPEMix, bool nmssmIsIt)
Function that calculates the light Higgs decays to two vector bosons, both on-shell or one off-shell...
Definition: twoBodyDecays.cpp:4751
double neutralinoamplitudedecayneutralinoZboson(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int ineutralino, int fneutralino)
m2 must be Z mass, ineutralino is i in T&B whilst fneutralino is j
Definition: twoBodyDecays.cpp:1588
double higgsAamplitudedecayquarkantiquarkNMSSM(double m1, double m2, double beta, DoubleMatrix &CPOMix, int uord, int higgs, bool QCD, double alphas)
uord indicates if it&#39;s an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1805
bool outputPartialWidths
Default: don&#39;t output partial widths in decays.
Definition: def.cpp:61
double stau2amplitudedecaystau1phi(double m1, double m2, double m3, double g, double gp, double thetatau, double beta, double alpha, double mWboson, double mtau, double greekmu, double Atau, char phi)
Definition: twoBodyDecays.cpp:858
double neutralinoamplitudedecayneutralinoneutHiggs(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, double mixingangle, int ineutralino, int fneutralino, char phi)
m2 must be phi mass, ineutralino is i in T&B whilst fneutralino is j
Definition: twoBodyDecays.cpp:1615
double higgsCPevenamplitudedecayAANMSSM(double m1, double mA1, double mA2, double mWboson, double runmt, double runmb, double g, double gp, double beta, DoubleMatrix &CPEMix, DoubleMatrix &CPOMix, double lam, double kappa, double Alambda, double Akappa, double mueff, int higgs, int pseudoscalar1, int pseudoscalar2)
Definition: threeBodyDecays.cpp:4646
const double root2
used to flag diabolical problems
Definition: def.h:35
DoubleVector higgsmatrixelementgammagammaviasbottoms(double m1, double msbottom1, double msbottom2, double mbottom, double mtop, double mWboson, double thetab, double g, double gprime, double alpha, double beta, double Atop, double Abottom, double greekmu, char higgstype)
function that calculates the part of the matrix element that differs depending on the loop particles...
Definition: twoBodyDecays.cpp:4131
double higgsHamplitudedecayAA(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for H->AA
Definition: twoBodyDecays.cpp:2196
double gluinoamplitudedecaydgausscharginoqqpbarfirsttwogen(double mgluino, double mchargino, double mquark, double mquarkp, double msqL, double msqpL, double g, double thetaL, double thetaR, double alphas, int charg, bool onetothree)
Functions for 1->3 decays via dgauss:
Definition: threeBodyDecays.cpp:85
double higgshamplitudedecaystau2stau2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for h->stau2 antistau2
Definition: twoBodyDecays.cpp:3477
double stauamplitudedecaysnustauHminus(double m1, double m2, double m3, double g, double mWboson, double beta, double thetatau, double mtau, double greekmu, double Atau, int oneortwo)
Also does decay mode snustau to Hplus and stau1/2 just with m1, m2 and m3 permuted.
Definition: twoBodyDecays.cpp:766
double neutralinoamplitudedecaysfermionfermionfirst2genNMSSM(double mneut, double msf, double mf, double g, double gp, DoubleMatrix &mixNeut, int neut, char type, char LorR)
type indicates type of fermion, &#39;u&#39; for up type quark, &#39;d&#39; for down type quark, &#39;l&#39; for charged lepto...
Definition: threeBodyDecays.cpp:5322
double higgsHamplitudedecaysbottom1sbottom2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->sbottom1 antisbottom2 or sbottom2 antisbottom1 as they are the same ...
Definition: twoBodyDecays.cpp:3398
double higgsHamplitudedecaysbottom2sbottom2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtop, double mbottom, double greekmu, double Atop, double Abottom, double theta)
calculates partial width for H->sbottom2 antisbottom2
Definition: twoBodyDecays.cpp:3359
double higgsHamplitudedecay2squarksamehand(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mupq, double mdownq, int sq)
calculates partial width for H->squark squark with no mixing and squarks of same handedness, therefore for first two generations, int sq tells the function which squarks it&#39;s decaying into - uL, dL, uR, dR for sq = 1,2,3,4 respectively
Definition: twoBodyDecays.cpp:2661
double stausneutrinoamplitudedecaytauneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: twoBodyDecays.cpp:661
double higgsphiamplitudedecaysamecharginoNMSSM(double m1, double m2, double g, double thetaL, double thetaR, double lam, DoubleMatrix &CPEMix, int chargino, int higgs)
Definition: twoBodyDecays.cpp:1986
int theta(double a)
Standard theta function: 1 is a>0, 0 otherwise.
Definition: utils.cpp:25
double higgsHamplitudedecay2squarksamehandNMSSM(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mq, DoubleMatrix &CPEMix, int sq)
Definition: twoBodyDecays.cpp:2704
double higgsAamplitudedecayHpmWboson(double m1, double mWboson, double mHpm, double g, double thetaA, int pseudoscalar, bool nmssmIsIt)
Does A/A2 -> H+W- (note multiply by 2 if want to include H-W+ as well as H+W- in this) ...
Definition: threeBodyDecays.cpp:4606
double higgsHamplitudedecay2sleptondiffhand(double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl)
calculates partial width for H->slepton slepton with no mixing and sleptons of different handedness...
Definition: twoBodyDecays.cpp:2934
double neutralinoamplitudedecaycharginoHplus(double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix &mixNeut, int neutralino, int chargino)
m2 must be Hplus mass, neutralino is i in T&B whilst chargino is j
Definition: twoBodyDecays.cpp:1547
double chargino2amplitudedecaychargino1Zboson(double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR)
m2 must be Zboson mass
Definition: twoBodyDecays.cpp:1196
double higgshamplitudedecaystau1stau2(double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mtau, double greekmu, double Atau, double theta)
calculates partial width for h->stau1 antistau2 or h->stau2 antistau1 as these widths have the same v...
Definition: twoBodyDecays.cpp:3516
double higgsphiamplitudedecayneutralinoneutralino(double m1, double m2, double m3, double g, double tanthetaW, double mixingangle, DoubleMatrix &mixNeut, int ineutralino, int fneutralino, char phi)
phi tells it whether a "h", "H" or "A" is decaying
Definition: twoBodyDecays.cpp:1843