*=======================================================================
      function s_correct_cal_pos(x,y)
*=======================================================================
*-
*-      Purpose: Returns the impact point correction factor. This
*-               factor is to be applied to the energy depositions.
*-               The final energy is the ADC value TIMES the correction factor.
*-
*-      Input Parameters: x,y - impact point coordinates
*-
*-      Created 15 Mar 1994      Tsolak A. Amatuni
*
* $Log: s_correct_cal_pos.f,v $
* Revision 1.4  2003/03/21 22:58:02  jones
* Subroutines had arguments with abort,errmsg . But these arguments were not
* used when the subroutine was called. Also abort ,errmsg were not used in the
* subroutines. So eliminate abort,errmsg. (E. Brash)
*
* Revision 1.3  1999/06/10 17:04:19  csa
* (JRA) Changed s_correct_cal_pos calculation
*
* Revision 1.2  1999/02/25 20:18:40  saw
* Vardan Tadevosyan shower code updates
*
* Revision 1.1  1999/01/29 17:34:57  saw
* Add variables for second tubes on shower counter
*
*
*-----------------------------------------------------------------------
*
      implicit none
      save
*
*      logical abort
*      character*(*) errmsg
      character*17 here
      parameter (here='S_CORRECT_CAL_POS')
      real a,b,c	! Fit parameters.
      parameter (a=2.3926,b=-0.371375,c=-0.25401)
*
      real*4 x,y         !Impact point coordinates
      real*4 s_correct_cal_pos
      real*4 d,al	! Auxiliary variables.
*
      include 'sos_data_structures.cmn'
      include 'sos_calorimeter.cmn'
*
*     Fit to MC data in the range of y [-30,+30].
*
!      d=y-35.		!need to insure d is never less than -35!!!
!      al=alog(d)
!      s_correct_cal_pos=1./(a+b*al+c/al)

*
*      Fit to data (run23121)
*
      s_correct_cal_pos=exp(-y/210.7) !~200 cm atten length.
      s_correct_cal_pos=s_correct_cal_pos/(1.+y*y/22000.)

      return
      end