USCInter/TankResponse.h

Go to the documentation of this file.

#ifndef _USCInterTankResponseNS_USCInterTankResponse_h_
#define _USCInterTankResponseNS_USCInterTankResponse_h_

#include <tls/VTankResponse.h>

#include <utl/SplineInterpolator.h>
#include <utl/AugerUnits.h>


namespace utl {
  class Branch;
}

namespace USCInterTankResponseNS {

  typedef utl::Spline::VectorInterpolator2D USCTankModel;
  typedef utl::Spline::Interpolator2D USCAvgTankModel;
  typedef utl::Spline::Interpolator2D USCSigmaTankModel;


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

    virtual double GetRadiusMin() const { return pow(10, GetLogrMin())*utl::meter; }
    virtual double GetRadiusMax() const { return pow(10, GetLogrMax())*utl::meter; }

  private:
    TankResponse() = default;
    TankResponse(const utl::Branch branch);

    virtual ~TankResponse() = default;
    TankResponse& operator=(const TankResponse &);

  private:
    unsigned long GetMaxMuonNumber() const { return 8; }

    double GetThetaMinDeg() const { return 60; }

    double GetThetaMaxDeg() const { return 88; }

    double GetLogrMin() const { return 2; }

    double GetLogrMax() const { return 3.75; }

    static const int NBIN = 2860;

    bool fSimulationMode = true;
    double fDS = 0;
    double fTriggerParameters[2][2] = { { 0 } };
    double fProbNoTriggerSilent[2][16][30] = { { { 0 } } };
    USCTankModel fUSCTankModel[8]; // 8 is the number of muons!
    USCAvgTankModel fUSCAvgTankModel;
    USCSigmaTankModel fUSCSigmaTankModel;

  };

}


#endif