Units.H

Go to the documentation of this file.

 
#ifndef _UNITS_H_
#define _UNITS_H_
 
#include <gsl/gsl_const_cgsm.h>
#include <gsl/gsl_const_mksa.h>
#include "Types.H"
#include "../Definitions.H"
 
namespace criptic {
 
  namespace units {
    
#ifdef CRIPTIC_UNITS_CGS
 
    // CGS units
    static constexpr Real amu
    = GSL_CONST_CGSM_UNIFIED_ATOMIC_MASS;    
    static constexpr Real eV
    = GSL_CONST_CGSM_ELECTRON_VOLT;          
    static constexpr Real pc
    = GSL_CONST_CGSM_PARSEC;                 
    static constexpr Real lyr
    = GSL_CONST_CGSM_LIGHT_YEAR;             
    static constexpr Real mb = 1e-27;        
    static constexpr Real cm = 1.0;          
    static constexpr Real m = 100.0;         
    static constexpr Real kg = 1000.0;       
    static constexpr Real g = 1.0;           
#elif defined(CRIPTIC_UNITS_MKS)
    
    // MKS units
    static constexpr Real amu
    = GSL_CONST_MKSA_UNIFIED_ATOMIC_MASS;    
    static constexpr Real eV
    = GSL_CONST_MKSA_ELECTRON_VOLT;          
    static constexpr Real pc
    = GSL_CONST_MKSA_PARSEC;                 
    static constexpr Real lyr
    = GSL_CONST_MKSA_LIGHT_YEAR;             
    static constexpr Real mb = 1e-31;        
    static constexpr Real cm = 0.01;          
    static constexpr Real m = 1.0;         
    static constexpr Real kg = 1.0;          
    static constexpr Real g = 0.001;         
#endif
 
    // Additional derived units
    static constexpr Real keV = 1e3*eV;        
    static constexpr Real MeV = 1e6*eV;        
    static constexpr Real GeV = 1e9*eV;        
    static constexpr Real TeV = 1e12*eV;       
    static constexpr Real PeV = 1e15*eV;       
    static constexpr Real hr = 3600.0;         
    static constexpr Real day = 24.0 * hr;     
    static constexpr Real yr = 365.25 * day;   
    static constexpr Real kyr = 1e3*yr;        
    static constexpr Real Myr = 1e6*yr;        
  }
}
#endif
// _UNITS_H_