OpDetGeo.cxx

Go to the documentation of this file.

// class header
#include "larcorealg/Geometry/OpDetGeo.h"

// LArSoft libraries
#include "larcorealg/Geometry/geo_vectors_utils.h"                // geo::vect::makeFromCoords()
#include "larcoreobj/SimpleTypesAndConstants/PhysicalConstants.h" // util::pi()

// ROOT libraries
#include "TGeoManager.h"
#include "TGeoNode.h"
#include "TGeoTube.h"

// C/C++ standard libraries
#include <cmath>

namespace geo {

  //-----------------------------------------
  OpDetGeo::OpDetGeo(TGeoNode const* node, TransformationMatrix&& trans)
    : fTrans(std::move(trans)), fOpDetNode{node}, fCenter{toWorldCoords(origin<LocalPoint_t>())}
  {}

  //......................................................................

  double OpDetGeo::RMax() const
  {
    if (TGeoSphere const* sphere = asSphere(); sphere) { return sphere->GetRmax(); }
    if (TGeoTube const* tube = asTube(); tube) { return tube->GetRmax(); }
    throw std::bad_cast{};
  }

  //......................................................................

  double OpDetGeo::HalfL() const
  {
    TGeoBBox const* pBox = asBox();
    return pBox ? pBox->GetDZ() : 0.0;
  }

  //......................................................................

  double OpDetGeo::HalfW() const
  {
    TGeoBBox const* pBox = asBox();
    return pBox ? pBox->GetDX() : 0.0;
  }

  //......................................................................

  double OpDetGeo::HalfH() const
  {
    TGeoBBox const* pBox = asBox();
    return pBox ? pBox->GetDY() : 0.0;
  }

  //......................................................................

  double OpDetGeo::RMin() const
  {
    if (TGeoSphere const* sphere = asSphere(); sphere) { return sphere->GetRmin(); }
    if (TGeoTube const* tube = asTube(); tube) { return tube->GetRmin(); }
    throw std::bad_cast{};
  }

  //......................................................................
  double OpDetGeo::ThetaZ() const
  {
    auto const& center = GetCenter();
    auto const& end = toWorldCoords(LocalPoint_t{0.0, 0.0, HalfL()});

    // TODO change this into something generic
    //either y or x will be 0, so adding both will always catch the right
    //one
    double angle = (end.Y() - center.Y() + end.X() - center.X()) /
                   std::abs(end.Y() - center.Y() + center.X() - end.X()) *
                   std::acos((end.Z() - center.Z()) / HalfL());
    if (angle < 0) angle += util::pi();
    return angle;
  }

  //......................................................................
  double OpDetGeo::ThetaZ(bool degree) const
  {
    return degree ? util::RadiansToDegrees(ThetaZ()) : ThetaZ();
  }

  //......................................................................
  double OpDetGeo::DistanceToPoint(Point_t const& point) const
  {
    return (point - GetCenter()).R();
  }

  //......................................................................
  std::string OpDetGeo::OpDetInfo(std::string indent /* = "" */,
                                  unsigned int verbosity /* = 0 */) const
  {
    std::ostringstream sstr;
    PrintOpDetInfo(sstr, indent, verbosity);
    return sstr.str();
  }

  //......................................................................
  double OpDetGeo::CosThetaFromNormal(Point_t const& point) const
  {
    auto const& local = toLocalCoords(point);
    return local.Z() / local.R();
  }

  //......................................................................
  void OpDetGeo::UpdateAfterSorting(OpDetID opdetid)
  {
    fID = opdetid;
  }

}