PhysicalConstants.h

Go to the documentation of this file.

#ifndef UTIL_PHYSICALCONSTANTS_H
#define UTIL_PHYSICALCONSTANTS_H

namespace util {

  constexpr double kRecombA = 0.800;  
  constexpr double kRecombk = 0.0486; 

  constexpr double kModBoxA = 0.930; 
  constexpr double kModBoxB = 0.212; 
  // constexpr double kWion          = 23.6e-6;  ///< ionization potenial in LAr, 23.6 eV = 1e, Wion in MeV/e

  // Conversion for energy deposited in GeV to number of ionization electrons produced
  constexpr double kGeVToElectrons = 4.237e7; 

  // More constants
  constexpr double kc = 29.9792458;
  //static double khbar = 6.58211899e-22;  ///< MeVs

  // Conversion factors
  constexpr double kMeterToCentimeter = 1.e2; 
  constexpr double kCentimeterToMeter = 1. / kMeterToCentimeter;
  constexpr double kMeterToKilometer = 1.e-3; 
  constexpr double kKilometerToMeter = 1. / kMeterToKilometer;

  constexpr double keVToMeV = 1.e-6; 
  constexpr double kMeVToeV = 1. / keVToMeV;

  constexpr double kBogusD = -999.; 
  constexpr int kBogusI = -999;     
  constexpr float kBogusF = -999.;  

  constexpr double quietCompiler =
    kBogusD * kBogusI * kBogusF * kRecombA * kRecombk * kGeVToElectrons;

  template <typename T = double>
  inline constexpr T pi()
  {
    return 3.14159265358979323846264338327950288L;
  }

  template <typename T>
  inline constexpr T DegreesToRadians(T angle)
  {
    return angle / 180 * pi<T>();
  }

  template <typename T>
  inline constexpr T RadiansToDegrees(T angle)
  {
    return angle / pi<T>() * 180;
  }

} // namespace util

#endif //UTIL_PHYSICALCONSTANTS_H