RadioSimulation.cc

Go to the documentation of this file.

#include <evt/RadioSimulation.h>
#include <utl/Vector.h>
#include <det/Detector.h>
#include <utl/LameShadowPtr_imp.h>

using namespace evt;
using namespace std;


LAMESHADOWPTR_INST(evt::RadioSimulation);


utl::CoordinateSystemPtr
RadioSimulation::GetLocalCoordinateSystem()
  const
{
  if (!fLocalCoordinateSystem.IsValid()) {
    throw utl::AugerException("Local coordinate system was not set!");
  }
  return fLocalCoordinateSystem.Get();
}


void
RadioSimulation::SetCoreCoordinates(const double x, const double y, const double z)
{
  const utl::CoordinateSystemPtr& cs = det::Detector::GetInstance().GetReferenceCoordinateSystem();
  fCorePosition = utl::Point(x, y, z, cs);
  fLocalCoordinateSystem = fwk::LocalCoordinateSystem::Create(fCorePosition.Get());
}


void
RadioSimulation::SetCorePosition(const utl::Point& core)
{
  fCorePosition = core;
  fLocalCoordinateSystem = fwk::LocalCoordinateSystem::Create(core);
}


utl::Point
RadioSimulation::GetCorePosition()
  const
{
  if (!fCorePosition.IsValid()) {
    throw utl::AugerException("Core position was not set!");
  }
  return fCorePosition.Get();
}


bool
RadioSimulation::GoToFirstSimRadioPulse()
{
  if (fSimRadioPulseVector.empty())
    return false;
  else {
    fPulseIter = fSimRadioPulseVector.begin();
    return true;
  }
}


const SimRadioPulse&
RadioSimulation::GetNextSimRadioPulse(bool& ok)
{
  // if fPulseIter has never been initialized, initialize it now
  if (!fPulseIter.IsValid())
    fPulseIter = fSimRadioPulseVector.begin();

  // if no more elements in vector, return illegal simpulse and ok=false
  if (fPulseIter.Get() == fSimRadioPulseVector.end()) {
    ok = false;
    fPulseIter.SetValid(false);
    return *fPulseIter.Get();
  }
  // otherwise increment iterator for next call and return the current element
  ok = true;
  const SimRadioPulse& srp = *fPulseIter.Get()++;

  return srp;
}


//const SimRadioPulse&
//RadioSimulation::GetSimRadioPulse() 
//{
//  const SimRadioPulse& srp = *fPulseIter.Get();
//  srp.SetLocalCoordinateSystem(GetLocalCoordinateSystem());
//  return srp;
//}


const SimRadioPulse&
RadioSimulation::GetSimPulseByIndex(const int index) 
{
  const SimRadioPulse& srp = fSimRadioPulseVector.at(index);

  return srp;
}