Tools/InclinedShowers/TankResponse/Aachen/TankResponse.cc

Go to the documentation of this file.

#include "TankResponse.h"

#include <utl/AugerException.h>
#include <utl/Reader.h>
#include <utl/Branch.h>
#include <utl/AugerUnits.h>
#include <utl/NormalDistribution.h>

#include <cmath>
#include <string>
#include <sstream>

using namespace AachenTankResponseNS;
using namespace utl;
using namespace std;
using namespace tls;

TankResponse& TankResponse::GetInstance(const Branch topBranch)
{
  static TankResponse singleton(topBranch);
  return singleton;
}

TankResponse::TankResponse(const Branch branch)
{
  string dataFilePath;
  branch.GetChild("DataFile").GetData(dataFilePath);

  fTable.OpenFile(dataFilePath);
}

double
TankResponse::Mean(const double theta, const double /*r*/, const ulong nmuon)
  const
{
  if (nmuon == 0) return 0;
  if (nmuon < ulong(fTable.GetTableLimit()))
    return fTable.Moment(theta,nmuon,1);
  return fTable.GetFittedMean(theta,nmuon);
}

double
TankResponse::StDev(const double theta, const double /*r*/, const ulong nmuon)
  const
{
  if (nmuon == 0) return 0;
  if (nmuon < ulong(fTable.GetTableLimit())) {
    const double m = fTable.Moment(theta,nmuon,1);
    return sqrt(fTable.Moment(theta,nmuon,2)-m*m);
  }
  return fTable.GetFittedSigma(theta,nmuon);
}

double
TankResponse::PDF(const double signal, const double theta, const double /*r*/, const ulong nmuon)
  const
{
  if (nmuon == 0) return signal == 0;
  if (nmuon < (unsigned long)(fTable.GetTableLimit()))
    return fTable.PDF(signal,theta,nmuon);
  else {
    const double m = fTable.GetFittedMean(theta,nmuon);
    const double s = fTable.GetFittedSigma(theta,nmuon);
    return NormalPDF(signal,m,s);
  }
}

double
TankResponse::CDF(const double signalmin, const double theta, const double /*r*/, const ulong nmuon)
  const
{
  if (nmuon == 0) return 0;
  if (nmuon < (unsigned long)(fTable.GetTableLimit()))
    return fTable.CDF(signalmin,theta,nmuon);
  else {
    const double m = fTable.GetFittedMean(theta,nmuon);
    const double s = fTable.GetFittedSigma(theta,nmuon);
    return NormalCDF(signalmin,m,s);
  }
}