softsusy is hosted by Hepforge, IPPP Durham
SOFTSUSY  4.1
twoBodyDecays.h
Go to the documentation of this file.
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 //double charginoToNeutralino21pion(const MssmSoftsusy * m);
39 
40 double gluinoamplitudedecay (double m1, double m2, double m3,
41  double alphastrong);
42 double gluinoamplitudedecaymix (double m1, double m2, double m3,
43  double alphastrong, double squarkmix,
44  double theta);
45 double squarkamplitudedecaygluino (double m1, double m2, double m3,
46  double alphastrong);
47 double squarkamplitudedecaygluinomix (double m1, double m2, double m3,
48  double alphastrong, double squarkmix,
49  double theta);
50 double squarkamplitudedecaycharginoW1 (double m1, double m2, double m3,
51  double g, double gamma);
52 double squarkamplitudedecaycharginoW2 (double m1, double m2, double m3,
53  double g, double gamma);
55 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
56  double theta, double beta, double mWboson, double runmt, double runmb,
57  double torb);
59 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
60  double theta, double beta, double mWboson, double runmt, double runmb,
61  double torb);
63 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
64  double theta, double beta, double mWboson, double mup, double mdown,
65  double torb);
67 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
68  double theta, double beta, double mWboson, double mup, double mdown,
69  double torb);
71 (double m1, double m2, double m3, double g, double gprime,
72  DoubleMatrix & mixNeut, int neutralino, int uord );
74 (double m1, double m2, double m3, double g, double gprime,
75  DoubleMatrix & mixNeut, int neutralino, int uord );
77 (double m1, double m2, double m3, double mWboson, double theta, double beta,
78  DoubleMatrix & mixNeut, double g, double gp, double runmt, int squark,
79  int oneortwo, int neutralino);
81 (double m1, double m2, double m3, double g, double thetat, double thetab,
82  int m1torb, int m1oneortwo, int m3torb, int m3oneortwo);
84 (double m1, double m2, double m3, double g, double mWboson, double beta,
85  double thetat, double thetab, double greekmu, double At, double Ab,
86  double mt, double mb, int t1or2, int b1or2);
87  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);
88  double squark32amplitudedecaysquark3Zboson (double m1, double m2, double m3, double g, double gp, double theta);
89  double sleptonamplitudedecayleptonneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, char LorR, int neutralino);
90  double sneutrinoamplitudedecayneutrinoneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
91  double sleptonamplitudedecaychargino (double m1, double m2, double m3, double g, double theta, int chargino);
92  double stauamplitudedecaytauneutralino (double m1, double m2, double m3, double g, double gp, double mWboson, DoubleMatrix & mixNeut, double theta, double beta, int oneortwo, int neutralino);
93  double stausneutrinoamplitudedecaytauneutrinoneutralino (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
94  double stauamplitudedecaynutauchargino (double m1, double m2, double m3, double g, double mWboson, double theta, double thetaL, double beta, double mtau, int oneortwo, int chargino);
95  double stausneutrinoamplitudedecaytauchargino (double m1, double m2, double m3, double g, double mWboson, double beta, double thetaL, double thetaR, int chargino);
96  double stauamplitudedecaysnustauHminus (double m1, double m2, double m3, double g, double mWboson, double beta, double thetatau, double mtau, double mu, double Atau, int oneortwo);
97  double stauamplitudedecaysnustauWboson (double m1, double m2, double m3, double g, double thetatau, int oneortwo);
98  double stau2amplitudedecaystau1Zboson (double m1, double m2, double m3, double g, double gp, double thetatau);
99  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);
100  double charginoamplitudedecayquarksquarkL (double m1, double m2, double m3, double g, double theta, int chargino);
101  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);
102  double charginoamplitudedecayleptonsleptonL (double m1, double m2, double m3, double g, double thetaLorR, int chargino);
103  double charginoamplitudedecaysnutautau (double m1, double m2, double m3, double g, double thetaL, double thetaR, double beta, double mWboson, int chargino);
104  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);
105  double charginoamplitudedecayWbosonneutralino (double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int chargino, int neutralino);
106  double charginoamplitudedecayHminusneutralino (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, DoubleMatrix & mixNeut, int chargino, int neutralino);
107  double chargino2amplitudedecaychargino1Zboson (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR);
108  double chargino2amplitudedecaychargino1neutHiggs (double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR, double beta, double alpha, char phi);
109  double neutralinoamplitudedecayquarksquarkLorR (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int uordtype , char LorR, int neutralino);
110  double neutralinoamplitudedecayleptonsleptonLorR (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, char LorR, int neutralino);
111  double neutralinoamplitudedecayneutrinosneutrinoL (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
112  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);
113  double neutralinoamplitudedecaystautau (double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix & mixNeut, double g, double gp, int oneortwo, int neutralino);
114  double neutralinoamplitudedecaycharginoWboson (double m1, double m2, double m3, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
115  double neutralinoamplitudedecaycharginoHplus (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
116  double neutralinoamplitudedecayneutralinoZboson (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, int ineutralino, int fneutralino);
117  double neutralinoamplitudedecayneutralinoneutHiggs (double m1, double m2, double m3, double g, double gp, DoubleMatrix & mixNeut, double mixingangle, int ineutralino, int fneutralino, char phi);
118 
119  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);
120  double higgsAamplitudedecayquarkantiquark (double m1, double m2, double g, double beta, double mWboson, int uord, bool QCD, double alphas);
121  double higgsAamplitudedecayquarkantiquarkNMSSM (double m1, double m2, double beta, DoubleMatrix & CPOMix, int uord, int higgs, bool QCD, double alphas);
122  double higgsphiamplitudedecayneutralinoneutralino (double m1, double m2, double m3, double g, double tanthetaW, double mixingangle, DoubleMatrix & mixNeut, int ineutralino, int fneutralino, char phi);
123  double higgsphiamplitudedecaysamechargino (double m1, double m2, double g, double thetaL, double thetaR, double alpha, double beta, int chargino, char phi);
124  double higgsphiamplitudedecaydifchargino (double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double beta, char phi);
125  double higgshamplitudedecayAA (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
126  double higgsHamplitudedecayhh (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
127  double higgsHamplitudedecayAA (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
128  double higgsHamplitudedecayHplusHminus (double m1, double m2, double g, double gp, double alpha, double beta, double mWboson);
129  double higgshamplitudedecayhiggsAZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
130  double higgsHamplitudedecayhiggsAZboson (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 higgsAamplitudedecayhiggsHZboson (double m1, double m2, double m3, double g, double gp, double alpha, double beta);
133  double higgshamplitudedecay2squarksamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mupq, double mdownq, int sq);
134  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);
135  double higgsHamplitudedecay2squarksamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mupq, double mdownq, int sq);
136  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);
137  double higgshamplitudedecay2sleptonsamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl);
138  double higgshamplitudedecay2sleptondiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl);
139  double higgsHamplitudedecay2sleptonsamehand (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson, double mel, int sl);
140  double higgsHamplitudedecay2sleptondiffhand (double m1, double m2, double m3, double g, double alpha, double beta, double mWboson, double mel, double greekmu, double Aelectron, int sl);
141  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);
142  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);
143  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);
144  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);
145  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);
146  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);
147  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);
148  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);
149  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);
150  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);
151  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);
152  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);
153  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);
154  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);
155  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);
156  double higgshtestamplitudedecaystau1stau1 (double m1, double m2, double thetatau, double g, double gp, double mWboson, double alpha, double beta, double mtau, double greekmu, double Atau);
157  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);
158  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);
159  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);
160  double higgsAamplitudedecaysfermions (double m1, double m2, double m3, double g, double mWboson, double mf, double greekmu, double Asf, double beta, char uord);
161  double higgsHplusamplitudedecayquarkantiquark (double m1, double m2, double m3, double g, double mWboson, double beta, DoubleMatrix & VCKM, int quark, int antiquark);
162  double higgsHplusamplitudedecayneutralinochargino (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
163  double higgsHplusamplitudedecayneutralinocharginosusyhitway (double m1, double m2, double m3, double g, double gp, double beta, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
164  double higgsHplusamplitudedecayWbosonhiggsh (double m1, double m2, double m3, double g, double alpha, double beta);
165  DoubleVector higgsHplusamplitudedecaysquarksquark (double m1, double m2, double m3, double g, double beta, double mWboson, double mup, double mdown, double greekmu, double Aup, double Adown);
166  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);
167  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);
168  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);
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 higgsHamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix & CPEMix, bool nmssmIsIt);
172  DoubleVector higgsH3amplitudedecayVVNMSSM(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix & CPEMix, bool nmssmIsIt);
173  double higgsAamplitudedecayHpmWboson(double m1, double mWboson, double mHpm, double g, double thetaA, int pseudoscalar, bool nmssmIsIt);
174  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);
175  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);
176  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);
177  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);
178  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);
179  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);
180  double chargino2amplitudedecaychargino1CPevenhiggsNMSSM (double mchar2, double mchar1, double mh, double g, double lam, double thetaL, double thetaR, DoubleMatrix & CPEMix, int higgs);
181  double chargino2amplitudedecaychargino1CPoddhiggsNMSSM (double mchar2, double mchar1, double mA, double g, double lam, double thetaL, double thetaR, DoubleMatrix & CPOMix, int higgs);
182  double neutralinoamplitudedecaycharginoWNMSSM (double mneut, double mchar, double mWboson, double g, double thetaL, double thetaR, DoubleMatrix & mixNeut, int neutralino, int chargino);
183  double neutralinoamplitudedecayneutralinoZNMSSM (double mneuti, double mneutj, double mZboson, double g, double gp, DoubleMatrix & mixNeut, int neutralinoi, int neutralinoj);
184  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);
185  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);
186  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);
187  double neutralinoamplitudedecaysfermionfermionfirst2genNMSSM (double mneut, double msf, double mf, double g, double gp, DoubleMatrix & mixNeut, int neut, char type, char LorR);
188  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);
189  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);
190  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);
191  double neutralinoamplitudestauneutrinotauneutrinoNMSSM (double mneut, double mstaunu, double mtaunu, double g, double gp, DoubleMatrix & mixNeut, int neut);
192 
193  double squarkamplitudedecayquarkneutralinoNMSSM (double m1, double mq, double mneut, double g, double gp, DoubleMatrix & mixNeut, char uord, char LorR, int neut);
194  double sleptonamplitudedecayleptonneutralinoNMSSM (double m1, double ml, double mneut, double g, double gp, DoubleMatrix & mixNeut, char uord, char LorR, int neut);
195  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);
196  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);
197  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);
198  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);
199  double charginoiamplitudedecayneutralinojWNMSSM (double mchar, double mneut, double mWboson, double g, double gp, double thetaL, double thetaR, DoubleMatrix & mixNeut, int chargino, int neut);
200  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);
201  double snutauamplitudedecaynutauneutralinoNMSSM (double m1, double mneut, double g, double gp, DoubleMatrix & mixNeut, int neutralino);
202 
203  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);
204 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);
206 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);
207  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);
208  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);
209  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);
210  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);
211  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);
212  double gluinoamplitudedecaygravitino (double m1, double mgrav, double MPlreduced, int gravonoff, int gluNLSP);
213  double squarkamplitudedecaygravitino(double m1, double mgrav, double mquark, double MPlreduced, int gravonoff, int squNLSP);
214  double neutralinoamplitudedecayphotongravitino(double m1, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double g, double gp, int neutralino, int gravonoff, int neutNLSP);
215  double neutralinoamplitudedecayZgravitino(double m1, double mZ, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double g, double gp, double beta, int neutralino, int gravonoff, int neutNLSP);
216  double neutralinoamplitudedecayphigravitino(double m1, double mphi, double mgrav, double MPlreduced, DoubleMatrix & mixNeut, double alpha, double beta, int neutralino, int gravonoff, char phi, int neutNLSP);
217 
219  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);
220  double higgsAamplitudedecaysamecharginoNMSSM (double m1, double m2, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix & CPOMix, int chargino, int pseudoscalar);
221  double higgsAamplitudedecaydifcharginoNMSSM (double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix & CPOMix, int pseudoscalar);
222  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);
223  double higgsAamplitudedecayhiggshorHZbosonNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double thetaA, DoubleMatrix & CPEMix, int pseudoscalar, int higgs);
224  double higgslHamplitudedecayquarkantiquarkNMSSM (double MSSMamplitude, DoubleMatrix & CPEMix, double alpha, int higgs, char uord);
225  double higgshamplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
226  double higgshamplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
227  double higgsHamplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
228  double higgsHamplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
229  double higgsH3amplitudedecay2squarksamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sq);
230  double higgsH3amplitudedecay2sleptonsamehandNMSSM (double m1, double m2, double m3, double g, double gp, double alpha, double beta, double mWboson,double mq, DoubleMatrix & CPEMix, int sl);
231  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);
232  double higgsphiamplitudedecaysamecharginoNMSSM (double m1, double m2, double g, double thetaL, double thetaR, double lam, DoubleMatrix & CPEMix, int chargino, int higgs);
233  double higgsphiamplitudedecaydiffcharginoNMSSM (double m1, double m2, double m3, double g, double thetaL, double thetaR, double lam, DoubleMatrix & CPEMix, int higgs);
234  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);
235  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);
236  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);
237  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);
238  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);
239  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);
240  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);
241  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);
242  double higgsCPevenamplitudedecaypseudoscalarZNMSSM (double m1, double mA, double mZboson, double g, double gp, double beta, DoubleMatrix & CPEMix, DoubleMatrix & CPOMix, int higgs, int pseudoscalar);
243  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);
244  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);
245  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);
246  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);
247  double higgsCPevenamplitudedecayWHpmNMSSM (double m1, double mWboson, double mHpm, double beta, double g, DoubleMatrix & CPEMix, int higgs);
248  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);
249  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);
250  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);
251  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);
252  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);
253  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);
254  double higgsCPevenamplitudedecaysnusnuNMSSM (double m1, double msnu, double g, double gp, double mWboson, double beta, DoubleMatrix & CPEMix, int higgs);
255 
256  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);
258 void slhaDecays(ostream & fout, vector<Particle> & decayTable, bool outputPartialWidths);
259 DoubleVector hggQCDcorrections(double amplitudeW, double alphas, int Nf, char higgs, double prefactor, double SMtotr, double SMtoti, double sqtotr, double sqtoti);
260  DoubleVector higgsmatrixelementgammagammaviatops (double m1, double mtop, double alpha, double beta, char higgstype);
261  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);
262  DoubleVector higgsmatrixelementgammagammaviabottoms (double m1, double mbottom, double alpha, double beta, char higgstype);
263  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);
264  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);
265  DoubleVector higgsmatrixelementgammagammaviaWbosons (double m1, double mWboson, double alpha, double beta, double g, double gprime, char higgstype);
266  DoubleVector higgsmatrixelementgammagammaviaHpms (double m1, double mHpm, double mWboson, double alpha, double beta, double g, double gprime, char higgstype);
267  DoubleVector higgsmatrixelementgammagammaviachargino1s (double m1, double mchargino1, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype);
268  DoubleVector higgsmatrixelementgammagammaviachargino2s (double m1, double mchargino2, double mWboson, double alpha, double beta, double thetaL, double thetaR, char higgstype);
269  DoubleVector higgsmatrixelementgammagammaviacharms (double m1, double mcharm, double alpha, double beta, char higgstype);
270  DoubleVector higgsmatrixelementgammagammaviataus (double m1, double mtau, double alpha, double beta, char higgstype);
271 #endif
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:214
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:74
double m1
decay global variable declarations
Definition: decays.cpp:14
Code calculates decay modes and prints out an SLHA format file with them in. For R-parity conserving ...
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...
bool outputPartialWidths
Default: don't output partial widths in decays.
Definition: def.cpp:61
const double root2
used to flag diabolical problems
Definition: def.h:39
NmssmSoftsusy object contains all the NMSSM functionality: SUSY breaking and preserving parameters,...
Definitions of particles and container for decay widths.
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...
Soft SUSY breaking parameters.
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
MssmSusy contains all SUSY couplings and tan beta, as well as their beta functions.
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:5232
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:2342
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:3680
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:437
double squark2amplitudedecaycharginoW1mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double mup, double mdown, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:225
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:2831
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:1554
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:2065
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:671
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:4308
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:4627
double squark3amplitudedecayneutralino(double m1, double m2, double m3, double mWboson, double theta, double beta, DoubleMatrix &mixNeut, double g, double gp, double runmt, int squark, int oneortwo, int neutralino)
Definition: twoBodyDecays.cpp:299
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:3326
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:4863
double sleptonamplitudedecayleptonneutralinoNMSSM(double m1, double ml, double mneut, double g, double gp, DoubleMatrix &mixNeut, char uord, char LorR, int neut)
Definition: threeBodyDecays.cpp:5723
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:2266
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:4001
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:640
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:4708
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:190
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:3053
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:1199
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:3785
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:1660
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:903
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:4042
DoubleVector higgshamplitudedecayVV(double m1, double mWboson, double mZboson, double g, double gp, double alpha, double beta, char Vtype, DoubleMatrix &CPEMix, bool nmssmIsIt)
double higgsamplitudedecayVVstar (double m1, double mboson, double g, double gp, double beta,...
Definition: twoBodyDecays.cpp:4796
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:2316
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:2521
double charginoamplitudedecayquarksquarkL(double m1, double m2, double m3, double g, double theta, int chargino)
quark mass is m2
Definition: twoBodyDecays.cpp:942
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:1368
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:3365
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:260
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:2634
double higgsphiamplitudedecaydifchargino(double m1, double m2, double m3, double g, double thetaL, double thetaR, double alpha, double beta, char phi)
this function calls the function higgsphicharginocouplings to calculate the couplings for it
Definition: twoBodyDecays.cpp:2121
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:2749
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)
NMSSM functions (where MSSM functions haven't been recycled)
Definition: twoBodyDecays.cpp:3751
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:1438
double stauamplitudedecaysnustauHminus(double m1, double m2, double m3, double g, double mWboson, double beta, double thetatau, double mtau, double mu, double Atau, int oneortwo)
Also does decay mode snustau to Hplus and stau1/2 just with m1, m2 and m3 permuted.
Definition: twoBodyDecays.cpp:811
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:3640
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:495
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:4772
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:3482
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:2191
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:4930
double charginoToNeutralino1pion(const MssmSoftsusy *m)
Two body partial width calculations.
Definition: twoBodyDecays.cpp:14
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)
BEN.
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:3561
double squark2amplitudedecaycharginoW2mix(double m1, double m2, double m3, double g, double gammaL, double gammaR, double theta, double beta, double mWboson, double mup, double mdown, double torb)
the variable torb depends on if it is stop (torb=1) or sbottom (torb =2) decaying and changes Aprimeu...
Definition: twoBodyDecays.cpp:241
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:4531
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:3812
double gluinoamplitudedecay(double m1, double m2, double m3, double alphastrong)
First, do hadronic decays.
Definition: twoBodyDecays.cpp:87
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:2597
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:94
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:4635
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:1270
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:4518
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:5265
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:976
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:5248
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:2369
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:3720
double stauamplitudedecaysnustauWboson(double m1, double m2, double m3, double g, double thetatau, int oneortwo)
m3 must be mw here
Definition: twoBodyDecays.cpp:847
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:2216
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:4306
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:1409
double chargino2amplitudedecaychargino1Zboson(double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR)
m2 must be Zboson mass
Definition: twoBodyDecays.cpp:1241
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:4348
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:727
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:3248
double stau2amplitudedecaystau1Zboson(double m1, double m2, double m3, double g, double gp, double thetatau)
m3 must be mz here
Definition: twoBodyDecays.cpp:879
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:1510
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:2907
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:2394
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:4732
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:1031
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:3953
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:207
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:3853
double stausneutrinoamplitudedecaytauneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: twoBodyDecays.cpp:706
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:1633
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:2559
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:1065
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:2790
double gluinoamplitudedecaygravitino(double m1, double mgrav, double MPlreduced, int gravonoff, int gluNLSP)
Decays to gravitinos.
Definition: twoBodyDecays.cpp:5206
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:3600
double higgsphiamplitudedecaysamecharginoNMSSM(double m1, double m2, double g, double thetaL, double thetaR, double lam, DoubleMatrix &CPEMix, int chargino, int higgs)
Definition: twoBodyDecays.cpp:2031
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:2241
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:2979
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:4439
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:4476
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)
uord indicates if it's an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1704
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:2940
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:3911
double stausneutrinoamplitudedecaytauchargino(double m1, double m2, double m3, double g, double mWboson, double beta, double thetaL, double thetaR, int chargino)
Definition: twoBodyDecays.cpp:776
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:3131
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:4997
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:1162
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,...
Definition: twoBodyDecays.cpp:2439
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:3287
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:1976
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:2481
double sneutrinoamplitudedecayneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: twoBodyDecays.cpp:617
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:1888
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'...
Definition: twoBodyDecays.cpp:363
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:1102
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:1317
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:2868
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:4176
double higgsAamplitudedecayquarkantiquarkNMSSM(double m1, double m2, double beta, DoubleMatrix &CPOMix, int uord, int higgs, bool QCD, double alphas)
uord indicates if it's an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1850
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:3404
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:3013
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:4448
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:3879
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:4348
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, 'u' for up type quark, 'd' for down type quark, 'l' for charged lepto...
Definition: threeBodyDecays.cpp:5448
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:4267
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:4135
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:3092
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:4401
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,...
Definition: twoBodyDecays.cpp:2706
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:2672
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:3443
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:4570
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:2291
void slhaDecays(ostream &fout, vector< Particle > &decayTable, bool outputPartialWidths)
outputs a decay table in SLHA format
Definition: twoBodyDecays.cpp:5288
double snutauamplitudedecaynutauneutralinoNMSSM(double m1, double mneut, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: threeBodyDecays.cpp:6020
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:3209
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:3522
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:1592
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:4706
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,...
Definition: twoBodyDecays.cpp:5219
double higgsAamplitudedecayquarkantiquark(double m1, double m2, double g, double beta, double mWboson, int uord, bool QCD, double alphas)
uord indicates if it's an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: twoBodyDecays.cpp:1804
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:3170
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:277
int theta(double a)
Standard theta function: 1 is a>0, 0 otherwise.
Definition: utils.cpp:25