Module physicalconstants_mod

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


Description of Variables

LightSpeed

real (kind=r8), public, PARAMETER :: LightSpeed = 29.979245800
 cm/ns

PiMass

real (kind=r8), public, PARAMETER :: PiMass = 139.57018
 MeV/c^2

MuMass

real (kind=r8), public, PARAMETER :: MuMass = 105.6583568
 MeV/c^2

ElMass

real (kind=r8), public, PARAMETER :: ElMass = 0.510998902
 MeV/c^2

ElMass_Gev

real (kind=r8), public, PARAMETER :: ElMass_Gev = ElMass / 1000.
 GeV/c^2

MuMass_GeV

real (kind=r8), public, PARAMETER :: MuMass_GeV = MuMass / 1000.
 GeV/c^2

Pi

real (kind=R8), public, PARAMETER :: Pi = 3.1415926535897932384626433832795029
 Pi

TwoPi

real (kind=R8), public, PARAMETER :: TwoPi = 2. * Pi
 2 * Pi

raddeg

real (kind=R8), public, PARAMETER :: raddeg = 180. / Pi
 radians to degrees

degrad

real (kind=R8), public, PARAMETER :: degrad = 1. / raddeg
 degrees to radians

rt2

real (kind=R8), public, PARAMETER :: rt2 = 1.41421356237309504880
 sqrt(2.)

Clight

real (kind=R8), public, PARAMETER :: Clight = lightspeed / 10.
 Mev/c /cm Tesla

Clight_GeV

real (kind=R8), public, PARAMETER :: Clight_GeV = lightspeed / 100.
 Gev/c /m Tesla

Emax

real (kind=R8), public, PARAMETER :: Emax = Mumass_gev / 2. * (1. + (ElMass / MuMass) ** 2)
Kinematic limit in GeV