VEyeRecObj.cc

Go to the documentation of this file.

#include <cmath>
#include <iostream>
#include <sstream>
#include <utl/AugerUnits.h>
#include "ValidationUtil.h"
#include "VEyeRecObj.h"

using namespace std;
using namespace utl;


ClassImp(VEyeRecObj);


VEyeRecObj::VEyeRecObj() :
  VEye_Id(0),
  VThetaSDP(-999.),
  VPhiSDP(-999.),
  VThetaAxis(-999.),
  VPhiAxis(-999.),
  VXMax(-999.),
  VEem(-999.),
  VEtot(-999.)
{
}


void
VEyeRecObj::SetRecord(const Int_t neye, const Double_t thsdp, const Double_t phsdp,
                      const Double_t thaxi, const Double_t phaxi, const Double_t xmax,
                      const Double_t eem, const Double_t etot)
{
  VEye_Id = neye;
  VThetaSDP = thsdp;
  VPhiSDP = phsdp;
  VThetaAxis = thaxi;
  VPhiAxis = phaxi;
  VXMax = xmax;
  VEem = eem;
  VEtot = etot;
}


void
VEyeRecObj::Dump()
{
  cout << "EyeId " << GetEyeId() << '\n'
       << "ThetaSDP " << GetThetaSDP() << '\n'
       << "PhiSDP " << GetPhiSDP() << '\n'
       << "ThetaAxis " << GetThetaAxis() << '\n'
       << "PhiAxis " << GetPhiAxis() << '\n'
       << "Xmax " << GetXMax() << '\n'
       << "Eem " <<GetEem() << '\n'
       << "Etot " << GetEtot() << endl;
}


bool
VEyeRecObj::Match(const VEyeRecObj& e)
{
  if (!GetXMax() || !GetEem() || !GetEtot())
    return false;

  if (e.GetEyeId() != GetEyeId())
    return false;

  const double angleSDP =
    Angle(e.GetThetaSDP()*degree, e.GetPhiSDP()*degree,
          GetThetaSDP()*degree, GetPhiSDP()*degree);

  if (angleSDP > 1e-2*degree)
    return false;

  const double angleAxis =
    Angle(e.GetThetaAxis()*degree, e.GetPhiAxis()*degree,
          GetThetaAxis()*degree, GetPhiAxis()*degree);
  if (angleAxis > 1e-2*degree)
    return false;

  if (RelativeError(e.GetXMax(), GetXMax()) > 5e-3)
    return false;

  if (RelativeError(e.GetEem(), GetEem()) > 5e-3)
    return false;

  if (RelativeError(e.GetEtot(), GetEtot()) > 5e-3)
    return false;

  return true;
}


template<typename T>
string
OutputEqual(const T& x, const T& y)
{
  ostringstream os;
  os << x << (x == y ? " == " : " != ") << y;
  return os.str();
}


bool
operator==(const VEyeRecObj& ex, const VEyeRecObj& ey)
{
  cout << "EyeId "     << OutputEqual(ex.GetEyeId(), ey.GetEyeId()) << '\n'
       << "ThetaSDP "  << OutputEqual(ex.GetThetaSDP(), ey.GetThetaSDP()) << '\n'
       << "PhiSDP "    << OutputEqual(ex.GetPhiSDP(), ey.GetPhiSDP()) << '\n'
       << "ThetaAxis " << OutputEqual(ex.GetThetaAxis(), ey.GetThetaAxis()) << '\n'
       << "PhiAxis "   << OutputEqual(ex.GetPhiAxis(), ey.GetPhiAxis()) << '\n'
       << "Xmax "      << OutputEqual(ex.GetXMax(), ey.GetXMax()) << '\n'
       << "Eem "       << OutputEqual(ex.GetEem(), ey.GetEem()) << '\n'
       << "Etot "      << OutputEqual(ex.GetEtot(), ey.GetEtot()) << endl;
  return (ex.GetEyeId() == ey.GetEyeId() &&
          ex.GetThetaSDP() == ey.GetThetaSDP() &&
          ex.GetPhiSDP() == ey.GetPhiSDP() &&
          ex.GetThetaAxis() == ey.GetThetaAxis() &&
          ex.GetPhiAxis() == ey.GetPhiAxis() &&
          ex.GetXMax() == ey.GetXMax() &&
          ex.GetEem() == ey.GetEem() &&
          ex.GetEtot() == ey.GetEtot());
}