Beams.hh

// -*- C++ -*-
#ifndef RIVET_Beams_HH
#define RIVET_Beams_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Math/LorentzTrans.hh"
 
namespace Rivet {
 


  ParticlePair beams(const Event& e);
 

  bool validBeams(const ParticlePair& beams);

  inline bool validBeams(const Event& e) {
    return validBeams(beams(e));
  }
 

  inline PdgIdPair beamIDs(const ParticlePair& beams) {
    return pids(beams);
  }

  inline PdgIdPair beamIDs(const Event& e) {
    return beamIDs(beams(e));
  }
 

  inline PdgIdPair beamEnergies(const ParticlePair& beams) {
    return energies(beams);
  }

  inline PdgIdPair beamEnergies(const Event& e) {
    return beamEnergies(beams(e));
  }
 

  double sqrtS(double ea, double eb);

  inline double sqrtS(const pair<double, double>& energies) {
    return sqrtS(energies.first, energies.second);
  }

  double sqrtS(const FourMomentum& pa, const FourMomentum& pb);

  inline double sqrtS(const ParticlePair& beams) {
    return sqrtS(beams.first.momentum(), beams.second.momentum());
  }

  inline double sqrtS(const Event& e) {
    return sqrtS(beams(e));
  }
 

  double asqrtS(const FourMomentum& pa, const FourMomentum& pb);

  double asqrtS(const ParticlePair& beams);

  inline double asqrtS(const Event& e) {
    return asqrtS(beams(e));
  }
 

  inline FourMomentum cmsBoostVec(const FourMomentum& pa, const FourMomentum& pb) {
    return pa + pb;
  }

  inline FourMomentum cmsBoostVec(const ParticlePair& beams) {
    return cmsBoostVec(beams.first, beams.second);
  }

  FourMomentum acmsBoostVec(const FourMomentum& pa, const FourMomentum& pb);

  FourMomentum acmsBoostVec(const ParticlePair& beams);
 

  Vector3 cmsBetaVec(const FourMomentum& pa, const FourMomentum& pb);

  inline Vector3 cmsBetaVec(const ParticlePair& beams) {
    return cmsBetaVec(beams.first, beams.second);
  }
 

  Vector3 acmsBetaVec(const FourMomentum& pa, const FourMomentum& pb);

  Vector3 acmsBetaVec(const ParticlePair& beams);
 

  Vector3 cmsGammaVec(const FourMomentum& pa, const FourMomentum& pb);

  inline Vector3 cmsGammaVec(const ParticlePair& beams) {
    return cmsGammaVec(beams.first, beams.second);
  }
 

  Vector3 acmsGammaVec(const FourMomentum& pa, const FourMomentum& pb);

  Vector3 acmsGammaVec(const ParticlePair& beams);
 

  LorentzTransform cmsTransform(const FourMomentum& pa, const FourMomentum& pb);

  inline LorentzTransform cmsTransform(const ParticlePair& beams) {
    return cmsTransform(beams.first, beams.second);
  }
 

  LorentzTransform acmsTransform(const FourMomentum& pa, const FourMomentum& pb);

  LorentzTransform acmsTransform(const ParticlePair& beams);

 
 
}
 
#endif