test/pythia8test/chi.c

0001 #include <stdio.h>
0002 #include <math.h>
0003
0004 int main(void)
0005 {
0006   double y, dy, y0, dy0;
0007   double chi;
0008
0009   scanf("%lE%lE%lE%lE", &y, &dy, &y0, &dy0);
0010
0011   if (dy == 0. && dy0 == 0.)
0012     chi = -1.;
0013   else
0014     chi = (y - y0) / sqrt(dy*dy + dy0*dy0);
0015
0016   printf("%lE %lE  %lE %lE  %lf ", y, dy, y0, dy0, chi);
0017
0018   if ( (fabs(y0) < 5*dy0 || fabs(y) <= 4*dy) && fabs(chi) <= 3.)
0019     printf("[BADSTAT]\n");
0020   else
0021     if (fabs(chi) > 3. || fabs(chi) <= 3.) { // this if is to deal with NaN
0022       if (fabs(chi) > 3.) printf("[DEVIATION]\n");
0023       if (fabs(chi) <= 3.) printf("[OK]\n");
0024     }
0025     else
0026       printf("[DEVIATION]\n");
0027
0028   return 0;
0029 }