GammaGammaKinematics.hh

// -*- C++ -*-
#ifndef RIVET_GammaGammaKinematics_HH
#define RIVET_GammaGammaKinematics_HH
 
#include "Rivet/Event.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Projection.hh"
#include "Rivet/Projections/Beam.hh"
#include "Rivet/Projections/GammaGammaLeptons.hh"
 
namespace Rivet {
 

  class GammaGammaKinematics : public Projection {
  public:

    GammaGammaKinematics(const GammaGammaLeptons& lepton = GammaGammaLeptons())
        : _theQ2(make_pair(-1.0, -1.0)), _theW2(-1.0) {
      setName("GammaGammaKinematics");
      //addPdgIdPair(ANY, hadid);
      declare(Beam(), "Beam");
      declare(lepton, "Lepton");
    }

    RIVET_DEFAULT_PROJ_CLONE(GammaGammaKinematics);
 

    using Projection::operator=;
 
 
  protected:

    virtual void project(const Event& e);

    virtual CmpState compare(const Projection& p) const;
 
 
  public:

    pair<double, double> Q2() const {
      return _theQ2;
    }

    double W2() const {
      return _theW2;
    }

    const ParticlePair& beamLeptons() const {
      return _inLepton;
    }

    const ParticlePair& scatteredLeptons() const {
      return _outLepton;
    }
 
 
  protected:

    pair<double, double> _theQ2;

    double _theW2;

    ParticlePair _inLepton, _outLepton;
  };
 
 
}
 
#endif