DISLepton.hh

// -*- C++ -*-
#ifndef RIVET_DISLepton_HH
#define RIVET_DISLepton_HH

#include "Rivet/Projections/Beam.hh"
#include "Rivet/Projections/PromptFinalState.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"

namespace Rivet {


  class DISLepton : public Projection {
  public:



    DISLepton(){
      setName("DISLepton");
      addProjection(Beam(), "Beam");
      addProjection(PromptFinalState(), "FS");
    }

    DEFAULT_RIVET_PROJ_CLONE(DISLepton);



  protected:

    virtual void project(const Event& e);

    virtual int compare(const Projection& p) const;


  public:

    const Particle& in() const { return _incoming; }

    const Particle& out() const { return _outgoing; }

    int pzSign() const { return sign(_incoming.pz()); }


  private:

    Particle _incoming;

    Particle _outgoing;

    // /// The charge sign of the DIS current
    // double _charge;

  };

}


#endif