52public:
54
55#ifdef USE_GOLEM_MODE
56# define EVALUNDEF return std::numeric_limits<double>::signaling_NaN();
59 virtual ncomplex A(
int ep,
int i,
int j ) { EVALUNDEF }
60 virtual ncomplex A(
int ep,
int i,
int j,
int k ) { EVALUNDEF }
61 virtual ncomplex A(
int ep,
int i,
int j,
int k,
int l ) { EVALUNDEF }
62 virtual ncomplex A(
int ep,
int i,
int j,
int k,
int l,
int m ) { EVALUNDEF }
65 virtual ncomplex B(
int ep,
int i,
int j ) { EVALUNDEF }
66 virtual ncomplex B(
int ep,
int i,
int j,
int k ) { EVALUNDEF }
69# undef EVALUNDEF
70#endif
71
72
73 inline static int ns(
int i,
int j )
CONST {
74 return ( i <= j ? ( i - 1 ) + ( ( j - 1 ) * j ) / 2 : ( j - 1 ) + ( ( i - 1 ) * i ) / 2 );
75 }
76
77 inline static int nss(
int i,
int j )
CONST
78 {
79 return ( i - 1 ) + ( ( j - 1 ) * j ) / 2;
80 }
81
82
83 inline static int is(
int i,
int j )
CONST {
84 return ( i <= j ? i + j * ( j + 1 ) / 2 : j + i * ( i + 1 ) / 2 );
85 }
86
87 inline static int is(
int i,
int j,
int k )
CONST {
88 if ( i <= j ) {
return ( j <= k ? i + ti2[j] + ti3[k] :
is( i, k ) + ti3[j] ); }
89 else {
return ( i > k ?
is( j, k ) + ti3[i] : j + ti2[i] + ti3[k] ); }
90 }
91
92 inline static int is(
int i,
int j,
int k,
int l )
CONST {
93 if ( i <= j )
94 {
95 if ( j <= k )
96 {
return ( k <= l ? i + ti2[j] + ti3[k] + ti4[l] :
is( i, j, l ) + ti4[k] ); }
97 else {
return ( j > l ?
is( i, k, l ) + ti4[j] :
is( i, k ) + ti3[j] + ti4[l] ); }
98 }
99 else
100 {
101 if ( i > k )
102 {
return ( i > l ?
is( j, k, l ) + ti4[i] :
is( j, k ) + ti3[i] + ti4[l] ); }
103 else {
return ( k <= l ? j + ti2[i] + ti3[k] + ti4[l] :
is( i, j, l ) + ti4[k] ); }
104 }
105 }
std::complex< double > ncomplex
***************************************************************************************Pseudo Class RRes *****************************************************************************************Parameters and physical constants **Maarten sept ************************************************************************DOUBLE PRECISION xsmu **************************************************************************PARTICLE DATA all others are from PDG *Only resonances with known widths into electron pairs are sept ************************************************************************C Declarations C
static int nss(int i, int j) CONST