CellGeo.cxx

Go to the documentation of this file.

#include <iostream>
#include "Geometry/CellGeo.h"
#include "TGeoManager.h"
#include "TGeoBBox.h"
#include "TGeoMatrix.h"
#include "TGeoNode.h"
#include "Geometry/CellUniqueId.h"
using namespace geo;

CellGeo::CellGeo(std::vector<const TGeoNode*>& n, int depth)
{
  fCellNode = n[depth];
  
  // Build the matrix that takes us to the top world frame
  fGeoMatrix = new TGeoHMatrix(*n[0]->GetMatrix());
  for (int i=1; i<=depth; ++i) {
    fGeoMatrix->Multiply(n[i]->GetMatrix());
  }

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

//......................................................................

void CellGeo::LocalToWorld(const double* local, double* world) const
{
  fGeoMatrix->LocalToMaster(local,world);
}

//......................................................................    

void CellGeo::LocalToWorldVect(const double* local, double* world) const
{
  fGeoMatrix->LocalToMasterVect(local,world);
}
    
//......................................................................

void CellGeo::WorldToLocal(const double* local, double* world) const
{
  fGeoMatrix->MasterToLocal(local,world);
}

//......................................................................

void CellGeo::WorldToLocalVect(const double* local, double* world) const
{
  fGeoMatrix->MasterToLocalVect(local,world);
}

//......................................................................

void CellGeo::GetCenter(double* xyz, double localz) const
{
  static double xyzLocal[3] = {0,0,0};
  xyzLocal[2] = localz;
  this->LocalToWorld(xyzLocal, xyz);
}

//......................................................................

double CellGeo::DistToReadOut(double localz) const 
{
  return this->HalfL()-localz;
}

//......................................................................

double CellGeo::HalfW() const 
{
  return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDX();
}
    
//......................................................................

double CellGeo::HalfL() const 
{
  return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDZ();
}

//......................................................................

double CellGeo::HalfD() const 
{
  return ((TGeoBBox*)fCellNode->GetVolume()->GetShape())->GetDY();
}

---