Corrector.cc

Go to the documentation of this file.

#include <cmath>
#include <sstream>

#include <boost/lexical_cast.hpp>

#include <fdet/Corrector.h>

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


#include <fdet/FDetector.h>
#include <fdet/Eye.h>
#include <fdet/Telescope.h>
#include <det/VManager.h>
#include <det/Detector.h>

#include <utl/Math.h>

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


Corrector::Corrector(const unsigned int eyeId, const unsigned int telId,
                     const std::string& physEyeIdString,
                     const std::string& physTelIdString) :
  fEyeId(eyeId),
  fTelescopeId(telId),
  fPhysicalEyeIdString(physEyeIdString),
  fPhysicalTelescopeIdString(physTelIdString),
  fInnerRadius(0),
  fOuterRadius(0),
  fLensThickness(0),
  fSigmaNormal(0),
  fTransmittance(0),
  fRefractiveIndex(0)
{
}


Corrector::~Corrector()
{
  delete fInnerRadius;
  delete fOuterRadius;
  delete fLensThickness;
  delete fSigmaNormal;
  delete fTransmittance;
  delete fRefractiveIndex;
}

void Corrector::Update() {

  //std::list<CorrectorSegment*> fSegment;

  delete fInnerRadius;
  fInnerRadius = 0;

  delete fOuterRadius;
  fOuterRadius = 0;

  delete fLensThickness;
  fLensThickness = 0;

  delete fSigmaNormal;
  fSigmaNormal = 0;
  
  delete fTransmittance;
  fTransmittance = 0;
  
  delete fRefractiveIndex;
  fRefractiveIndex = 0;
}


double
Corrector::GetInnerRadius()
  const
{
  return GetCorrectorData(fInnerRadius, "innerRadius", "corrector", "corrector ring inner radius");
}


double
Corrector::GetOuterRadius()
  const
{
  return GetCorrectorData(fOuterRadius, "outerRadius", "corrector", "corrector ring outer radius");
}


double
Corrector::GetMeanLensThickness()
  const
{
  return GetCorrectorData(fLensThickness, "lensThickness", "corrector", "corrector lens mean thickness");
}


double
Corrector::GetSigmaNormal()
  const
{
  return GetCorrectorData(fSigmaNormal, "sigmaNormal", "corrector", "corrector lens sigma normal");
}


const TabulatedFunction&
Corrector::GetTransmittance()
  const
{
  return GetCorrectorData(fTransmittance, "transmittance", "corrector", "corrector lens transmittance");
}


const TabulatedFunction&
Corrector::GetRefractiveIndex()
  const
{
  // -------------------
  // Tilo's KA prototyp:
  // return 1.476 + pow ((69.*nanometer/lambda), 2);
  // -------------------

  // Schott BK7 glass
  // input: wave length lambda of transmittance

  // data from 
  // http://www.schott.com/optics_devices/german/download/schott-11-01-2006.zip

  if (!fRefractiveIndex) {

    fRefractiveIndex = new TabulatedFunction;
    const TabulatedFunction& trans = GetTransmittance();

    const double b1 = 1.039612120E+00;
    const double c1 = 6.000698670E-03;
    const double b2 = 2.317923440E-01;
    const double c2 = 2.001791440E-02;
    const double b3 = 1.010469450E+00;
    const double c3 = 1.035606530E+02;

    // dispersion relation from 
    // http://www.schott.com/optics_devices/german/download/katalog-optisches-glas-infos-2003.pdf

    const TabulatedFunction::ConstIterator transEnd = trans.End();
    for (TabulatedFunction::ConstIterator iWl = trans.Begin();
         iWl != transEnd; ++iWl) {

      const double lambda = iWl->X() / nanometer;
      const double l2 = Sqr(lambda)*1.e-6;
      const double n = sqrt(b1*l2/(l2-c1) + b2*l2/(l2-c2) + b3*l2/(l2-c3) + 1);
      fRefractiveIndex->PushBack(lambda*nanometer, n);

    }

  }

  return *fRefractiveIndex;
}


template<typename T>
inline
const T&
Corrector::GetCorrectorData(T*& requestedData,
                            const string& property,
                            const string& component,
                            const string& errorMsg)
  const
{  
  if (!requestedData) {
    
    requestedData = new T;
    
    const VManager& manager = Detector::GetInstance().GetFManagerRegister();

    VManager::IndexMap indexMap;
    indexMap["eyeId"] = fPhysicalEyeIdString;
    indexMap["telescopeId"] = fPhysicalTelescopeIdString;

    const VManager::Status status =
      manager.GetData(*requestedData, property, component, indexMap);

    if (status == VManager::eNotFound) {

      ostringstream err;
      err << "Did not find requested component for " << errorMsg << "; "
          << VManager::QueryInfoMessage(property, component, indexMap);
      ERROR(err);
      throw NonExistentComponentException(err.str());
    }
  }
  
  return *requestedData;
}


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