SurfXYZTest.cc File Reference

#include "boost/test/unit_test.hpp"
#include <cmath>
#include "lardata/RecoObjects/SurfXYZPlane.h"
#include "lardata/RecoObjects/KalmanLinearAlgebra.h"
#include "cetlib_except/exception.h"

Go to the source code of this file.

Classes
struct	SurfXYZTestFixture

Macros
#define	BOOST_TEST_MODULE   ( SurfXYZTest )

Functions
	BOOST_AUTO_TEST_CASE (Equality)
	BOOST_AUTO_TEST_CASE (Parallel)
	BOOST_AUTO_TEST_CASE (Transformation)
	BOOST_AUTO_TEST_CASE (Separation)
	BOOST_AUTO_TEST_CASE (NotParallel)
	BOOST_AUTO_TEST_CASE (TrackParameters)

Variables
auto const	tol = 1.e-6% tolerance()

Macro Definition Documentation

#define BOOST_TEST_MODULE   ( SurfXYZTest )

Definition at line 1 of file SurfXYZTest.cc.

Function Documentation

BOOST_AUTO_TEST_CASE

(

Equality

)

Definition at line 35 of file SurfXYZTest.cc.

                               {
  BOOST_TEST(surf1.isEqual(surf2));
  BOOST_TEST(!surf1.isEqual(surf3));
  BOOST_TEST(!surf1.isEqual(surf4));
  BOOST_TEST(!surf2.isEqual(surf3));
  BOOST_TEST(!surf2.isEqual(surf4));
  BOOST_TEST(!surf3.isEqual(surf4));
}

BOOST_AUTO_TEST_CASE

(

Parallel

)

Definition at line 46 of file SurfXYZTest.cc.

                               {
  BOOST_TEST(surf1.isParallel(surf2));
  BOOST_TEST(surf1.isParallel(surf3));
  BOOST_TEST(!surf1.isParallel(surf4));
  BOOST_TEST(surf2.isParallel(surf3));
  BOOST_TEST(!surf2.isParallel(surf4));
  BOOST_TEST(!surf3.isParallel(surf4));
}

BOOST_AUTO_TEST_CASE

(

Transformation

)

Definition at line 57 of file SurfXYZTest.cc.

                                     {
  double xyz1[3] = {1., 2., 3.};
  double xyz2[3];
  double uvw[3];
  surf4.toLocal(xyz1, uvw);
  surf4.toGlobal(uvw, xyz2);
  for(int i=0; i<3; ++i)
    BOOST_TEST(xyz1[i] == xyz2[i], tol);
}

BOOST_AUTO_TEST_CASE

(

Separation

)

Definition at line 69 of file SurfXYZTest.cc.

                                 {
  BOOST_TEST(surf1.distanceTo(surf2) == 0.);
  BOOST_TEST(surf1.distanceTo(surf3) == 1.);
  BOOST_TEST(surf3.distanceTo(surf1) == -1.);
}

BOOST_AUTO_TEST_CASE

(

NotParallel

)

Definition at line 77 of file SurfXYZTest.cc.

                                  {
  BOOST_CHECK_EXCEPTION( surf1.distanceTo(surf4), cet::exception, \
                         [](cet::exception const & e)             \
                         {                                        \
                           return e.category() == "SurfXYZPlane";  \
                         } );
}

BOOST_AUTO_TEST_CASE

(

TrackParameters

)

Definition at line 87 of file SurfXYZTest.cc.

                                      {
  trkf::TrackVector v(5);
  v(0) = 0.1;   // u.
  v(1) = 0.2;   // v.
  v(2) = 2.;    // du/dw.
  v(3) = 3.;    // dv/dw.
  v(4) = 0.5;   // p = 2 GeV.

  // For this vector, the direction cosines are.
  // du/ds = 2./sqrt(14.);
  // dv/ds = 3./sqrt(14.);
  // dw/ds = 1./sqrt(14.);

  double xyz[3];
  double mom[3];
  surf1.getPosition(v, xyz);
  BOOST_TEST(xyz[0] == 0.1, tol);
  BOOST_TEST(xyz[1] == 0.2, tol);
  BOOST_TEST(xyz[2] == 0.0, tol);
  surf3.getPosition(v, xyz);
  BOOST_TEST(xyz[0] == 1.1, tol);
  BOOST_TEST(xyz[1] == 1.2, tol);
  BOOST_TEST(xyz[2] == 1.0, tol);
  surf1.getMomentum(v, mom, trkf::Surface::FORWARD);
  BOOST_TEST(mom[0] == 4./std::sqrt(14.), tol);
  BOOST_TEST(mom[1] == 6./std::sqrt(14.), tol);
  BOOST_TEST(mom[2] == 2./std::sqrt(14.), tol);
  surf1.getMomentum(v, mom, trkf::Surface::BACKWARD);
  BOOST_TEST(mom[0] == -4./std::sqrt(14.), tol);
  BOOST_TEST(mom[1] == -6./std::sqrt(14.), tol);
  BOOST_TEST(mom[2] == -2./std::sqrt(14.), tol);
  surf4.getMomentum(v, mom, trkf::Surface::FORWARD);
  BOOST_TEST(mom[0] == (4.*std::cos(0.1) + 2.*std::sin(0.1))/std::sqrt(14.), tol);
  BOOST_TEST(mom[1] == (4.*std::sin(0.1)*std::sin(1.) + 6.*std::cos(1.) - 2.*std::cos(0.1)*std::sin(1.))/std::sqrt(14.), tol);
  BOOST_TEST(mom[2] == (-4.*std::sin(0.1)*std::cos(1.) + 6.*std::sin(1.) + 2.*std::cos(0.1)*std::cos(1.))/std::sqrt(14.), tol);

  // Should throw exception (no direction).

  BOOST_CHECK_EXCEPTION( surf1.getMomentum(v, mom), cet::exception, \
                         [](cet::exception const & e)               \
                         {                                          \
                           return e.category() == "SurfXYZPlane";    \
                         } );
}

Variable Documentation

auto const tol = 1.e-6% tolerance()

Definition at line 16 of file SurfXYZTest.cc.