LorentzVector.hh

Go to the documentation of this file.

 
#pragma once
 
#include <complex>
#include <utility>
 
#include "CLHEP/Vector/LorentzVector.h"
 
#include "HEJ/Particle.hh"
 
namespace HEJ {
  inline
  auto dot(
    CLHEP::HepLorentzVector const & pi,
    CLHEP::HepLorentzVector const & pj
  ) {
    return pi.dot(pj);
  }
 
  std::complex<double> angle(
    CLHEP::HepLorentzVector const & pi,
    CLHEP::HepLorentzVector const & pj
  );
 
  std::complex<double> square(
    CLHEP::HepLorentzVector const & pi,
    CLHEP::HepLorentzVector const & pj
  );
 
  inline
  auto m2(CLHEP::HepLorentzVector const & h1) {
    return h1.m2();
  }
 
  inline
  auto plus(CLHEP::HepLorentzVector const & h1) {
    return h1.plus();
  }
 
  inline
  auto minus(CLHEP::HepLorentzVector const & h1) {
    return h1.minus();
  }
 
  inline
  auto perphat(CLHEP::HepLorentzVector const & h1) {
    const auto perp = std::complex<double>{h1.x(), h1.y()};
    return perp/std::abs(perp);
  }
 
 
  std::pair<CLHEP::HepLorentzVector, CLHEP::HepLorentzVector>
  split_into_lightlike(CLHEP::HepLorentzVector const & P);
 
  inline
  CLHEP::HepLorentzVector to_HepLorentzVector(fastjet::PseudoJet const & mom){
    return {mom.px(), mom.py(), mom.pz(), mom.e()};
  }
 
  inline
  CLHEP::HepLorentzVector to_HepLorentzVector(Particle const & particle){
    return to_HepLorentzVector(particle.p);
  }
 
  inline
  fastjet::PseudoJet to_PseudoJet(CLHEP::HepLorentzVector const & mom){
    return {mom.px(), mom.py(), mom.pz(), mom.e()};
  }
} // namespace HEJ