RdSdSimCheck/Validation.cc

Go to the documentation of this file.

#include <evt/Event.h>
#include <evt/ShowerRecData.h>
#include <evt/ShowerSRecData.h>

#include <sevt/SEvent.h>
#include <sevt/SortCriteria.h>

#include <utl/ErrorLogger.h>
#include <utl/Branch.h>
#include <utl/AugerUnits.h>

#include <tst/Validatrix.h>

#include <fwk/CentralConfig.h>
#include <fwk/LocalCoordinateSystem.h>

#include "Validation.h"

using namespace std;
using namespace evt;
using namespace sevt;
using namespace fevt;
using namespace fwk;
using namespace utl;
using namespace tst::Validatrix;


namespace ValidationNS {

  VModule::ResultFlag
  Validation::Init()
  {
    const Branch topB = CentralConfig::GetInstance()->GetTopBranch("Validation");
    topB.GetChild("RefFilename").GetData(fRefFileName);

    fResults = new ofstream(fRefFileName + ".new");

    return eSuccess;
  }


  VModule::ResultFlag
  Validation::Run(evt::Event& event)
  {
    // SD
    if (!event.HasSEvent()) {
      ERROR("Event should contain an SEvent.");
      return eFailure;
    }

    event.GetSEvent().SortStations(ByIncreasingId());
    const SEvent& sEvent = event.GetSEvent();

    const string eventId = event.GetHeader().GetId();

    *fResults << BeLabel("SD_rec");

    for (const auto& s : sEvent.StationsRange()) {
      if (s.HasRecData()) {
        ostringstream label;
        label << "E" << eventId << "_S" << s.GetId();
        *fResults << BeLabel(label.str());
        *fResults << BeCloseRel(s.GetRecData().GetTotalSignal(), /*tolerance*/0.05);
      }
    }

    *fResults << BeLabel("Shower_rec");

    // Shower-level rec data
    if (!event.HasRecShower())
      ERROR("Event does not contain a RecShower");

    const ShowerRecData& showerRec = event.GetRecShower();

    // SD rec (shower level)

    if (!showerRec.HasSRecShower())
      ERROR("RecShower does not have SRecShower");

    *fResults << BeLabel("Shower_S_rec");

    const ShowerSRecData& showerSRec = showerRec.GetSRecShower();

    const Point& corePos = showerSRec.GetCorePosition();
    const Point& baryPos = showerSRec.GetBarycenter();

    const CoordinateSystemPtr coreCS = LocalCoordinateSystem::Create(corePos);

    const CoordinateSystemPtr baryCS = LocalCoordinateSystem::Create(baryPos);

    const double tolerance = 1e-3;

    ostringstream label;

    label << eventId << "_ShowerSize";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetShowerSize(), tolerance);
    label.str("");

    label << eventId << "_ShowerSizeError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetShowerSizeError(), tolerance);
    label.str("");

    label << eventId << "_LDFChi2";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetLDFChi2(), tolerance);
    label.str("");

    label << eventId << "_LDFNdof";
    *fResults << BeLabel(label.str());
    *fResults << BeEqual(showerSRec.GetLDFNdof());
    label.str("");

    label << eventId << "_Theta";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetAxis().GetTheta(coreCS), tolerance);
    label.str("");

    label << eventId << "_Phi";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetAxis().GetPhi(coreCS), tolerance);
    label.str("");

    label << eventId << "_ThetaError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetThetaError(), tolerance);
    label.str("");

    label << eventId << "_PhiError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetPhiError(), tolerance);
    label.str("");

    label << eventId << "_CoreX";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseAbs(corePos.GetX(baryCS), tolerance);
    label.str("");

    label << eventId << "_CoreY";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseAbs(corePos.GetY(baryCS), tolerance);
    label.str("");

    label << eventId << "_CoreZ";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseAbs(corePos.GetZ(baryCS), tolerance);
    label.str("");

    label << eventId << "_CorrelationXY";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseAbs(showerSRec.GetCorrelationXY(), tolerance);
    label.str("");

    label << eventId << "_CoreXVariance";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetParameterCovariance(sevt::eShowerAxisX, sevt::eShowerAxisX), tolerance);
    label.str("");

    label << eventId << "_CoreYVariance";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetParameterCovariance(sevt::eShowerAxisY, sevt::eShowerAxisY), tolerance);
    label.str("");

    label << eventId << "_CoreXYCovariance";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetParameterCovariance(sevt::eShowerAxisX, sevt::eShowerAxisY), tolerance);
    label.str("");

    label << eventId << "_CoreTSecond";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetCoreTime().GetGPSSecond(), tolerance);
    label.str("");

    label << eventId << "_CoreTNanoSecond";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetCoreTime().GetGPSNanoSecond(), tolerance);
    label.str("");

    label << eventId << "_CoreXError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetCoreError().GetX(baryCS), tolerance);
    label.str("");

    label << eventId << "_CoreYError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetCoreError().GetY(baryCS), tolerance);
    label.str("");

    label << eventId << "_CoreTError";
    *fResults << BeLabel(label.str());
    *fResults << BeCloseRel(showerSRec.GetCoreTimeError().GetInterval(), tolerance);
    label.str("");

    return eSuccess;
  }


  VModule::ResultFlag
  Validation::Finish()
  {
    delete fResults;
    fResults = nullptr;

    const bool comp = Compare(fRefFileName, fRefFileName + ".new", /*failOnFirst=*/false);

    return comp ? eSuccess : eFailure;
  }

}