lwbb.F

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

C ********************************************************************

      SUBROUTINE lwbb(ENU)

      IMPLICIT NONE
      REAL emean,slope,emin,emax,enu 
      REAL a1,a2,e
      REAL rlu

C...Gives energy (ENU) of a (anti-)neutrino chosen from a simple
C...parametrization of a wide band beam.

      DATA emean,slope,emin,emax/30.,0.02,12.,300./
      a1=1./(emean-12.)
      a2=exp(emean*slope)
  100 enu=300.*rlu(0)
      IF(enu.LT.emean)THEN
      e=a1*(enu-12.)
      ELSE
      e=a2*exp(-enu*slope)
      ENDIF
      IF(enu.LT.emin.OR.enu.GT.emax) goto 100
      IF(e.LT.rlu(0)) goto 100
      RETURN
      END