testGaisserHillas4Parameter.cc

Go to the documentation of this file.

#include <evt/GaisserHillas4Parameter.h>
#include <Shower/testGaisserHillas4Parameter.h>

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

#include <tst/Verify.h>
#include <cppunit/extensions/HelperMacros.h>
#include <vector>

#include <iostream>

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


void
testGaisserHillas4Parameter::CheckClassicToWidthConversion()
{

  const double xMax = 700;
  const double xMaxError = 30;
  const double nMax = 1e6;
  const double nMaxError = 1e5;
  const double x0 = -150;
  const double x0Error = 2;
  const double lambda = 60;
  const double lambdaError = 6;

  GaisserHillas4Parameter ghClassic(gh::eClassic);

  ghClassic.SetXMax(xMax, xMaxError);
  ghClassic.SetNMax(nMax, nMaxError);
  ghClassic.SetShapeParameter(gh::eX0, x0, x0Error);
  ghClassic.SetShapeParameter(gh::eLambda, lambda,lambdaError);
  //ghClassic.SetChiSquare(0.1, 1);

  GaisserHillas4Parameter ghWidth(gh::eWidth);

  ghWidth.SetXMax( ghClassic.GetXMax(), ghClassic.GetXMaxError() );
  ghWidth.SetNMax( ghClassic.GetNMax(), ghClassic.GetXMaxError() );
  ghWidth.SetShapeParameter(gh::eFWHM,
                            ghClassic.GetShapeParameter(gh::eFWHM),
                            ghClassic.GetShapeParameterError(gh::eFWHM));
  ghWidth.SetShapeParameter(gh::eAsym,
                            ghClassic.GetShapeParameter(gh::eAsym),
                            ghClassic.GetShapeParameterError(gh::eAsym));

  CPPUNIT_ASSERT(Verify<Equal>(ghWidth.GetXMax(), xMax));
  CPPUNIT_ASSERT(Verify<Equal>(ghWidth.GetNMax(), nMax));
  CPPUNIT_ASSERT(Verify<CloseTo>(ghWidth.GetShapeParameter(gh::eX0), x0));
  CPPUNIT_ASSERT(Verify<CloseTo>(ghWidth.GetShapeParameter(gh::eLambda), lambda));

}

void
testGaisserHillas4Parameter::CheckWidthToClassicConversion()
{

  const double xMax = 700;
  const double xMaxError = 30;
  const double nMax = 1e8;
  const double nMaxError = 1e7;
  const double fwhm = 550;
  const double fwhmError = 50;
  const double asym = 0.45;
  const double asymError = 0.01;

  GaisserHillas4Parameter ghWidth(gh::eWidth);

  ghWidth.SetXMax(xMax,xMaxError);
  ghWidth.SetNMax(nMax,nMaxError);
  ghWidth.SetShapeParameter(gh::eFWHM, fwhm, fwhmError);
  ghWidth.SetShapeParameter(gh::eAsym, asym, asymError);

  GaisserHillas4Parameter ghClassic(gh::eClassic);

  ghClassic.SetXMax(ghWidth.GetXMax(), ghWidth.GetXMaxError());
  ghClassic.SetNMax(ghWidth.GetNMax(), ghWidth.GetXMaxError());
  ghClassic.SetShapeParameter(gh::eX0,
                              ghWidth.GetShapeParameter(gh::eX0),
                              ghWidth.GetShapeParameterError(gh::eX0));
  ghClassic.SetShapeParameter(gh::eLambda,
                              ghWidth.GetShapeParameter(gh::eLambda),
                              ghWidth.GetShapeParameterError(gh::eLambda));

  CPPUNIT_ASSERT(Verify<Equal>(ghClassic.GetXMax(), xMax));
  CPPUNIT_ASSERT(Verify<Equal>(ghClassic.GetNMax(), nMax));
  CPPUNIT_ASSERT(Verify<CloseTo>(ghClassic.GetShapeParameter(gh::eFWHM), fwhm));
  CPPUNIT_ASSERT(Verify<CloseTo>(ghClassic.GetShapeParameter(gh::eAsym), asym));

}

void
testGaisserHillas4Parameter::CheckWrongSetterWidth()
{
  GaisserHillas4Parameter ghWidth(gh::eWidth);
  ghWidth.SetShapeParameter(gh::eX0, 10, 10);
}

void
testGaisserHillas4Parameter::CheckWrongSetterClassic()
{
  GaisserHillas4Parameter ghClassic(gh::eClassic);
  ghClassic.SetShapeParameter(gh::eFWHM, 10, 10);
}