driftvolumes_test.cxx File Reference

Test for neighbourhood discovery in simple geometries. More...

#include "test/Geometry/geometry_unit_test_base.h"
#include "larcorealg/Geometry/ChannelMapStandardAlg.h"
#include "larcorealg/Geometry/GeometryCore.h"
#include "larcorealg/Geometry/DriftPartitions.h"
#include "larcorealg/CoreUtils/DumpUtils.h"
#include "messagefacility/MessageLogger/MessageLogger.h"

Go to the source code of this file.

Typedefs
using	StandardGeometryConfiguration = testing::BasicGeometryEnvironmentConfiguration< geo::ChannelMapStandardAlg >
using	StandardGeometryTestEnvironment = testing::GeometryTesterEnvironment< StandardGeometryConfiguration >

Functions
int	main (int argc, char const **argv)
	Runs the test. More...

Detailed Description

Test for neighbourhood discovery in simple geometries.

Author

Gianluca Petrillo (petri.nosp@m.llo@.nosp@m.fnal..nosp@m.gov)

Date

July 13, 2017

Usage:

driftvolumes_test configuration.fcl

The configuration file must contain a geometry service provider configuration under services.Geometry, That geometry configuration must use the standard channel mapping.

Definition in file driftvolumes_test.cxx.

Typedef Documentation

using StandardGeometryConfiguration = testing::BasicGeometryEnvironmentConfiguration<geo::ChannelMapStandardAlg>

Definition at line 35 of file driftvolumes_test.cxx.

using StandardGeometryTestEnvironment = testing::GeometryTesterEnvironment<StandardGeometryConfiguration>

Definition at line 45 of file driftvolumes_test.cxx.

Function Documentation

int main	(	int	argc,
		char const **	argv
	)

Runs the test.

---

Parameters

argc	number of arguments in argv
argv	arguments to the function

Returns

number of detected errors (0 on success)

Exceptions

cet::exception

most of error situations throw

The arguments in argv are: 0. name of the executable ("Geometry_test")

1. path to the FHiCL configuration file
2. FHiCL path to the configuration of the geometry (default: services.Geometry)

Definition at line 67 of file driftvolumes_test.cxx.

                                      {

  StandardGeometryConfiguration config("driftvolumes_test");

  //
  // parameter parsing
  //
  int iParam = 0;

  // first argument: configuration file (mandatory)
  if (++iParam < argc) config.SetConfigurationPath(argv[iParam]);

  // second argument: path of the parameter set for geometry configuration
  // (optional; default: "services.Geometry" from the inherited object)
  if (++iParam < argc) config.SetGeometryParameterSetPath(argv[iParam]);

  //
  // testing environment setup
  //
  StandardGeometryTestEnvironment TestEnvironment(config);
  auto const& geom = *(TestEnvironment.Provider<geo::GeometryCore>());

  //
  // run the test algorithm
  //

  unsigned int nErrors = 0;
  for (auto const& cryo: geom.IterateCryostats()) {
    auto partition = geo::buildDriftVolumes(cryo);
    mf::LogVerbatim("driftvolumes_test")
      << "Partition for cryostat " << cryo.ID() << ":";
    partition.print(mf::LogVerbatim("driftvolumes_test"));

    //
    // test that the partition topology is correct
    //
    for (geo::TPCGeo const& TPC: geom.IterateTPCs(cryo.ID())) {

      auto const& center = TPC.GetCenter<geo::Point_t>();

      auto where = partition.TPCat(center);
      if (!where) {
        mf::LogProblem("driftvolumes_test")
          << "Center of TPC " << TPC.ID() << " " << lar::dump::vector3D(center)
          << " not assigned to any TPC!";
        ++nErrors;
      }
      else if (where->ID() != TPC.ID()) {
        mf::LogProblem log("driftvolumes_test");
        log
          << "Center of TPC " << TPC.ID() << " " << lar::dump::vector3D(center)
          << " assigned to TPC " << where->ID() << ":\n";
        where->PrintTPCInfo(log, "  ", /* verbosity */ 5);
        ++nErrors;
      }

      //
      // test ownership of a uniform distribution of points inside the TPC
      //
      constexpr int nXsteps = 5, nYsteps = 5, nZsteps = 5;
      const double xstep = TPC.HalfSizeX() / (nXsteps + 1);
      const double ystep = TPC.HalfSizeY() / (nYsteps + 1);
      const double zstep = TPC.HalfSizeZ() / (nZsteps + 1);
      for (int xs = -nXsteps; xs <= nXsteps; ++xs) {
        double const x = center.X() + xs * xstep;
        for (int ys = -nYsteps; ys <= nYsteps; ++ys) {
          double const y = center.Y() + ys * ystep;
          for (int zs = -nZsteps; zs <= nZsteps; ++zs) {
            double const z = center.Z() + zs * zstep;

            auto where = partition.TPCat({ x, y, z});
            if (!where) {
              mf::LogProblem("driftvolumes_test")
                << "Point " << lar::dump::vector3D(center) << " within TPC "
                << TPC.ID() << " (" << xs << "," << ys << "," << zs
                << ") not assigned to any TPC!";
              ++nErrors;
            }
            else if (where->ID() != TPC.ID()) {
              mf::LogProblem log("driftvolumes_test");
              log
                << "Point " << lar::dump::vector3D(center) << " within TPC "
                << TPC.ID() << " (" << xs << "," << ys << "," << zs
                << ") assigned to TPC " << where->ID() << ":\n";
              where->PrintTPCInfo(log, "  ", /* verbosity */ 5);
              ++nErrors;
            }

          } // for zs
        } // for ys
      } // for xs

    } // for TPCs

  } // for

  // 4. And finally we cross fingers.
  if (nErrors > 0) {
    mf::LogError("geometry_test") << nErrors << " errors detected!";
  }

  return nErrors;
} // main()