AugerUnits.h

Go to the documentation of this file.

// Auger system of units

#ifndef _utl_AugerUnits_h_
#define _utl_AugerUnits_h_

#include <utl/SymbolTable.h>


namespace utl {

  // Prefixes
  constexpr double yocto = 1e-24;
  constexpr double zepto = 1e-21;
  constexpr double atto  = 1e-18;
  constexpr double femto = 1e-15;
  constexpr double pico  = 1e-12;
  constexpr double nano  = 1e-9;
  constexpr double micro = 1e-6;
  constexpr double milli = 1e-3;
  constexpr double centi = 1e-2;
  constexpr double deci  = 1e-1;
  constexpr double deka  = 1e+1;
  constexpr double hecto = 1e+2;
  constexpr double kilo  = 1e+3;
  constexpr double mega  = 1e+6;
  constexpr double giga  = 1e+9;
  constexpr double tera  = 1e+12;
  constexpr double peta  = 1e+15;
  constexpr double exa   = 1e+18;
  constexpr double zetta = 1e+21;
  constexpr double yotta = 1e+24;

  // Length [L]
  constexpr double meter  = 1;
  constexpr double meter2 = meter*meter;
  constexpr double meter3 = meter*meter*meter;

  constexpr double millimeter  = milli*meter;
  constexpr double millimeter2 = millimeter*millimeter;
  constexpr double millimeter3 = millimeter*millimeter*millimeter;

  constexpr double centimeter  = centi*meter;
  constexpr double centimeter2 = centimeter*centimeter;
  constexpr double centimeter3 = centimeter*centimeter*centimeter;

  constexpr double kilometer  = kilo*meter;
  constexpr double kilometer2 = kilometer*kilometer;
  constexpr double kilometer3 = kilometer*kilometer*kilometer;

  constexpr double parsec = 3.0856775807e+16 * meter;
  constexpr double kiloParsec = kilo * parsec;
  constexpr double megaParsec = mega * parsec;

  constexpr double micrometer = micro*meter;
  constexpr double nanometer  = nano*meter;
  constexpr double angstrom   = 1e-10*meter;
  constexpr double fermi      = femto*meter;

  constexpr double barn      = 1e-28*meter2;
  constexpr double millibarn = milli*barn;
  constexpr double microbarn = micro*barn;
  constexpr double nanobarn  = nano*barn;
  constexpr double picobarn  = pico*barn;

  // symbols
  constexpr double mm  = millimeter;
  constexpr double mm2 = millimeter2;
  constexpr double mm3 = millimeter3;

  constexpr double cm  = centimeter;
  constexpr double cm2 = centimeter2;
  constexpr double cm3 = centimeter3;

  constexpr double m  = meter;
  constexpr double m2 = meter2;
  constexpr double m3 = meter3;

  constexpr double km  = kilometer;
  constexpr double km2 = kilometer2;
  constexpr double km3 = kilometer3;

  // Angle
  constexpr double radian      = 1;
  constexpr double milliradian = milli*radian;
  constexpr double degree      = (3.14159265358979323846/180)*radian;

  constexpr double steradian   = 1;

  // symbols
  constexpr double rad  = radian;
  constexpr double mrad = milliradian;
  constexpr double sr   = steradian;
  constexpr double deg  = degree;

  // Time [T]
  constexpr double nanosecond  = 1;
  constexpr double nanosecond2  = nanosecond*nanosecond;
  constexpr double second      = giga*nanosecond;
  constexpr double millisecond = milli*second;
  constexpr double microsecond = micro*second;
  constexpr double picosecond  = pico*second;
  constexpr double minute      = 60*second;
  constexpr double hour        = 60*minute;
  constexpr double day         = 24*hour;

  constexpr double hertz = 1/second;
  constexpr double kilohertz = kilo*hertz;
  constexpr double megahertz = mega*hertz;

  constexpr double Hz = hertz;
  constexpr double kHz = kilohertz;
  constexpr double MHz = megahertz;

  // symbols
  constexpr double ns = nanosecond;
  constexpr double s  = second;
  constexpr double ms = millisecond;

  // Electric charge [Q]
  constexpr double eplus   = 1;                // positron charge
  constexpr double eSI     = 1.602176462e-19;  // positron charge in coulomb
  constexpr double coulomb = eplus/eSI;        // coulomb = 6.24150 e+18*eplus

  // Energy [E]
  constexpr double electronvolt      = 1;
  constexpr double megaelectronvolt  = mega*electronvolt;
  constexpr double kiloelectronvolt  = kilo*electronvolt;
  constexpr double gigaelectronvolt  = giga*electronvolt;
  constexpr double teraelectronvolt  = tera*electronvolt;
  constexpr double petaelectronvolt  = peta*electronvolt;
  constexpr double exaelectronvolt   = exa*electronvolt;
  constexpr double zettaelectronvolt = zetta*electronvolt;

  constexpr double joule = electronvolt/eSI;  // joule = 6.24150 e+12 * MeV

  // symbols
  constexpr double MeV = megaelectronvolt;
  constexpr double eV  = electronvolt;
  constexpr double keV = kiloelectronvolt;
  constexpr double GeV = gigaelectronvolt;
  constexpr double TeV = teraelectronvolt;
  constexpr double PeV = petaelectronvolt;
  constexpr double EeV = exaelectronvolt;
  constexpr double ZeV = zettaelectronvolt;

  // Mass [E][T^2][L^-2]
  constexpr double kilogram  = joule*second*second/(meter*meter);
  constexpr double gram      = milli*kilogram;
  constexpr double milligram = milli*gram;

  // symbols
  constexpr double kg = kilogram;
  constexpr double g  = gram;
  constexpr double mg = milligram;

  // Power [E][T^-1]
  constexpr double watt = joule/second;  // watt = 6.24150 e+3 * MeV/ns

  // Force [E][L^-1]
  constexpr double newton = joule/meter;  // newton = 6.24150 e+9 * MeV/mm

  // Pressure [E][L^-3]
  //constexpr double hep_pascal = newton/m2;    // pascal = 6.24150 e+3 * MeV/mm3
  //constexpr double pascal     = hep_pascal;   // pascal = 6.24150 e+3 * MeV/mm3
  constexpr double pascal     = newton/m2;      // pascal = 6.24150 e+3 * MeV/mm3
  constexpr double bar        = 100000*pascal;  // bar    = 6.24150 e+8 * MeV/mm3
  constexpr double millibar   = milli*bar;
  constexpr double atmosphere = 101325*pascal;  // atm    = 6.32420 e+8 * MeV/mm3

  // symbols
  constexpr double hPa = hecto*pascal;

  // Electric current [Q][T^-1]
  constexpr double ampere      = coulomb/second;  // ampere = 6.24150 e+9 * eplus/ns
  constexpr double milliampere = milli*ampere;
  constexpr double microampere = micro*ampere;
  constexpr double nanoampere  = nano*ampere;

  // Electric potential [E][Q^-1]
  constexpr double megavolt = megaelectronvolt/eplus;
  constexpr double kilovolt = milli*megavolt;
  constexpr double volt     = micro*megavolt;
  constexpr double millivolt = milli*volt;
  constexpr double microvolt = micro*volt;

  constexpr double V = volt;

  // Electric resistance [E][T][Q^-2]
  constexpr double ohm = volt/ampere;  // ohm = 1.60217e-16*(MeV/eplus)/(eplus/ns)

  // Electric capacitance [Q^2][E^-1]
  constexpr double farad      = coulomb/volt;  // farad = 6.24150e+24 * eplus/Megavolt
  constexpr double millifarad = milli*farad;
  constexpr double microfarad = micro*farad;
  constexpr double nanofarad  = nano*farad;
  constexpr double picofarad  = pico*farad;

  // Magnetic Flux [T][E][Q^-1]
  constexpr double weber = volt*second;  // weber = 1000*megavolt*ns

  // Magnetic Field [T][E][Q^-1][L^-2]
  constexpr double tesla      = volt*second/meter2;  // tesla =0.001*megavolt*ns/mm2
  constexpr double microtesla = micro*tesla;

  constexpr double gauss     = 1e-4*tesla;
  constexpr double kilogauss = deci*tesla;

  // Inductance [T^2][E][Q^-2]
  constexpr double henry = weber/ampere;  // henry = 1.60217e-7*MeV*(ns/eplus)**2

  // Temperature
  constexpr double kelvin = 1;

  // Amount of substance
  constexpr double mole = 1;

  // Activity [T^-1]
  constexpr double becquerel = 1/second;
  constexpr double curie = 3.7e+10 * becquerel;

  // Absorbed dose [L^2][T^-2]
  constexpr double gray = joule/kilogram;

  // Luminous intensity [I]
  constexpr double candela = 1;

  // Luminous flux [I]
  constexpr double lumen = candela*steradian;

  // Illuminance [I][L^-2]
  constexpr double lux = lumen/meter2;

  // Miscellaneous
  constexpr double fraction    = 1;
  constexpr double perCent     = 0.01;
  constexpr double percent     = perCent;
  constexpr double perThousand = 0.001;
  constexpr double permil      = perThousand;
  constexpr double perMillion  = 0.000001;

  // System of unit settings for the CLHEP
  // expression Evaluator.  This puts the Evaluator
  // into the same units system defined in this file.
  // (note Evaluator is in SI by default)


  namespace Evaluator {
    constexpr double kLength = meter;
    constexpr double kMass = kilogram;
    constexpr double kTime = second;
    constexpr double kCurrent = ampere;
    constexpr double kTemperature = kelvin;
    constexpr double kAmount = mole;
    constexpr double kLuminousIntensity = candela;
  }


  class AugerUnits {

  public:
    AugerUnits() { FillMap(); }

    static const SymbolTable& GetSymbolMap() { FillMap(); return fgUnitMap; }

  private:
    static void FillMap();

    static SymbolTable fgUnitMap;
  };

}


#endif