*
* $Id$
*
* $Log$
* Revision 1.1  2000/06/19 20:00:34  eugenio
* Initial revision
*
* Revision 1.1.1.1  1994/10/08  02:21:28  zfiles
* first version of qqlib in CVS
*
*
#include "sys/CLEO_machine.h"
#include "pilot.h"
*CMZ :  1.01/00 03/11/90  16.54.51  by  Paul Avery
*CMZ :  1.00/00 04/06/90  18.55.48  by  Jorge L. Rodriguez
*-- Author :

      DOUBLE PRECISION FUNCTION PIHINT(S)
C...HADRONIC VACUUM POLARIZATION BY LOGARITHMIC INTERPOLATION
C...THIS IS A GOOD APPROXIMATION FAR FROM RESONANCES
C...THE VALUES OF SS ARE CHOSEN TO BE LOGARITHMICALLY
C...DISTRIBUTED ; THE VALUES OF REPIH ARE CALCULATED
C...USING PIH (GAUSSIAN INTEGRATION + ANALYTIC CONTRIBUTIONS OF RES.)
      IMPLICIT DOUBLE PRECISION (A-H,K,O-Z)
      DIMENSION SS(20),REPIH(20)
      DATA LUNTTO/6/
      DATA SS/ 0.100D-01,0.196D-01,0.384D-01,0.753D-01,
     *   0.148D+00,0.289D+00,0.567D+00,0.111D+01,
     *   0.218D+01,0.427D+01,0.837D+01,0.164D+02,
     *   0.321D+02,0.630D+02,0.123D+03,0.242D+03,
     *   0.474D+03,0.930D+03,0.182D+04,0.357D+04/
      DATA REPIH/ -0.101D-03,-0.194D-03,-0.367D-03,-0.673D-03,
     *   -0.118D-02,-0.195D-02,-0.299D-02,-0.425D-02,
     *   -0.569D-02,-0.728D-02,-0.905D-02,-0.110D-01,
     *   -0.131D-01,-0.153D-01,-0.175D-01,-0.198D-01,
     *   -0.221D-01,-0.245D-01,-0.268D-01,-0.292D-01/
      DATA DELTA/0.672944D0/
      COMMON/QED000/IMODDD,LOUT
      IF((S.GT.(0.04D0)).AND.(S.LT.(120.D0))) WRITE (LOUT,5) S
      IF((S.GT.(0.04D0)).AND.(S.LT.(120.D0))) WRITE(LUNTTO,5)S
    5 FORMAT(/,10X,'S =',F7.3,' GEV**2;',/,
     Z 10X,'FOR THIS VALUE PIH MAY BE BADLY APPROXIMATED DUE',/,
     Z 10X,'TO THE PRESENCE OF MANY RESONANCES UNDER W=10 GEV')
    6 S=DABS(S)
      IF(S.LE.SS(1)) GO TO 3
      DO 1 J=2,20
      J2=J
    1 IF(S.LE.SS(J)) GO TO 2
    2 PIHINT=REPIH(J2-1)
     Z     +(REPIH(J2)-REPIH(J2-1))*DLOG(S/SS(J2-1))/DELTA
      RETURN
    3 PIHINT=0.D0
      RETURN
      END