TEvtGen/EvtGenBase/EvtDiLog.hh

   1 // $Id: EvtDiLog.hh,v 1.3 2009-03-16 16:44:33 robbep Exp $
   2 #ifndef EVTGENBASE_EVTDILOG_HH
   3 #define EVTGENBASE_EVTDILOG_HH 1
   4
   5 // Include files
   6
   7 /** @namespace EvtDiLog EvtDiLog.hh EvtGenBase/EvtDiLog.hh
   8  *  Dilogarithm function (replaces CERNLIB DDILOG)
   9  *
  10  *  @author Patrick Robbe
  11  *  @date   2007-01-23
  12  */
  13 namespace EvtDiLog {
  14   double DiLog( double x ) ;
  15
  16   // constants for computation
  17   static const double Z1 = 1. ;
  18   static const double HF = Z1 / 2. ;
  19   static const double PI = 3.14159265358979324 ;
  20   static const double PI3 = PI * PI / 3. ;
  21   static const double PI6 = PI * PI / 6. ;
  22   static const double PI12 = PI * PI / 12. ;
  23   static const double C[ 20 ]  = {
  24     0.42996693560813697 ,
  25     0.40975987533077105 ,
  26     -0.01858843665014592 ,
  27     0.00145751084062268 ,
  28     -0.00014304184442340 ,
  29     0.00001588415541880 ,
  30     -0.00000190784959387 ,
  31     0.00000024195180854 ,
  32     -0.00000003193341274 ,
  33     0.00000000434545063 ,
  34     -0.00000000060578480 ,
  35     0.00000000008612098 ,
  36     -0.00000000001244332 ,
  37     0.00000000000182256 ,
  38     -0.00000000000027007 ,
  39     0.00000000000004042 ,
  40     -0.00000000000000610 ,
  41     0.00000000000000093 ,
  42     -0.00000000000000014 ,
  43     0.00000000000000002 } ;
  44 }
  45 #endif // EVTGENBASE_EVTDILOG_HH