AachenTable.h

Go to the documentation of this file.

#ifndef _sdet_AachenTable_h_
#define _sdet_AachenTable_h_

#include <vector>
#include <string>
#include <algorithm>
#include <cmath>

#include "TabulatedPDF.h"

class AachenTable {

public:
  AachenTable();
  AachenTable(const std::string& fileName);

  ~AachenTable();

public:

  void OpenFile(const std::string& fileName);

  double PDF(double signal, double theta, int nmu) const;

  double CDF(double signalmin, double theta, int nmu) const;

  double Moment(double theta, int nmu, int order) const;

  double GetFittedMean(double theta, int nmu) const;

  double GetFittedSigma(double theta, int nmu) const;

  inline
  int GetTableLimit(){return 10;}

  inline
  int GetThetaMinDeg(){return 60;} // warning, in degree!

  inline
  int GetThetaMaxDeg(){return 88;} // warning, in degree!

private:

  class TableItem {
  public :
    double theta;
    double mean[9];
    double parMean;
    double parSigma[2];
  private :
    tls::TabulatedPDF tpf[9];
  public :
    operator double() const { return theta; } // TableItem can be casted to double, this makes sorting possible

    inline
    tls::TabulatedPDF& GetTPF(int nmu){ return tpf[nmu-1]; }
    inline
    const tls::TabulatedPDF& GetTPF(int nmu) const { return tpf[nmu-1]; }

    inline
    double& GetMean(int nmu){ return mean[nmu-1]; }
    inline
    const double& GetMean(int nmu) const { return mean[nmu-1]; }

    inline
      double GetFittedMean(int nmu) const { return parMean*nmu; }
    inline
      double GetFittedSigma(int nmu) const { return std::sqrt(parSigma[0]+parSigma[1]*nmu); }
  };

  typedef std::vector< TableItem > InternalTable;
  InternalTable fTable;

  inline
  const InternalTable::const_iterator UpperLimit(double theta) const
  { return std::upper_bound(fTable.begin(),fTable.end(),theta); }
};

#endif