TSIOValidationObj.cc

Go to the documentation of this file.

#include <cmath>
#include <iostream>
#include <iomanip>
#include <utl/AugerUnits.h>
#include <TVector3.h>
#include "ValidationUtil.h"
#include "TSIOValidationObj.h"

using namespace std;


ClassImp(TSIOValidationObj);


TSIOValidationObj::TSIOValidationObj() :
  TValidationObj(),
  VEventId(0),
  VTheta(-999.),
  VPhi(-999.),
  VS1000(-999.)
{
}


void
TSIOValidationObj::SetRecord(const Int_t nev, const Double_t th, const Double_t ph, const Double_t s1000)
{
  VEventId = nev;
  VTheta = th;
  VPhi = ph;
  VS1000 = s1000;
}


bool
TSIOValidationObj::Validate(const TValidationObj* const vo)
{
  using namespace utl;

  const TSIOValidationObj* const o = dynamic_cast<const TSIOValidationObj*>(vo);
  if (!o) {
    cerr << "Fault in casting TSIOValidationObj!!!!" << endl;
    return false;
  }

  if (o->VEventId != VEventId)
    return false;

  if (o->VS1000 <= 0)
    return false;
  if (RelativeError(o->VS1000, VS1000) > 5e-3)
    return false;

  const double angle =
    Angle(o->VTheta*degree, o->VPhi*degree, VTheta*degree, VPhi*degree);
  const double angleTolerance = 1e-2*degree;
  if (angle > angleTolerance)
    return false;

  return true;
}


void
TSIOValidationObj::Dump()
{
  cout << "Event ID " << VEventId << '\n'
       << "Theta  "<< setprecision(16) << VTheta << '\n'
       << "Phi " << setprecision(16) << VPhi << '\n'
       << "S(1000) " << setprecision(16) << VS1000 << endl;
}