In This Package:

Element.h

Go to the documentation of this file.

// $Id: Element.h,v 1.16 2008/10/28 12:04:36 cattanem Exp $
#ifndef DETDESC_ELEMENT_H
#define DETDESC_ELEMENT_H 1 
#include <vector>
#include "DetDesc/Material.h"
#include "DetDesc/CLIDElement.h"

class Isotope;

class Element : virtual public Material
{
public:
  typedef  std::pair<double,SmartRef<Isotope> >  Entry    ;
  typedef  std::vector<Entry>                    Isotopes ;  
public:
  Element( const std::string&  name     = ""   , 
           const std::string&  symb     = "??" ,
           const double        a        =  0   , 
           const double        z        =  0   , 
           const double        i        =  0   ,
           const double        density  =  0   ,
           const double        rl       =  0   , 
           const double        al       =  0   ,  
           const double        temp     =  Gaudi::Units::STP_Temperature,
           const double        press    =  Gaudi::Units::STP_Pressure,
           const eState        s        =  stateUndefined );
  // Destructor
  virtual ~Element();
  virtual inline const CLID& clID    () const { return Element::classID() ; };
  static         const CLID& classID ()       { return CLID_Element       ; }; 
  inline const std::string &      symbol  () const { return m_symb ; }
  inline void                  setSymbol (const std::string& S ) { m_symb = S; } 
  int                      nOfIsotopes() const;
  const SmartRef<Isotope>& isotope        ( unsigned int i ) const;
        SmartRef<Isotope>& isotope        ( unsigned int i )      ;
  const Isotopes&          isotopes       () const ;
        Isotopes&          isotopes       ()       ;
  double                   isotopeFraction( unsigned int i ) const;
  
  void addIsotope ( const SmartRef<Isotope>& iPtr  , 
                    const double             Fract , const bool comp = false );
  void addIsotope ( const Entry&             iPtr  , const bool comp = false );
  
  void removeIsotope ( const SmartRef<Isotope>& iPtr , const bool comp = false );

  virtual inline double    A() const;
  virtual inline void   setA( const double value ) ;
  virtual inline double    Z() const;
  virtual inline void   setZ( const double value );
  virtual inline double    I() const;
  virtual inline void   setI( const double value );
  virtual inline double    C() const;
  virtual inline void   setC( const double value);
  virtual inline double    a() const;
  virtual inline void   seta( const double value);
  virtual inline double    m() const;
  virtual inline void   setm( const double value);
  virtual inline double    X0() const;
  virtual inline void   setX0( const double value);
  virtual inline double    X1() const;
  virtual inline void   setX1( const double value);
  virtual inline double    N() const;
  inline double  coulombFactor() const ;
  inline double  tsaiFactor   () const ;

  void    compute();
  void    ComputeCoulombFactor  ();
  void    ComputeLradTsaiFactor ();

  void ComputeInteractionLength();

  void ComputeRadiationLength();

  void ComputeMeanExcitationEnergy();

  void ComputeDensityEffect();


  virtual std::ostream& fillStream ( std::ostream& s ) const ;
  virtual MsgStream&    fillStream ( MsgStream&    s ) const ;
private:  
  double                  m_Aeff;
  double                  m_Zeff;
  // Effective mean escitation energy
  double                  m_Ieff;
  // Parameters for density effect correction
  double                  m_C;
  double                  m_a;
  double                  m_m;
  double                  m_X0;
  double                  m_X1;
  Isotopes                m_isotopes;  
  double                  m_coulomb;
  double                  m_tsai;
  std::string             m_symb;
};
#include "DetDesc/Element.icpp"

#endif // DETDESC_MATERIAL_ELEMENT_H

---