Tabular/TankResponse.h

Go to the documentation of this file.

#ifndef _TabularTankResponseNS_TabularTankResponse_h_
#define _TabularTankResponseNS_TabularTankResponse_h_

#include <tls/VTankResponse.h>
#include <vector>
#include <algorithm>
#include <boost/multi_array.hpp>
#include <utl/NormalDistribution.h>

namespace utl {
  class Branch;
}

namespace TabularTankResponseNS {

  class TankResponse: public tls::VTankResponse {
  public:

    static TankResponse& GetInstance(const utl::Branch branch);

    virtual
    double
    PDF(const double signal,
        const double theta,
        const double r,
        const ulong muons) const;

    virtual
    double
    CDF(const double threshold,
        const double theta,
        const double r,
        const ulong muons) const;

    virtual
    double
    Mean(const double theta,
         const double r,
         const ulong muons) const;

    virtual
    double
    StDev(const double theta,
          const double r,
          const ulong muons) const;

    virtual double GetThetaMin() const { return fX; }
    virtual double GetThetaMax() const { return fX+fDX; }

    virtual double GetRadiusMin() const { return pow(10,fY); }
    virtual double GetRadiusMax() const { return pow(10,fY+fDY); }

  private:
    TankResponse() {}
    TankResponse(const utl::Branch branch);

    virtual ~TankResponse() {}
    TankResponse& operator=(const TankResponse &);

  private:

    inline
    double
    EvaluateCdf(const std::vector<double>& v, const double signal, const double mean, const double sigma)
      const
    {
      if (signal <= 0)
        return 0.0;
      if (signal >= fDS*(v.size()-2))
        return utl::NormalCDF(signal, mean, sigma);
      double z = signal/fDS;
      const size_t i = static_cast<size_t>(z);
      z -= i;
      return (1.0-z)*v[i] + z*v[i+1];
    }

    size_t fNX;
    double fX, fDX;
    size_t fNY;
    double fY, fDY;
    double fDS;
    boost::multi_array<double, 2> fMean;
    boost::multi_array<double, 2> fSigma;
    boost::multi_array<std::vector<double>, 3> fCdf;
    double fSmooth;
  };

} // NS TabularTankResponseNS

#endif