Garfield::GeometryRoot Class Reference

#include <GeometryRoot.hh>

Inheritance diagram for Garfield::GeometryRoot:

Classes
struct	material

Public Member Functions
	GeometryRoot ()
	~GeometryRoot ()
void	SetGeometry (TGeoManager *geoman)
Medium *	GetMedium (const double x, const double y, const double z) const
int	GetNumberOfMaterials ()
TGeoMaterial *	GetMaterial (const int i)
TGeoMaterial *	GetMaterial (const char *name)
void	SetMedium (const int imat, Medium *med)
void	SetMedium (const char *mat, Medium *med)
bool	IsInside (const double x, const double y, const double z) const
bool	GetBoundingBox (double &xmin, double &ymin, double &zmin, double &xmax, double &ymax, double &zmax)
void	EnableDebugging ()
void	DisableDebugging ()
Public Member Functions inherited from Garfield::GeometryBase
	GeometryBase ()
virtual	~GeometryBase ()
virtual Medium *	GetMedium (const double x, const double y, const double z) const =0
virtual bool	IsInside (const double x, const double y, const double z) const =0
virtual bool	GetBoundingBox (double &xmin, double &ymin, double &zmin, double &xmax, double &ymax, double &zmax)=0

Protected Attributes
TGeoManager *	m_geoManager
int	m_nMaterials
std::vector< material >	m_materials
bool	m_debug
Protected Attributes inherited from Garfield::GeometryBase
std::string	m_className

Detailed Description

Definition at line 15 of file GeometryRoot.hh.

Constructor & Destructor Documentation

GeometryRoot()

Garfield::GeometryRoot::GeometryRoot

(

)

Definition at line 12 of file GeometryRoot.cc.

    : m_geoManager(NULL), m_nMaterials(0), m_debug(false) {
 
  m_className = "GeometryRoot";
}

~GeometryRoot()

Garfield::GeometryRoot::~GeometryRoot ( )

inline

Definition at line 21 of file GeometryRoot.hh.

{}

Member Function Documentation

DisableDebugging()

void Garfield::GeometryRoot::DisableDebugging ( )

inline

Definition at line 54 of file GeometryRoot.hh.

{ m_debug = false; }

EnableDebugging()

void Garfield::GeometryRoot::EnableDebugging ( )

inline

Definition at line 53 of file GeometryRoot.hh.

{ m_debug = true; }

GetBoundingBox()

bool Garfield::GeometryRoot::GetBoundingBox ( double &  xmin, double &  ymin, double &  zmin, double &  xmax, double &  ymax, double &  zmax  )

virtual

Implements Garfield::GeometryBase.

Definition at line 168 of file GeometryRoot.cc.

                                                                            {
 
  if (!m_geoManager) return false;
  TGeoBBox* box = (TGeoBBox*)m_geoManager->GetTopVolume()->GetShape();
  const double dx = box->GetDX();
  const double dy = box->GetDY();
  const double dz = box->GetDZ();
  const double ox = box->GetOrigin()[0];
  const double oy = box->GetOrigin()[1];
  const double oz = box->GetOrigin()[2];
  xmin = ox - dx;
  xmax = ox + dx;
  ymin = oy - dy;
  ymax = oy + dy;
  zmin = oz - dz;
  zmax = oz + dz;
  return true;
}

GetMaterial() [1/2]

TGeoMaterial * Garfield::GeometryRoot::GetMaterial

(

const char *

name

)

Definition at line 72 of file GeometryRoot.cc.

                                                        {
 
  if (!m_geoManager) {
    std::cerr << "GeometryRoot::GetMaterial:\n";
    std::cerr << "    ROOT geometry is not defined.\n";
    std::cerr << "    Call SetGeometry first.\n";
    return NULL;
  }
 
  return m_geoManager->GetMaterial(name);
}

GetMaterial() [2/2]

TGeoMaterial * Garfield::GeometryRoot::GetMaterial

(

const int

i

)

Definition at line 60 of file GeometryRoot.cc.

                                                   {
 
  if (!m_geoManager) {
    std::cerr << "GeometryRoot::GetMaterial:\n";
    std::cerr << "    ROOT geometry is not defined.\n";
    std::cerr << "    Call SetGeometry first.\n";
    return NULL;
  }
 
  return m_geoManager->GetMaterial(i);
}

GetMedium()

Medium * Garfield::GeometryRoot::GetMedium ( const double  x, const double  y, const double  z  ) const

virtual

Implements Garfield::GeometryBase.

Definition at line 31 of file GeometryRoot.cc.

                                                      {
 
  if (!m_geoManager) return NULL; 
  m_geoManager->SetCurrentPoint(x, y, z);
  if (m_geoManager->IsOutside()) return NULL;
  TGeoNode* cnode = m_geoManager->GetCurrentNode();
  std::string name(cnode->GetMedium()->GetMaterial()->GetName());
 
  for (int i = m_nMaterials; i--;) {
    if (m_materials[i].name == name) {
      return m_materials[i].medium;
    }
  }
  return NULL;
}

GetNumberOfMaterials()

int Garfield::GeometryRoot::GetNumberOfMaterials

(

)

Definition at line 48 of file GeometryRoot.cc.

                                       {
 
  if (!m_geoManager) {
    std::cerr << "GeometryRoot::GetNumberOfMaterials:\n";
    std::cerr << "    ROOT geometry is not defined.\n";
    std::cerr << "    Call SetGeometry first.\n";
    return 0;
  }
 
  return m_geoManager->GetListOfMaterials()->GetEntries();
}

IsInside()

bool Garfield::GeometryRoot::IsInside ( const double  x, const double  y, const double  z  ) const

inlinevirtual

Implements Garfield::GeometryBase.

Definition at line 39 of file GeometryRoot.hh.

                                                                      {
 
    if (m_geoManager) {
      m_geoManager->SetCurrentPoint(x, y, z);
      return !m_geoManager->IsOutside();
    }
    return false;
  }

SetGeometry()

void Garfield::GeometryRoot::SetGeometry

(

TGeoManager *

geoman

)

Definition at line 18 of file GeometryRoot.cc.

                                                  {
 
  if (geoman == NULL) {
    std::cerr << "GeometryRoot::SetGeometry:\n";
    std::cerr << "    Pointer to TGeoManager is null.\n";
    return;
  }
 
  m_geoManager = geoman;
  m_materials.clear();
  m_nMaterials = 0;
}

SetMedium() [1/2]

void Garfield::GeometryRoot::SetMedium	(	const char *	mat,
		Medium *	med
	)

Definition at line 142 of file GeometryRoot.cc.

                                                          {
 
  if (!m_geoManager) {
    std::cerr << "GeometryRoot::SetMedium:\n";
    std::cerr << "    ROOT geometry is not defined.\n";
    std::cerr << "    Call SetGeometry first.\n";
    return;
  }
 
  if (!med) {
    std::cerr << "GeometryRoot::SetMedium:\n";
    std::cerr << "    Medium pointer is null.\n";
    return;
  }
 
  const int imat = m_geoManager->GetMaterialIndex(name);
  if (imat < 0) {
    std::cerr << "GeometryRoot::SetMedium:" << std::endl;
    std::cerr << "    ROOT material " << name << " does not exist."
              << std::endl;
    return;
  }
 
  SetMedium(imat, med);
}

SetMedium() [2/2]

void Garfield::GeometryRoot::SetMedium	(	const int	imat,
		Medium *	med
	)

Definition at line 84 of file GeometryRoot.cc.

                                                        {
 
  if (!m_geoManager) {
    std::cerr << "GeometryRoot::SetMedium:\n";
    std::cerr << "    ROOT geometry is not defined.\n";
    std::cerr << "    Call SetGeometry first.\n";
    return;
  }
 
  if (!med) {
    std::cerr << "GeometryRoot::SetMedium:\n";
    std::cerr << "    Medium pointer is null.\n";
    return;
  }
 
  TGeoMaterial* mat = m_geoManager->GetMaterial(imat);
  if (!mat) {
    std::cerr << "GeometryRoot::SetMedium:\n";
    std::cerr << "    ROOT material with index " << imat
              << " does not exist.\n";
    return;
  }
 
  std::string name(mat->GetName());
  bool isNew = true;
  // Check if this material has already been associated with a medium
  for (int i = m_nMaterials; i--;) {
    if (name == m_materials[i].name) {
      std::cout << "GeometryRoot::SetMedium:\n";
      std::cout << "    Current association of material " << name
                << " with medium " << med->GetName() << " is overwritten.\n";
      m_materials[i].medium = med;
      isNew = false;
      break;
    }
  }
 
  if (isNew) {
    material newMaterial;
    newMaterial.name = name;
    newMaterial.medium = med;
    m_materials.push_back(newMaterial);
    ++m_nMaterials;
  }
 
  // Check if material properties match
  const double rho1 = mat->GetDensity();
  const double rho2 = med->GetMassDensity();
  std::cout << "GeometryROOT::SetMedium:\n";
  std::cout << "    ROOT material: " << name << "\n";
  std::cout << "      Density: " << rho1 << " g / cm3\n";
  std::cout << "    Medium: " << med->GetName() << "\n";
  std::cout << "      Density: " << rho2 << " g / cm3\n";
  if (rho1 > 0 && fabs(rho1 - rho2) / rho1 > 0.01) {
    std::cout << "    WARNING: Densities differ by > 1%.\n";
  }
}

Referenced by SetMedium().

Member Data Documentation

m_debug

bool Garfield::GeometryRoot::m_debug

protected

Definition at line 69 of file GeometryRoot.hh.

Referenced by DisableDebugging(), and EnableDebugging().

m_geoManager

TGeoManager* Garfield::GeometryRoot::m_geoManager

protected

Definition at line 58 of file GeometryRoot.hh.

Referenced by GetBoundingBox(), GetMaterial(), GetMedium(), GetNumberOfMaterials(), IsInside(), SetGeometry(), and SetMedium().

m_materials

std::vector<material> Garfield::GeometryRoot::m_materials

protected

Definition at line 66 of file GeometryRoot.hh.

Referenced by GetMedium(), SetGeometry(), and SetMedium().

m_nMaterials

int Garfield::GeometryRoot::m_nMaterials

protected

Definition at line 61 of file GeometryRoot.hh.

Referenced by GetMedium(), SetGeometry(), and SetMedium().

---

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

• GeometryRoot.hh
• GeometryRoot.cc