testPhysicalFunctions.cc

Go to the documentation of this file.

#include <cppunit/extensions/HelperMacros.h>

#include <utl/ErrorLogger.h>
#include <utl/MathConstants.h>
#include <utl/PhysicalConstants.h>
#include <utl/PhysicalFunctions.h>
#include <utl/AugerUnits.h>
#include <tst/Verify.h>
#include <utl/AugerException.h>
#include <string>

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


class testPhysicalFunctions : public CppUnit::TestFixture {

  CPPUNIT_TEST_SUITE(testPhysicalFunctions);
  CPPUNIT_TEST(testGaisserHillas);
  CPPUNIT_TEST(testEnergy);
  CPPUNIT_TEST(testBeta);
  CPPUNIT_TEST(testRefractionIndex);
  CPPUNIT_TEST(testCherenkovThreshold);
  CPPUNIT_TEST(testShowerAge);
  CPPUNIT_TEST(testEnergyDeposit);
  CPPUNIT_TEST_SUITE_END();

public:

  void testGaisserHillas() {

    double x0      = 12.89*(g/cm/cm);
    double nMax    = 6.245e+9;
    double xMax    = 700.0*(g/cm/cm);
    double l       = kLambdaGH;

    CPPUNIT_ASSERT(Verify<CloseTo>(GaisserHillas(xMax, x0, xMax, nMax, l),
                                 nMax));

  }

  void testEnergy() {

    double energy = 2.121*MeV;

    double beta = Beta(energy);
    CPPUNIT_ASSERT(Verify<CloseTo>(Energy(beta), energy));

  }

  void testBeta() {

    double energy = 12.123*MeV;
    double beta   = Beta(energy);
    CPPUNIT_ASSERT(Verify<CloseTo>(Beta(energy), beta));

  }

  void testRefractionIndex() {

    double depth  = 756.778*(g/cm/cm);

    CPPUNIT_ASSERT(Verify<CloseTo>(RefractionIndex::LorentzLorentz(depth), 1.0002145));

  }

  void testCherenkovThreshold() {

    double n = 1.5111;
    CPPUNIT_ASSERT(Verify<CloseTo>(CherenkovThreshold(1.5111),
                                      CherenkovThreshold(n)));

  }

  void testShowerAge() {
    double age = 1.0;
    CPPUNIT_ASSERT(Verify<CloseTo>(ShowerAge(300.3,300.3), age));
  }

  void testEnergyDeposit() {
    double age = ShowerAge(300*g/cm/cm,600*g/cm/cm);
    CPPUNIT_ASSERT(Verify<CloseTo>(EnergyDeposit(age,1*MeV),
                                      EnergyDeposit(300*g/cm/cm,600*g/cm/cm,1*MeV)));
  }


};

CPPUNIT_TEST_SUITE_REGISTRATION(testPhysicalFunctions);

// Configure (x)emacs for this file ...
// Local Variables:
// mode:c++
// compile-command: "make -C .. -k testTime && ../testTime"
// End: