module physicalconstants_mod
        ! Uses
    use precision_mod
        ! Variables
    real (kind=r8), public, PARAMETER :: LightSpeed = 29.979245800
    real (kind=r8), public, PARAMETER :: PiMass = 139.57018
    real (kind=r8), public, PARAMETER :: MuMass = 105.6583568
    real (kind=r8), public, PARAMETER :: ElMass = 0.510998902
    real (kind=r8), public, PARAMETER :: ElMass_Gev = ElMass / 1000.
    real (kind=r8), public, PARAMETER :: MuMass_GeV = MuMass / 1000.
    real (kind=R8), public, PARAMETER :: Pi = 3.1415926535897932384626433832795029
    real (kind=R8), public, PARAMETER :: TwoPi = 2. * Pi
    real (kind=R8), public, PARAMETER :: raddeg = 180. / Pi
    real (kind=R8), public, PARAMETER :: degrad = 1. / raddeg
    real (kind=R8), public, PARAMETER :: rt2 = 1.41421356237309504880
    real (kind=R8), public, PARAMETER :: Clight = lightspeed / 10.
    real (kind=R8), public, PARAMETER :: Clight_GeV = lightspeed / 100.
    real (kind=R8), public, PARAMETER :: Emax = Mumass_gev / 2. * (1. + (ElMass / MuMass) ** 2)
end module physicalconstants_mod
Long past time that we collected these in a module and used consistant values
Author: Art Olin
Version: June 03
real (kind=r8), public, PARAMETER :: LightSpeed = 29.979245800 cm/ns
real (kind=r8), public, PARAMETER :: PiMass = 139.57018 MeV/c^2
real (kind=r8), public, PARAMETER :: MuMass = 105.6583568 MeV/c^2
real (kind=r8), public, PARAMETER :: ElMass = 0.510998902 MeV/c^2
real (kind=r8), public, PARAMETER :: ElMass_Gev = ElMass / 1000. GeV/c^2
real (kind=r8), public, PARAMETER :: MuMass_GeV = MuMass / 1000. GeV/c^2
real (kind=R8), public, PARAMETER :: Pi = 3.1415926535897932384626433832795029 Pi
real (kind=R8), public, PARAMETER :: TwoPi = 2. * Pi 2 * Pi
real (kind=R8), public, PARAMETER :: raddeg = 180. / Pi radians to degrees
real (kind=R8), public, PARAMETER :: degrad = 1. / raddeg degrees to radians
real (kind=R8), public, PARAMETER :: rt2 = 1.41421356237309504880 sqrt(2.)
real (kind=R8), public, PARAMETER :: Clight = lightspeed / 10. Mev/c /cm Tesla
real (kind=R8), public, PARAMETER :: Clight_GeV = lightspeed / 100. Gev/c /m Tesla
real (kind=R8), public, PARAMETER :: Emax = Mumass_gev / 2. * (1. + (ElMass / MuMass) ** 2) Kinematic limit in GeV