pyhiw2au.f

Go to the documentation of this file. Or view the newest version in sPHENIX GitHub for file pyhiw2au.f

    
C***********************************************************************    
    
      FUNCTION pyhiw2au(EPS,IREIM)    
    
C...Calculates real and imaginary parts of the auxiliary function W2;   
C...see R. K. Ellis, I. Hinchliffe, M. Soldate and J. J. van der Bij,   
C...FERMILAB-Pub-87/100-T, LBL-23504, June, 1987    
      common/ludat1/mstu(200),paru(200),mstj(200),parj(200) 
      SAVE /ludat1/ 
    
      asinh(x)=log(x+sqrt(x**2+1.)) 
      acosh(x)=log(x+sqrt(x**2-1.)) 
    
      IF(eps.LT.0.) THEN    
        w2re=4.*(asinh(sqrt(-1./eps)))**2   
        w2im=0. 
      ELSEIF(eps.LT.1.) THEN    
        w2re=4.*(acosh(sqrt(1./eps)))**2-paru(1)**2 
        w2im=-4.*paru(1)*acosh(sqrt(1./eps))    
      ELSE  
        w2re=-4.*(asin(sqrt(1./eps)))**2    
        w2im=0. 
      ENDIF 
    
      IF(ireim.EQ.1) pyhiw2au=w2re    
      IF(ireim.EQ.2) pyhiw2au=w2im    
    
      RETURN    
      END