Unitsystem.h

// -*- C++ -*-
//
// Unitsystem.h is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2019 Leif Lonnblad, David Grellscheid
//
// ThePEG is licenced under version 3 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
#ifndef ThePEG_Units_H
#define ThePEG_Units_H
 
#include "ThePEG/Vectors/Lorentz5Vector.fh"
#include "ThePEG/Vectors/LorentzVector.fh"
#include "ThePEG/Vectors/ThreeVector.fh"
#include "ThePEG/Vectors/Transverse.fh"
 
#include "PhysicalQty.h"
#include "PhysicalQtyOps.h"
#include "PhysicalQtyComplex.h"
 
 
namespace ThePEG {
 
namespace Units {
 
template<long int L, long int E, long int Q, long int DL=1, long int DE=1, long int DQ=1>
using Qty = ThePEG::Qty<std::ratio<L,DL>, std::ratio<E,DE>, std::ratio<Q,DQ>>;
 
typedef Qty<0,1,0> Energy;
 
typedef Energy Mass;
 
typedef Qty<1,0,0> Length;
 
typedef Length Time;
 
typedef Qty<-1,0,0> InvLength;
 
typedef double Velocity;
 
typedef Qty<0,0,1> Charge;
 
typedef Qty<1,1,0> AngularMomentum;
 
typedef Qty<-1,1,0> Tension;
 
typedef Qty<2,0,0> Area;
 
typedef Qty<-2,0,0> InvArea;
 
typedef Area CrossSection;
 
typedef Qty<0, 2, 0> Energy2;
typedef Qty<0, 3, 0> Energy3;
typedef Qty<0, 4, 0> Energy4;
typedef Qty<0, 5, 0> Energy5;
typedef Qty<0, 6, 0> Energy6;
typedef Qty<0, 7, 0> Energy7;
typedef Qty<0, 8, 0> Energy8;
typedef Qty<0, 9, 0> Energy9;
typedef Qty<0,10, 0> Energy10;
typedef Qty<0,11, 0> Energy11;
typedef Qty<0,12, 0> Energy12;
 
typedef Qty<0, 1,0, 1,2,1> SqrtEnergy;
typedef Qty<0,-1,0, 1,2,1> InvSqrtEnergy;
 
typedef Qty<0, -1, 0> InvEnergy;
typedef Qty<0, -2, 0> InvEnergy2;
typedef Qty<0, -3, 0> InvEnergy3;
typedef Qty<0, -4, 0> InvEnergy4;
typedef Qty<0, -5, 0> InvEnergy5;
typedef Qty<0, -6, 0> InvEnergy6;
typedef Qty<0, -7, 0> InvEnergy7;
typedef Qty<0, -8, 0> InvEnergy8;
typedef Qty<0, -9, 0> InvEnergy9;
typedef Qty<0,-10, 0> InvEnergy10;
typedef Qty<0,-11, 0> InvEnergy11;
typedef Qty<0,-12, 0> InvEnergy12;
 
typedef Qty<2,2,0> Energy2XSec;
 
typedef Qty<2,-2,0> DiffXSec;
 
typedef Qty<2,-4,0> Diff2XSec;
 
typedef Qty<2,-6,0> Diff3XSec;
 
typedef Energy2 Scale;
 
typedef ThreeVector<Length> Point;
 
typedef ThreeVector<Length> Distance;
 
typedef ThreeVector<double> Axis;
 
typedef ThreeVector<Energy> Momentum3;
 
typedef ThreeVector<double> Boost;
 
typedef LorentzVector<Length> LorentzDistance;
 
typedef Lorentz5Vector<Length> Lorentz5Distance;
 
typedef LorentzVector<Length> LorentzPoint;
 
typedef LorentzVector<Energy> LorentzMomentum;
 
typedef Lorentz5Vector<Energy> Lorentz5Momentum;
 
typedef Transverse<Energy> TransverseMomentum;
 
 
constexpr Length operator "" _mm( long double x ) {
  return Length{Length::baseunit(), static_cast<double>(x)};
}
constexpr Length operator "" _mm( unsigned long long x ) {
  return Length{Length::baseunit(), static_cast<double>(x)};
}
 
constexpr Length meter      = 1.0e+3_mm;
constexpr Length millimeter = 1_mm;
constexpr Length mm         = 1_mm;
constexpr Length centimeter = 10_mm;
constexpr Length micrometer = 1.0e-3_mm;
constexpr Length nanometer  = 1.0e-6_mm;
constexpr Length picometer  = 1.0e-9_mm;
constexpr Length femtometer = 1.0e-12_mm;
 
constexpr Energy operator "" _MeV( long double x ) {
  return Energy{Energy::baseunit(), static_cast<double>(x)};
}
constexpr Energy operator "" _MeV( unsigned long long x ) {
  return Energy{Energy::baseunit(), static_cast<double>(x)};
}
 
constexpr Energy operator "" _GeV( long double x ) {
  return Energy{1000_MeV, static_cast<double>(x)};
}
constexpr Energy operator "" _GeV( unsigned long long x ) {
  return Energy{1000_MeV, static_cast<double>(x)};
}
 
constexpr Energy operator "" _TeV( long double x ) {
  return Energy{1000_GeV, static_cast<double>(x)};
}
constexpr Energy operator "" _TeV( unsigned long long x ) {
  return Energy{1000_GeV, static_cast<double>(x)};
}
 
constexpr Energy keV = 1.0e-3_MeV;
constexpr Energy MeV = 1_MeV;
constexpr Energy GeV = 1_GeV;
constexpr Energy TeV = 1_TeV;
 
 
 
constexpr Energy2 operator "" _MeV2( long double x ) {
  return Energy2{Energy2::baseunit(),   static_cast<double>(x)};
}
constexpr Energy2 operator "" _MeV2( unsigned long long x ) {
  return Energy2{Energy2::baseunit(),   static_cast<double>(x)};
}
 
constexpr Energy2 operator "" _GeV2( long double x ) {
  return Energy2{1.0e+6_MeV2, static_cast<double>(x)};
}
constexpr Energy2 operator "" _GeV2( unsigned long long x ) {
  return Energy2{1.0e+6_MeV2, static_cast<double>(x)};
}
 
constexpr Energy2 MeV2 = 1_MeV2;
constexpr Energy2 GeV2 = 1_GeV2;
 
constexpr InvEnergy InvGeV = 1/GeV;
 
 
constexpr Area operator "" _pb( long double x ) {
  return Area{1.0e-34 * Area::baseunit(), static_cast<double>(x)};
}
constexpr Area operator "" _pb( unsigned long long x ) {
  return Area{1.0e-34 * Area::baseunit(), static_cast<double>(x)};
}
 
constexpr Area femtobarn = 1.0e-03_pb;
constexpr Area picobarn  = 1_pb;
constexpr Area nanobarn  = 1.0e+03_pb;
constexpr Area microbarn = 1.0e+06_pb;
constexpr Area millibarn = 1.0e+09_pb;
constexpr Area barn      = 1.0e+12_pb; 
 
constexpr Charge eplus = Charge::baseunit();
 
constexpr Qty<1,1,0> hbarc = 197.326968e-15 * MeV * meter;
constexpr Qty<1,1,0> hbar_Planck = hbarc / 1.0; // c is one
}
 
namespace UnitRemoval {
 
  constexpr Units::Energy E = Units::Energy::baseunit();
 
  constexpr Units::Energy2 E2 = E*E;
  constexpr Units::Energy3 E3 = E*E2;
  constexpr Units::Energy4 E4 = E2*E2;
 
  constexpr Units::InvEnergy InvE = 1.0/E;
  constexpr Units::InvEnergy2 InvE2 = 1.0/E2;
  constexpr Units::InvEnergy3 InvE3 = 1.0/E3;
  constexpr Units::InvEnergy4 InvE4 = 1.0/E4;
 
  constexpr Units::SqrtEnergy SqrtE = Units::SqrtEnergy::baseunit();
  constexpr Units::InvSqrtEnergy InvSqrtE = Units::InvSqrtEnergy::baseunit();
}
 
}
 
#endif /* ThePEG_Units_H */