geo::Geometry Class Reference

The geometry of one entire detector (near, far, ipnd). More...

#include <Geometry.h>

List of all members.

Public Member Functions
const PlaneGeo &	Plane (unsigned int i) const
const std::set< unsigned int > &	GetPlanesByView (View_t v=kXorY) const
const unsigned int	NextPlaneInView (unsigned int p, int d=+1) const
const unsigned int	NextPlaneOtherView (unsigned int p, int d=+1) const
void	CellInfo (unsigned int ip, unsigned int ic, View_t *view=0, double *pos=0, double *dpos=0) const
const CellGeo &	IdToCell (const CellUniqueId &id, int *iplane=0, int *icell=0) const
unsigned int	NPlanes () const
double	DetHalfWidth () const
double	DetHalfHeight () const
double	DetLength () const
double	SurfaceY () const
void	WorldBox (double *xlo, double *xhi, double *ylo, double *yhi, double *zlo, double *zhi) const
int	CurrentCell (int *ip, int *ic) const
const CellUniqueId	CurrentCellId () const
const TGeoMaterial *	Material (double x, double y, double z) const
double	TotalMass () const
double	MassBetweenPoints (double *p1, double *p2) const
TGeoManager *	ROOTGeoManager () const
Static Public Member Functions
Geometry &	Instance (int runnum, const char *fname)
Private Types
typedef std::pair< unsigned short, unsigned short >	PlaneCellPair
typedef std::map< CellUniqueId, PlaneCellPair >	UIDMap
Private Member Functions
	Geometry (const char *gdmlfile)
	~Geometry ()
void	FindPlanes (std::vector< const TGeoNode * > &n, unsigned int depth)
void	MakePlane (std::vector< const TGeoNode * > &n, unsigned int depth)
void	SetDrawOptions ()
void	BuildMaps ()
void	SetDetectorSize ()
Private Attributes
int	fRunNumber
	Run number of configuration.
std::string	fGDMLFile
	gdml file holding the geometry
std::vector< PlaneGeo * >	fPlanes
	The detector planes.
UIDMap	fIdMap
	Unique ID -> Plane,Cell.
std::set< unsigned int >	fAllPlanes
	List of all planes.
std::set< unsigned int >	fXplanes
	List of X measuring planes.
std::set< unsigned int >	fYplanes
	List of Y measuring planes.
double	fDetLength
	Detector 1/2 length (cm).
double	fDetHalfHeight
	Detector 1/2 height (cm).
double	fDetHalfWidth
	Detector 1/2 width (cm).
Static Private Attributes
Geometry *	fInstance = 0
	Instance of geometry.

---

Detailed Description

The geometry of one entire detector (near, far, ipnd).

Definition at line 31 of file Geometry.h.

---

Member Typedef Documentation

typedef std::pair<unsigned short, unsigned short> geo::Geometry::PlaneCellPair [private]

Definition at line 77 of file Geometry.h. Referenced by BuildMaps().

typedef std::map<CellUniqueId, PlaneCellPair> geo::Geometry::UIDMap [private]

Definition at line 78 of file Geometry.h.

---

Constructor & Destructor Documentation

Geometry::Geometry (  const char *  gdmlfile  )  [private]

Definition at line 86 of file Geometry.cxx. References BuildMaps(), FindPlanes(), fPlanes, SetDetectorSize(), and SetDrawOptions(). Referenced by Instance(). { std::string fname; // Test if the gdml file can be openned using exactly the name we've // been passed while (1) { struct stat sb; fname = gdmlfile; if (stat(fname.c_str(), &sb)==0) break; const char* srt_private = getenv("SRT_PRIVATE_CONTEXT"); if (srt_private!=0) { fname = srt_private; fname += "/Geometry/gdml/"; fname += gdmlfile; if (stat(fname.c_str(), &sb)==0) break; } const char* srt_public = getenv("SRT_PUBLIC_CONTEXT"); if (srt_public!=0) { fname = srt_public; fname += "/Geometry/gdml/"; fname += gdmlfile; if (stat(fname.c_str(), &sb)==0) break; } // Failed to resolve the file name throw Exception(__FILE__, __LINE__, Exception::FILE_NOT_FOUND, "GDML file not found."); } // Open the GDML file TGeoManager::Import(fname.c_str()); this->SetDrawOptions(); std::vector<const TGeoNode*> n(16); n[0] = gGeoManager->GetTopNode(); this->FindPlanes(n, 0); sort(fPlanes.begin(), fPlanes.end(), plane_sort); this->SetDetectorSize(); this->BuildMaps(); std::cout << "Detector geometry loaded from " << fname << std::endl; }

Geometry::~Geometry (   )  [private]

Definition at line 142 of file Geometry.cxx. References fPlanes. { for (unsigned int i=0; i<fPlanes.size(); ++i) { if (fPlanes[i]) { delete fPlanes[i]; fPlanes[i] = 0; } } }

---

Member Function Documentation

void Geometry::BuildMaps (   )  [private]

Build maps used for quick access to the geometry Definition at line 463 of file Geometry.cxx. References geo::PlaneGeo::Cell(), fAllPlanes, fIdMap, fPlanes, fXplanes, fYplanes, geo::CellGeo::Id(), geo::PlaneGeo::Ncells(), Plane(), and PlaneCellPair. Referenced by Geometry(). { // Sets which list planes by view for (unsigned int i=0; i<fPlanes.size(); ++i) { fAllPlanes.insert(i); if (fPlanes[i]->View()==kX) fXplanes.insert(i); else if (fPlanes[i]->View()==kY) fYplanes.insert(i); else { throw Exception(__FILE__,__LINE__,Exception::BAD_GEO_CONFIG, "Bad plane view!"); } } // Unique Id -> cell map for (unsigned int i=0; i<fPlanes.size(); ++i) { const PlaneGeo* p = fPlanes[i]; for (unsigned int j=0; j<p->Ncells(); ++j) { // Assert that the cell ids are in fact unique if (fIdMap.find(this->Plane(i).Cell(j).Id())!=fIdMap.end()) { throw Exception(__FILE__,__LINE__,Exception::BAD_CELL_UNIQUE_ID, "Duplicate Cell UIDs"); } PlaneCellPair p(i,j); fIdMap[this->Plane(i).Cell(j).Id()] = p; } } }

void Geometry::CellInfo (  unsigned int  ip, unsigned int  ic, View_t *  view = 0, double *  pos = 0, double *  dpos = 0 )  const

Given a plane and cell number, return which view the cell measures, its nominal center position, and nominal size about that center position. Parameters: ip  - plane number ic  - cell number view  - on retrun, the view this cell measures in pos  - the nominal (xyz) center position of this cell (cm) dpos  - half-dimensions of the cell in xyz (cm) Exceptions: geo::Exception  Definition at line 339 of file Geometry.cxx. References geo::PlaneGeo::Cell(), geo::CellGeo::GetCenter(), geo::CellGeo::HalfD(), geo::CellGeo::HalfL(), geo::CellGeo::HalfW(), geo::PlaneGeo::Ncells(), NPlanes(), Plane(), and geo::PlaneGeo::View(). Referenced by evd::Display3D::DrawRawData(), and evd::RecoBaseDrawer::GetClusterOutlines(). { if (ip<0 || ip>=this->NPlanes()) { throw Exception(__FILE__,__LINE__,Exception::BAD_PLANE_NUMBER); } const PlaneGeo& p = this->Plane(ip); if (view) *view = p.View(); if (ic<0 || ic>=p.Ncells()) { throw Exception(__FILE__,__LINE__,Exception::BAD_CELL_NUMBER); } const::CellGeo& c = p.Cell(ic); if (pos) c.GetCenter(pos); if (dpos) { if (p.View() == kX) { dpos[0] = c.HalfW(); dpos[1] = c.HalfL(); dpos[2] = c.HalfD(); } else if (p.View() == kY) { dpos[0] = c.HalfL(); dpos[1] = c.HalfW(); dpos[2] = c.HalfD(); } } }

int Geometry::CurrentCell (  int *  planeid, int *  cellid )  const

Return the plane and cell number for the current tracking position held by the gGeoManager Parameters: planeid  : pointer to plane number, >=0 on success, -1 on failure cellid  : pointer to cell number, >=0 on success, -1 on failure Returns: >0 if cell is found, <0 if not found. Definition at line 158 of file Geometry.cxx. References IdToCell(). Referenced by testFindCell(). { try { this->IdToCell(this->CurrentCellId(), planeid, cellid); } catch (Exception e) { *planeid = -1; *cellid = -1; return -1; } return 1; }

const CellUniqueId Geometry::CurrentCellId (   )  const

Return the unique cell Id# (well, #'s) for the current cell Definition at line 175 of file Geometry.cxx. References geo::CellUniqueId, and geo::PathToUniqueId(). { return PathToUniqueId(gGeoManager->GetPath()); }

double Geometry::DetHalfHeight (   )  const

Definition at line 197 of file Geometry.cxx. Referenced by mcchk::NeutrinoAna::Ana(), evd::DetectorView::DetectorView(), evd::Display3D::Draw(), main(), and evd::TZProjPad::TZProjPad(). { return fDetHalfHeight; }

double Geometry::DetHalfWidth (   )  const

Definition at line 190 of file Geometry.cxx. Referenced by mcchk::NeutrinoAna::Ana(), evd::DetectorView::DetectorView(), evd::Display3D::Draw(), main(), and evd::TZProjPad::TZProjPad(). { return fDetHalfWidth; }

double Geometry::DetLength (   )  const

Definition at line 204 of file Geometry.cxx. Referenced by mcchk::NeutrinoAna::Ana(), mcchk::CosmicAna::Ana(), clusterss::CompareClusters::Ana(), clusterss::ClusterCheck::Ana(), evd::DetectorView::DetectorView(), evd::Display3D::Draw(), and evd::TZProjPad::TZProjPad(). { return fDetLength; }

void Geometry::FindPlanes (  std::vector< const TGeoNode * > &  n, unsigned int  depth )  [private]

Recursively search for planes in the geometry Parameters: n  - Array holding the path of nodes tranversed depth  - Depth of the search to this point Exceptions: geo::Exception  is maximum depth is exceeded. Definition at line 278 of file Geometry.cxx. References MakePlane(). Referenced by Geometry(). { const char* nm = n[depth]->GetName(); if (strncmp(nm,"vPlane", 6)==0) { this->MakePlane(n, depth); return; } // Explore the next layer down unsigned int deeper = depth+1; if (deeper>=n.size()) { throw Exception(__FILE__,__LINE__,Exception::BAD_NODE, "Exceeded maximum depth."); } const TGeoVolume* v = n[depth]->GetVolume(); int nd = v->GetNdaughters(); for (int i=0; i<nd; ++i) { n[deeper] = v->GetNode(i); this->FindPlanes(n, deeper); } }

const std::set< unsigned int > & Geometry::GetPlanesByView (  View_t  v = kXorY  )  const

Return list of planes which measure the requested projection Parameters: v  : X/Y or All (see View_t typedef in PlaneGeo.h) Returns: : a set of plane indicies satifying the request Exceptions: :  geo::Exception if initialization fails Definition at line 400 of file Geometry.cxx. Referenced by mcchk::DetAna::Ana(), mcchk::CosmicAna::Ana(), vali::Validator::BookPT(), clusterss::MakeClusterSS::ClusterWindows(), clusterss::MakeClusterSS::FindLongWindow(), subshower::RecoSubShower::FindMaxWindow(), clusterss::MakeClusterSS::FindPlaneFromShwMax(), and testFindPlanes(). { // Return the user's choice, fall back on all planes if (v==kX) return fXplanes; else if (v==kY) return fYplanes; return fAllPlanes; }

const CellGeo & Geometry::IdToCell (  const CellUniqueId &  id, int *  iplane = 0, int *  icell = 0 )  const

Given a unique cell identifier, look up the plane and cell numbers. Parameters: id  : Cell identifier iplane  : optional, pointer to return plane number icell  : optional, pointer to return cell number Returns: cell geometry description Exceptions: geo::Exception  if id is not found Definition at line 380 of file Geometry.cxx. References geo::PlaneGeo::Cell(), fIdMap, and Plane(). Referenced by CurrentCell(), novamc::MCApplication::Stepping(), and testUniqueId(). { UIDMap::const_iterator itr = fIdMap.find(id); if (itr == fIdMap.end()) { throw Exception(__FILE__,__LINE__,Exception::BAD_CELL_UNIQUE_ID); } *iplane = (int)itr->second.first; *icell = (int)itr->second.second; return this->Plane(*iplane).Cell(*icell); }

Geometry & Geometry::Instance (  int  runnum, const char *  gdmlfile )  [static]

Build a detector geometry Parameters: runnum  : Run number to configure geometry for detId  : Detector ID (0=near, 1=far, 2=ipnd) Exceptions: geo::Exception  if requested detector/run combination cannot be satisfied Definition at line 49 of file Geometry.cxx. References fGDMLFile, fInstance, fRunNumber, and Geometry(). Referenced by clusterss::CompareClusters::Ana(), clusterss::ClusterCheck::Ana(), util::EDMUtils::GetGeometry(), main(), clusterss::MakeClusterSS::Reco(), evgen::CosmicsGen::Reco(), and cellhitmerge::CellHitMerge::Reco(). { // If our current geometry doesn't match request, get rid of it if (fInstance!=0) { if (fInstance->fRunNumber != runnum || fInstance->fGDMLFile.compare(gdmlfile) != 0) { delete fInstance; fInstance = 0; gGeoManager->Clear(); } } // Build a new geometry if we need to if (fInstance==0) { std::string gdml; if (gdmlfile!=0) { // User has specficied the file name to load gdml = gdmlfile; } else { char msg[256]; sprintf(msg,"no geometry file specified"); throw Exception(__FILE__,__LINE__,Exception::BAD_GEO_CONFIG,msg); } if (runnum<0) { char msg[256]; std::sprintf(msg,"Bad run number: %d",runnum); throw Exception(__FILE__,__LINE__,Exception::BAD_GEO_CONFIG,msg); } fInstance = new Geometry(gdml.c_str()); fInstance->fGDMLFile = gdml; fInstance->fRunNumber = runnum; }//end if no instance return *fInstance; }

void Geometry::MakePlane (  std::vector< const TGeoNode * > &  n, unsigned int  depth )  [private]

Definition at line 303 of file Geometry.cxx. References fPlanes. Referenced by FindPlanes(). { fPlanes.push_back(new PlaneGeo(n, depth)); }

double Geometry::MassBetweenPoints (  double *  p1, double *  p2 )  const

Return the column density between 2 points in gm/cm^2 Parameters: p1  : pointer to array holding xyz of first point in world coordinates p2  : pointer to array holding xyz of second point in world coorinates Definition at line 531 of file Geometry.cxx. { // The purpose of this method is to determine the column density // between the two points given. Do that by starting at p1 and // stepping until you get to the node of p2. calculate the distance // between the point just inside that node and p2 to get the last // bit of column density double columnD = 0.; // first initialize a track - get the direction cosines double length = TMath::Sqrt(TMath::Power(p2[0]-p1[0], 2.) + TMath::Power(p2[1]-p1[1], 2.) + TMath::Power(p2[2]-p1[2], 2.)); double dxyz[3] = {(p2[0]-p1[0])/length, (p2[1]-p1[1])/length, (p2[2]-p1[2])/length}; gGeoManager->InitTrack(p1,dxyz); // might be helpful to have a point to a TGeoNode TGeoNode *node = gGeoManager->GetCurrentNode(); // check that the points are not in the same volume already. // if they are in different volumes, keep stepping until you // are in the same volume as the second point while(!gGeoManager->IsSameLocation(p2[0], p2[1], p2[2])){ gGeoManager->FindNextBoundary(); columnD += gGeoManager->GetStep() * node->GetMedium()->GetMaterial()->GetDensity(); // the act of stepping puts you in the next node and returns that // node node = gGeoManager->Step(); } //end loop to get to volume of second point // now you are in the same volume as the last point, but not at that // point. get the distance between the current point and the last // one const double *current = gGeoManager->GetCurrentPoint(); length = TMath::Sqrt(TMath::Power(p2[0]-current[0], 2.) + TMath::Power(p2[1]-current[1], 2.) + TMath::Power(p2[2]-current[2], 2.)); columnD += length*node->GetMedium()->GetMaterial()->GetDensity(); return columnD; }

const TGeoMaterial * Geometry::Material (  double  x, double  y, double  z )  const

Definition at line 182 of file Geometry.cxx. { const TGeoNode* node = gGeoManager->FindNode(x,y,z); return node->GetMedium()->GetMaterial(); }

const unsigned int Geometry::NextPlaneInView (  unsigned int  p1, int  d = +1 )  const

Return the index of the next plane in a particular view Parameters: p1  : Index of current plane d  : Direction to step in (>0 next down stream, <0 upstream) Returns: Index of next plane in same view. If no plane can be found (eg. reach end of detector) return kPLANE_NOT_FOUND. Exceptions: :  Asserts that p1 is in valid range Definition at line 419 of file Geometry.cxx. References fPlanes. Referenced by clusterss::MakeClusterSS::ClusterWindows(), subshower::RecoSubShower::LinkedBlobs(), subshower::RecoSubShower::NumXYPlanes(), calhit::CalHit::RoughReco(), clusterss::MakeClusterSS::StupidHoughNoWork(), testFindPlanes(), and subshower::RecoSubShower::TransCluster(). { assert(p1<fPlanes.size()); // March along from p1 in the direction d until we find a match int s = (d>0 ? 1 : -1); if (s<0 && p1==0) return kPLANE_NOT_FOUND; if (s>0 && p1==fPlanes.size()-1) return kPLANE_NOT_FOUND; for (unsigned int p2=p1+d; 1; p2+=s) { if (fPlanes[p1]->View() == fPlanes[p2]->View()) return p2; if (p2==fPlanes.size()-1) return kPLANE_NOT_FOUND; if (p2==0) return kPLANE_NOT_FOUND; } }

const unsigned int Geometry::NextPlaneOtherView (  unsigned int  p1, int  d = +1 )  const

Return the index of the next plane in a particular view Parameters: p1  : Index of current plane d  : Direction to step in (>0 next down stream, <0 upstream) Returns: index of next plane in same view. Exceptions: :  Asserts that p1 is in valid range Definition at line 444 of file Geometry.cxx. References fPlanes. Referenced by calhit::CalHit::RoughReco(), and testFindPlanes(). { assert(p1<fPlanes.size()); // March along from p1 in the direction d until we find a match int s = (d>0 ? 1 : -1); if (s<0 && p1==0) return kPLANE_NOT_FOUND; if (s>0 && p1==fPlanes.size()-1) return kPLANE_NOT_FOUND; for (unsigned int p2=p1+d; 1; p2+=s) { if (fPlanes[p1]->View() != fPlanes[p2]->View()) return p2; if (p2==fPlanes.size()) return kPLANE_NOT_FOUND; if (p2==0) return kPLANE_NOT_FOUND; } }

unsigned int geo::Geometry::NPlanes (   )  const [inline]

Definition at line 49 of file Geometry.h. Referenced by CellInfo(), calib::Calibrator::GetXYZD(), and cellhitmerge::CellHitMerge::Reco(). { return fPlanes.size(); }

const PlaneGeo & Geometry::Plane (  unsigned int  i  )  const

Return the geometry description of the ith plane in the detector. Parameters: i  : input plane number, starting from 0 Returns: plane geometry for ith plane Exceptions: geo::Exception  if iplane is outside allowed range Definition at line 318 of file Geometry.cxx. References fPlanes. Referenced by spider::Filament::AddHit(), mcchk::DetAna::Ana(), mcchk::CosmicAna::Ana(), subshower::RecoSubShower::BestHough(), clusterss::MakeClusterSS::BestHough(), vali::Validator::BookPT(), BuildMaps(), rpr::TrackReco::CalcChi2Ndof(), subshower::RecoSubShower::CalculateEnergyVertexAngle(), evd::RecoBaseDrawer::CellHit2D(), CellInfo(), clusterss::MakeClusterSS::ClusterWindows(), rpr::TrackReco::Connection(), evd::Display3D::DrawMCHits(), evd::Display3D::DrawRawData(), clusterss::MakeClusterSS::FindAllTransverseWindows(), spider::Filament::FitDist(), evd::SimulationDrawer::FLSHit2D(), calib::Calibrator::GetXYZD(), subshower::RecoSubShower::HoughTransCluster(), clusterss::MakeClusterSS::HoughTransCluster(), IdToCell(), main(), calib::Calibrator::MakeCellHit(), ctrk::CosmicTrack::MakeClusters(), ctrk::CosmicTrack::MakeTracks(), vali::Validator::PlotCell(), vali::Validator::PlotPT(), raw_digit_sort(), evd::RawDataDrawer::RawDigit2D(), rpr::TrackReco::Reco(), spider::SpiderWeb::Reco(), photrans::SimpleTransport::Reco(), rpr::FindTrackSeg::Reco(), trk::DemoShell::Reco(), cellhitmerge::CellHitMerge::Reco(), calhit::CalHit::RoughReco(), SetDetectorSize(), clusterss::MakeClusterSS::StupidHoughNoWork(), testCellPos(), testFindCell(), testFindPlanes(), testStepping(), testUniqueId(), subshower::RecoSubShower::TransCluster(), spider::Filament::Update(), and cluster::Clusterer::UpdateClusterStat(). { if (i>=fPlanes.size()) throw Exception(__FILE__,__LINE__, Exception::BAD_PLANE_NUMBER, "Plane index out of range"); return *fPlanes[i]; }

TGeoManager * Geometry::ROOTGeoManager (   )  const

Definition at line 135 of file Geometry.cxx. Referenced by photrans::SimpleTransport::Reco(). { return gGeoManager; }

void Geometry::SetDetectorSize (   )  [private]

Calcualte the detector size based on the geomtry. Requires a little bit of work, so do this once when the geometry is updated Definition at line 586 of file Geometry.cxx. References geo::PlaneGeo::Cell(), fDetHalfHeight, fDetHalfWidth, fDetLength, fPlanes, geo::CellGeo::GetCenter(), and Plane(). Referenced by Geometry(). { double dummy; // Compute the height and width based on the sizes of the vertical // and horizontal planes. Remember, in the local plane frame z goes // along the length of the modules. Hence the use of axis=3 to get // the detector height (from the vertical planes) and width (from // the horizontal planes). const TGeoVolume* vvp = gGeoManager->GetVolume("vPlaneV"); const TGeoVolume* hvp = gGeoManager->GetVolume("vPlaneH"); // Handle the case of the IPND where the volumes are names slightly // different accounting for the different widths of the modules. if (vvp==0) vvp = gGeoManager->GetVolume("vPlaneVIPND"); if (hvp==0) hvp = gGeoManager->GetVolume("vPlaneHIPND"); if (vvp==0 && hvp==0) { vvp = gGeoManager->GetVolume("vPlaneVIPND"); hvp = gGeoManager->GetVolume("vPlaneHIPND"); } if (vvp==0 || hvp==0) { throw Exception(__FILE__,__LINE__, Exception::BAD_GEO_CONFIG, "Unable to find shapes to set detector size"); } const TGeoShape* vp = vvp->GetShape(); const TGeoShape* hp = hvp->GetShape(); vp->GetAxisRange(3,dummy,fDetHalfHeight); hp->GetAxisRange(3,dummy,fDetHalfWidth); // Remember -- in the local plane coordinate frame, the "depth" is // along y, hence axis=2. In the world coordinate this is along z. double vplanez1, vplanez2; vp->GetAxisRange(2,vplanez1,vplanez2); // Compute the detector length based on the center positions of // cells in the first and last plane. Adjust to account for the // plane widths. This ajustment works if the vertical and horizontal // planes have the same depth, or if the detector begins and ends // with vertical planes. As of the TDR both of these conditions was // true for all the NOvA detectors. double xyz1[3] = {0,0,0}; double xyz2[3] = {0,0,0}; this->Plane(0). Cell(0).GetCenter(xyz1); this->Plane(fPlanes.size()-1).Cell(0).GetCenter(xyz2); fDetLength = (xyz2[2]-xyz1[2])+(vplanez2-vplanez1); }

void Geometry::SetDrawOptions (   )  [private]

Definition at line 267 of file Geometry.cxx. Referenced by Geometry(). { }

double Geometry::SurfaceY (   )  const

A typical y-position value at the surface (where earth meets air) for this detector site Returns: typical y position at surface in units of cm Exceptions: geo::Exception  if detector ID is not set properly Todo: : Only far det number is reasonable accurate. Near and IPND need some thought Definition at line 222 of file Geometry.cxx. References fGDMLFile. { if(fGDMLFile.find("near") != std::string::npos) return 130.0E2; // near: surface ~130m above det. center else if(fGDMLFile.find("far") != std::string::npos) return 9.94E2; // far : surface ~10m above det. center else if(fGDMLFile.find("ipnd") != std::string::npos) return -1.966E2;// IPND: surface ~2m below det. center else throw Exception(__FILE__,__LINE__,Exception::BAD_GEO_CONFIG); return 0.0; }

double Geometry::TotalMass (   )  const

Return the total mass of the detector Definition at line 496 of file Geometry.cxx. Referenced by main(). { double mass = 0.; // the detector resides in the vDetEnclosure volume so start there const TGeoVolume *enclosure = gGeoManager->FindVolumeFast("vDetEnclosure"); // // loop over the daughters in the volume and let ROOT calculate the // mass in kg for you for the daughter nodes that have vBlock in // their names this only goes one level down, so no worries of // double counting non-superblock blocks // for(int d = 0; d < enclosure->GetNdaughters(); ++d){ const TGeoNode* node = enclosure->GetNode(d); const char* nm = node->GetName(); if(strncmp(nm,"vSuperBlock",11) == 0){ std::string name = nm; name.erase(name.find('_')); mass += gGeoManager->FindVolumeFast(name.c_str())->Weight(); } } return mass; }

void Geometry::WorldBox (  double *  xlo, double *  xhi, double *  ylo, double *  yhi, double *  zlo, double *  zhi )  const

Return the ranges of x,y and z for the "world volume" that the entire geometry lives in. If any pointers are 0, then those coordinates are ignored. Parameters: xlo  : On return, lower bound on x positions xhi  : On return, upper bound on x positions ylo  : On return, lower bound on y positions yhi  : On return, upper bound on y positions zlo  : On return, lower bound on z positions zhi  : On return, upper bound on z positions Definition at line 244 of file Geometry.cxx. Referenced by testProject(). { const TGeoShape* s = gGeoManager->GetVolume("vWorld")->GetShape(); assert(s); if (xlo || xhi) { double x1, x2; s->GetAxisRange(1,x1,x2); if (xlo) *xlo = x1; if (xhi) *xhi = x2; } if (ylo || yhi) { double y1, y2; s->GetAxisRange(2,y1,y2); if (ylo) *ylo = y1; if (yhi) *yhi = y2; } if (zlo || zhi) { double z1, z2; s->GetAxisRange(3,z1,z2); if (zlo) *zlo = z1; if (zhi) *zhi = z2; } }

---

Member Data Documentation

std::set<unsigned int> geo::Geometry::fAllPlanes [private]

List of all planes. Definition at line 86 of file Geometry.h. Referenced by BuildMaps().

double geo::Geometry::fDetHalfHeight [private]

Detector 1/2 height (cm). Definition at line 91 of file Geometry.h. Referenced by SetDetectorSize().

double geo::Geometry::fDetHalfWidth [private]

Detector 1/2 width (cm). Definition at line 92 of file Geometry.h. Referenced by SetDetectorSize().

double geo::Geometry::fDetLength [private]

Detector 1/2 length (cm). Definition at line 90 of file Geometry.h. Referenced by SetDetectorSize().

std::string geo::Geometry::fGDMLFile [private]

gdml file holding the geometry Definition at line 82 of file Geometry.h. Referenced by Instance(), and SurfaceY().

UIDMap geo::Geometry::fIdMap [private]

Unique ID -> Plane,Cell. Definition at line 84 of file Geometry.h. Referenced by BuildMaps(), and IdToCell().

Geometry * Geometry::fInstance = 0 [static, private]

Instance of geometry. Definition at line 28 of file Geometry.cxx. Referenced by Instance().

std::vector<PlaneGeo*> geo::Geometry::fPlanes [private]

The detector planes. Definition at line 83 of file Geometry.h. Referenced by BuildMaps(), Geometry(), MakePlane(), NextPlaneInView(), NextPlaneOtherView(), Plane(), SetDetectorSize(), and ~Geometry().

int geo::Geometry::fRunNumber [private]

Run number of configuration. Definition at line 81 of file Geometry.h. Referenced by Instance().

std::set<unsigned int> geo::Geometry::fXplanes [private]

List of X measuring planes. Definition at line 87 of file Geometry.h. Referenced by BuildMaps().

std::set<unsigned int> geo::Geometry::fYplanes [private]

List of Y measuring planes. Definition at line 88 of file Geometry.h. Referenced by BuildMaps().

---

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

• Geometry.h
• Geometry.cxx

---