Simple/TankResponse.h

Go to the documentation of this file.

#ifndef _SimpleTankResponseNS_TankResponse_h_
#define _SimpleTankResponseNS_TankResponse_h_

#include <tls/VTankResponse.h>
#include <utl/AugerUnits.h>
#include <vector>
#include <limits>

namespace utl {
  class Branch;
}

namespace SimpleTankResponseNS {

  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 60*utl::degree; }
    virtual double GetThetaMax() const { return 88*utl::degree; }

    virtual double GetRadiusMin() const { return 0.0; }
    virtual double GetRadiusMax() const { return std::numeric_limits<double>::max(); }

  private:
    // Singleton pattern: Keep constructor private
    TankResponse(const utl::Branch branch);
    virtual ~TankResponse() {}

    const double fVerticalTankTrackLength;
    
    std::vector<double> fParMean;
    std::vector<double> fParSigma;

    void CalculateAverageAndSigma(double& mean, double& sigma,
                                  const double theta, const ulong muons) const;
  };

} // SimpleTankResponseNS

#endif