GaisserHillas6Parameter.cc

Go to the documentation of this file.

#include <evt/GaisserHillas6Parameter.h>
#include <utl/ErrorLogger.h>
#include <utl/AugerUnits.h>
#include <utl/Math.h>

#include <cmath>

using namespace evt;
using namespace utl;
using namespace std;


GaisserHillas6Parameter::GaisserHillas6Parameter() :
  fXZero(0),
  fXZeroError(0),
  fA(0),
  fAError(0),
  fB(0),
  fBError(0),
  fC(0),
  fCError(0)
{ }


double
GaisserHillas6Parameter::Eval(const double depth)
  const
{
  if (depth <= fXZero)
    return 0;

  const double dx = fXMax - fXZero;
  const double lambda = fA + depth*fB + Sqr(depth)*fC;

  return fNMax * pow((depth - fXZero) / dx, dx / lambda) *
    exp((fXMax - depth) / lambda);
}


void
GaisserHillas6Parameter::Dump(ostream& os)
  const
{
  os << "GH function with 6 parameters\n"
        "  X_0   = " << fXZero/(g/cm2) << " g/cm2 (sigma = " << fXZeroError/(g/cm2) << ")\n"
        "  X_max = " << fXMax/(g/cm2) << " g/cm2 (sigma = " << fXMaxError/(g/cm2) << ")\n"
        "  N_max = " << fNMax << " (sigma = " << fNMaxError << ")\n"
        "  a     = " << fA/(g/cm2) << " (sigma = " << fAError/(g/cm2) << ")\n"
        "  b     = " << fB << " (sigma = " << fBError << ")\n"
        "  c     = " << fC*(g/cm2) << " (sigma = " << fCError*(g/cm2) << ")\n";
}


double
GaisserHillas6Parameter::GetIntegral()
  const
{
  WARNING("not implemented!");
  return 0;
}


double
GaisserHillas6Parameter::GetIntegralError()
  const
{
  WARNING("not implemented!");
  return 0;
}