geo::CellGeo Class Reference

Encapsulate the cell geometry. More...

#include "/cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-18/GeometryObjects/CellGeo.h"

Public Member Functions
	CellGeo (const std::vector< const TGeoNode * > &n, int depth)
	Construct a cell geometry. More...
	~CellGeo ()
void	GetCenter (double *xyz, double localz=0.0) const
void	GetCenter (TVector3 &xyz, double localz=0.0) const
double	DistToReadOut (double localz) const
double	HalfW () const
double	HalfL () const
double	HalfD () const
void	LocalToWorld (const double *local, double *world) const
void	LocalToWorldVect (const double *local, double *world) const
void	WorldToLocal (const double *local, double *world) const
void	WorldToLocalVect (const double *local, double *world) const
void	TranslateCell (double dx, double dy=0, double dz=0)
void	TranslateCell (TGeoCombiTrans *comb)
void	RestoreOriginal ()
	Restore geometry to original. More...
const TGeoNode *	Node () const
const CellUniqueId &	Id () const

Static Public Member Functions
static void	clear ()

Private Member Functions
TGeoHMatrix	GetGeoMatrix () const
	Transformation matrix to world frame. More...

Private Attributes
CellUniqueId	fId
	Node ID number. More...
const TGeoNode *	fCellNode
	Pointer to the cell node. More...
TGeoHMatrix *	fTranslatedMatrix
	Adjusted matrix, if TranslateCell() has been called. More...
int	fParentMatrixIdx
	Index into fgParentMatrices. More...

Static Private Attributes
static std::vector< TGeoHMatrix >	fgParentMatrices
	Collection of unique parent matrices seen. More...

Detailed Description

Encapsulate the cell geometry.

A note on the cell geometry: Cells are constructed such that, in their local frame, their profile occupies the x-y plane with their long dimension running along the z-axis. The dimension transverse to the beam runs along the local x and is called the width. The dimension parallel to the beam runs along the local y and is called the depth. The length of the cell runs along the local z-axis with z increasing as one moves toward the readout end of the cell.

Definition at line 25 of file CellGeo.h.

Constructor & Destructor Documentation

geo::CellGeo::CellGeo	(	const std::vector< const TGeoNode * > &	n,
		int	depth
	)

Construct a cell geometry.

Definition at line 23 of file CellGeo.cxx.

References om::cerr, allTimeWatchdog::endl, fCellNode, fgParentMatrices, fId, fParentMatrixIdx, MECModelEnuComparisons::i, and geo::NodesToUniqueId().

    : fTranslatedMatrix(0)
  {
    fCellNode = n[depth];

    // The current split between parentMatrix and the final step is a
    // compromise between speed and memory usage. This way, we can find our
    // transformation matrix from the world frame with only one multiply, but
    // we use something like 16x less storage than if we kept the final matrix
    // for every single cell.

    TGeoHMatrix parentMatrix = *n.front()->GetMatrix();
    // This loop is omitting the last entry in 'n'. That'll be available
    // through fCellNode.
    for(int i = 1; i < depth; ++i){
      parentMatrix.Multiply(n[i]->GetMatrix());
    }

    // We know that the geometry construction traverses depth-first, so if this
    // parent matrix has seen before it will actually have been immediately
    // prior.
    if(fgParentMatrices.empty() ||
       !(parentMatrix == fgParentMatrices.back())){
      fgParentMatrices.push_back(parentMatrix);

      // Every module in the detector is made of two halves, so count
      // 2*#modules to get the expected number of parent matrices.
      if(fgParentMatrices.size() == 64*14*12*2+1){
        std::cerr << "Warning: CellGeo::fgParentMatrices size exceeds "
                  << 64*14*12*2 << " expected." << std::endl;
      }

      // If all is well, this is the last matrix, so don't waste any additional
      // space in the vector.
      if(fgParentMatrices.size() == 64*14*12*2){
        fgParentMatrices.shrink_to_fit();
      }
    }
    // Either it matches the previous matrix, or we just added it to the list,
    // either way it's the last index.
    fParentMatrixIdx = fgParentMatrices.size()-1;

    // Construct the unique id for this cell
    fId = NodesToUniqueId(n, depth);
  }

geo::CellGeo::~CellGeo

(

)

Definition at line 70 of file CellGeo.cxx.

References fTranslatedMatrix.

  {
    delete fTranslatedMatrix;
  }

Member Function Documentation

static void geo::CellGeo::clear ( )

inlinestatic

Definition at line 57 of file CellGeo.h.

References fgParentMatrices.

Referenced by geo::GeometryBase::LoadGeometryFile().

{fgParentMatrices.clear();}

double geo::CellGeo::DistToReadOut

(

double

localz

)

const

Return the distance from a point in the cell to the read out side of the cell

Parameters

localz

: Z position in cm in local cell frame

Returns

Distance in cm

Definition at line 184 of file CellGeo.cxx.

References HalfL().

  {
    return this->HalfL()-localz;
  }

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

Return the center position of a cell.

Parameters

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

Definition at line 159 of file CellGeo.cxx.

References LocalToWorld().

Referenced by geo::GeometryTest::analyze(), slicer::S4DParamCalc::analyze(), mcchk::RockAna::analyze(), align::Alignment::analyze(), sn::SNSlicerAna::analyze(), calib::HitEfficiency::analyze(), trk::KalmanTrackAna::analyze(), remid::ReMIdDedxRock::analyze(), remid::ReMIdDedxFD::analyze(), remid::ReMIdDedxStudies::analyze(), calib::CalibAna::ApproximateHitPos(), trk::KalmanTrackMerge::CanJoinTracks(), geo::GeometryBase::CellInfo(), dt::Chunk::CellMeasurements(), geo::GeometryBase::CellTpos(), geo::GeometryBase::ClosestApproach(), dt::Cand::ClosestToEndCell(), geo::GeometryBase::CountCellsOnLineFast(), murem::TrackCleanUpAlg::DeDxInPlane(), evd::GeometryDrawer::DrawCells2D(), evd::GeometryDrawer::DrawHighlightCell(), geo::LiveGeometry::FillBadBoxes(), bpfit::BreakPoint::FillHitList(), trk::KalmanGeoHelper::fillprivategeo(), bsf::BremShowerFilter::filter(), trk::KalmanTrack::FindTracks(), trident::DimuonSliceAna::FindVertex(), bpfit::DimuonFitter::FindVertexZ(), trk::CosmicTrackUtilities::FindZBoundaries(), calib::FindZBoundaries(), bpfit::DimuonFitter::FitView(), trk::WindowTrackingAlg::FitWindow(), evd::GeometryDrawer::GetBox(), jmshower::RecoJMShower::GetCellDistToPoint(), jmshower::RecoJMShower::GetCellDistToTrk(), slid::ParticleIDAlg::GetCellNodePos(), jmshower::RecoJMShower::GetCellNodePos(), GetCenter(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), geo::LiveGeometry::GetDetectorEdges(), skim::ParametersNumu::GetDetectorEdges(), airshower::AirSlicer::GetHitPos(), ndreco::NDRecoFxs::getProngActivity(), ndreco::NDRecoFxs::getProngDedx(), ndreco::NDRecoFxs::getTrackActivity(), ndreco::NDRecoFxs::getTrackDedx(), jmshower::JMShower::GetTransHits(), jmshower::RecoJMShower::GetTrkHitPath(), jmshower::RecoJMShower::GetTrkHitPos(), jmshower::RecoJMShower::GetTrkPlaneCell(), jmshower::RecoJMShower::GetTrkPlaneDistToEdge(), jmshower::RecoJMShower::GetTrkPlanePos(), calib::GetXYZD(), rb::HitList::HitList(), dt::Chunk::HitsOnLine(), nerd::NERDProng::HitToHitDistance(), me::MEFinder::HitToHitDistance(), bsf::BremShowerFilter::inFiducial(), dif::DiFShowerFinder::inFiducial(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), lem::FindLEMMatches::MatchToVertex(), lem::LEM::MatchToVertex(), zcl::SMMCluster::MissC(), numue::NumuEAlg::NumuEAlg(), trk::CosmicTrackUtilities::PathLengthInCell(), calib::PathLengthInCell(), geo::plane_sort(), rb::Track::PlaneDirMap(), slid::ParticleIDAlg::PlaneHitCell(), slid::ParticleIDAlg::PlaneHitXYZ(), hv::HoughVertex::produce(), trk::KalmanTrack::produce(), trk::KalmanTrackMerge::produce(), jmshower::JMTrackMerge::produce(), zcl::SMMCluster::produce(), vdt::VertexDT::produce(), cvn::RegCVNMapper::produce(), cheat::MCCheater::produce(), remid::RecoMuon::produce(), evd::RawDataDrawer::RawDigit2D(), jmshower::RecoJMShower::RecoShowers(), trk::RLFit::RLFit(), hough::MultiHoughT::Scrub(), earms::ElasticArmsHS::Scrub(), vdt::VertexDT::Scrub(), fuzz::FuzzyKVertex::Scrub(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), trk::WindowTrackingAlg::ShortTrackExtraPlane(), geo::sort_hori(), geo::sort_vert(), murem::TrackCleanUpAlg::SortByDistFromTrack(), geo::GeometryTest::testCellId(), geo::GeometryTest::testCellIdFromPos(), geo::GeometryTest::testCellPos(), slid::NuEEnergyAlg::VertexEnergy(), and lem::VertexToPlaneAndCell().

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

void geo::CellGeo::GetCenter	(	TVector3 &	xyz,
		double	localz = 0.0
	)		const

Return the center position of a cell.

Parameters

xyz	:a Tvector3 object. The returned location.
localz	: Distance along the length of the cell (cm). Default is center of cell

Definition at line 170 of file CellGeo.cxx.

References GetCenter().

  {
    double xyzArray[3];
    this->GetCenter(xyzArray, localz);
    xyz.SetXYZ(xyzArray[0], xyzArray[1], xyzArray[2]);
  }

TGeoHMatrix geo::CellGeo::GetGeoMatrix ( ) const

private

Transformation matrix to world frame.

Definition at line 211 of file CellGeo.cxx.

References fCellNode, fgParentMatrices, fParentMatrixIdx, fTranslatedMatrix, and runNovaSAM::ret.

Referenced by LocalToWorld(), LocalToWorldVect(), TranslateCell(), WorldToLocal(), and WorldToLocalVect().

  {
    if(fTranslatedMatrix) return *fTranslatedMatrix;

    TGeoHMatrix ret = fgParentMatrices[fParentMatrixIdx];
    ret.Multiply(fCellNode->GetMatrix());
    return ret;
  }

double geo::CellGeo::HalfD

(

)

const

Definition at line 205 of file CellGeo.cxx.

References fCellNode.

Referenced by dif::DiFShowerFinder::adjustPlane(), dif::DiFShowerFinder::adjustPlane_end(), geo::GeometryTest::analyze(), remid::ReMIdDedxRock::analyze(), remid::ReMIdDedxFD::analyze(), remid::ReMIdDedxStudies::analyze(), trk::CosmicTrackAna::analyze(), trk::CosmicTrackAna::beginJob(), calib::BestPathEstimates(), slid::ParticleIDAlg::CalcTrkHitPath(), trk::KalmanTrackMerge::CanJoinTracks(), geo::GeometryBase::CellInfo(), dt::Chunk::CellMeasurements(), geo::GeometryBase::CountCellsOnLineFast(), murem::TrackCleanUpAlg::DeDxInPlane(), evd::GeometryDrawer::DrawCells2D(), evd::GeometryDrawer::DrawHighlightCell(), calib::BetheBlochFit::endJob(), dif::DiFShowerFinder::eparm(), calib::StopperThreshold::FillHist(), bpfit::BreakPoint::FillHitList(), trk::KalmanGeoHelper::fillprivategeo(), calib::StopperThreshold::FillTree(), bsf::BremShowerFilter::filter(), trk::WindowTrackingAlg::FindEndPoint(), bsf::BremShowerFilter::findShowerByReco(), dif::DiFShowerFinder::findShowerByReco(), trk::KalmanTrack::FindTracks(), trk::CosmicTrackUtilities::FindZBoundaries(), calib::FindZBoundaries(), evd::GeometryDrawer::GetBox(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), nuesand::FillNueSandbox::GetdEdx(), jmshower::RecoJMShower::GetTrkHitPath(), rb::HitList::HitList(), dt::Chunk::HitsOnLine(), bsf::BremShowerFilter::inFiducial(), dif::DiFShowerFinder::inFiducial(), skim::EvaluatorNue::KeepSlice(), trk::CosmicTrackAlg::MakeTrack(), trk::WindowTrackingAlg::MakeTrack(), calib::MakeZBoundaryMap(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), dif::DiFShowerFinder::muonstub(), trk::CosmicTrackUtilities::PathLengthInCell(), calib::PathLengthInCell(), hv::HoughVertex::produce(), trk::KalmanTrack::produce(), trk::KalmanTrackMerge::produce(), remid::RecoMuon::produce(), caf::CAFMaker::produce(), evd::RawDataDrawer::RawDigit2D(), jmshower::RecoJMShower::RecoShowers(), trk::KalmanTrackMerge::ShiftInterpolationPoints(), calib::StopperThreshold::testPath(), trk::KalmanTrackMerge::ViewMergeTracks(), and calib::ZBounds().

  {
    return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDY();
  }

double geo::CellGeo::HalfL

(

)

const

Definition at line 198 of file CellGeo.cxx.

References fCellNode.

Referenced by geo::GeometryTest::analyze(), LightLevels::analyze(), trk::CosmicTrackAna::beginJob(), geo::GeometryBase::CellInfo(), geo::GeometryBase::ClosestApproach(), LightLevels::CreateHistos(), geo::GeometryBase::DistToElectronics(), DistToReadOut(), calib::Calibrator::GetAttenCurve(), align::SplitTracks::GetCellEndpoints(), calib::CosmicCalib::GetChannelHists(), calib::MakeAttenuationProfiles::GetChannelHists(), photrans::ImprovedTransport::GetPosCorr(), jmshower::RecoJMShower::GetTrkHitPath(), calib::AttenCache::LoadFromDatabase(), trk::CosmicTrackUtilities::PathLengthInCell(), calib::PathLengthInCell(), trk::KalmanTrack::produce(), trk::KalmanTrackMerge::produce(), photrans::PhotonTransport::produce(), and photrans::ImprovedTransport::StepAlongHit().

  {
    return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDZ();
  }

double geo::CellGeo::HalfW

(

)

const

Definition at line 191 of file CellGeo.cxx.

References fCellNode.

Referenced by geo::GeometryTest::analyze(), calib::HitEfficiency::analyze(), remid::ReMIdDedxRock::analyze(), remid::ReMIdDedxFD::analyze(), remid::ReMIdDedxStudies::analyze(), trk::CosmicTrackAna::beginJob(), calib::BestPathEstimates(), trk::KalmanTrackMerge::CanJoinTracks(), geo::GeometryBase::CellInfo(), dt::Chunk::CellMeasurements(), geo::GeometryBase::CountCellsOnLineFast(), trk::KalmanGeoHelper::CountMissedCellsOnLine(), evd::GeometryDrawer::DrawCells2D(), evd::GeometryDrawer::DrawHighlightCell(), bpfit::BreakPoint::FillHitList(), trk::CosmicTrackAna::FillTrackHistograms(), calib::FindBelowThresholdCalibCandidates(), trk::KalmanTrack::FindTracks(), bpfit::DimuonFitter::FitView(), evd::GeometryDrawer::GetBox(), evd::RecoBaseDrawer::GetClusterOutlinesConvexHull(), jmshower::JMShower::GetTransHits(), jmshower::RecoJMShower::GetTrkPlaneCell(), rb::HitList::HitList(), dt::Chunk::HitsOnLine(), skim::EvaluatorNue::KeepSlice(), fuzz::ViewMatchAlg::Matching(), nerd::ViewMatchAlg::Matching(), trk::KalmanGeoHelper::MatchTrajectoryToPlane(), trk::CosmicTrackUtilities::PathLengthInCell(), calib::PathLengthInCell(), slid::ParticleIDAlg::PlaneHitCell(), hv::HoughVertex::produce(), trk::KalmanTrack::produce(), trk::KalmanTrackMerge::produce(), slid::Recluster::produce(), remid::RecoMuon::produce(), caf::CAFMaker::produce(), and evd::RawDataDrawer::RawDigit2D().

  {
    return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDX();
  }

const CellUniqueId& geo::CellGeo::Id ( ) const

inline

Definition at line 55 of file CellGeo.h.

References fId.

Referenced by geo::GeometryBase::BuildMaps(), geo::GeometryBase::DistToElectronics(), calib::MakeAttenuationProfiles::event(), trk::CosmicTrackAna::FillTrackHistograms(), calib::CosmicCalib::filter(), photrans::PhotonTransport::produce(), photrans::ImprovedTransport::StepAlongHit(), geo::GeometryTest::testCellId(), geo::GeometryTest::testCellIdFromPos(), geo::GeometryTest::testUniqueId(), and trk::CosmicTrackUtilities::TruePathLengthInCell().

{ return fId; }

void geo::CellGeo::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 80 of file CellGeo.cxx.

References GetGeoMatrix().

Referenced by zcl::SMMTriggerAna::analyze(), mcchk::CosmicAna::analyze(), calib::AssessCalib::analyze(), zcl::FmmTrackerValidation::analyze(), htk::HoughTrack::analyze(), mcchk::ShowerAnaCheck::EMMoliereRadius(), mcchk::ShowerAnaCheck::EMShowerWidth(), align::SplitTracks::GeoEstimatePos(), align::SplitTracks::GetCellEndpoints(), GetCenter(), zcl::SMMCluster::produce(), cheat::MCCheater::produce(), geo::GeometryTest::testFindCell(), and geo::GeometryTest::testStepping().

  {
    GetGeoMatrix().LocalToMaster(local,world);
  }

void geo::CellGeo::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 90 of file CellGeo.cxx.

References GetGeoMatrix().

Referenced by geo::GeometryTest::testFindCell(), and geo::GeometryTest::testStepping().

  {
    GetGeoMatrix().LocalToMasterVect(local,world);
  }

const TGeoNode* geo::CellGeo::Node ( ) const

inline

Definition at line 54 of file CellGeo.h.

References fCellNode.

{ return fCellNode; }

void geo::CellGeo::RestoreOriginal

(

)

Restore geometry to original.

Definition at line 148 of file CellGeo.cxx.

References fTranslatedMatrix.

  {
    delete fTranslatedMatrix;
  }

void geo::CellGeo::TranslateCell	(	double	dx,
		double	dy = 0,
		double	dz = 0
	)

Translate plane and all cells in plane by dx, dy and dz New transform will be stored in a separate TGeoHMatrix

Set rotation angles to 0 for now

Build the translation and rotation matrices phi, theta and psi are euler angles

dx, dy and dz are the translations we want

Definition at line 117 of file CellGeo.cxx.

References chisquared::theta.

  {

    /// Set rotation angles to 0 for now
    double phi = 0;
    double theta = 0;
    double psi = 0;

    /// Build the translation and rotation matrices
    /// phi, theta and psi are euler angles
    TGeoRotation *rot = new TGeoRotation("alignmentRot",phi,theta,psi);
    /// dx, dy and dz are the translations we want
    TGeoTranslation *trans = new TGeoTranslation("",dx,dy,dz);
    TGeoCombiTrans *comb = new TGeoCombiTrans(*trans,*rot);

    this->TranslateCell(comb);

  }

void geo::CellGeo::TranslateCell

(

TGeoCombiTrans *

comb

)

Translate plane using a previously defined TGeoCombiTrans New transform will be stored in a separate TGeoHMatrix

Definition at line 138 of file CellGeo.cxx.

References fTranslatedMatrix, and GetGeoMatrix().

  {
    if(fTranslatedMatrix) delete fTranslatedMatrix;

    fTranslatedMatrix = new TGeoHMatrix(GetGeoMatrix());
    fTranslatedMatrix->MultiplyLeft(comb);
  }

void geo::CellGeo::WorldToLocal	(	const double *	local,
		double *	world
	)		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 100 of file CellGeo.cxx.

References GetGeoMatrix().

Referenced by calib::GetXYZD(), and g4n::FLSHitListAction::SteppingAction().

  {
    GetGeoMatrix().MasterToLocal(local,world);
  }

void geo::CellGeo::WorldToLocalVect	(	const double *	local,
		double *	world
	)		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 110 of file CellGeo.cxx.

References GetGeoMatrix().

  {
    GetGeoMatrix().MasterToLocalVect(local,world);
  }

Member Data Documentation

const TGeoNode* geo::CellGeo::fCellNode

private

Pointer to the cell node.

Definition at line 60 of file CellGeo.h.

Referenced by CellGeo(), GetGeoMatrix(), HalfD(), HalfL(), HalfW(), and Node().

std::vector< TGeoHMatrix > geo::CellGeo::fgParentMatrices

staticprivate

Collection of unique parent matrices seen.

Definition at line 70 of file CellGeo.h.

Referenced by CellGeo(), clear(), and GetGeoMatrix().

CellUniqueId geo::CellGeo::fId

private

Node ID number.

Definition at line 59 of file CellGeo.h.

Referenced by CellGeo(), and Id().

int geo::CellGeo::fParentMatrixIdx

private

Index into fgParentMatrices.

Definition at line 68 of file CellGeo.h.

Referenced by CellGeo(), and GetGeoMatrix().

TGeoHMatrix* geo::CellGeo::fTranslatedMatrix

private

Adjusted matrix, if TranslateCell() has been called.

Definition at line 65 of file CellGeo.h.

Referenced by GetGeoMatrix(), RestoreOriginal(), TranslateCell(), and ~CellGeo().

---

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

• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-18/GeometryObjects/CellGeo.h
• /cvmfs/nova-development.opensciencegrid.org/novasoft/releases/N21-01-18/GeometryObjects/CellGeo.cxx