EventDis.cxx

Go to the documentation of this file.

#include "eicsmear/erhic/EventDis.h"

#include "eicsmear/erhic/Kinematics.h"

namespace erhic {

EventDis::~EventDis() { }

EventDis::EventDis()
: x(NAN)
, QSquared(NAN)
, y(NAN)
, WSquared(NAN)
, nu(NAN)
, yJB(NAN)
, QSquaredJB(NAN)
, xJB(NAN)
, WSquaredJB(NAN)
, yDA(NAN)
, QSquaredDA(NAN)
, xDA(NAN)
, WSquaredDA(NAN) {
}

EventDis::EventDis(const EventDis& that)
: VirtualEvent(that) {
  CopyKinematics(that);
}

EventDis& EventDis::operator=(const EventDis& that) {
  if (this != &that) {  // Protect against self-assignment
    CopyKinematics(that);
  }  // if
  return *this;
}

void EventDis::CopyKinematics(const EventDis& that) {
  SetLeptonKinematics(DisKinematics(that.x, that.y, that.nu,
                                    that.QSquared, that.WSquared));
  SetJacquetBlondelKinematics(DisKinematics(that.xJB, that.yJB, 0.,
                                            that.QSquaredJB, that.WSquaredJB));
  SetDoubleAngleKinematics(DisKinematics(that.xDA, that.yDA, 0.,
                                         that.QSquaredDA, that.WSquaredDA));
}

void EventDis::SetLeptonKinematics(const DisKinematics& kin) {
  x = kin.mX;
  QSquared = kin.mQ2;
  WSquared = kin.mW2;
  nu = kin.mNu;
  y = kin.mY;
}

void EventDis::SetJacquetBlondelKinematics(const DisKinematics& kin) {
  xJB = kin.mX;
  QSquaredJB = kin.mQ2;
  WSquaredJB = kin.mW2;
  yJB = kin.mY;
}

void EventDis::SetDoubleAngleKinematics(const DisKinematics& kin) {
  xDA = kin.mX;
  QSquaredDA = kin.mQ2;
  WSquaredDA = kin.mW2;
  yDA = kin.mY;
}

}  // namespace erhic