Aachen/TankResponse.h

Go to the documentation of this file.

#ifndef _sdet_AachenTankResponse_h_
#define _sdet_AachenTankResponse_h_

#include <tls/VTankResponse.h>

#include <utl/AugerUnits.h>

#include "AachenTable.h"

namespace utl {
  class Branch;
}

namespace AachenTankResponseNS {

 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; }
    virtual double GetRadiusMax() const { return 5; }

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

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

  private:
    double fMuonEnergy;
    
    mutable AachenTable fTable;
  };
}

#endif