DISRapidityGap.hh

// -*- C++ -*-
#ifndef RIVET_DISRapidityGap_HH
#define RIVET_DISRapidityGap_HH
 
#include "Rivet/Projections/DISKinematics.hh"
#include "Rivet/Projections/DISFinalState.hh"
#include "Rivet/Particle.hh"
#include "Rivet/Event.hh"
 
namespace Rivet {
 
  
  class DISRapidityGap : public Projection {
 
  public:
 
    DISRapidityGap() {
      setName("DISRapidityGap");
      declare(DISKinematics(), "DISKIN");
      declare(DISFinalState(DISFrame::HCM), "DISFS");
    }
 
    RIVET_DEFAULT_PROJ_CLONE(DISRapidityGap);
 
    using Projection::operator =;
 
    double t() const { return _t; }
 
    double gap() const { return _gap; }
 
    double gapUpp() const { return _gapUpp; }
 
    double gapLow() const { return _gapLow; }
 
    double EpPzX(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _ePpzX_LAB;
      else if (f == DISFrame::XCM) return _ePpzX_XCM;
      else return _ePpzX_HCM;
    }
 
    double emPzX(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _eMpzX_LAB;
      else if (f == DISFrame::XCM) return _eMpzX_XCM;
      else return _eMpzX_HCM;
    }
 
 
    const Particles& systemX(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _pX_LAB;
      else if (f == DISFrame::XCM) return _pX_XCM;
      else return _pX_HCM;
    }
 
    const Particles& systemY(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _pY_LAB;
      else if (f == DISFrame::XCM) return _pY_XCM;
      else return _pY_HCM;
    }
 
    const FourMomentum pX(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _momX_LAB;
      else if (f == DISFrame::XCM) return _momX_XCM;
      else return _momX_HCM;
    }
 
    const FourMomentum pY(const DISFrame& f) const {
      if (f == DISFrame::LAB) return _momY_LAB;
      else if (f == DISFrame::XCM) return _momY_XCM;
      else return _momY_HCM;
    }
 
    double m2X() const { return _M2X; }
 
    double m2Y() const { return _M2Y; }
 
 
  protected:
 
    virtual CmpState compare(const Projection& p) const;
 
    virtual void project(const Event& e);
 
    void clear();
 
    void findGap(const Particles& particles, const DISKinematics& diskin);
 
 
  protected:
 
    double _M2X, _M2Y, _t;
    double _gap, _gapUpp, _gapLow;
    double _ePpzX_LAB, _eMpzX_LAB;
    double _ePpzX_HCM, _eMpzX_HCM;
    double _ePpzX_XCM, _eMpzX_XCM;
    FourMomentum _momX_HCM, _momY_HCM;
    FourMomentum _momX_LAB, _momY_LAB;
    FourMomentum _momX_XCM, _momY_XCM;
    Particles _pX_HCM, _pY_HCM, _pX_LAB, _pY_LAB, _pX_XCM, _pY_XCM;
 
  };
 
  
}
 
#endif