#include <AuxDetGeo.h>

Public Member Functions
	AuxDetGeo (std::vector< const TGeoNode * > &path, int depth)

	~AuxDetGeo ()

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

std::string	Name () const

size_t	FindSensitiveVolume (double const worldLoc[3]) const

AuxDetSensitiveGeo const &	PositionToSensitiveVolume (double const worldLoc[3], size_t &sv) const

AuxDetSensitiveGeo const &	SensitiveVolume (size_t sv) const

size_t	NSensitiveVolume () const

void	SortSubVolumes (geo::AuxDetGeoObjectSorter const &sorter)

Private Member Functions
void	FindAuxDetSensitive (std::vector< const TGeoNode * > &path, unsigned int depth)

void	MakeAuxDetSensitive (std::vector< const TGeoNode * > &path, int depth)

Private Attributes
const TGeoVolume *	fTotalVolume
	Total volume of AuxDet, called vol*. More...

TGeoHMatrix *	fGeoMatrix
	Transformation matrix to world frame. 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...

std::vector< AuxDetSensitiveGeo * >	fSensitive
	sensitive volumes in the detector More...

Detailed Description

Definition at line 24 of file AuxDetGeo.h.

Constructor & Destructor Documentation

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

Definition at line 30 of file AuxDetGeo.cxx.

     : fTotalVolume(0)
     {
       
       TGeoVolume *vc = path[depth]->GetVolume();
       if(vc){
         fTotalVolume = vc;
         if(!fTotalVolume)
           throw cet::exception("AuxDetGeo") << "cannot find AuxDet volume\n";
         
       }// end if found volume
       
       MF_LOG_DEBUG("Geometry") << "detector 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());
       }
       
       // look for sensitive volumes - if there are none then this aux det
       // could be from an older gdml file than the introduction of AuxDetSensitiveGeo
       // in that case assume the full AuxDetGeo is sensitive and copy its information
       // into a single AuxDetSensitive
       this->FindAuxDetSensitive(path, depth);
       if( fSensitive.size() < 1) fSensitive.push_back(new AuxDetSensitiveGeo(fTotalVolume, fGeoMatrix));
       
       // 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::AuxDetGeo::~AuxDetGeo ( )

Definition at line 86 of file AuxDetGeo.cxx.

     {
       if(fGeoMatrix)  delete fGeoMatrix;
       return;
     }

Member Function Documentation

double gar::geo::AuxDetGeo::DistanceToPoint ( double * xyz ) const

Definition at line 181 of file AuxDetGeo.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::AuxDetGeo::FindAuxDetSensitive	(	std::vector< const TGeoNode * > &	path,
		unsigned int	depth
	)

private

Definition at line 93 of file AuxDetGeo.cxx.

     {
       // Check if the current node is a senstive volume - we already know
       // we are in an Auxiliary detector
       std::string pathName(path[depth]->GetName());
       if( pathName.find("Sensitive") != std::string::npos){
         this->MakeAuxDetSensitive(path, depth);
         return;
       }
       
       // Explore the next layer down
       unsigned int deeper = depth+1;
       if (deeper>=path.size()) {
         throw cet::exception("ExceededMaxDepth") << "Exceeded maximum depth\n";
       }
       const TGeoVolume* v = path[depth]->GetVolume();
       int nd = v->GetNdaughters();
       for (int i=0; i<nd; ++i) {
         path[deeper] = v->GetNode(i);
         this->FindAuxDetSensitive(path, deeper);
       }
     }

size_t gar::geo::AuxDetGeo::FindSensitiveVolume ( double const worldLoc[3] ) const

Definition at line 191 of file AuxDetGeo.cxx.

     {
       double local[3] = {0.};
       for(unsigned int a = 0; a < fSensitive.size(); ++a) {
         
         this->SensitiveVolume(a).WorldToLocal(worldPos, local);
         double HalfCenterWidth = (this->SensitiveVolume(a).HalfWidth1() + this->SensitiveVolume(a).HalfWidth2()) / 2;
         
         if( local[2] >= - this->SensitiveVolume(a).Length()/2       &&
            local[2] <=   this->SensitiveVolume(a).Length()/2       &&
            local[1] >= - this->SensitiveVolume(a).HalfHeight()     &&
            local[1] <=   this->SensitiveVolume(a).HalfHeight()     &&
            // if SensitiveVolume a is a box, then HalfSmallWidth = HalfWidth
            local[0] >= - HalfCenterWidth + local[2]*(HalfCenterWidth-this->SensitiveVolume(a).HalfWidth2())/(this->SensitiveVolume(a).Length()/2) &&
            local[0] <=   HalfCenterWidth - local[2]*(HalfCenterWidth-this->SensitiveVolume(a).HalfWidth2())/(this->SensitiveVolume(a).Length()/2)
            )  return a;
         
       }// for loop over AuxDetSensitive a
       
       throw cet::exception("Geometry") << "Can't find AuxDetSensitive for position ("
       << worldPos[0] << ","
       << worldPos[1] << ","
       << worldPos[2] << ")\n";
       
       return UINT_MAX;
     } 

void gar::geo::AuxDetGeo::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 164 of file AuxDetGeo.cxx.

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

void gar::geo::AuxDetGeo::GetNormalVector ( double * xyzDir ) const

Definition at line 172 of file AuxDetGeo.cxx.

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

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

inline

Definition at line 38 of file AuxDetGeo.h.

38 { return fHalfHeight; }

gar::geo::AuxDetGeo::fHalfHeight

double fHalfHeight

half height of volume

Definition: AuxDetGeo.h:72

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

inline

Definition at line 36 of file AuxDetGeo.h.

36 { return fHalfWidth1; }

gar::geo::AuxDetGeo::fHalfWidth1

double fHalfWidth1

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

Definition: AuxDetGeo.h:70

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

inline

Definition at line 37 of file AuxDetGeo.h.

37 { return fHalfWidth2; }

gar::geo::AuxDetGeo::fHalfWidth2

double fHalfWidth2

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

Definition: AuxDetGeo.h:71

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

inline

Definition at line 35 of file AuxDetGeo.h.

35 { return fLength; }

gar::geo::AuxDetGeo::fLength

double fLength

length of volume, along z direction in local

Definition: AuxDetGeo.h:69

void gar::geo::AuxDetGeo::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 127 of file AuxDetGeo.cxx.

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

void gar::geo::AuxDetGeo::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 136 of file AuxDetGeo.cxx.

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

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

private

Definition at line 118 of file AuxDetGeo.cxx.

     {
       fSensitive.push_back(new AuxDetSensitiveGeo(path, depth));
     }

std::string gar::geo::AuxDetGeo::Name ( ) const

inline

Definition at line 48 of file AuxDetGeo.h.

48 { return fTotalVolume->GetName(); }

gar::geo::AuxDetGeo::fTotalVolume

const TGeoVolume * fTotalVolume

Total volume of AuxDet, called vol*.

Definition: AuxDetGeo.h:67

size_t gar::geo::AuxDetGeo::NSensitiveVolume ( ) const

inline

Definition at line 56 of file AuxDetGeo.h.

56 { return fSensitive.size(); }

gar::geo::AuxDetGeo::fSensitive

std::vector< AuxDetSensitiveGeo * > fSensitive

sensitive volumes in the detector

Definition: AuxDetGeo.h:73

AuxDetSensitiveGeo const & gar::geo::AuxDetGeo::PositionToSensitiveVolume	(	double const	worldLoc[3],
		size_t &	sv
	)		const

Definition at line 219 of file AuxDetGeo.cxx.

     {
       sv = this->FindSensitiveVolume(worldLoc);
       return this->SensitiveVolume(sv);
     }

AuxDetSensitiveGeo const& gar::geo::AuxDetGeo::SensitiveVolume ( size_t sv ) const

inline

Definition at line 55 of file AuxDetGeo.h.

55 { return *fSensitive[sv]; }

gar::geo::AuxDetGeo::fSensitive

std::vector< AuxDetSensitiveGeo * > fSensitive

sensitive volumes in the detector

Definition: AuxDetGeo.h:73

void gar::geo::AuxDetGeo::SortSubVolumes ( geo::AuxDetGeoObjectSorter const & sorter )

Definition at line 227 of file AuxDetGeo.cxx.

     {
       sorter.SortAuxDetSensitive(fSensitive);
       
       return;
     }

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

inline

Definition at line 39 of file AuxDetGeo.h.

39 { return fTotalVolume; }

gar::geo::AuxDetGeo::fTotalVolume

const TGeoVolume * fTotalVolume

Total volume of AuxDet, called vol*.

Definition: AuxDetGeo.h:67

void gar::geo::AuxDetGeo::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 145 of file AuxDetGeo.cxx.

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

void gar::geo::AuxDetGeo::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 154 of file AuxDetGeo.cxx.

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

Member Data Documentation

TGeoHMatrix* gar::geo::AuxDetGeo::fGeoMatrix

private

Transformation matrix to world frame.

Definition at line 68 of file AuxDetGeo.h.

double gar::geo::AuxDetGeo::fHalfHeight

private

half height of volume

Definition at line 72 of file AuxDetGeo.h.

double gar::geo::AuxDetGeo::fHalfWidth1

private

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

Definition at line 70 of file AuxDetGeo.h.

double gar::geo::AuxDetGeo::fHalfWidth2

private

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

Definition at line 71 of file AuxDetGeo.h.

double gar::geo::AuxDetGeo::fLength

private

length of volume, along z direction in local

Definition at line 69 of file AuxDetGeo.h.

std::vector<AuxDetSensitiveGeo*> gar::geo::AuxDetGeo::fSensitive

private

sensitive volumes in the detector

Definition at line 73 of file AuxDetGeo.h.

const TGeoVolume* gar::geo::AuxDetGeo::fTotalVolume

private

Total volume of AuxDet, called vol*.

Definition at line 67 of file AuxDetGeo.h.

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

garsoft/Geometry/AuxDetGeo.h
garsoft/Geometry/AuxDetGeo.cxx

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation