PixelToApertureLightConverter.cc

Go to the documentation of this file.

#include "PixelToApertureLightConverter.h"

#include <string>

#include <evt/Event.h>
#include <fevt/Eye.h>
#include <fevt/Telescope.h>
#include <fevt/TelescopeSimData.h>
#include <fevt/FdComponentSelector.h>

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



using namespace FdLightCollectionEfficiencyKG;
using namespace std;


/*************************************************************************/
PixelToApertureLightConverter::PixelToApertureLightConverter(const fdet::FDetector& theFDet,
                                                             const fevt::Eye& eye)
  : fReferenceWl(theFDet.GetReferenceLambda()),
    fDetEye(theFDet.GetEye(eye.GetId())),
    fEye(eye),
    fCurTelId(0),
    fCurConfigSignature(NULL), fCurTelescope(NULL), fCurPixConverter(NULL)
{
}


/*************************************************************************/
double
PixelToApertureLightConverter::GetConversionConstant(const unsigned int telId,
                                                     const unsigned int pixId)
{
  if (telId != fCurTelId) { // currently at the wrong id, switch telescope
    fCurTelId           = telId;
    fCurTelescope       = &fDetEye.GetTelescope(telId);
    fCurPixConverter    = &fConverterMap[telId];
    const fevt::Telescope& tel = fEye.GetTelescope(telId, fevt::ComponentSelector::eExists);
    fCurConfigSignature = &tel.GetSimData().GetConfigSignature();
  }
  
  // Return calculated conversion factor if known
  PixelConversionConstantMap::const_iterator pixConvIter = fCurPixConverter->find(pixId);
  if (pixConvIter != fCurPixConverter->end())
    return pixConvIter->second;

  // Calculate and cache new conversion factor
  const fdet::Pixel& detPixel = fCurTelescope->GetPixel(pixId);
  const double pixQEff = detPixel.GetQEfficiency().Y(fReferenceWl);
  const double elecGain = fCurTelescope->GetChannel(detPixel).GetElectronicsGain();
  const double simCalibConst = detPixel.GetSimulatedEndToEndCalibration(*fCurConfigSignature);
  const double pixConversion = pixQEff * elecGain * simCalibConst;
  (*fCurPixConverter)[pixId] = pixConversion;

  return pixConversion;
}

// Configure (x)emacs for this file ...
// Local Variables:
// mode: c++
// End: