25 const double RELATIVE_PRECISION = 0.00001;
28 const double ABSOLUTE_PRECISION = 0.0;
33 const double MIN_MASS = 0.1;
34 const double ZERO_MASS = 0.000;
35 const double SCANSTEP = 0.1;
44 void set_momenta(
double *pa0,
double *pb0,
double* pmiss0);
45 void set_mn(
double mn);
55 int nsols(
double Dsq);
56 int nsols_massless(
double Dsq);
57 inline int signchange_n(
long double t1,
long double t2,
long double t3,
long double t4,
long double t5);
58 inline int signchange_p(
long double t1,
long double t2,
long double t3,
long double t4,
long double t5);
59 int scan_high(
double &Deltasq_high);
60 int find_high(
double &Deltasq_high);
62 double pax, pay, ma, Ea;
63 double pmissx, pmissy;
64 double pbx, pby, mb, Eb;
65 double mn, mn_unscale;
70 double pmissxsq, pmissysq;
74 double a1, b1, c1, a2, b2, c2, d1, e1, f1, d2, e2, f2;
75 double d11, e11, f12, f10, d21, d20, e21, e20, f22, f21, f20;