gar::geo::AuxDetSensitiveGeo Class Reference

#include <AuxDetSensitiveGeo.h>

Public Member Functions
	AuxDetSensitiveGeo (std::vector< const TGeoNode * > &path, int depth)
	AuxDetSensitiveGeo (const TGeoVolume *volume, TGeoHMatrix *rotation)
	~AuxDetSensitiveGeo ()
void	GetCenter (double *xyz, double localz=0.0) const
void	GetNormalVector (double *xyzDir) const
double	Length () const
double	HalfWidth1 () const
double	HalfWidth2 () const
double	HalfHeight () const
const TGeoVolume *	TotalVolume () const
double	DistanceToPoint (double *xyz) const
void	LocalToWorld (const double *local, double *world) const
void	LocalToWorldVect (const double *local, double *world) const
void	WorldToLocal (const double *world, double *local) const
void	WorldToLocalVect (const double *world, double *local) const

Private Attributes
TGeoHMatrix *	fGeoMatrix
	Transformation matrix to world frame. More...
const TGeoVolume *	fTotalVolume
	Total volume of AuxDet, called vol*. More...
double	fLength
	length of volume, along z direction in local More...
double	fHalfWidth1
	1st half width of volume, at -z/2 in local coordinates More...
double	fHalfWidth2
	2nd half width (width1==width2 for boxes), at +z/2 More...
double	fHalfHeight
	half height of volume More...

Detailed Description

Definition at line 20 of file AuxDetSensitiveGeo.h.

Constructor & Destructor Documentation

gar::geo::AuxDetSensitiveGeo::AuxDetSensitiveGeo	(	std::vector< const TGeoNode * > &	path,
		int	depth
	)

Definition at line 29 of file AuxDetSensitiveGeo.cxx.

    : fTotalVolume(0)
    {
      
      TGeoVolume *vc = path[depth]->GetVolume();
      if(vc){
        fTotalVolume = vc;
        if(!fTotalVolume)
          throw cet::exception("AuxDetSensitiveGeo") << "cannot find AuxDetSensitive volume\n";
        
      }// end if found volume
      
      MF_LOG_DEBUG("Geometry") << "detector sensitive total  volume is " << fTotalVolume->GetName();
      
      // Build the matrix that takes us to the top world frame
      // build a matrix to take us from the local to the world coordinates
      // in one step
      fGeoMatrix = new TGeoHMatrix(*path[0]->GetMatrix());
      for(int i = 1; i <= depth; ++i){
        fGeoMatrix->Multiply(path[i]->GetMatrix());
      }
      
      // 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;
      }
      
    }

gar::geo::AuxDetSensitiveGeo::AuxDetSensitiveGeo	(	const TGeoVolume *	volume,
		TGeoHMatrix *	rotation
	)

Definition at line 77 of file AuxDetSensitiveGeo.cxx.

    : fGeoMatrix(rotation)
    , fTotalVolume(volume)
    {
      
      MF_LOG_DEBUG("Geometry") << "detector sensitive total  volume is " << fTotalVolume->GetName();
      
      // 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;
      }
      
    }

gar::geo::AuxDetSensitiveGeo::~AuxDetSensitiveGeo

(

)

Definition at line 111 of file AuxDetSensitiveGeo.cxx.

    {
      if(fGeoMatrix)  delete fGeoMatrix;
      return;
    }

Member Function Documentation

double gar::geo::AuxDetSensitiveGeo::DistanceToPoint

(

double *

xyz

)

const

Definition at line 174 of file AuxDetSensitiveGeo.cxx.

    {
      double Center[3];
      GetCenter(Center);
      return std::sqrt((Center[0]-xyz[0])*(Center[0]-xyz[0]) +
                       (Center[1]-xyz[1])*(Center[1]-xyz[1]) +
                       (Center[2]-xyz[2])*(Center[2]-xyz[2]));
    }

void gar::geo::AuxDetSensitiveGeo::GetCenter	(	double *	xyz,
		double	localz = 0.0
	)		const

Return the center position of an AuxDet

Parameters

xyz	: 3-D array. The returned location.
localz	: Distance along the length of the volume (z) (cm). Default is 0.

Definition at line 158 of file AuxDetSensitiveGeo.cxx.

    {
      double xyzLocal[3] = {0.,0.,localz};
      this->LocalToWorld(xyzLocal, xyz);
    }

void gar::geo::AuxDetSensitiveGeo::GetNormalVector

(

double *

xyzDir

)

const

Definition at line 166 of file AuxDetSensitiveGeo.cxx.

    {
      double normal[3]={0.,0.,1.};
      this->LocalToWorldVect(normal,xyzDir);
    }

double gar::geo::AuxDetSensitiveGeo::HalfHeight ( ) const

inline

Definition at line 36 of file AuxDetSensitiveGeo.h.

{ return fHalfHeight;  }

double gar::geo::AuxDetSensitiveGeo::HalfWidth1 ( ) const

inline

Definition at line 34 of file AuxDetSensitiveGeo.h.

{ return fHalfWidth1;  }

double gar::geo::AuxDetSensitiveGeo::HalfWidth2 ( ) const

inline

Definition at line 35 of file AuxDetSensitiveGeo.h.

{ return fHalfWidth2;  }

double gar::geo::AuxDetSensitiveGeo::Length ( ) const

inline

Definition at line 33 of file AuxDetSensitiveGeo.h.

{ return fLength;      }

void gar::geo::AuxDetSensitiveGeo::LocalToWorld	(	const double *	local,
		double *	world
	)		const

Transform a position from local frame to world frame

Parameters

local	: 3D array. Position in the local frame Input.
world	: 3D array. Position in the world frame. Returned.

Definition at line 121 of file AuxDetSensitiveGeo.cxx.

    {
      fGeoMatrix->LocalToMaster(local,world);
    }

void gar::geo::AuxDetSensitiveGeo::LocalToWorldVect	(	const double *	local,
		double *	world
	)		const

Transform a 3-vector from local frame to world frame

Parameters

local	: 3D array. Position in the local frame Input.
world	: 3D array. Position in the world frame. Returned.

Definition at line 130 of file AuxDetSensitiveGeo.cxx.

    {
      fGeoMatrix->LocalToMasterVect(local,world);
    }

const TGeoVolume* gar::geo::AuxDetSensitiveGeo::TotalVolume ( ) const

inline

Definition at line 37 of file AuxDetSensitiveGeo.h.

{ return fTotalVolume; }

void gar::geo::AuxDetSensitiveGeo::WorldToLocal	(	const double *	world,
		double *	local
	)		const

Transform a position from world frame to local frame

Parameters

world	: 3D array. Position in the world frame. Input.
local	: 3D array. Position in the local frame Returned.

Definition at line 139 of file AuxDetSensitiveGeo.cxx.

    {
      fGeoMatrix->MasterToLocal(world,local);
    }

void gar::geo::AuxDetSensitiveGeo::WorldToLocalVect	(	const double *	world,
		double *	local
	)		const

Transform a 3-vector from world frame to local frame

Parameters

world	: 3D array. Position in the world frame. Input.
local	: 3D array. Position in the local frame Returned.

Definition at line 148 of file AuxDetSensitiveGeo.cxx.

    {
      fGeoMatrix->MasterToLocalVect(world,local);
    }

Member Data Documentation

TGeoHMatrix* gar::geo::AuxDetSensitiveGeo::fGeoMatrix

private

Transformation matrix to world frame.

Definition at line 48 of file AuxDetSensitiveGeo.h.

double gar::geo::AuxDetSensitiveGeo::fHalfHeight

private

half height of volume

Definition at line 53 of file AuxDetSensitiveGeo.h.

double gar::geo::AuxDetSensitiveGeo::fHalfWidth1

private

1st half width of volume, at -z/2 in local coordinates

Definition at line 51 of file AuxDetSensitiveGeo.h.

double gar::geo::AuxDetSensitiveGeo::fHalfWidth2

private

2nd half width (width1==width2 for boxes), at +z/2

Definition at line 52 of file AuxDetSensitiveGeo.h.

double gar::geo::AuxDetSensitiveGeo::fLength

private

length of volume, along z direction in local

Definition at line 50 of file AuxDetSensitiveGeo.h.

const TGeoVolume* gar::geo::AuxDetSensitiveGeo::fTotalVolume

private

Total volume of AuxDet, called vol*.

Definition at line 49 of file AuxDetSensitiveGeo.h.

---

The documentation for this class was generated from the following files:

• garsoft/Geometry/AuxDetSensitiveGeo.h
• garsoft/Geometry/AuxDetSensitiveGeo.cxx