AuxDetSensitiveGeo.cxx

Go to the documentation of this file.

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

// LArSoft libraries
#include "larcorealg/Geometry/geo_vectors_utils.h" // geo::vect namespace

#include "TGeoTrd2.h"
#include "TGeoBBox.h"

// Framework includes
#include "messagefacility/MessageLogger/MessageLogger.h"
#include "cetlib_except/exception.h"


namespace geo{
  
  //-----------------------------------------
  AuxDetSensitiveGeo::AuxDetSensitiveGeo(std::vector<const TGeoNode*> const& path, int depth)
    : fTrans(path, depth)
  {
    
    fTotalVolume = path[depth]->GetVolume();
    if(!fTotalVolume){
      throw cet::exception("AuxDetSensitiveGeo") << "cannot find AuxDetSensitive volume\n";
    }
    
    LOG_DEBUG("Geometry") << "detector sensitive total  volume is " << fTotalVolume->GetName();
    
    InitShapeSize();

  }

  //-----------------------------------------
  AuxDetSensitiveGeo::AuxDetSensitiveGeo(const TGeoVolume* volume, 
                                         TGeoHMatrix const& rotation)
    : fTrans(rotation)
    , fTotalVolume(volume)
  {
    assert(fTotalVolume);
    LOG_DEBUG("Geometry") << "detector sensitive total  volume is " << fTotalVolume->GetName();
    
    InitShapeSize();

  }
  
  //-----------------------------------------
  AuxDetSensitiveGeo::AuxDetSensitiveGeo(const TGeoVolume* volume, 
                                         TGeoHMatrix&& rotation)
    : fTrans(std::move(rotation))
    , fTotalVolume(volume)
  {
    assert(fTotalVolume);
    LOG_DEBUG("Geometry") << "detector sensitive total  volume is " << fTotalVolume->GetName();
    
    InitShapeSize();

  }
  
  //......................................................................
  geo::Point_t AuxDetSensitiveGeo::GetCenter(double localz /* = 0.0 */) const
    { return toWorldCoords(LocalPoint_t{ 0.0, 0.0, localz }); }
  
  //......................................................................
  void AuxDetSensitiveGeo::GetCenter(double* xyz, double localz) const {
    auto const& center = GetCenter(localz);
    xyz[0] = center.X();
    xyz[1] = center.Y();
    xyz[2] = center.Z();
  } // AuxDetSensitiveGeo::GetCenter(double*)
  
  //......................................................................
  
  // Return the unit normal vector (0,0,1) in local coordinates to global coordinates
  geo::Vector_t AuxDetSensitiveGeo::GetNormalVector() const
    { return toWorldCoords(geo::Zaxis<LocalVector_t>()); }
  
  //......................................................................
  
  // Return the unit normal vector (0,0,1) in local coordinates to global coordinates
  void AuxDetSensitiveGeo::GetNormalVector(double* xyzDir) const {
    auto const& norm = GetNormalVector();
    xyzDir[0] = norm.X();
    xyzDir[1] = norm.Y();
    xyzDir[2] = norm.Z();
  } // AuxDetSensitiveGeo::GetNormalVector(double*)
  
  
  //......................................................................
  geo::Length_t AuxDetSensitiveGeo::DistanceToPoint(double const* point) const
    { return DistanceToPoint(geo::vect::makePointFromCoords(point)); }
  
  
  //......................................................................
  void AuxDetSensitiveGeo::InitShapeSize() {
    // set the ends depending on whether the shape is a box or trapezoid
    std::string volName(fTotalVolume->GetName());
    if( volName.find("Trap") != std::string::npos ) {

      //       Small Width
      //          ____          Height is the thickness
      //         /    \     T     of the trapezoid
      //        /      \    |
      //       /        \   | Length
      //      /__________\  _ 
      //         Width 
      fHalfHeight      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDy1(); // same as Dy2()
      fLength          = 2.0*((TGeoTrd2*)fTotalVolume->GetShape())->GetDz();
      fHalfWidth1      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDx1(); // at -Dz
      fHalfWidth2      =     ((TGeoTrd2*)fTotalVolume->GetShape())->GetDx2(); // at +Dz
    } 
    else {
      fHalfWidth1      =     ((TGeoBBox*)fTotalVolume->GetShape())->GetDX();
      fHalfHeight      =     ((TGeoBBox*)fTotalVolume->GetShape())->GetDY();
      fLength          = 2.0*((TGeoBBox*)fTotalVolume->GetShape())->GetDZ();
      fHalfWidth2      = fHalfWidth1;
    }
  } // AuxDetSensitiveGeo::InitShapeSize()
}