SIOValidationObj.cc

Go to the documentation of this file.

#include "SIOValidationObj.h"
#include "TSIOValidationObj.h"

#include <utl/ErrorLogger.h>

#include <det/Detector.h>

#include <utl/TimeStamp.h>

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

#include <sevt/SEvent.h>
#include <sevt/Header.h>
#include <sevt/Station.h>
#include <sevt/PMT.h>
#include <sevt/PMTRecData.h>
#include <sevt/StationRecData.h>
#include <sevt/SortCriteria.h>
#include <fwk/LocalCoordinateSystem.h>

#include <utl/Trace.h>
#include <utl/AugerUnits.h>
#include <utl/MathConstants.h>
#include <utl/PhysicalConstants.h>
#include <utl/TabulatedFunctionErrors.h>
#include <utl/UTMPoint.h>
#include <fwk/CentralConfig.h>

using namespace det;
using namespace evt;
using namespace sevt;
using namespace std;
using namespace fwk;


SIOValidationObj::SIOValidationObj() :
  ValidationObj("SIOValidationObj")
{
  _myTObject = new TSIOValidationObj();
  _myTObjectName = _myTObject->GetName();
}


bool 
SIOValidationObj::Fill(const evt::Event& event)
{
  TSIOValidationObj* eSIO =
    dynamic_cast<TSIOValidationObj*>(_myTObject);
  if (!eSIO) {
    cout << "Fault in casting TSIOValidationObj!!!!" << endl;
    return false;
  }
  if (event.HasSEvent() && event.HasRecShower()) {
        
    const sevt::SEvent& sevent = event.GetSEvent();
    const ShowerSRecData& showersrec = event.GetRecShower().GetSRecShower();
    const utl::Vector& axis = showersrec.GetAxis();
    const utl::CoordinateSystemPtr localCS =
      LocalCoordinateSystem::Create(showersrec.GetCorePosition()); 

    eSIO->SetEventId(sevent.GetHeader().GetId());
    eSIO->SetTheta(axis.GetTheta(localCS)/utl::degree);
    eSIO->SetPhi(axis.GetPhi(localCS)/utl::degree);
    eSIO->SetShowerSize(showersrec.GetShowerSize());

    cout << "SIOValidationObj filled!" << endl;
    return true;

  }     
  return false;
}