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 valuesAuthor: 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