AuxDetSensitiveGeo.h

Go to the documentation of this file.

#ifndef LARCOREALG_GEOMETRY_AUXDETSENSITIVEGEO_H
#define LARCOREALG_GEOMETRY_AUXDETSENSITIVEGEO_H

// LArSoft libraries
#include "larcorealg/CoreUtils/RealComparisons.h"
#include "larcorealg/Geometry/GeoVectorLocalTransformation.h"
#include "larcorealg/Geometry/LocalTransformationGeo.h"
#include "larcorealg/Geometry/TransformationMatrix.h"
#include "larcoreobj/SimpleTypesAndConstants/geo_vectors.h"

// ROOT libraries
#include "Math/GenVector/DisplacementVector3D.h"
#include "Math/GenVector/PositionVector3D.h"
#include "Math/GenVector/Transform3D.h"

// C/C++ standard libraries
#include <string>
#include <type_traits>

class TGeoNode;
class TGeoVolume;

namespace geo {

  class AuxDetSensitiveGeo {
  public:

    struct AuxDetSensitiveGeoCoordinatesTag {};

    using LocalPoint_t = geo::Point3DBase_t<AuxDetSensitiveGeoCoordinatesTag>;

    using LocalVector_t = geo::Vector3DBase_t<AuxDetSensitiveGeoCoordinatesTag>;


    AuxDetSensitiveGeo(TGeoNode const& node, geo::TransformationMatrix&& trans);

    geo::Point_t GetCenter(double localz = 0.0) const;

    geo::Vector_t GetNormalVector() const;

    //box geometry
    double Length() const { return fLength; }
    double HalfLength() const { return Length() / 2.0; }
    double HalfWidth1() const { return fHalfWidth1; }
    double HalfWidth2() const { return fHalfWidth2; }
    double HalfCenterWidth() const { return (HalfWidth1() + HalfWidth2()) / 2.0; }
    double HalfHeight() const { return fHalfHeight; }
    const TGeoVolume* TotalVolume() const { return fTotalVolume; }

    geo::Length_t DistanceToPoint(geo::Point_t const& point) const
    {
      return (point - GetCenter()).R();
    }
    geo::Length_t DistanceToPoint(double const* point) const;


    geo::Point_t toWorldCoords(LocalPoint_t const& local) const
    {
      return fTrans.toWorldCoords(local);
    }

    geo::Vector_t toWorldCoords(LocalVector_t const& local) const
    {
      return fTrans.toWorldCoords(local);
    }

    LocalPoint_t toLocalCoords(geo::Point_t const& world) const
    {
      return fTrans.toLocalCoords(world);
    }

    LocalVector_t toLocalCoords(geo::Vector_t const& world) const
    {
      return fTrans.toLocalCoords(world);
    }


    template <typename Stream>
    void PrintAuxDetInfo(Stream&& out, std::string indent = "", unsigned int verbosity = 0) const;

    std::string AuxDetInfo(std::string indent = "", unsigned int verbosity = 1) const;

    static constexpr unsigned int MaxVerbosity = 2;

  private:
    using LocalTransformation_t =
      geo::LocalTransformationGeo<ROOT::Math::Transform3D, LocalPoint_t, LocalVector_t>;

    LocalTransformation_t fTrans;   
    const TGeoVolume* fTotalVolume; 
    double fLength;                 
    double fHalfWidth1;             
    double fHalfWidth2;             
    double fHalfHeight;             

    void InitShapeSize();

  }; // class AuxDetSensitiveGeo

  static_assert(std::is_move_assignable_v<geo::AuxDetSensitiveGeo>);
  static_assert(std::is_move_constructible_v<geo::AuxDetSensitiveGeo>);

} // namespace geo

//------------------------------------------------------------------------------
//--- template implementation
//---
template <typename Stream>
void geo::AuxDetSensitiveGeo::PrintAuxDetInfo(Stream&& out,
                                              std::string indent /* = "" */,
                                              unsigned int verbosity /* = 0 */
) const
{

  //----------------------------------------------------------------------------
  out << "centered at " << GetCenter() << " cm";

  if (verbosity-- <= 0) return; // 0

  //----------------------------------------------------------------------------
  lar::util::RealComparisons<double> coordIs(1e-4);
  out << ", size ( " << (2.0 * HalfWidth1());
  if (coordIs.nonEqual(HalfWidth1(), HalfWidth2())) out << "/" << (2.0 * HalfWidth2());
  out << " x " << (2.0 * HalfHeight()) << " x " << Length() << " ) cm";

  if (verbosity-- <= 0) return; // 1

  //----------------------------------------------------------------------------
  out << ", normal facing " << GetNormalVector();

  //  if (verbosity-- <= 0) return; // 2

  //----------------------------------------------------------------------------

} // geo::AuxDetSensitiveGeo::PrintAuxDetInfo()

//------------------------------------------------------------------------------

#endif // LARCOREALG_GEOMETRY_AUXDETSENSITIVEGEO_H