Beam.hh

// -*- C++ -*-
#ifndef RIVET_Beam_HH
#define RIVET_Beam_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Math/LorentzTrans.hh"
#include "Rivet/Tools/Beams.hh"
 
namespace Rivet {
 

  class Beam : public Projection {
  public:

    using Projection::operator=;

    Beam() {
      setName("Beam");
    }

    RIVET_DEFAULT_PROJ_CLONE(Beam);
 


    const ParticlePair& beams() const {
      return _theBeams;
    }

    double sqrtS() const {
      return Rivet::sqrtS(beams());
    }

    FourMomentum cmsBoostVec() const {
      return Rivet::cmsBoostVec(beams());
    }

    LorentzTransform cmsTransform() const {
      return Rivet::cmsTransform(beams());
    }

    Vector3 cmsBetaVec() const {
      return Rivet::cmsBetaVec(beams());
    }

    Vector3 cmsGammaVec() const {
      return Rivet::cmsGammaVec(beams());
    }

 


    double asqrtS() const {
      return Rivet::asqrtS(beams());
    }

    Vector3 acmsBetaVec() const {
      return Rivet::acmsBetaVec(beams());
    }

    Vector3 acmsGammaVec() const {
      return Rivet::acmsGammaVec(beams());
    }

    LorentzTransform acmsTransform() const {
      return Rivet::acmsTransform(beams());
    }

 

    FourVector pv() const;
 

    virtual void project(const Event& e);
 
 
  protected:

    ParticlePair _theBeams;
 
 
  protected:

    virtual CmpState compare(const Projection&) const {
      return CmpState::EQ;
    }
  };
 
 
}
 
#endif