GeoMaterial.h

/*
  Copyright (C) 2002-2017 CERN for the benefit of the ATLAS collaboration
*/

#ifndef GEOMODELKERNEL_GEOMATERIAL_H
#define GEOMODELKERNEL_GEOMATERIAL_H

#include "GeoModelKernel/RCBase.h"
#include "GeoModelKernel/GeoElement.h"
#include <vector>

class GeoMaterial : public RCBase
{
 public:
  GeoMaterial (const std::string &Name, double Density);
  
  //    Add an element to the material.
  void add (const GeoElement* element, double fraction = 1.0);
  
  //    Add another material to the material (this copies all of
  //    the element from the other material into this one).
  //    Fraction is by mass.
  void add (const GeoMaterial* material, double fraction);
  
  //    Lock the material against the addition of other
  //    materials or elements.
  void lock ();
  
  //    Constant dEdx term.
  double getDeDxConstant () const;
  
  //    I0 term (ionization potential).  From tables. Units:  eV.
  double getDeDxI0 () const;
  
  //    Get dEdx_min.  dEdxConstant*11.528 Paul Avery, CBX-92-39.
  double getDeDxMin () const;
  
  //    Returns the radiation length, computed from the density
  //    and the list of constituents, and their properties.
  double getRadLength () const;
  
  //    Return the nuclear interaction length, computed from the
  //    density, the list of constituents, and their properties.
  double getIntLength () const;
  
  unsigned int getNumElements () const;
  
  //    Gets the ith element.
  const GeoElement* getElement (unsigned int i) const;
  
  //    Gets the fraction by weight of the ith element
  double getFraction (int i) const;
  
  //    The name of the material.
  const std::string& getName () const;
  
  //    The density of the material.
  const double& getDensity () const;
  
  //    Gives an integral identifier for the material.  For
  //    convenience.
  const unsigned int& getID () const;
  
 protected:
  virtual ~GeoMaterial();
  
 private:
  GeoMaterial(const GeoMaterial &right);
  GeoMaterial & operator=(const GeoMaterial &right);
  
  std::string m_name;
  double m_density;
  unsigned int m_iD;

  //    A list of the fractional composition of each material.
  //    Fraction is by mass.
  std::vector<double> m_fraction;

  //    The radiation length of the material.
  double m_radLength;
  
  //    The interaction length of the material.
  double m_intLength;

  //    A static used to assign unique identifiers to new
  //    materials.
  static unsigned int s_lastID;

  //    The constant term in the formula governing dE/dx.
  double m_dedDxConst;

  //    The ionization potential in the formula governing dE/dx.
  double m_deDxI0;

  //    A flag used to lock the material from further addition
  //    of elements or other constituents.
  bool m_locked;

  static const double s_ionizationPotential[93];

 private:
  //    The list of GeoElements composing a GeoMaterial.
  std::vector<const GeoElement *> m_element;  
};

inline const std::string& GeoMaterial::getName () const
{
  return m_name;
}

inline const double& GeoMaterial::getDensity () const
{
  return m_density;
}

inline const unsigned int& GeoMaterial::getID () const
{
  return m_iD;
}

#endif