*=======================================================================
      function s_correct_cal(x,y)
*=======================================================================
*-
*-      Purpose: Returns the impact point correction factor. This
*-               factor is to be applied to the energy depositions.
*-
*-      Input Parameters: x,y - impact point coordinates
*-
*-      Created 15 Mar 1994      Tsolak A. Amatuni
* $Log: s_correct_cal.f,v $
* Revision 1.7  2003/04/03 00:45:01  jones
* Update to calorimeter calibration (V. Tadevosyan)
*
* Revision 1.6  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.5  1999/01/29 17:34:57  saw
* Add variables for second tubes on shower counter
*
* Revision 1.4  1995/05/22 19:45:34  cdaq
* (SAW) Split gen_data_data_structures into gen, hms, sos, and coin parts"
*
* Revision 1.3  1995/05/11  21:02:59  cdaq
* (???) Tweak hardwired attenuation length
*
* Revision 1.2  1994/11/22  21:09:22  cdaq
* (SPB) Recopied from hms file and modified names for SOS
*
* Revision 1.1  1994/04/13  18:10:02  cdaq
* Initial revision
*
*-----------------------------------------------------------------------
*
      implicit none
      save
*
*      logical abort
*      character*(*) errmsg
      character*13 here
      parameter (here='S_CORRECT_CAL')
*
      real*4 x,y         !Impact point coordinates
      real*4 s_correct_cal
*
      include 'sos_data_structures.cmn'
      include 'sos_calorimeter.cmn'
*
*

c     Check calorimeter boundaries.

      if(y.lt.scal_ymin) y=scal_ymin
      if(y.gt.scal_ymax) y=scal_ymax
*
      s_correct_cal=exp(-y/400.)  !400 cm atten length.
      s_correct_cal=s_correct_cal/(1. + y*y/12000)
*
      return
      end