TelescopeData.h

Go to the documentation of this file.

#ifndef _TelescopeData_h_
#define _TelescopeData_h_

#include <vector>
#include <map>
#include <fevt/FdConstants.h>
#include <utl/Point.h>
#include <utl/TabulatedFunction.h>
#include "ZetaPixel.h"

namespace FdEnergyDepositFinderKG {

  //--------- variables of a time-bin relevant for CFM
  class TelescopeDataBin {

  public:
    TelescopeDataBin(const double minDepth, const double maxDepth,
                     const double hXmean,
                     const utl::Point& pFirst, const utl::Point& pLast,
                     const utl::Point& pXMean,
                     const double t, const double dt,
                     const double sig, const double sigErr, const double vbg,
                     const std::map<fevt::FdConstants::LightSource, double>&
                     lightEfficiencies, const int parentBin);
    ~TelescopeDataBin();

    double fMinDepth;
    double fMaxDepth;
    double fHeight;
    utl::Point fFirstPoint;
    utl::Point fLastPoint;
    utl::Point fMeanDepthPoint;
    double fTime;
    double fTimeBinHalfWidth;
    double fSignal;
    double fSignalErr;
    double fBackgroundVariance;
    int fParentBin;

    void AddZetaPixel(const ZetaPixel& pixel);

    double GetMaxZeta() const;
    const std::vector<ZetaPixel>& GetZetaPixels() const;
    const utl::TabulatedFunction& GetArcLengthFunction() const;

    double GetMeanDepth() const;

    bool HasEfficiency(fevt::FdConstants::LightSource source) const;
    double GetEfficiency(fevt::FdConstants::LightSource source) const;

  private:
    std::map<fevt::FdConstants::LightSource, double> fLightEfficiencies;
    std::vector<ZetaPixel> fZetaPixels;
    mutable utl::TabulatedFunction fArcLengthFunction;
    double fMaxZeta;
  };

  inline bool operator<(const TelescopeDataBin& lhs,
                        const TelescopeDataBin& rhs)
  {
    if (lhs.fMinDepth == 0. || rhs.fMinDepth == 0.) {
      if (lhs.fTime  < rhs.fTime)
        return true;
      else
        return false;
    } else {
      if (lhs.fMinDepth  < rhs.fMinDepth)
        return true;
      else
        return false;
    }
  }


  //--------- collection of telescope data relevant for CFM
  class TelescopeData {

  public:

    enum ETelDataType {
      eAll = 0,
      eOutsideFOV,
      eInsideFOV
    };


  public:

    TelescopeData(int telId, ETelDataType type = eInsideFOV);

    void AddTelescopeDataBin(const TelescopeDataBin& telDataBin);

    unsigned int GetId() const;
    ETelDataType GetType() const;
    double GetMinDepth() const;
    double GetMaxDepth() const;

    double GetPhotonToPhotoElectron() const
    { return fPhotonToPhotoElectron; }
    double GetGainVariance() const
    { return fGainVariance; }
    double GetDiaphragmArea() const
    { return fDiaphragmArea; }

    double GetZeta() const
    { return fZeta; }

    const std::vector<TelescopeDataBin>& GetTelescopeDataBins() const;
    std::vector<TelescopeDataBin>::const_iterator TelDataBinsBegin() const;
    std::vector<TelescopeDataBin>::const_iterator TelDataBinsEnd() const;

    void SortBins();
    bool DepthInRange(const double depth) const;

    void SetTelescopeParameters(double peFactor, double gainVariance,
                                double diaArea);
    void SetZeta(double val)
    { fZeta=val; }


  private:
    TelescopeData();
    int fTelescopeId;
    ETelDataType fTelDataType;
    std::vector<TelescopeDataBin> fTelescopeDataBins;
    double fPhotonToPhotoElectron;
    double fGainVariance;
    double fDiaphragmArea;
    double fZeta;
  };


  // --- for sorting

  inline bool operator<(const TelescopeData& lhs, const TelescopeData& rhs)
  {
    if (lhs.GetMinDepth() < rhs.GetMinDepth())
      return true;
    else
      return false;
  }

  // --- predicate for iterating
  class TelescopeFilter {

  public:


  public:
    TelescopeFilter(const TelescopeData::ETelDataType type =
                    TelescopeData::eAll) :
      fFilterCriterion(type) {}

    bool operator()(const TelescopeData& tel) const
    {
      if (fFilterCriterion == TelescopeData::eAll)
        return true;
      else
        return tel.GetType() == fFilterCriterion;
    }

  private:
    TelescopeData::ETelDataType fFilterCriterion;

  };
}

#endif