gar::geo::AuxDetGeometryCore Class Reference

Description of geometry of one set of auxiliary detectors. More...

#include <AuxDetGeometryCore.h>

Public Types
using	AuxDetList_t = AuxDetGeometryData_t::AuxDetList_t
	Type of list of auxiliary detectors. More...

Public Member Functions
	AuxDetGeometryCore (fhicl::ParameterSet const &pset)
	Initialize geometry from a given configuration. More...
	~AuxDetGeometryCore ()
	Destructor. More...
	AuxDetGeometryCore (AuxDetGeometryCore const &)=delete
	AuxDetGeometryCore (AuxDetGeometryCore &&)=delete
AuxDetGeometryCore &	operator= (AuxDetGeometryCore const &)=delete
AuxDetGeometryCore &	operator= (AuxDetGeometryCore &&)=delete
std::string	ROOTFile () const
	Returns the full directory path to the geometry file source. More...
std::string	GDMLFile () const
	Returns the full directory path to the GDML file source. More...
std::string	DetectorName () const
	Returns a string with the name of the detector, as configured. More...
unsigned int	NAuxDets () const
	Returns the number of auxiliary detectors. More...
unsigned int	NAuxDetSensitive (size_t const &aid) const
	Returns the number of sensitive components of auxiliary detector. More...
std::vector< AuxDetGeo * > const &	AuxDetGeoVec () const
	Returns the full list of pointer to the auxiliary detectors. More...
AuxDetGeo const &	AuxDet (unsigned int const ad=0) const
	Returns the specified auxiliary detector. More...
unsigned int	FindAuxDetAtPosition (double const worldLoc[3]) const
	Returns the index of the auxiliary detector at specified location. More...
void	FindAuxDetSensitiveAtPosition (double const worldLoc[3], size_t &adg, size_t &sv) const
	Fills the indices of the sensitive auxiliary detector at location. More...
AuxDetGeo const &	PositionToAuxDet (double const worldLoc[3], unsigned int &ad) const
	Returns the auxiliary detector at specified location. More...
const AuxDetSensitiveGeo &	PositionToAuxDetSensitive (double const worldLoc[3], size_t &ad, size_t &sv) const
	Returns the auxiliary detector at specified location. More...
uint32_t	PositionToAuxDetChannel (double const worldLoc[3], size_t &ad, size_t &sv) const
const TVector3	AuxDetChannelToPosition (uint32_t const &channel, std::string const &auxDetName) const
const AuxDetGeo &	ChannelToAuxDet (std::string const &auxDetName, uint32_t const &channel) const
const AuxDetSensitiveGeo &	ChannelToAuxDetSensitive (std::string const &auxDetName, uint32_t const &channel) const
Geometry initialization
void	LoadGeometryFile (std::string gdmlfile, std::string rootfile)
	Loads the geometry information from the specified files. More...
bool	hasAuxDetChannelMap () const
	Returns whether we have a channel map. More...
void	ApplyChannelMap (std::shared_ptr< geo::seg::AuxDetChannelMapAlg > pChannelMap)
	Initializes the geometry to work with this channel map. More...

Protected Member Functions
geo::seg::AuxDetChannelMapAlg const *	AuxDetChannelMap () const
	Returns the object handling the channel map. More...
AuxDetList_t &	AuxDets ()
	Return the internal auxiliary detectors list. More...
AuxDetList_t const &	AuxDets () const

Private Member Functions
void	FindAuxDet (std::vector< const TGeoNode * > &path, unsigned int depth)
void	MakeAuxDet (std::vector< const TGeoNode * > &path, int depth)
void	ClearGeometry ()
	Deletes the detector geometry structures. More...

Private Attributes
AuxDetGeometryData_t	fGeoData
	The detector description data. More...
std::string	fDetectorName
	Name of the detector. More...
std::string	fGDMLfile
	path to geometry file used for Geant4 simulation More...
std::string	fROOTfile
	path to geometry file for geometry in GeometryCore More...
std::shared_ptr< const geo::seg::AuxDetChannelMapAlg >	fChannelMapAlg
	Object containing the channel to wire mapping. More...

Detailed Description

Description of geometry of one set of auxiliary detectors.

---

Note

All lengths are specified in centimetres

How to correctly instantiate a GeometryCore object

Instantiation is a multi-step procedure:

1. construct a GeometryCore object (the "service provider"), with the full configuration; at this step, configuration is just stored
2. load a geometry with GeometryCore::LoadGeometryFile(); this loads the detector geometry information
3. prepare a channel map algorithm object (might use for example GeometryCore::DetectorName() or the detector geometry from the newly created object, but any use of channel mapping related functions is forbidden and it would yield undefined behaviour (expected to be catastrophic)
4. acquire the channel mapping algorithm with GeometryCore::ApplyChannelMap(); at this point, the ChannelMapAlg object is asked to initialize itself and to perform whatever modifications to the geometry provider is needed.

Step 3 (creation of the channel mapping algorithm object) can be performed at any time before step 4, provided that no GeometryCore instance is needed for it.

Configuration parameters

• Name (string; mandatory): string identifying the detector; it can be different from the base name of the file used to initialize the geometry; standard names are recommended by each experiment. This name can be used, for example, to select which channel mapping algorithm to use.
• SurfaceY (real; mandatory): depth of the detector, in centimetrs; see SurfaceY() for details
• MinWireZDist (real; default: 3)
• PositionEpsilon (real; default: 0.01%) set the default tolerance (see DefaultWiggle())

Definition at line 98 of file AuxDetGeometryCore.h.

Member Typedef Documentation

using gar::geo::AuxDetGeometryCore::AuxDetList_t = AuxDetGeometryData_t::AuxDetList_t

Type of list of auxiliary detectors.

Definition at line 102 of file AuxDetGeometryCore.h.

Constructor & Destructor Documentation

gar::geo::AuxDetGeometryCore::AuxDetGeometryCore

(

fhicl::ParameterSet const &

pset

)

Initialize geometry from a given configuration.

Parameters

pset	configuration parameters

This constructor does not load any geometry description. The next step is to do exactly that, by GeometryCore::LoadGeometryFile().

Definition at line 42 of file AuxDetGeometryCore.cxx.

    : fDetectorName(pset.get< std::string >("Name"))
    {
      std::transform(fDetectorName.begin(), fDetectorName.end(), fDetectorName.begin(), ::tolower);
    } // AuxDetGeometryCore::AuxDetGeometryCore()

gar::geo::AuxDetGeometryCore::~AuxDetGeometryCore

(

)

Destructor.

Definition at line 50 of file AuxDetGeometryCore.cxx.

    {
      ClearGeometry();
    } // AuxDetGeometryCore::~AuxDetGeometryCore()

gar::geo::AuxDetGeometryCore::AuxDetGeometryCore ( AuxDetGeometryCore const &  )

delete

gar::geo::AuxDetGeometryCore::AuxDetGeometryCore ( AuxDetGeometryCore &&  )

delete

Member Function Documentation

void gar::geo::AuxDetGeometryCore::ApplyChannelMap

(

std::shared_ptr< geo::seg::AuxDetChannelMapAlg >

pChannelMap

)

Initializes the geometry to work with this channel map.

Parameters

pChannelMap

a pointer to the channel mapping algorithm to be used

See also

LoadGeometryFile()

The specified channel mapping is used with this geometry. The algorithm object is asked and allowed to make the necessary modifications to the geometry description. These modifications typically involve some resorting of the objects.

The ownership of the algorithm object is shared, usually with a calling framework: we maintain it alive as long as we need it (and no other code can delete it), and we delete it only if no other code is sharing the ownership.

This method needs to be called after LoadGeometryFile() to complete the geometry initialization.

Definition at line 57 of file AuxDetGeometryCore.cxx.

    {
      pChannelMap->Initialize(fGeoData);
      fChannelMapAlg = pChannelMap;
    } // AuxDetGeometryCore::ApplyChannelMap()

const AuxDetGeo & gar::geo::AuxDetGeometryCore::AuxDet

(

unsigned int const

ad = 0

)

const

Returns the specified auxiliary detector.

Parameters

ad	the auxiliary detector index

Returns

a constant reference to the specified auxiliary detector

Todo:

what happens if it does not exist?

remove the default parameter?

Definition at line 134 of file AuxDetGeometryCore.cxx.

    {
      if(ad >= NAuxDets())
        throw cet::exception("AuxDetGeometryCore") << "AuxDet "
        << ad
        << " does not exist\n";

      return *(AuxDets()[ad]);
    }

geo::seg::AuxDetChannelMapAlg const* gar::geo::AuxDetGeometryCore::AuxDetChannelMap ( ) const

inlineprotected

Returns the object handling the channel map.

Definition at line 305 of file AuxDetGeometryCore.h.

{ return fChannelMapAlg.get(); }

const TVector3 gar::geo::AuxDetGeometryCore::AuxDetChannelToPosition	(	uint32_t const &	channel,
		std::string const &	auxDetName
	)		const

Definition at line 191 of file AuxDetGeometryCore.cxx.

    {
      return fChannelMapAlg->AuxDetChannelToPosition(channel, auxDetName, AuxDets());
    }

std::vector<AuxDetGeo*> const& gar::geo::AuxDetGeometryCore::AuxDetGeoVec ( ) const

inline

Returns the full list of pointer to the auxiliary detectors.

Definition at line 178 of file AuxDetGeometryCore.h.

{ return AuxDets(); }

AuxDetList_t& gar::geo::AuxDetGeometryCore::AuxDets ( )

inlineprotected

Return the internal auxiliary detectors list.

Definition at line 309 of file AuxDetGeometryCore.h.

{ return fGeoData.auxDets; }

AuxDetList_t const& gar::geo::AuxDetGeometryCore::AuxDets ( ) const

inlineprotected

Definition at line 310 of file AuxDetGeometryCore.h.

{ return fGeoData.auxDets; }

const AuxDetGeo & gar::geo::AuxDetGeometryCore::ChannelToAuxDet	(	std::string const &	auxDetName,
		uint32_t const &	channel
	)		const

Definition at line 198 of file AuxDetGeometryCore.cxx.

    {
      size_t adIdx = fChannelMapAlg->ChannelToAuxDet(AuxDets(), auxDetName, channel);
      return this->AuxDet(adIdx);
    }

const AuxDetSensitiveGeo & gar::geo::AuxDetGeometryCore::ChannelToAuxDetSensitive	(	std::string const &	auxDetName,
		uint32_t const &	channel
	)		const

Definition at line 206 of file AuxDetGeometryCore.cxx.

    {
      auto idx = fChannelMapAlg->ChannelToSensitiveAuxDet(AuxDets(), auxDetName, channel);
      return this->AuxDet(idx.first).SensitiveVolume(idx.second);
    }

void gar::geo::AuxDetGeometryCore::ClearGeometry ( )

private

Deletes the detector geometry structures.

Definition at line 103 of file AuxDetGeometryCore.cxx.

    {
      // auxiliary detectors
      std::for_each(AuxDets().begin(), AuxDets().end(), std::default_delete<AuxDetGeo>());
      AuxDets().clear();

    } // AuxDetGeometryCore::ClearGeometry()

std::string gar::geo::AuxDetGeometryCore::DetectorName ( ) const

inline

Returns a string with the name of the detector, as configured.

Definition at line 143 of file AuxDetGeometryCore.h.

{ return fDetectorName; }

void gar::geo::AuxDetGeometryCore::FindAuxDet ( std::vector< const TGeoNode * > &  path, unsigned int  depth  )

private

Definition at line 214 of file AuxDetGeometryCore.cxx.

    {
      const char* nm = path[depth]->GetName();
      if( (strncmp(nm, "volAuxDet", 9) == 0) ){
        this->MakeAuxDet(path, depth);
        return;
      }

      //explore the next layer down
      unsigned int deeper = depth+1;
      if(deeper >= path.size()){
        throw cet::exception("AuxDetGeometryCore")
        << "exceeded maximum TGeoNode 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->FindAuxDet(path, deeper);
      }

    }

unsigned int gar::geo::AuxDetGeometryCore::FindAuxDetAtPosition

(

double const

worldLoc[3]

)

const

Returns the index of the auxiliary detector at specified location.

Parameters

worldLoc

3D coordinates of the point (world reference frame)

Returns

the index of the detector, or UINT_MAX if no detector is there

Todo:

replace with numeric_limits<>?

Definition at line 146 of file AuxDetGeometryCore.cxx.

    {
      return fChannelMapAlg->NearestAuxDet(worldPos, AuxDets());
    } // AuxDetGeometryCore::FindAuxDetAtPosition()

void gar::geo::AuxDetGeometryCore::FindAuxDetSensitiveAtPosition	(	double const	worldLoc[3],
		size_t &	adg,
		size_t &	sv
	)		const

Fills the indices of the sensitive auxiliary detector at location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
adg	(output) auxiliary detector index
sv	(output) sensitive volume index

Definition at line 162 of file AuxDetGeometryCore.cxx.

    {
      adg = this->FindAuxDetAtPosition(worldPos);
      sv  = fChannelMapAlg->NearestSensitiveAuxDet(worldPos, AuxDets(), adg);

      return;
    } // AuxDetGeometryCore::FindAuxDetAtPosition()

std::string gar::geo::AuxDetGeometryCore::GDMLFile ( ) const

inline

Returns the full directory path to the GDML file source.

Returns

the full directory path to the GDML file source

This is the full path of the source of the detector geometry handed to the detector simulation (GEANT).

Definition at line 139 of file AuxDetGeometryCore.h.

{ return fGDMLfile; }

bool gar::geo::AuxDetGeometryCore::hasAuxDetChannelMap ( ) const

inline

Returns whether we have a channel map.

Definition at line 277 of file AuxDetGeometryCore.h.

{ return bool(fChannelMapAlg); }

void gar::geo::AuxDetGeometryCore::LoadGeometryFile	(	std::string	gdmlfile,
		std::string	rootfile
	)

Loads the geometry information from the specified files.

Parameters

gdmlfile	path to file to be used for Geant4 simulation
rootfile	path to file for internal geometry representation

See also

ApplyChannelMap()

Both paths must directly resolve to an available file, as no search is performed for them.

The gdmlfile parameter does not have to necessarily be in GDML format, as long as it's something supported by Geant4. This file is not used by the geometry, but its path is provided on request by the simulation modules (see GArSoft LArG4 module). The rootfile also does not need to be a ROOT file, but just anything that TGeoManager::Import() supports. This file is parsed immediately and the internal geometry representation is built out of it.

Note

After calling this method, the detector geometry information can be considered complete, but the geometry service provider is not fully initialized yet, since it's still necessary to provide or update the channel mapping.

Definition at line 64 of file AuxDetGeometryCore.cxx.

    {

      if (gdmlfile.empty()) {
        throw cet::exception("AuxDetGeometryCore")
        << "No GDML Geometry file specified!\n";
      }

      if (rootfile.empty()) {
        throw cet::exception("AuxDetGeometryCore")
        << "No ROOT Geometry file specified!\n";
      }

      ClearGeometry();

      // Open the GDML file, and convert it into ROOT TGeoManager format.
      // try to be efficient - if the GeometryCore object already imported
      // the file, then the gGeoManager will be non-null.  If not, import it.
      // Then lock the gGeoManager to prevent future imports.
      if( !gGeoManager ){
        TGeoManager::Import(rootfile.c_str());
        gGeoManager->LockGeometry();
      }

      std::vector<const TGeoNode*> path(8);
      path[0] = gGeoManager->GetTopNode();
      FindAuxDet(path, 0);

      fGDMLfile = gdmlfile;
      fROOTfile = rootfile;

      MF_LOG_INFO("AuxDetGeometryCore")
      << "New detector geometry loaded from "
      << "\n\t" << fROOTfile
      << "\n\t" << fGDMLfile << "\n";

    } // AuxDetGeometryCore::LoadGeometryFile()

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

private

Definition at line 240 of file AuxDetGeometryCore.cxx.

    {
      AuxDets().push_back(new AuxDetGeo(path, depth));
    }

unsigned int gar::geo::AuxDetGeometryCore::NAuxDets ( ) const

inline

Returns the number of auxiliary detectors.

Todo:

use a AutDetID_t instead of unsigned int?

This method returns the total number of scintillator paddles (Auxiliary Detectors aka AuxDet) outside of the cryostat

Todo:

Change return type to size_t

Definition at line 163 of file AuxDetGeometryCore.h.

{ return AuxDets().size(); }

unsigned int gar::geo::AuxDetGeometryCore::NAuxDetSensitive

(

size_t const &

aid

)

const

Returns the number of sensitive components of auxiliary detector.

Parameters

aid	ID of the auxiliary detector

Returns

number of sensitive components in the auxiliary detector aid cet::exception (category "Geometry") if aid does not exist

Definition at line 113 of file AuxDetGeometryCore.cxx.

    {
      if( aid > NAuxDets() - 1)
        throw cet::exception("Geometry")
        << "Requested AuxDet index "
        << aid
        << " is out of range: "
        << NAuxDets();

      return AuxDets()[aid]->NSensitiveVolume();
    }

AuxDetGeometryCore& gar::geo::AuxDetGeometryCore::operator= ( AuxDetGeometryCore const &  )

delete

AuxDetGeometryCore& gar::geo::AuxDetGeometryCore::operator= ( AuxDetGeometryCore &&  )

delete

const AuxDetGeo & gar::geo::AuxDetGeometryCore::PositionToAuxDet	(	double const	worldLoc[3],
		unsigned int &	ad
	)		const

Returns the auxiliary detector at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
ad	(output) the auxiliary detector index

Returns

constant reference to AuxDetGeo object of the auxiliary detector

Todo:

what happens if it does not exist?

Definition at line 152 of file AuxDetGeometryCore.cxx.

    {
      // locate the desired Auxiliary Detector
      ad = this->FindAuxDetAtPosition(worldLoc);

      return this->AuxDet(ad);
    }

uint32_t gar::geo::AuxDetGeometryCore::PositionToAuxDetChannel	(	double const	worldLoc[3],
		size_t &	ad,
		size_t &	sv
	)		const

Definition at line 183 of file AuxDetGeometryCore.cxx.

    {
      return fChannelMapAlg->PositionToAuxDetChannel(worldLoc, AuxDets(), ad, sv);
    }

const AuxDetSensitiveGeo & gar::geo::AuxDetGeometryCore::PositionToAuxDetSensitive	(	double const	worldLoc[3],
		size_t &	ad,
		size_t &	sv
	)		const

Returns the auxiliary detector at specified location.

Parameters

worldLoc	3D coordinates of the point (world reference frame)
ad	(output) the auxiliary detector index
sv	(output) the auxiliary detector sensitive volume index

Returns

reference to AuxDetSensitiveGeo object of the auxiliary detector

Todo:

what happens if it does not exist?

Definition at line 173 of file AuxDetGeometryCore.cxx.

    {
      // locate the desired Auxiliary Detector
      this->FindAuxDetSensitiveAtPosition(worldLoc, ad, sv);
      return this->AuxDet(ad).SensitiveVolume(sv);
    }

std::string gar::geo::AuxDetGeometryCore::ROOTFile ( ) const

inline

Returns the full directory path to the geometry file source.

Returns

the full directory path to the geometry file source

This is the full path of the source of the detector geometry GeometryCore relies on.

Definition at line 130 of file AuxDetGeometryCore.h.

{ return fROOTfile; }

Member Data Documentation

std::shared_ptr<const geo::seg::AuxDetChannelMapAlg> gar::geo::AuxDetGeometryCore::fChannelMapAlg

private

Object containing the channel to wire mapping.

Definition at line 327 of file AuxDetGeometryCore.h.

std::string gar::geo::AuxDetGeometryCore::fDetectorName

private

Name of the detector.

Definition at line 324 of file AuxDetGeometryCore.h.

std::string gar::geo::AuxDetGeometryCore::fGDMLfile

private

path to geometry file used for Geant4 simulation

Definition at line 325 of file AuxDetGeometryCore.h.

AuxDetGeometryData_t gar::geo::AuxDetGeometryCore::fGeoData

private

The detector description data.

Definition at line 322 of file AuxDetGeometryCore.h.

std::string gar::geo::AuxDetGeometryCore::fROOTfile

private

path to geometry file for geometry in GeometryCore

Definition at line 326 of file AuxDetGeometryCore.h.

---

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

• garsoft/Geometry/AuxDetGeometryCore.h
• garsoft/Geometry/AuxDetGeometryCore.cxx