LArVoxelCalculator.cc

Go to the documentation of this file.

// Framework includes
#include "larsim/Simulation/LArVoxelCalculator.h"

#include "fhiclcpp/ParameterSet.h"

#include <cmath>

namespace sim {

  //----------------------------------------------------------------------------
  LArVoxelCalculator::LArVoxelCalculator(fhicl::ParameterSet const& pset)
    : m_voxelSize{pset.get<double>("VoxelSizeX"),
                  pset.get<double>("VoxelSizeY"),
                  pset.get<double>("VoxelSizeZ"),
                  pset.get<double>("VoxelSizeT")}
    , m_voxelOffset{pset.get<double>("VoxelOffsetX"),
                    pset.get<double>("VoxelOffsetY"),
                    pset.get<double>("VoxelOffsetZ"),
                    pset.get<double>("VoxelOffsetT")}
    , m_energyCut{pset.get<double>("VoxelEnergyCut")}
  {}

  //----------------------------------------------------------------------------
  double LArVoxelCalculator::SuggestedStepSize() const
  {
    return std::min(m_voxelSize[0], std::min(m_voxelSize[1], m_voxelSize[2]));
  }

  //----------------------------------------------------------------------------
  int LArVoxelCalculator::AxisToBin(const int axis, const double coord) const
  {
    // We have to be careful of how to handle the case when coord -
    // offset < 0.  The standard floor() function rounds the number in
    // the correct direction.
    return static_cast<int>(floor((coord - m_voxelOffset[axis]) / m_voxelSize[axis]));
  }

  //----------------------------------------------------------------------------
  double LArVoxelCalculator::BinToAxis(const int axis, const int bin) const
  {
    return (static_cast<double>(bin) + 0.5) * m_voxelSize[axis] + m_voxelOffset[axis];
  }

} // namespace sim