LArVoxelCalculator_service.cc

Go to the documentation of this file.

// Framework includes

#include "larsim/Simulation/LArVoxelCalculator.h"

#include <cmath>

namespace sim {

  //----------------------------------------------------------------------------
  LArVoxelCalculator::LArVoxelCalculator(fhicl::ParameterSet const& pset, 
                                         art::ActivityRegistry & /* reg */) 
  {

    this->reconfigure(pset);

  }

  //----------------------------------------------------------------------------
  LArVoxelCalculator::~LArVoxelCalculator()
  {
  }

  //----------------------------------------------------------------------------
  void LArVoxelCalculator::reconfigure(fhicl::ParameterSet const& pset)
  {
    m_energyCut = pset.get<double>("VoxelEnergyCut");

    double voxelSizeX   = pset.get<double>("VoxelSizeX");
    double voxelSizeY   = pset.get<double>("VoxelSizeY");
    double voxelSizeZ   = pset.get<double>("VoxelSizeZ");
    double voxelSizeT   = pset.get<double>("VoxelSizeT");
    double voxelOffsetX = pset.get<double>("VoxelOffsetX");
    double voxelOffsetY = pset.get<double>("VoxelOffsetY");
    double voxelOffsetZ = pset.get<double>("VoxelOffsetZ");
    double voxelOffsetT = pset.get<double>("VoxelOffsetT");
    
    m_voxelSize.clear();
    m_voxelSize.push_back( voxelSizeX );
    m_voxelSize.push_back( voxelSizeY );
    m_voxelSize.push_back( voxelSizeZ );
    m_voxelSize.push_back( voxelSizeT );
    m_voxelOffset.clear();
    m_voxelOffset.push_back( voxelOffsetX );
    m_voxelOffset.push_back( voxelOffsetY );
    m_voxelOffset.push_back( voxelOffsetZ );
    m_voxelOffset.push_back( voxelOffsetT );

    return;
  }

  //----------------------------------------------------------------------------
  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
 
namespace sim {

  DEFINE_ART_SERVICE(LArVoxelCalculator)


} // namespace sim