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