NKGBetaOnlyLDF.h

Go to the documentation of this file.

#ifndef _LDFFinderKG_NKGBetaOnlyLDF_h_
#define _LDFFinderKG_NKGBetaOnlyLDF_h_


#include "NKGLDF.h"
#include <cmath>
#include <vector>
#include <utl/Math.h>
#include <utl/AugerUnits.h>


namespace LDFFinderKG {

  class NKGBetaOnlyLDF : public NKGLDF {

  public:
    NKGBetaOnlyLDF(const double refDistance) : NKGLDF(refDistance, 6) { }

    std::vector<double>
    ShapeModel(const double cosTheta, const double showerSize)
      const
    {
      const double lgSRef = std::log10(showerSize);
      const double secTheta = 1 / cosTheta;

      const double a0 = fShapeModelVector[0];
      const double a1 = fShapeModelVector[1];
      const double b0 = fShapeModelVector[2];
      const double b1 = fShapeModelVector[3];
      const double c0 = fShapeModelVector[4];
      const double c1 = fShapeModelVector[5];

      std::vector<double> shape(2);
      double& beta = shape[0];
      double& gamma = shape[1];

      beta = a0 + a1*lgSRef + secTheta*(b0 + b1*lgSRef + secTheta*(c0 + c1*lgSRef));

      gamma = 0;

      return shape;
    }

  };

}


#endif