FDetector.cc

Go to the documentation of this file.

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

#include <det/VManager.h>
#include <det/Detector.h>

#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/Pixel.h>
#include <fevt/Channel.h>

#include <utl/AugerException.h>
#include <utl/ErrorLogger.h>

using namespace fdet;
using namespace std;
using namespace det;
using namespace utl;


unsigned int
FDetector::GetFirstEyeId()
  const
{
  if (!fFullEyeList)
    return 0;
  return *std::min_element(fFullEyeList->begin(), fFullEyeList->end());
}


unsigned int
FDetector::GetLastEyeId()
  const
{
  if (!fFullEyeList)
    return 0;
  return *std::max_element(fFullEyeList->begin(), fFullEyeList->end());
}


FDetector::~FDetector()
{
  for (InternalEyeIterator it = fCommissionedEyeList.begin();
       it != fCommissionedEyeList.end(); ++it)
    delete *it;
  delete fFullEyeList;
  delete fReferenceLambda;
}


void
FDetector::Update()
{
  // Get all eyes indicated in the XML file (uses lazy evaluation)
  //
  GetFullEyeList();

  for (InternalEyeIterator it = fCommissionedEyeList.begin();
       it != fCommissionedEyeList.end(); ++it)
    (*it)->Update();
}


const Eye&
FDetector::GetEye(const unsigned int eyeId)
  const
{
  for (AllEyeIterator it = AllEyesBegin(); it != AllEyesEnd(); ++it)
    if (eyeId == it->GetId())
      return *it;

  ostringstream err;
  err << "Eye with Id = " << eyeId << " does not exist, "
         "or else detector description was not properly configured. "
         "detector time = " << Detector::GetInstance().GetTime();
  ERROR(err);
  throw utl::NonExistentComponentException(err.str());

}


const Eye&
FDetector::GetEye(const string& eyeName)
  const
{
  for (AllEyeIterator it = AllEyesBegin(); it != AllEyesEnd(); ++it)
    if (eyeName == it->GetName())
      return *it;

  ostringstream err;
  err << "Eye with name = " << eyeName << " does not exist, "
         "or else detector description was not properly configured. "
         "detector time = " << Detector::GetInstance().GetTime();
  ERROR(err);
  throw utl::NonExistentComponentException(err.str());
}


const FDetector::InternalEyeIdCollection&
FDetector::GetFullEyeList()
  const
{
  if (!fFullEyeList) {

    const VManager& manager =
      Detector::GetInstance().GetFManagerRegister();

    fFullEyeList = new InternalEyeIdCollection;

    VManager::IndexMap indexMap;
    const VManager::Status foundFlag =
      manager.GetData(*fFullEyeList, "eyeList", "", indexMap);

    if (foundFlag == VManager::eNotFound) {
      DEBUGLOG("Could not find a list of Eye Id's. "
               "No eye information will be available!");
    }

    for (InternalEyeIdCollection::iterator it = fFullEyeList->begin();
         it != fFullEyeList->end(); ++it) {
      Eye* const eye = new Eye(*it);
      fCommissionedEyeList.push_back(eye);
    }

  }
  return *fFullEyeList;
}


const Eye&
FDetector::GetEye(const fevt::Eye& eventEye)
  const
{
  return GetEye(eventEye.GetId());
}


const Telescope&
FDetector::GetTelescope(const fevt::Telescope& eventTel)
  const
{
  const unsigned int eyeid = eventTel.GetEyeId();
  const unsigned int telid = eventTel.GetId();

  return GetEye(eyeid).GetTelescope(telid);
}


const Channel&
FDetector::GetChannel(const fevt::Channel& eventChannel)
  const
{
  const unsigned int eyeid = eventChannel.GetEyeId();
  const unsigned int telid = eventChannel.GetTelescopeId();
  const unsigned int chid  = eventChannel.GetId();

  return GetEye(eyeid).GetTelescope(telid).GetChannel(chid);
}


const Channel&
FDetector::GetChannel(const fdet::Pixel& p)
  const
{
  const unsigned int eyeid = p.GetEyeId();
  const unsigned int telid = p.GetTelescopeId();
  const unsigned int chid  = p.GetChannelId();

  return GetEye(eyeid).GetTelescope(telid).GetChannel(chid);
}


const Channel&
FDetector::GetChannel(const fevt::Pixel& p)
  const
{
  const unsigned int eyeid = p.GetEyeId();
  const unsigned int telid = p.GetTelescopeId();
  const unsigned int chid =
    det::Detector::GetInstance().GetFDetector().GetEye(eyeid).GetTelescope(telid).GetPixel(p.GetId()).GetChannelId();

  return GetEye(eyeid).GetTelescope(telid).GetChannel(chid);
}


const Pixel&
FDetector::GetPixel(const fevt::Pixel& eventPixel)
  const
{
  const unsigned int eyeid = eventPixel.GetEyeId();
  const unsigned int telid = eventPixel.GetTelescopeId();
  const unsigned int pixid = eventPixel.GetId();

  return GetEye(eyeid).GetTelescope(telid).GetPixel(pixid);
}


const Pixel&
FDetector::GetPixel(const fevt::Channel& c)
  const
{
  const unsigned int eyeid = c.GetEyeId();
  const unsigned int telid = c.GetTelescopeId();
  const unsigned int pixid =
    det::Detector::GetInstance().GetFDetector().GetEye(eyeid).GetTelescope(telid).GetChannel(c.GetId()).GetPixelId();

  return GetEye(eyeid).GetTelescope(telid).GetPixel(pixid);
}


const Pixel&
FDetector::GetPixel(const fdet::Channel& c)
  const
{
  const unsigned int eyeid = c.GetEyeId();
  const unsigned int telid = c.GetTelescopeId();
  const unsigned int pixid = c.GetPixelId();

  return GetEye(eyeid).GetTelescope(telid).GetPixel(pixid);
}


double
FDetector::GetCDASUpTimeFraction()
  const
{
  UpdateFdUpTime();
  return fUpTimeFraction;
}


bool
FDetector::IsCDASInAquisition()
  const
{
  UpdateFdUpTime();
  return fUpTimeFraction > 0;
}


int
FDetector::GetCDASDAQStatus()
  const
{
  UpdateFdUpTime();
  return fStatus;
}


bool
FDetector::HasFdUpTime()
  const
{
  UpdateFdUpTime();
  return fHasFdUpTime;
}


void
FDetector::UpdateFdUpTime()
  const
{
  if (fUpTimeValidityStamp.IsValid())
    return;

  GetFDetectorData(fHasFdUpTime, "has_fduptime", "CDAS", "set events up time bit");
  int startTime;
  GetFDetectorData(startTime, "uptime_gpsStart", "CDAS", "start validity CDAS up time");
  int endTime;
  GetFDetectorData(endTime, "uptime_gpsStop", "CDAS", "end validity CDAS up time");
  GetFDetectorData(fUpTimeFraction, "uptime_fraction", "CDAS", "CDAS up time fraction");
  //GetFDetectorData(fVetoFraction, "veto_fraction", "CDAS", "CDAS Veto fraction");
  GetFDetectorData(fStatus, "status", "CDAS", "CDAS DAQ status");
  fUpTimeValidityStamp.SetValidityInterval(startTime, endTime);
}


// returns the maximum minimal-wavelength validity of all FDet-Components
double
FDetector::GetModelMinWavelength()
  const
{
  double minWavelength = 0;
  bool first = true;
  for (InternalEyeIterator it = fCommissionedEyeList.begin();
       it != fCommissionedEyeList.end(); ++it) {
    if (first) {
      first = false;
      minWavelength = (*it)->GetModelMinWavelength();
    }
    else {
      minWavelength = max(minWavelength, (*it)->GetModelMinWavelength());
    }
  }
  return minWavelength;
}


// returns the minimum maximal-wavelength validity of all FDet-Components
double
FDetector::GetModelMaxWavelength()
  const
{
  double maxWavelength = 0;
  bool first = true;
  for (InternalEyeIterator it = fCommissionedEyeList.begin();
       it != fCommissionedEyeList.end(); ++it) {
    if (first) {
      first = false;
      maxWavelength = (*it)->GetModelMaxWavelength();
    }
    else {
      maxWavelength = min(maxWavelength, (*it)->GetModelMaxWavelength());
    }
  }
  return maxWavelength;
}


double
FDetector::GetReferenceLambda()
  const
{
  if (!fReferenceLambda)
    GetFDetectorData(fReferenceLambda, "referenceLambda", "telescopeList", "reference lambda");

  return *fReferenceLambda;
}


template<typename T>
inline
const T&
FDetector::GetFDetectorData(T*& requestedData,
                            const string& property,
                            const string& component,
                            const string& errorMsg)
  const
{
  if (!requestedData) {
    requestedData = new T;
    GetFDetectorData(*requestedData, property, component, errorMsg);
  }
  return *requestedData;
}


template<typename T>
inline
void
FDetector::GetFDetectorData(T& requestedData,
                            const string& property,
                            const string& component,
                            const string& errorMsg)
  const
{
  const VManager& manager = Detector::GetInstance().GetFManagerRegister();
  VManager::IndexMap indexMap; // empty (there is only one CDAS !)
  const VManager::Status status =
    manager.GetData(requestedData, property, component, indexMap);

  // ERROR
  if (status == VManager::eNotFound) {
    ostringstream err;
    err << "Did not find requested component: '" << errorMsg << "' for CDAS";
    ERROR(err);
    throw NonExistentComponentException(err.str());
  }
}


// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// compile-command: "make -C .. -k"
// End: