SoftSUSY is hosted by Hepforge, IPPP Durham
SOFTSUSY  4.1
threeBodyDecays.h
Go to the documentation of this file.
1 
12 #ifndef THREEBODYDECAYS_H
13 #define THREEBODYDECAYS_H
14 
15 #include "nmssmsoftsusy.h"
16 #include "softsusy.h"
17 #include "physpars.h"
18 #include "lowe.h"
19 #include "def.h"
20 #include "softpars.h"
21 #include "softsusy.h"
22 #include "flavoursoft.h"
23 #include "susy.h"
24 #include "particle.h"
25 #include "twoBodyDecays.h"
26 #include <iostream>
27 #include <cstring>
28 #include <stdlib.h>
29 #include <cmath>
30 #include <stdlib.h>
31 #include <vector>
32 #include <iomanip>
33 #include <complex>
34 
35 using namespace std;
36 
38 const double accuracy = 0.01;
39 
42 const double mcpole = 1.40;
45 const double mspole = 0.19;
48 const double hadronicScale = 1.5;
49 
50 void printRowPDG(ostream & cout, double x);
51 void OutputNoPWs(ostream & cout, Particle & P);
52 void OutputYesPWs(ostream & cout, Particle & P);
53 
55  double Zsfintegralsum(double m1, double m2, double msf, double mf, double mz, double min, double max, double Nsteps, int adaptive, double approx);
56  double Aintegralsum(double m1, double m2, double mz, double mA, double mf, double min, double max, double Nsteps, int adaptive, double approx);
57  double G3integralsum(double m1, double m2, double m3, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
58  double G2integralsum(double m1, double m2, double m3, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
59 double chiprimeintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
60  double Yintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
61  double Xintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
62  double zetaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
63  double chiintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
64  double rhointegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
65  double xsiintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
66  double phitildaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
67  double psitildaintegralsum(double m1, double m2, double m3, double m4, double mq, double min, double max, double Nsteps, int adaptive, double approx);
68 
69  double Zintegralsum(double m1, double m2, double mz, double mf, double min, double max, double Nsteps, int adaptive, double approx);
70  double G1integralsum(double m1, double m2, double m3, double mq, double min, double max, double Nsteps, int adaptive, double approx);
71 
72  double G4integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
73  double G5integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
74  double G6integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
75  double G7integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
76  double G8integralsum(double m1, double m2, double m3, double m4, double mt, double mb, double min, double max, double Nsteps, int adaptive, double approx);
77 
78  double Jintegralsum(double m1, double m2, double msf, double mphi, double mf, double min, double max, double Nsteps, double approx, int adaptive, int AorhorH);
79 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);
80 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);
82  //only done via usual method (not via a separate integrating function compute_areai) as can't pass a DoubleMatrix & in a function pointer
83 
84 
85 
87  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);
88  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);
89  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);
90  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);
91  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);
92 
93 
96 double fdgauss(double x);
97 double gpsitildadgauss(double Et);
98 double gphitildadgauss(double Et);
99 double gxsidgauss (double Et);
100 double grhodgauss (double Et);
101 double gchidgauss (double Et);
102 double gzetadgauss (double Et);
103 double gXdgauss (double Et);
104 double gYdgauss (double Et);
105 double gchiprimedgauss (double Et);
106 double gG1dgauss (double Et);
107 double gG4dgauss (double Et);
108 double gG5dgauss(double Et);
109 double gG6dgauss(double Et);
110 double gG7dgauss(double Et);
111 double gG8dgauss(double Et);
112 double gG2dgauss(double Eb);
113 double gG3dgauss(double Eb);
114 double gZdgauss(double E);
115 double ghHdgauss (double E);
116 double gAdgauss (double E);
117 double gZsfdgauss (double s);
118 double gJdgauss(double s);
119 
120 double gneutineutjffZ1dgauss(double s);
121 double gneutineutjffZ2dgauss(double s);
122 double gneutineutjffZ3dgauss(double s);
123 double gneutineutjffZ4dgauss(double s);
124 
125 double gintegralhdgauss(double E);
126 double gintegralHdgauss(double E);
127 double gintegralh1dgauss(double E);
128 double gintegralh2dgauss(double E);
129 double gintegralh3dgauss(double E);
130 double gintegralh4dgauss(double E);
131 double gintegralH1dgauss(double E);
132 double gintegralH2dgauss(double E);
133  double gintegralH3dgauss(double E);
134  double gintegralH4dgauss(double E);
135  double gintegralhH1dgauss(double E);
136  double gintegralhH2dgauss(double E);
137  double gintegralhH3dgauss(double E);
138  double gintegralhH4dgauss(double E);
139  double gintegralA1dgauss(double E);
140  double gintegralA2dgauss(double E);
141  double gintegralA3dgauss(double E);
142  double gintegralA4dgauss(double E);
143  double gintegral1Zsfdgauss(double E);
144  double gintegral2Zsfdgauss(double E);
145  double gintegral3Zsfdgauss(double E);
146  double gintegral4Zsfdgauss(double E);
147  double gintegral5Zsfdgauss(double E);
148  double gintegral6Zsfdgauss(double E);
149  double gintegral7Zsfdgauss(double E);
150  double gintegral8Zsfdgauss(double E);
151  double gintegral1hsfdgauss(double E);
152  double gintegral2hsfdgauss(double E);
153  double gintegral3hsfdgauss(double E);
154  double gintegral4hsfdgauss(double E);
155  double gintegral5hsfdgauss(double E);
156  double gintegral6hsfdgauss(double E);
157  double gintegral7hsfdgauss(double E);
158  double gintegral8hsfdgauss(double E);
159  double gintegral1Hsfdgauss(double E);
160  double gintegral2Hsfdgauss(double E);
161  double gintegral3Hsfdgauss(double E);
162  double gintegral4Hsfdgauss(double E);
163  double gintegral5Hsfdgauss(double E);
164  double gintegral6Hsfdgauss(double E);
165  double gintegral7Hsfdgauss(double E);
166  double gintegral8Hsfdgauss(double E);
167  double gintegral1ZAdgauss(double E);
168  double gintegral2ZAdgauss(double E);
169  double gintegral3ZAdgauss(double E);
170  double gintegral4ZAdgauss(double E);
171  double gneutineutjffgA1dgauss(double E);
172  double gneutineutjffgA2dgauss(double E);
173  double gneutineutjffgA3dgauss(double E);
174  double gneutineutjffgA4dgauss(double E);
175 
176  double gneuticharjffpW1dgauss(double E);
177  double gneuticharjffpW2dgauss(double E);
178  double gneuticharjffpHpm1dgauss(double E);
179  double gneuticharjffpHpm2dgauss(double E);
180  double gneuticharjffpHpm3dgauss(double E);
181  double gneuticharjffpHpm4dgauss(double E);
182  double gneuticharjffp1sf1sf2dgauss(double E);
183  double gneuticharjffp2sf1sf2dgauss(double E);
184  double gneuticharjffp3sf1sf2dgauss(double E);
185  double gneuticharjffp4sf1sf2dgauss(double E);
186  double gneuticharjffp1sfp1sf2dgauss(double E);
187  double gneuticharjffp2sfp1sf2dgauss(double E);
188  double gneuticharjffp3sfp1sf2dgauss(double E);
189  double gneuticharjffp4sfp1sf2dgauss(double E);
190  double gneuticharjffp5sfp1sf2dgauss(double E);
191  double gneuticharjffp6sfp1sf2dgauss(double E);
192  double gneuticharjffp7sfp1sf2dgauss(double E);
193  double gneuticharjffp8sfp1sf2dgauss(double E);
194  double gneuticharjffp1WHpmdgauss(double E);
195  double gneuticharjffp2WHpmdgauss(double E);
196  double gneuticharjffp3WHpmdgauss(double E);
197  double gneuticharjffp4WHpmdgauss(double E);
198  double gneuticharjffpW1Sfpdgauss(double E);
199  double gneuticharjffpW2Sfpdgauss(double E);
200  double gneuticharjffpW3Sfpdgauss(double E);
201  double gneuticharjffpW4Sfpdgauss(double E);
202  double gneuticharjffpW5Sfpdgauss(double E);
203  double gneuticharjffpW6Sfpdgauss(double E);
204  double gneuticharjffpW7Sfpdgauss(double E);
205  double gneuticharjffpW8Sfpdgauss(double E);
206  double gneuticharjffpHg1dgauss(double E);
207  double gneuticharjffpHg2dgauss(double E);
208  double gneuticharjffpHg3dgauss(double E);
209  double gneuticharjffpHg4dgauss(double E);
210  double gneuticharjffp1gsfpdgauss(double E);
211  double gneuticharjffp2gsfpdgauss(double E);
212  double gneuticharjffp3gsfpdgauss(double E);
213  double gneuticharjffp4gsfpdgauss(double E);
214  double gneuticharjffp5gsfpdgauss(double E);
215  double gneuticharjffp6gsfpdgauss(double E);
216  double gneuticharjffp7gsfpdgauss(double E);
217  double gneuticharjffp8gsfpdgauss(double E);
218  double gneuticharjffp1sfpsfpdgauss(double E);
219  double gneuticharjffp21sfpsfpdgauss(double E);
220  double gneuticharjffp31sfpsfpdgauss(double E);
221  double gneuticharjffp41sfpsfpdgauss(double E);
222 
223  DoubleVector squarkmixcharginocouplings (double g, double theta, double beta, double gammaL, double gammaR, double runmt, double runmb, double mWboson, int torb);
224  DoubleVector higgsphisamecharginocouplings(double alpha, double beta, double thetaL, double thetaR);
225  DoubleVector higgsphidifcharginocouplings (double alpha, double beta, double thetaL, double thetaR);
226  DoubleVector higgshsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
227  DoubleVector higgshsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
228  DoubleVector higgsHsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
229  DoubleVector higgsHsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
230  DoubleVector higgshsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel);
231  DoubleVector higgshsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
232  DoubleVector higgsHsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel) ;
233  DoubleVector higgsHsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
234  DoubleVector higgsHplussquarkcouplings (double mWboson, double g, double beta, double mup, double mdown, double greekmu, double Au, double Ad);
235  DoubleVector squarkmixcharginocouplings (double g, double theta, double beta, double gammaL, double gammaR, double runmt, double runmb, double mWboson, double mch1, double mch2, int torb);
236 DoubleVector higgsphisamecharginocouplings (double alpha, double beta, double thetaL, double thetaR);
237 DoubleVector higgsphidifcharginocouplings (double alpha, double beta, double thetaL, double thetaR);
238 DoubleVector higgshsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
239 DoubleVector higgshsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
240 DoubleVector higgsHsquarksamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mupq, double mdownq);
241 DoubleVector higgsHsquarkdiffhandcouplings (double mWboson, double g, double alpha, double beta, double mupq, double mdownq, double greekmu, double Aup, double Adown);
242 DoubleVector higgshsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel);
243 DoubleVector higgshsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
244 DoubleVector higgsHsleptonsamehandcouplings (double mWboson, double g, double gp, double alpha, double beta, double mel);
245 DoubleVector higgsHsleptondiffhandcouplings (double mWboson, double g, double alpha, double beta, double mel, double greekmu, double Ae);
246 DoubleVector higgsHplussquarkcouplings (double mWboson, double g, double beta, double mup, double mdown, double greekmu, double Au, double Ad);
248 DoubleVector foftau(double mpart, double mcomp);
249 DoubleVector goftau(double mpart, double mcomp);
250 double chToN2piInt(double qSq, const DoubleVector & v);
251 double charginoToNeutralino2pion(const MssmSoftsusy * m);
252 Complex fofqsq(double qSq);
253 Complex bw(double mSq, double gamma, double qSq);
254 
255 #endif
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: threeBodyDecays.cpp:6143
double gintegralhH1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: threeBodyDecays.cpp:692
double gZdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt, MZboson = mz
Definition: threeBodyDecays.cpp:459
double gpsitildadgauss(double Et)
Definition: threeBodyDecays.cpp:217
double ghHdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt, MZboson = mz
Definition: threeBodyDecays.cpp:473
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: threeBodyDecays.cpp:6158
double gneuticharjffp2sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: threeBodyDecays.cpp:1428
Two loop corrections to the bottom mass. This file has been generated at Fri 7 Apr 2017 21:00:53 with...
double gneutineutjffZ2dgauss(double s)
Definition: threeBodyDecays.cpp:551
DoubleVector is of variable length, and contains double precision.
Definition: linalg.h:35
double gintegral3ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: threeBodyDecays.cpp:1214
double gneutineutjffgA1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: threeBodyDecays.cpp:1244
double gneuticharjffpHpm1dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm; //sometimes mass order differs as c...
Definition: threeBodyDecays.cpp:1341
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: threeBodyDecays.cpp:794
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: threeBodyDecays.cpp:6071
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: threeBodyDecays.cpp:1076
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: threeBodyDecays.cpp:6056
double gneuticharjffpHpm3dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm; //sometimes mass order differs as c...
Definition: threeBodyDecays.cpp:1375
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: threeBodyDecays.cpp:6085
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: threeBodyDecays.cpp:875
const double hadronicScale
Definition: threeBodyDecays.h:48
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: threeBodyDecays.cpp:1167
switches (options) and parameters such as default fermion masses, required accuracy etc ...
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: threeBodyDecays.cpp:2278
double gintegralHdgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: threeBodyDecays.cpp:612
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: threeBodyDecays.cpp:911
double gneuticharjffpW5Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1825
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: threeBodyDecays.cpp:3115
double gneuticharjffp1WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: threeBodyDecays.cpp:1661
QedQcd object contains Standard Model quark and lepton masses. It integrates them using 3 loop qcd x ...
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)
Hintegral.
Definition: threeBodyDecays.cpp:2224
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: threeBodyDecays.cpp:843
double gneutineutjffgA3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: threeBodyDecays.cpp:1274
double gintegralhH4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: threeBodyDecays.cpp:718
double gneuticharjffpW1Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1736
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: threeBodyDecays.cpp:6115
double gintegralh2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:630
double gintegralA2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: threeBodyDecays.cpp:735
double gneutineutjffZ1dgauss(double s)
m1 = mneuti, m4 = mneutj, mq = mf, MZboson = mZboson
Definition: threeBodyDecays.cpp:534
double gneuticharjffpHg3dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: threeBodyDecays.cpp:1945
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: threeBodyDecays.cpp:967
double gneuticharjffp1sfpsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = -mf, m4 = charginoj, m5 = msfp1, m6 = msfp2
Definition: threeBodyDecays.cpp:2154
DoubleVector goftau(double mpart, double mcomp)
g(tau) function for use in h->Z gamma
Definition: threeBodyDecays.cpp:6222
double gintegralh4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:649
MssmSusy contains all SUSY couplings and tan beta, as well as their beta functions.
double gintegral4ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: threeBodyDecays.cpp:1229
double gneuticharjffp8gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2132
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: threeBodyDecays.cpp:893
double gintegralhH2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: threeBodyDecays.cpp:700
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: threeBodyDecays.cpp:1057
double gintegralh3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:640
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: threeBodyDecays.cpp:859
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: threeBodyDecays.cpp:2444
double gneuticharjffp4sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: threeBodyDecays.cpp:1460
double gneutineutjffgA2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: threeBodyDecays.cpp:1259
DoubleVector foftau(double mpart, double mcomp)
Function Declarations used to calculate Partial Widths in decays.cpp.
Definition: threeBodyDecays.cpp:6192
double gAdgauss(double E)
m1 = mZi, m4 = mZj, m2 = mstop1, m3 = mstop2, mq = mt
Definition: threeBodyDecays.cpp:492
double gneuticharjffp2sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1501
double gintegralA3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: threeBodyDecays.cpp:744
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: threeBodyDecays.cpp:827
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: threeBodyDecays.cpp:811
double gintegralA1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: threeBodyDecays.cpp:726
double gneuticharjffpHg4dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: threeBodyDecays.cpp:1961
double gneuticharjffp1sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1477
double gG6dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: threeBodyDecays.cpp:392
double gneuticharjffp2gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2000
double gintegralh1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:622
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: threeBodyDecays.cpp:6101
Header file for RP conserving MSSM object including all (real) soft SUSY breaking parameters and (rea...
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: threeBodyDecays.cpp:930
double gneuticharjffp3WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: threeBodyDecays.cpp:1698
double gneuticharjffp1sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: threeBodyDecays.cpp:1412
double gZsfdgauss(double s)
m2 = msfi where i is the sfermion contribution interference with Z you are considering ...
Definition: threeBodyDecays.cpp:502
double gneutineutjffZ4dgauss(double s)
Definition: threeBodyDecays.cpp:585
double gG1dgauss(double Et)
m1 = mgluino, m2 = mstopi, m6 = mtop, m8 = mcharginoj
Definition: threeBodyDecays.cpp:359
Matrix from 1..rows, 1..cols of double values.
Definition: linalg.h:211
double gneuticharjffp1gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:1978
double fdgauss(double x)
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: threeBodyDecays.cpp:5922
const double mspole
Definition: threeBodyDecays.h:45
double gneuticharjffp8sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1638
double gneutineutjffZ3dgauss(double s)
Definition: threeBodyDecays.cpp:568
double gneuticharjffpW2dgauss(double E)
m1 = mZi, m2 = mWj, m3 = mf, m4 = mfp, MWboson = mw
Definition: threeBodyDecays.cpp:1323
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: threeBodyDecays.cpp:1095
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 gneuticharjffp2WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: threeBodyDecays.cpp:1679
Particle class definition for decays.
Definition: particle.h:46
double gneuticharjffp4gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2044
double gneuticharjffp3gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2022
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: threeBodyDecays.cpp:6037
NmssmSoftsusy object contains all the NMSSM functionality: SUSY breaking and preserving parameters...
double gG5dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: threeBodyDecays.cpp:381
Flags for potential problems in sProblem structure, and structure for containing physical MSSM parame...
double gneuticharjffpW3Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1780
double gintegralH1dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: threeBodyDecays.cpp:657
double gneuticharjffpW7Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1869
double gneuticharjffp4WHpmdgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHpm
Definition: threeBodyDecays.cpp:1717
double gneuticharjffpHg1dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: threeBodyDecays.cpp:1913
double gneuticharjffpHpm2dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm; //sometimes mass order differs as c...
Definition: threeBodyDecays.cpp:1358
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: threeBodyDecays.cpp:6130
double gneuticharjffpHg2dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mWboson, m6 = mHP
Definition: threeBodyDecays.cpp:1929
double gintegralH4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:683
double gneuticharjffpW8Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1891
double gintegralH2dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: threeBodyDecays.cpp:665
double gG4dgauss(double Et)
m1 = mgluino, m2 = mstopi, m4 = msbottomi, m6 = mtop, m7 = mbottom , m8 = mcharginoj ...
Definition: threeBodyDecays.cpp:369
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: threeBodyDecays.cpp:1003
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: threeBodyDecays.cpp:761
Soft SUSY breaking parameters.
double m1
decay global variable declarations
Definition: decays.cpp:14
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: threeBodyDecays.cpp:985
Contains all supersymmetric MSSM parameters, incorporating R_p MSSM.
Definition: softsusy.h:77
double gneuticharjffpW2Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1758
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: threeBodyDecays.cpp:1021
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: threeBodyDecays.cpp:6015
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: threeBodyDecays.cpp:6171
double gneuticharjffp7gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2110
double gneuticharjffp3sf1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msf1, m6 = msf2
Definition: threeBodyDecays.cpp:1444
double gneuticharjffp4sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1546
double gintegral1ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: threeBodyDecays.cpp:1185
void OutputNoPWs(ostream &cout, Particle &P)
Outputs the decay table into the leshouchesOutput file with no PWs (partial widths) given...
Definition: threeBodyDecays.cpp:166
Definitions of particles and container for decay widths.
const double mcpole
Definition: threeBodyDecays.h:42
double gneuticharjffpHpm4dgauss(double E)
m1 = mneutralinoi, m2 = mcharginoj, m3 = mf, m4 = mfp, m5 = mHpm; //sometimes mass order differs as c...
Definition: threeBodyDecays.cpp:1392
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: threeBodyDecays.cpp:1149
double gintegralA4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mA = mhiggsA
Definition: threeBodyDecays.cpp:753
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: threeBodyDecays.cpp:949
double gneuticharjffpW1dgauss(double E)
m1 = mZi, m2 = mWj, m3 = mf, m4 = mfp, MWboson = mw
Definition: threeBodyDecays.cpp:1306
void printRowPDG(ostream &cout, double x)
Definition: threeBodyDecays.cpp:152
double gneuticharjffp3sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1524
double gneutineutjffgA4dgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, MZboson = mz
Definition: threeBodyDecays.cpp:1289
double gneuticharjffp5sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1568
double gintegralhH3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl, mH = mhiggsH
Definition: threeBodyDecays.cpp:709
double gneuticharjffpW6Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1847
double gneuticharjffp6sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1590
const double accuracy
Approximate accuracy with which 3 body decays are calculated.
Definition: threeBodyDecays.h:38
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: threeBodyDecays.cpp:1131
double gintegral2ZAdgauss(double E)
m1 = mZi, m4 = mZj, mq = mf, mA = mhiggsA, m2 = msfi, MZboson = mz
Definition: threeBodyDecays.cpp:1199
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: threeBodyDecays.cpp:2570
double gintegralH3dgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mH = mhiggsH
Definition: threeBodyDecays.cpp:674
double gneuticharjffpW4Sfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp
Definition: threeBodyDecays.cpp:1803
double gneuticharjffp7sfp1sf2dgauss(double E)
m1 = mneutralinoi, m2 = mf, m3 = mcharginoj, m4 = mfp, m5 = msfp1, m6 = msf2
Definition: threeBodyDecays.cpp:1612
void OutputYesPWs(ostream &cout, Particle &P)
Outputs the decay table into the leshouchesOutput file with PWs (partial widths) given after the comm...
Definition: threeBodyDecays.cpp:190
int theta(double a)
Standard theta function: 1 is a>0, 0 otherwise.
Definition: utils.cpp:25
double gintegralhdgauss(double E)
m1 = mZi, m4 = mZj, mq = mt, MZboson = mz, mh = mhiggsl
Definition: threeBodyDecays.cpp:603
double gneuticharjffp5gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2066
drop-in replacement for the original home-grown Complex class
Definition: mycomplex.h:17
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: threeBodyDecays.cpp:778
double gneuticharjffp6gsfpdgauss(double E)
m1 = mneutralinoi, m2 = mfp, m3 = mf, m4 = mcharginoj, m5 = mWboson, m6 = msfp1
Definition: threeBodyDecays.cpp:2088
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: threeBodyDecays.cpp:1039
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: threeBodyDecays.cpp:1113