Analysis Software
Documentation for sPHENIX simulation software
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
pyhiw1au.f
Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyhiw1au.f
1 
2 C***********************************************************************
3 
4  FUNCTION pyhiw1au(EPS,IREIM)
5 
6 C...Calculates real and imaginary parts of the auxiliary function W1;
7 C...see R. K. Ellis, I. Hinchliffe, M. Soldate and J. J. van der Bij,
8 C...FERMILAB-Pub-87/100-T, LBL-23504, June, 1987
9  common/ludat1/mstu(200),paru(200),mstj(200),parj(200)
10  SAVE /ludat1/
11 
12  asinh(x)=log(x+sqrt(x**2+1.))
13  acosh(x)=log(x+sqrt(x**2-1.))
14 
15  IF(eps.LT.0.) THEN
16  w1re=2.*sqrt(1.-eps)*asinh(sqrt(-1./eps))
17  w1im=0.
18  ELSEIF(eps.LT.1.) THEN
19  w1re=2.*sqrt(1.-eps)*acosh(sqrt(1./eps))
20  w1im=-paru(1)*sqrt(1.-eps)
21  ELSE
22  w1re=2.*sqrt(eps-1.)*asin(sqrt(1./eps))
23  w1im=0.
24  ENDIF
25 
26  IF(ireim.EQ.1) pyhiw1au=w1re
27  IF(ireim.EQ.2) pyhiw1au=w1im
28 
29  RETURN
30  END