LArVoxelCalculator.h

Go to the documentation of this file.

#ifndef sim_LArVoxelCalculator_H
#define sim_LArVoxelCalculator_H

#include "art/Framework/Services/Registry/ServiceDeclarationMacros.h"

namespace fhicl {
  class ParameterSet;
}

#include <vector>

namespace sim {

  class LArVoxelCalculator {
  public:
    LArVoxelCalculator(fhicl::ParameterSet const& pset);

    double VoxelSizeX() const { return m_voxelSize[0]; }
    double VoxelSizeY() const { return m_voxelSize[1]; }
    double VoxelSizeZ() const { return m_voxelSize[2]; }
    double VoxelSizeT() const { return m_voxelSize[3]; }
    double VoxelOffsetX() const { return m_voxelOffset[0]; }
    double VoxelOffsetY() const { return m_voxelOffset[1]; }
    double VoxelOffsetZ() const { return m_voxelOffset[2]; }
    double VoxelOffsetT() const { return m_voxelOffset[3]; }

    double EnergyCut() const { return m_energyCut; }

    double SuggestedStepSize() const;

    int AxisToBin(const int, const double) const;

    int XAxisToBin(const double value) const { return AxisToBin(0, value); }
    int YAxisToBin(const double value) const { return AxisToBin(1, value); }
    int ZAxisToBin(const double value) const { return AxisToBin(2, value); }
    int TAxisToBin(const double value) const { return AxisToBin(3, value); }

    double BinToAxis(const int, const int) const;

    double XBinToAxis(const int value) const { return BinToAxis(0, value); }
    double YBinToAxis(const int value) const { return BinToAxis(1, value); }
    double ZBinToAxis(const int value) const { return BinToAxis(2, value); }
    double TBinToAxis(const int value) const { return BinToAxis(3, value); }

  private:
    typedef std::vector<double> vector_type;

    vector_type const m_voxelSize;

    vector_type const m_voxelOffset;

    double const m_energyCut;
  };

} // namespace sim

DECLARE_ART_SERVICE(sim::LArVoxelCalculator, SHARED)
#endif //  sim_LArVoxelCalculator_H