Changeset View
Changeset View
Standalone View
Standalone View
src/EvtGenModels/EvtWilsonCoefficients.cpp
Show All 25 Lines | |||||
#include <stdlib.h> | #include <stdlib.h> | ||||
double li2spence( double ); | double li2spence( double ); | ||||
EvtWilsonCoefficients::EvtWilsonCoefficients() | EvtWilsonCoefficients::EvtWilsonCoefficients() | ||||
{ | { | ||||
int i, j; | int i, j; | ||||
double tmpa[8] = {14. / 23., 16. / 23., 6. / 23., -12. / 23., | double tmpa[8] = { 14. / 23., 16. / 23., 6. / 23., -12. / 23., | ||||
0.4086, -0.4230, -0.8994, 0.1456}; | 0.4086, -0.4230, -0.8994, 0.1456 }; | ||||
double tmph[8] = {2.2996, -1.0880, -3. / 7., -1. / 14., | double tmph[8] = { 2.2996, -1.0880, -3. / 7., -1. / 14., | ||||
-0.6494, -0.0380, -0.0186, -0.0057}; | -0.6494, -0.0380, -0.0186, -0.0057 }; | ||||
double tmpp[8] = {0, 0, -80. / 203., 8. / 33., | double tmpp[8] = { 0, 0, -80. / 203., 8. / 33., | ||||
0.0433, 0.1384, 0.1648, -0.0073}; | 0.0433, 0.1384, 0.1648, -0.0073 }; | ||||
double tmps[8] = {0, 0, -0.2009, -0.3579, 0.0490, -0.3616, -0.3554, 0.0072}; | double tmps[8] = { 0, 0, -0.2009, -0.3579, | ||||
0.0490, -0.3616, -0.3554, 0.0072 }; | |||||
double tmpq[8] = {0, 0, 0, 0, 0.0318, 0.0918, -0.2700, 0.0059}; | double tmpq[8] = { 0, 0, 0, 0, 0.0318, 0.0918, -0.2700, 0.0059 }; | ||||
double tmpg[8] = {313063. / 363036., 0, 0, 0, | double tmpg[8] = { 313063. / 363036., 0, 0, 0, | ||||
-0.9135, 0.0873, -0.0571, -0.0209}; | -0.9135, 0.0873, -0.0571, -0.0209 }; | ||||
double tmpk[6][8] = { | double tmpk[6][8] = { | ||||
{0, 0, 1. / 2., -1. / 2., 0, 0, 0, 0}, | { 0, 0, 1. / 2., -1. / 2., 0, 0, 0, 0 }, | ||||
{0, 0, 1. / 2., +1. / 2., 0, 0, 0, 0}, | { 0, 0, 1. / 2., +1. / 2., 0, 0, 0, 0 }, | ||||
{0, 0, -1. / 14., +1. / 6., 0.0510, -0.1403, -0.0113, 0.0054}, | { 0, 0, -1. / 14., +1. / 6., 0.0510, -0.1403, -0.0113, 0.0054 }, | ||||
{0, 0, -1. / 14., -1. / 6., 0.0984, +0.1214, +0.0156, 0.0026}, | { 0, 0, -1. / 14., -1. / 6., 0.0984, +0.1214, +0.0156, 0.0026 }, | ||||
{0, 0, 0, 0, -0.0397, 0.0117, -0.0025, +0.0304}, | { 0, 0, 0, 0, -0.0397, 0.0117, -0.0025, +0.0304 }, | ||||
{0, 0, 0, 0, +0.0335, 0.0239, -0.0462, -0.0112}}; | { 0, 0, 0, 0, +0.0335, 0.0239, -0.0462, -0.0112 } }; | ||||
double tmpr[2][8] = { | double tmpr[2][8] = { | ||||
{0, 0, +0.8966, -0.1960, -0.2011, 0.1328, -0.0292, -0.1858}, | { 0, 0, +0.8966, -0.1960, -0.2011, 0.1328, -0.0292, -0.1858 }, | ||||
{0, 0, -0.1193, +0.1003, -0.0473, 0.2323, -0.0133, -0.1799}}; | { 0, 0, -0.1193, +0.1003, -0.0473, 0.2323, -0.0133, -0.1799 } }; | ||||
for ( i = 0; i < 8; i++ ) { | for ( i = 0; i < 8; i++ ) { | ||||
a[i] = tmpa[i]; | a[i] = tmpa[i]; | ||||
h[i] = tmph[i]; | h[i] = tmph[i]; | ||||
p[i] = tmpp[i]; | p[i] = tmpp[i]; | ||||
s[i] = tmps[i]; | s[i] = tmps[i]; | ||||
q[i] = tmpq[i]; | q[i] = tmpq[i]; | ||||
g[i] = tmpg[i]; | g[i] = tmpg[i]; | ||||
for ( j = 0; j < 6; j++ ) | for ( j = 0; j < 6; j++ ) | ||||
▲ Show 20 Lines • Show All 515 Lines • Show Last 20 Lines |