testParticle.cc

Go to the documentation of this file.

#include <string>
#include <cmath>

#include <utl/Particle.h>
#include <utl/CoordinateSystem.h>
#include <utl/TimeInterval.h>
#include <utl/AugerUnits.h>
#include <utl/PhysicalConstants.h>
#include <utl/Math.h>

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

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

class testParticle
  : public CppUnit::TestFixture
{
  CPPUNIT_TEST_SUITE(testParticle);
  CPPUNIT_TEST(testProperties);
  CPPUNIT_TEST(testMomentum);
  //  CPPUNIT_TEST(testFN);
  //  CPPUNIT_TEST_EXCEPTION(testFN, anException);
  CPPUNIT_TEST_SUITE_END();
public:
  // shared data for tests

  void setUp() { }

  void tearDown() { }

  void testProperties()
  {
    CoordinateSystemPtr cs = CoordinateSystem::GetRootCoordinateSystem();
    Point        pos(0, 0, 0, cs);
    Vector       dir(0, 0, 0, cs);
    TimeInterval t;

    Particle p1(Particle::eElectron, Particle::eShower, pos, dir, t, 0, 0);
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetType(), int(Particle::eElectron)));
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetName(), string("e-")));
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetMass(), kElectronMass));
  }

  void testMomentum()
  {
    CoordinateSystemPtr cs = CoordinateSystem::GetRootCoordinateSystem();
    Point        pos(0, 0, 0, cs);
    Vector       p(0, 1*MeV, 0, cs);
    Vector       dir(p);
    dir.Normalize();
    TimeInterval t;

    Particle p1(Particle::ePhoton, Particle::eShower, pos, p, t, 0);
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetType(), int(Particle::ePhoton)));
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetName(), string("gamma")));
    CPPUNIT_ASSERT(Verify<Equal>(p1.GetMass(), 0.));
    CPPUNIT_ASSERT(Verify<CloseTo>(p1.GetKineticEnergy(), 1.*MeV));
    CPPUNIT_ASSERT(Verify<CloseTo>(p1.GetTotalEnergy(), 1.*MeV));

    Particle p2(Particle::eProton, Particle::eShower, pos, p, t, 0);
    CPPUNIT_ASSERT(Verify<Equal>(p2.GetType(), int(Particle::eProton)));
    CPPUNIT_ASSERT(Verify<Equal>(p2.GetName(), string("proton")));
    CPPUNIT_ASSERT(Verify<Equal>(p2.GetMass(), kProtonMass));
    CPPUNIT_ASSERT(Verify<CloseTo>(p2.GetKineticEnergy(),
                                 sqrt(Sqr(1.*MeV) + Sqr(p2.GetMass()))-p2.GetMass()));
    CPPUNIT_ASSERT(Verify<CloseTo>(p2.GetTotalEnergy() - p2.GetKineticEnergy(),
                                 p2.GetMass()));

    Vector p1000(1e3*p);
    p2.SetMomentum(p1000);
    CPPUNIT_ASSERT(Verify<CloseTo>(p1000.GetCoordinates(cs),
                                      p2.GetMomentum().GetCoordinates(cs)));
    CPPUNIT_ASSERT(Verify<CloseTo>(dir.GetCoordinates(cs),
                                      p2.GetDirection().GetCoordinates(cs)));
    CPPUNIT_ASSERT(Verify<CloseTo>(p2.GetKineticEnergy(),
                                 sqrt(Sqr(1.*GeV) + Sqr(p2.GetMass()))-p2.GetMass()));
  }

};


CPPUNIT_TEST_SUITE_REGISTRATION(testParticle);


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