FIOValidationObj.cc

Go to the documentation of this file.

#include "FIOValidationObj.h"
#include "TFIOValidationObj.h"

#include <utl/ErrorLogger.h>

#include <det/Detector.h>

#include <fdet/Pixel.h>
#include <fdet/Telescope.h>
#include <fdet/FDetector.h>
#include <fdet/Channel.h>
#include <fdet/Eye.h>

#include <utl/TimeStamp.h>
#include <utl/TabulatedFunctionErrors.h>

#include <evt/Event.h>
#include <evt/ShowerRecData.h>
#include <evt/ShowerFRecData.h>
#include <evt/VGaisserHillasParameter.h>

#include <fevt/FEvent.h>
#include <fevt/Header.h>
#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/EyeHeader.h>
#include <fevt/EyeRecData.h>
#include <fevt/Pixel.h>
#include <fevt/PixelRecData.h>


#include <utl/Trace.h>
#include <utl/AugerUnits.h>
#include <utl/MathConstants.h>
#include <utl/PhysicalConstants.h>

#include <fwk/LocalCoordinateSystem.h>

using namespace det;
using namespace evt;
using namespace utl;
using namespace fevt;
using namespace std;
using namespace fwk;


FIOValidationObj::FIOValidationObj() :
  ValidationObj("FIOValidationObj")
{
  _myTObject = new TFIOValidationObj();
  _myTObjectName = _myTObject->GetName();
}


bool 
FIOValidationObj::Fill(const evt::Event& event)
{
  TFIOValidationObj* eFIO = 
    dynamic_cast<TFIOValidationObj*>(_myTObject);
  if (!eFIO) {
    cout << "Fault in casting TFIOValidationObj!!!!" << endl;
    return false;
  }
  VEyeRecObj eyearray[4] = { VEyeRecObj(), VEyeRecObj(), VEyeRecObj(), VEyeRecObj() };
  int EventNum = 0;
  if (event.HasFEvent()) {
    const fevt::FEvent& fevent = event.GetFEvent();
    // iterate over eyes
    for (FEvent::ConstEyeIterator eyeiter = fevent.EyesBegin(); 
         eyeiter != fevent.EyesEnd(); ++eyeiter) {
      const fevt::Eye& eye = *eyeiter; // dereferencing the eye iterator
      // gives you a reference to the eye
      const fevt::EyeRecData& eyerec = eye.GetRecData();
      const int neye = eye.GetId();
      const ShowerFRecData& frs = eyerec.GetFRecShower();
      EventNum = eye.GetHeader().GetEventNumber();

      const Vector sdp = eyerec.GetSDP(); // get Shower Detector Plane
      const Vector frsa = frs.GetAxis(); // get Shower Reconstructed Axis

      const Point eyeposition = 
        det::Detector::GetInstance().GetFDetector().GetEye(eye).GetPosition();

      const CoordinateSystemPtr eye2CS = fwk::LocalCoordinateSystem::Create(eyeposition);

      const Double_t theta = sdp.GetTheta(eye2CS) / degree; 
      const Double_t phi = sdp.GetPhi(eye2CS) / degree;
      const Double_t theta_sra = frsa.GetTheta(eye2CS) / degree; 
      const Double_t phi_sra = frsa.GetPhi(eye2CS) / degree; 
      const Double_t e_em = frs.GetEmEnergy();
      const Double_t e_tot = frs.GetTotalEnergy();
      const Double_t xmax = frs.HasGHParameters() ?
        frs.GetGHParameters().GetXMax()/(g/cm2) : 0;

      const int index = neye - 1;
      eyearray[index] = VEyeRecObj(neye, theta, phi, theta_sra, phi_sra, xmax, e_em, e_tot);
    }
    eFIO->SetEventId(EventNum);
    eFIO->SetEyeRec(eyearray);

    cout << "FIOValidationObj input end!" << endl;
    return true;
  }

  return false;
}