SoftSUSY is hosted by Hepforge, IPPP Durham
SOFTSUSY  4.0
decays.h
Go to the documentation of this file.
1 
11 #ifndef DECAYS_H
12 #define DECAYS_H
13 
14 #include "nmssmsoftsusy.h"
15 #include "softsusy.h"
16 #include "physpars.h"
17 #include "lowe.h"
18 #include "softpars.h"
19 #include "softsusy.h"
20 #include "flavoursoft.h"
21 #include "susy.h"
22 #include "particle.h"
23 #include <iostream>
24 #include <cstring>
25 #include <stdlib.h>
26 #include <cmath>
27 #include <stdlib.h>
28 #include <vector>
29 
30 using namespace std;
31 
33 const double accuracy = 0.01;
34 
36 int calculateDecays(ostream & out, MssmSoftsusy * r,
37  const NmssmSoftsusy & nmssm, bool nmssmIsIt);
39 double fdgauss(double x);
40 
43 const double mcpole = 1.40;
46 const double mspole = 0.19;
47 
49 double gluinoamplitudedecay (double m1, double m2, double m3,
50  double alphastrong);
51 double gluinoamplitudedecaymix (double m1, double m2, double m3,
52  double alphastrong, double squarkmix,
53  double theta);
54 double squarkamplitudedecaygluino (double m1, double m2, double m3,
55  double alphastrong);
56 double squarkamplitudedecaygluinomix (double m1, double m2, double m3,
57  double alphastrong, double squarkmix,
58  double theta);
59 double squarkamplitudedecaycharginoW1 (double m1, double m2, double m3,
60  double g, double gamma);
61 double squarkamplitudedecaycharginoW2 (double m1, double m2, double m3,
62  double g, double gamma);
64 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
65  double theta, double beta, double mWboson, double runmt, double runmb,
66  double torb);
68 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
69  double theta, double beta, double mWboson, double runmt, double runmb,
70  double torb);
72 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
73  double theta, double beta, double mWboson, double mup, double mdown,
74  double torb);
76 (double m1, double m2, double m3, double g, double gammaL, double gammaR,
77  double theta, double beta, double mWboson, double mup, double mdown,
78  double torb);
80 (double m1, double m2, double m3, double g, double gprime,
81  DoubleMatrix & mixNeut, int neutralino, int uord );
82  double squarkRamplitudedecayneutralino (double m1, double m2, double m3, double g, double gprime, DoubleMatrix & mixNeut, int neutralino, int uord );
83 
84  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);
85  double squark3amplitudedecaysquark3Wboson (double m1, double m2, double m3, double g, double thetat, double thetab, int m1torb, int m1oneortwo, int m3torb, int m3oneortwo);
86  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);
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);
257 
258  DoubleVector squarkmixcharginocouplings (double g, double theta, double beta, double gammaL, double gammaR, double runmt, double runmb, double mWboson, int torb);
259  DoubleVector higgsphisamecharginocouplings(double alpha, double beta, double thetaL, double thetaR);
260  DoubleVector higgsphidifcharginocouplings (double alpha, double beta, double thetaL, double thetaR);
261  DoubleVector higgshsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
262  DoubleVector higgshsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
263  DoubleVector higgsHsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
264  DoubleVector higgsHsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
265  DoubleVector higgshsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel);
266  DoubleVector higgshsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
267  DoubleVector higgsHsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel) ;
268  DoubleVector higgsHsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
269  DoubleVector higgsHplussquarkcouplings (double mWboson, double g, double beta, double mup, double mdown, double greekmu, double Au, double Ad);
270 
271 
272  DoubleVector foftau(double mpart, double mcomp);
273 
274 
275  void printRowPDG(ostream & cout, double x);
276  void OutputNoPWs(ostream & cout, Particle & P, double BRTol);
277  void OutputYesPWs(ostream & cout, Particle & P, double BRTol);
278 
280  double Zsfintegralsum(double m1, double m2, double msf, double mf, double mz, double min, double max, double Nsteps, int adaptive, double approx);
281  double Aintegralsum(double m1, double m2, double mz, double mA, double mf, double min, double max, double Nsteps, int adaptive, double approx);
282  double G3integralsum(double m1, double m2, double m3, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
283  double G2integralsum(double m1, double m2, double m3, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
284  double chiprimeintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
285  double Yintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
286  double Xintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
287  double zetaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
288  double chiintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
289  double rhointegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
290  double xsiintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
291  double phitildaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
292  double psitildaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
293 
294  double Zintegralsum(double m1, double m2, double mz, double mf, double min, double max, double Nsteps, int adaptive, double approx);
295  double G1integralsum(double m1, double m2, double m3, double mq, double min, double max, double Nsteps, int adaptive, double approx);
296 
297  double G4integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
298  double G5integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
299  double G6integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
300  double G7integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
301  double G8integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
302 
303  double Jintegralsum(double m1, double m2, double msf, double mphi, double mf, double min, double max, double Nsteps, double approx, int adaptive, int AorhorH);
304 double ghHintegral(double m1, double m2, double mf, double mh, double mH, double g, double gp, double alpha, int neutralinoi, int neutralinoj, DoubleMatrix & mixNeut, double E);
305 double hHintegral (double m1, double m2, double mf, double mh, double mH, double fromE, double toE, double stepE, double g, double gp, double alpha, DoubleMatrix & mixNeut, int neutralinoi, int neutralinoj);
308 
309 
310 
312  double gluinoamplitudedecaydgaussneutralinoqqpbarfirsttwogen (double mgluino, double mneutralino, double msqL, double msqR, double mquark, double g, double gp, DoubleMatrix & mixNeut, double alphas, char uord, int neut, bool onetothree);
313  double gluinoamplitudedecaydgaussneutralinottbar (double mgluino, double mst1, double mst2, double mneutralino, double mt, double mWboson, double g, double gp, double thetat, double beta, double alphas, DoubleMatrix & mixNeut, double runmt, int neutralino, bool onetothree, char torb);
314  // double gluinoamplitudedecaydgaussneutralinobbbarorig(double mgluino, double msb1, double msb2, double mneutralino, double mb, double mWboson, double g, double gp, double thetab, double beta, double alphas, DoubleMatrix & mixNeut, double runmb, int neutralino, bool onetothree);
315  double gluinoamplitudedecaydgausschartbbar (double mgluino, double mst1, double mst2, double msb1, double msb2, double mtop, double mbottom, double mchar, double alphas, double thetat, double thetab, double MWBoson, double g, double gp, double gammaL, double gammaR, double beta, double runmt, double runmb, int chargino, bool onetothree);
316  // double gluinoamplitudedecaydgausschartbbarorig (double mgluino, double mst1, double mst2, double msb1, double msb2, double mtop, double mbottom, double mchar, double alphas, double thetat, double thetab, double MWBoson, double g, double gp, double gammaL, double gammaR, double beta, double runmt, double runmb, int chargino, bool onetothree);
317  double neutralinoamplitudedecaydgaussneutralinoffbar (double mneutralinoi, double msf1, double msf2, double mZboson, double mhiggsl, double mhiggsH, double mhiggsA, double mneutralinoj, double mf, double alphas, double thetaq, double mWboson, double g, double gp, double alpha, double beta, double runmq, DoubleMatrix & mixNeut, int ineutralino, int jneutralino, bool onetothree, char uordornuorl);
318  double neutralinoamplitudedecaycharginoffprimebar (double mneutralinoi, double msfp1, double msfp2, double msf1, double msf2, double mWboson, double mHP, double mcharginoj, double mfp, double mf, double thetaq, double thetaqp, double g, double gp, double alpha, double beta, double thetaL2, double thetaR2, double runmqp, double runmq, DoubleMatrix & mixNeut, int ineutralino, int jchargino, bool onetothree, char qorl, char norc);
319 
320 
323 double fdgauss(double x);
324 double gpsitildadgauss(double Et);
325 double gphitildadgauss(double Et);
326 double gxsidgauss (double Et);
327 double grhodgauss (double Et);
328 double gchidgauss (double Et);
329 double gzetadgauss (double Et);
330 double gXdgauss (double Et);
331 double gYdgauss (double Et);
332 double gchiprimedgauss (double Et);
333 double gG1dgauss (double Et);
334 double gG4dgauss (double Et);
335 double gG5dgauss(double Et);
336 double gG6dgauss(double Et);
337 double gG7dgauss(double Et);
338 double gG8dgauss(double Et);
339 double gG2dgauss(double Eb);
340 double gG3dgauss(double Eb);
341 double gZdgauss(double E);
342 double ghHdgauss (double E);
343 double gAdgauss (double E);
344 double gZsfdgauss (double s);
345 double gJdgauss(double s);
346 
347 double gneutineutjffZ1dgauss(double s);
348 double gneutineutjffZ2dgauss(double s);
349 double gneutineutjffZ3dgauss(double s);
350 double gneutineutjffZ4dgauss(double s);
351 
352 double gintegralhdgauss(double E);
353 double gintegralHdgauss(double E);
354 double gintegralh1dgauss(double E);
355 double gintegralh2dgauss(double E);
356 double gintegralh3dgauss(double E);
357 double gintegralh4dgauss(double E);
358 double gintegralH1dgauss(double E);
359 double gintegralH2dgauss(double E);
360  double gintegralH3dgauss(double E);
361  double gintegralH4dgauss(double E);
362  double gintegralhH1dgauss(double E);
363  double gintegralhH2dgauss(double E);
364  double gintegralhH3dgauss(double E);
365  double gintegralhH4dgauss(double E);
366  double gintegralA1dgauss(double E);
367  double gintegralA2dgauss(double E);
368  double gintegralA3dgauss(double E);
369  double gintegralA4dgauss(double E);
370  double gintegral1Zsfdgauss(double E);
371  double gintegral2Zsfdgauss(double E);
372  double gintegral3Zsfdgauss(double E);
373  double gintegral4Zsfdgauss(double E);
374  double gintegral5Zsfdgauss(double E);
375  double gintegral6Zsfdgauss(double E);
376  double gintegral7Zsfdgauss(double E);
377  double gintegral8Zsfdgauss(double E);
378  double gintegral1hsfdgauss(double E);
379  double gintegral2hsfdgauss(double E);
380  double gintegral3hsfdgauss(double E);
381  double gintegral4hsfdgauss(double E);
382  double gintegral5hsfdgauss(double E);
383  double gintegral6hsfdgauss(double E);
384  double gintegral7hsfdgauss(double E);
385  double gintegral8hsfdgauss(double E);
386  double gintegral1Hsfdgauss(double E);
387  double gintegral2Hsfdgauss(double E);
388  double gintegral3Hsfdgauss(double E);
389  double gintegral4Hsfdgauss(double E);
390  double gintegral5Hsfdgauss(double E);
391  double gintegral6Hsfdgauss(double E);
392  double gintegral7Hsfdgauss(double E);
393  double gintegral8Hsfdgauss(double E);
394  double gintegral1ZAdgauss(double E);
395  double gintegral2ZAdgauss(double E);
396  double gintegral3ZAdgauss(double E);
397  double gintegral4ZAdgauss(double E);
398  double gneutineutjffgA1dgauss(double E);
399  double gneutineutjffgA2dgauss(double E);
400  double gneutineutjffgA3dgauss(double E);
401  double gneutineutjffgA4dgauss(double E);
402 
403  double gneuticharjffpW1dgauss(double E);
404  double gneuticharjffpW2dgauss(double E);
405  double gneuticharjffpHpm1dgauss(double E);
406  double gneuticharjffpHpm2dgauss(double E);
407  double gneuticharjffpHpm3dgauss(double E);
408  double gneuticharjffpHpm4dgauss(double E);
409  double gneuticharjffp1sf1sf2dgauss(double E);
410  double gneuticharjffp2sf1sf2dgauss(double E);
411  double gneuticharjffp3sf1sf2dgauss(double E);
412  double gneuticharjffp4sf1sf2dgauss(double E);
413  double gneuticharjffp1sfp1sf2dgauss(double E);
414  double gneuticharjffp2sfp1sf2dgauss(double E);
415  double gneuticharjffp3sfp1sf2dgauss(double E);
416  double gneuticharjffp4sfp1sf2dgauss(double E);
417  double gneuticharjffp5sfp1sf2dgauss(double E);
418  double gneuticharjffp6sfp1sf2dgauss(double E);
419  double gneuticharjffp7sfp1sf2dgauss(double E);
420  double gneuticharjffp8sfp1sf2dgauss(double E);
421  double gneuticharjffp1WHpmdgauss(double E);
422  double gneuticharjffp2WHpmdgauss(double E);
423  double gneuticharjffp3WHpmdgauss(double E);
424  double gneuticharjffp4WHpmdgauss(double E);
425  double gneuticharjffpW1Sfpdgauss(double E);
426  double gneuticharjffpW2Sfpdgauss(double E);
427  double gneuticharjffpW3Sfpdgauss(double E);
428  double gneuticharjffpW4Sfpdgauss(double E);
429  double gneuticharjffpW5Sfpdgauss(double E);
430  double gneuticharjffpW6Sfpdgauss(double E);
431  double gneuticharjffpW7Sfpdgauss(double E);
432  double gneuticharjffpW8Sfpdgauss(double E);
433  double gneuticharjffpHg1dgauss(double E);
434  double gneuticharjffpHg2dgauss(double E);
435  double gneuticharjffpHg3dgauss(double E);
436  double gneuticharjffpHg4dgauss(double E);
437  double gneuticharjffp1gsfpdgauss(double E);
438  double gneuticharjffp2gsfpdgauss(double E);
439  double gneuticharjffp3gsfpdgauss(double E);
440  double gneuticharjffp4gsfpdgauss(double E);
441  double gneuticharjffp5gsfpdgauss(double E);
442  double gneuticharjffp6gsfpdgauss(double E);
443  double gneuticharjffp7gsfpdgauss(double E);
444  double gneuticharjffp8gsfpdgauss(double E);
445  double gneuticharjffp1sfpsfpdgauss(double E);
446  double gneuticharjffp21sfpsfpdgauss(double E);
447  double gneuticharjffp31sfpsfpdgauss(double E);
448  double gneuticharjffp41sfpsfpdgauss(double E);
449 
450 
451 #endif
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: decays.cpp:11436
void OutputYesPWs(ostream &cout, Particle &P, double BRTol)
Outputs the decay table into the leshouchesOutput file with PWs (partial widths) given after the comm...
Definition: decays.cpp:12647
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: decays.cpp:16817
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: decays.cpp:10386
double neutralinoamplitudedecaydgaussneutralinoffbar(double mneutralinoi, double msf1, double msf2, double mZboson, double mhiggsl, double mhiggsH, double mhiggsA, double mneutralinoj, double mf, double alphas, double thetaq, double mWboson, double g, double gp, double alpha, double beta, double runmq, DoubleMatrix &mixNeut, int ineutralino, int jneutralino, bool onetothree, char uordornuorl)
Definition: decays.cpp:15035
double higgsHamplitudedecayAA(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for H->AA
Definition: decays.cpp:8869
double gG4dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: decays.cpp:12839
double gneuticharjffp8gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14604
double gintegralhH1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: decays.cpp:13175
double gneuticharjffp7sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14120
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: decays.cpp:9493
double gintegral2Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13262
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: decays.cpp:16672
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: decays.cpp:9911
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
double gintegral4hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13443
const double accuracy
Approximate accuracy with which 3 body decays are calculated.
Definition: decays.h:33
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35
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: decays.cpp:7539
DoubleVector squarkmixcharginocouplings(double g, double theta, double beta, double gammaL, double gammaR, double runmt, double runmb, double mWboson, double mch1, double mch2, int torb)
Functions that calculate the couplings:
Definition: decays.cpp:12130
double gneuticharjffpHg3dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: decays.cpp:14428
double gintegralH3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: decays.cpp:13157
double gneuticharjffpHpm3dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm;
Definition: decays.cpp:13890
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: decays.cpp:7365
double gintegralh3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13123
double gneuticharjffp2gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14478
double gintegral1hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13384
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: decays.cpp:10542
double gluinoamplitudedecay(double m1, double m2, double m3, double alphastrong)
Function Declarations used to calculate Partial Widths in decays.cpp.
Definition: decays.cpp:6514
double gneuticharjffpW5Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14316
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: decays.cpp:17107
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: decays.cpp:9163
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: decays.cpp:7901
Contains all supersymmetric NMSSM parameters, incorporating R_p NMSSM.
Definition: nmssmsoftsusy.h:45
double gneuticharjffp2sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14014
double gintegralhH4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: decays.cpp:13201
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: decays.cpp:8261
const double mcpole
Definition: decays.h:43
double gintegral7Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13348
double gintegralH1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: decays.cpp:13140
double gintegral1Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13538
double gneuticharjffpW2Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14252
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: decays.cpp:9825
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: decays.cpp:12058
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: decays.cpp:9028
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: decays.cpp:9243
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: decays.cpp:6877
double gneuticharjffp3gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14499
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: decays.cpp:10212
DoubleVector higgshsleptondiffhandcouplings(double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae)
calculates the couplings of light scalar higgs h to two sleptons of different handedness ...
Definition: decays.cpp:12342
double gneutineutjffZ4dgauss(double s)
Definition: decays.cpp:13067
double gintegralH4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13166
double neutralinoamplitudedecayleptonsleptonLorR(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, char LorR, int neutralino)
m2 must be lepton mass
Definition: decays.cpp:7955
double sleptonamplitudedecayleptonneutralinoNMSSM(double m1, double ml, double mneut, double g, double gp, DoubleMatrix &mixNeut, char uord, char LorR, int neut)
Definition: decays.cpp:18142
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: decays.cpp:8190
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: decays.cpp:7462
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: decays.cpp:6679
DoubleVector higgsHplusamplitudedecaysquarksquark(double m1, double m2, double m3, double g, double beta, double mWboson, double mup, double mdown, double greekmu, double Aup, double Adown)
calculates partial width for Hplus->squarki antisquarkj where i,j are each L/R so no mixing...
Definition: decays.cpp:10636
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: decays.cpp:10256
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: decays.cpp:7184
DoubleVector foftau(double mpart, double mcomp)
f(tau) function for use in h->gamma gamma or Z gamma
Definition: decays.cpp:12404
double gintegralA1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: decays.cpp:13209
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...
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: decays.cpp:7600
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: decays.cpp:17082
double sneutrinoamplitudedecayneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: decays.cpp:7160
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: decays.cpp:8411
double higgshamplitudedecayAA(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for h->AA
Definition: decays.cpp:8817
DoubleVector higgsphidifcharginocouplings(double alpha, double beta, double thetaL, double thetaR)
calculates the couplings of a netural higgs (h,H,A) to Wtildai Wtildaj - i.e. to two charginos of dif...
Definition: decays.cpp:12249
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: decays.cpp:10677
double higgsphiamplitudedecaysamecharginoNMSSM(double m1, double m2, double g, double thetaL, double thetaR, double lam, DoubleMatrix &CPEMix, int chargino, int higgs)
Definition: decays.cpp:8647
double gintegral2Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13557
double gneutineutjffZ1dgauss(double s)
m1 = mneuti, m4 = mneutj, mq = mf, MZboson = mZboson
Definition: decays.cpp:13013
DoubleVector higgsHsleptonsamehandcouplings(double mWboson, double g, double gp, double alpha, double beta, double mel)
calculates the couplings of heavy scalar higgs H to two sleptons of same handedness ...
Definition: decays.cpp:12355
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: decays.cpp:6805
double gluinoamplitudedecaydgaussneutralinottbar(double mgluino, double mst1, double mst2, double mneutralino, double mt, double mWboson, double g, double gp, double thetat, double beta, double alphas, DoubleMatrix &mixNeut, double runmt, int neutralino, bool onetothree, char torb)
calculates PW for gluino -> neutralino + q qbar pair where q are t
Definition: decays.cpp:14738
double gneuticharjffp6gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14562
double gneuticharjffpW8Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14379
double gpsitildadgauss(double Et)
Definition: decays.cpp:12677
double gneuticharjffp5gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14541
double gintegral4Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13297
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: decays.cpp:12541
double gneuticharjffp5sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14078
double fdgauss(double x)
Dgauss numerical integration integrand.
Definition: decays.cpp:12673
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: decays.cpp:11614
double gneuticharjffp6sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14099
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: decays.cpp:7852
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&#39;s an up type quark (1) or down type quark (0) to choose which trig functions to ...
Definition: decays.cpp:8307
DoubleVector higgshsleptonsamehandcouplings(double mWboson, double g, double gp, double alpha, double beta, double mel)
calculates the couplings of light scalar higgs h to two sleptons of same handedness ...
Definition: decays.cpp:12327
double ghHdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt, MZboson = mz
Definition: decays.cpp:12952
DoubleVector higgsHplussquarkcouplings(double mWboson, double g, double beta, double mup, double mdown, double greekmu, double Au, double Ad)
calculates the couplings of charged scalar higgs H+ to two squarks
Definition: decays.cpp:12383
double gneuticharjffpHpm1dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm;
Definition: decays.cpp:13860
double gneuticharjffpW2dgauss(double E)
m1 = mZi, m2 = mWj, m3 = mf, m4 = mfp, MWboson = mw
Definition: decays.cpp:13844
double gneuticharjffp7gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14583
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: decays.cpp:8458
MssmSusy contains all SUSY couplings and tan beta, as well as their beta functions.
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: decays.cpp:8921
double gneuticharjffp3sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: decays.cpp:13956
DoubleVector higgsphisamecharginocouplings(double alpha, double beta, double thetaL, double thetaR)
calculates the couplings of a neutral higgs (h,H,A) to Wtildai Wtildai - i.e. to two charginos of the...
Definition: decays.cpp:12227
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: decays.cpp:12104
double gintegral5Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13314
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: decays.cpp:12071
double chargino2amplitudedecaychargino1Zboson(double m1, double m2, double m3, double g, double gp, double thetaL, double thetaR)
m2 must be Zboson mass
Definition: decays.cpp:7822
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: decays.cpp:8497
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: decays.cpp:10570
double gintegralh2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13113
double stauamplitudedecaysnustauWboson(double m1, double m2, double m3, double g, double thetatau, int oneortwo)
m3 must be mw here
Definition: decays.cpp:7403
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: decays.cpp:9739
double gintegral3hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13423
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: decays.cpp:9954
double gAdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt
Definition: decays.cpp:12971
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: decays.cpp:9868
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: decays.cpp:9997
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: decays.cpp:9782
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: decays.cpp:6705
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: decays.cpp:7998
double higgsAamplitudedecaysamecharginoNMSSM(double m1, double m2, double g, double thetaL, double thetaR, double alpha, double lam, DoubleMatrix &CPOMix, int chargino, int pseudoscalar)
Definition: decays.cpp:8683
double gintegral6Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13635
double gneutineutjffgA1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: decays.cpp:13759
double gneuticharjffp8sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14142
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: decays.cpp:10612
double Zsfintegralsum(double m1, double m2, double msf, double mf, double mz, double min, double max, double Nsteps, int adaptive, double approx)
Integral calculating functions for 1->3 decays.
double gXdgauss(double Et)
Definition: decays.cpp:12777
int calculateDecays(ostream &out, MssmSoftsusy *r, const NmssmSoftsusy &nmssm, bool nmssmIsIt)
Calculate Decays does all decay table calculations and outputs.
Definition: decays.cpp:29
double gneuticharjffpW4Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14295
double gneutineutjffgA3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: decays.cpp:13793
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: decays.cpp:8150
DoubleVector higgshsquarksamehandcouplings(double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq)
calculates the couplings of light scalar higgs h to two squarks of same handedness ...
Definition: decays.cpp:12268
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
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: decays.cpp:8947
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: decays.cpp:10040
double gG6dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: decays.cpp:12865
double gintegral5hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13462
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: decays.cpp:9002
double gintegral6hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13481
double stausneutrinoamplitudedecaytauneutrinoneutralino(double m1, double m2, double m3, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: decays.cpp:7254
double gneuticharjffp4WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: decays.cpp:14214
double gintegralH2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: decays.cpp:13148
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:211
double gneuticharjffp4gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14520
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: decays.cpp:16715
double gintegralA4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: decays.cpp:13236
double gintegral8hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13519
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: decays.cpp:8974
double gneuticharjffpW1Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14231
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: decays.cpp:9695
double gneuticharjffp1WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: decays.cpp:14164
double gintegral3Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13577
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: decays.cpp:9075
double gneuticharjffpW3Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14273
Particle class definition.
Definition: particle.h:14
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: decays.cpp:9203
void printRowPDG(ostream &cout, double x)
Definition: decays.cpp:12609
double neutralinoamplitudedecaycharginoffprimebar(double mneutralinoi, double msfp1, double msfp2, double msf1, double msf2, double mWboson, double mHP, double mcharginoj, double mfp, double mf, double thetaq, double thetaqp, double g, double gp, double alpha, double beta, double thetaL2, double thetaR2, double runmqp, double runmq, DoubleMatrix &mixNeut, int ineutralino, int jchargino, bool onetothree, char qorl, char norc)
Definition: decays.cpp:15583
DoubleVector higgsHsquarksamehandcouplings(double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq)
calculates the couplings of heavy scalar higgs H to two squarks of same handedness ...
Definition: decays.cpp:12297
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: decays.cpp:17149
double gneuticharjffpW1dgauss(double E)
m1 = mZi, m2 = mWj, m3 = mf, m4 = mfp, MWboson = mw
Definition: decays.cpp:13829
NmssmSoftsusy object contains all the NMSSM functionality: SUSY breaking and preserving parameters...
double gZsfdgauss(double s)
m2 = msfi where i is the sfermion contribution interference with Z you are considering ...
Definition: decays.cpp:12981
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: decays.cpp:6653
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...
double gneuticharjffpHg2dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: decays.cpp:14414
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: decays.cpp:10299
double gintegralh4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13132
double gneuticharjffp4sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: decays.cpp:13973
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: decays.cpp:9577
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: decays.cpp:9407
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: decays.cpp:6775
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: decays.cpp:7636
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: decays.cpp:6964
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: decays.cpp:17006
const double mspole
Definition: decays.h:46
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: decays.cpp:8588
double gintegralA3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: decays.cpp:13227
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 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: decays.cpp:16769
double snutauamplitudedecaynutauneutralinoNMSSM(double m1, double mneut, double g, double gp, DoubleMatrix &mixNeut, int neutralino)
Definition: decays.cpp:18452
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: decays.cpp:7739
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: decays.cpp:8233
double gintegralA2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: decays.cpp:13218
double gneuticharjffp1sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:13991
double charginoamplitudedecayquarksquarkL(double m1, double m2, double m3, double g, double theta, int chargino)
quark mass is m2
Definition: decays.cpp:7503
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: decays.cpp:7778
Soft SUSY breaking parameters.
double gchidgauss(double Et)
Definition: decays.cpp:12741
double gneuticharjffp2sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: decays.cpp:13939
double higgsHamplitudedecayhh(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for h->hh
Definition: decays.cpp:8843
double gneuticharjffpW7Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14358
double gneutineutjffZ3dgauss(double s)
Definition: decays.cpp:13049
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: decays.cpp:9360
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: decays.cpp:11684
double gintegralHdgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: decays.cpp:13095
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: decays.cpp:16903
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:83
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: decays.cpp:11754
double gintegral1ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: decays.cpp:13692
double gintegral1Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13244
double gintegral5Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13616
double gintegral8Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13673
double gintegralhH2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: decays.cpp:13183
double gneuticharjffp1gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: decays.cpp:14457
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: decays.cpp:10721
double gneutineutjffZ2dgauss(double s)
Definition: decays.cpp:13031
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: decays.cpp:7675
double gintegral7hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13500
double gintegral2ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: decays.cpp:13708
double gintegral6Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13331
double gluinoamplitudedecaydgaussneutralinoqqpbarfirsttwogen(double mgluino, double mneutralino, double msqL, double msqR, double mquark, double g, double gp, DoubleMatrix &mixNeut, double alphas, char uord, int neut, bool onetothree)
only done via usual method (not via a separate integrating function compute_areai) as can&#39;t pass a Do...
Definition: decays.cpp:14683
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: decays.cpp:10083
double gneuticharjffp4sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14057
double stau2amplitudedecaystau1Zboson(double m1, double m2, double m3, double g, double gp, double thetatau)
m3 must be mz here
Definition: decays.cpp:7437
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: decays.cpp:10169
double gintegralhdgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13086
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: decays.cpp:8028
Definitions of particles and container for decay widths.
double gintegral4Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13597
double gneuticharjffp3sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: decays.cpp:14036
double gintegral8Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13365
double gneuticharjffpHg4dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: decays.cpp:14442
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&#39;t been recycled)
Definition: decays.cpp:10507
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: decays.cpp:9121
double gneutineutjffgA4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: decays.cpp:13810
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: decays.cpp:6730
double gintegral4ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: decays.cpp:13742
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: decays.cpp:9282
double gintegral3Zsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, MZboson = mz, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13279
double gG1dgauss(double Et)
m1 = mgluino, m2 = mstopi, m6 = mtop, m8 = mcharginoj
Definition: decays.cpp:12829
double stausneutrinoamplitudedecaytauchargino(double m1, double m2, double m3, double g, double mWboson, double beta, double thetaL, double thetaR, int chargino)
Definition: decays.cpp:7328
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: decays.cpp:9614
double gluinoamplitudedecaygravitino(double m1, double mgrav, double MPlreduced, int gravonoff, int gluNLSP)
Decays to gravitinos.
Definition: decays.cpp:12045
double gintegral2hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13403
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: decays.cpp:7217
double gneutineutjffgA2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: decays.cpp:13776
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: decays.cpp:17855
double gneuticharjffp1sfpsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = -mf, m4 = charginoj, m5 = msfp1, m6 = msfp2
Definition: decays.cpp:14625
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: decays.cpp:10430
double gneuticharjffpW6Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: decays.cpp:14337
double gintegral7Hsfdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mhiggsl = mh, m2 = msfi where i is the sfermion index you&#39;re considering...
Definition: decays.cpp:13654
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: decays.cpp:10126
DoubleVector higgshsquarkdiffhandcouplings(double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown)
calculates the couplings of light scalar higgs h to two squarks of different handedness ...
Definition: decays.cpp:12283
double gneuticharjffpHpm4dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm;
Definition: decays.cpp:13905
double gG5dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: decays.cpp:12853
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: decays.cpp:7276
double gintegral3ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: decays.cpp:13725
DoubleVector higgsHsleptondiffhandcouplings(double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae)
calculates the couplings of heavy scalar higgs H to two sleptons of different handedness ...
Definition: decays.cpp:12370
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: decays.cpp:10474
double gneuticharjffp1sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: decays.cpp:13922
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: decays.cpp:9657
double gZdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt, MZboson = mz
Definition: decays.cpp:12938
double gneuticharjffp3WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: decays.cpp:14197
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: decays.cpp:8742
double gintegralhH3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: decays.cpp:13192
DoubleVector higgsHsquarkdiffhandcouplings(double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown)
calculates the couplings of heavy scalar higgs H to two squarks of different handedness ...
Definition: decays.cpp:12313
double gintegralh1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: decays.cpp:13105
int theta(double a)
Standard theta function: 1 is a>0, 0 otherwise.
Definition: utils.cpp:25
double gluinoamplitudedecaydgausschartbbar(double mgluino, double mst1, double mst2, double msb1, double msb2, double mtop, double mbottom, double mchar, double alphas, double thetat, double thetab, double MWBoson, double g, double gp, double gammaL, double gammaR, double beta, double runmt, double runmb, int chargino, bool onetothree)
Definition: decays.cpp:14906
double gneuticharjffpHpm2dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm;
Definition: decays.cpp:13875
double gneuticharjffp2WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: decays.cpp:14180
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: decays.cpp:9450
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: decays.cpp:12087
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: decays.cpp:8104
double higgsHamplitudedecayHplusHminus(double m1, double m2, double g, double gp, double alpha, double beta, double mWboson)
calculates partial width for H->H+H-
Definition: decays.cpp:8895
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: decays.cpp:11824
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: decays.cpp:9324
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: decays.cpp:7032
double gneuticharjffpHg1dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: decays.cpp:14400
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: decays.cpp:10342
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: decays.cpp:9534
void OutputNoPWs(ostream &cout, Particle &P, double BRTol)
Outputs the decay table into the leshouchesOutput file with no PWs (partial widths) given...
Definition: decays.cpp:12623
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: decays.cpp:6754