Material.h

Go to the documentation of this file.

#ifndef DDRec_Material_H
#define DDRec_Material_H

#include "DD4hep/Objects.h"
#include "DDSurfaces/IMaterial.h"

#include <list>

namespace DD4hep {
  namespace DDRec {
    
    
    class MaterialData : public DDSurfaces::IMaterial{
      
    protected:
      std::string _name ;
      double _Z ;
      double _A ;
      double _rho ;
      double _x0 ;
      double _lambda ;

    public:

      MaterialData( Geometry::Material m ) : 

        _name("unknown"),
        _Z( -1. ),
        _A( 0. ),
        _rho( 0. ),
        _x0( 0. ),
        _lambda( 0.)  {

        if( m.isValid() ) {

          _name= m.name() ;
          _Z = m.Z() ;
          _A = m.A() ;
          _rho = m.density() ;
          _x0 = m.radLength() ;
          _lambda = m.intLength() ;

        }
      }
      
      MaterialData()  : _name("unknown"),
                        _Z( -1. ),
                        _A( 0. ),
                        _rho( 0. ),
                        _x0( 0. ),
                        _lambda( 0.) {}

      MaterialData( const std::string& nam, double Z_val, double A_val, double density_val, double radLength, double intLength )
        : _name( nam ),
          _Z( Z_val ),
          _A( A_val ),
          _rho( density_val ),
          _x0( radLength ),
          _lambda(  intLength ) {}
      
      MaterialData( const MaterialData& m )  : _name( m.name() ),
                                               _Z( m.Z() ),
                                               _A( m.A() ),
                                               _rho( m.density() ),
                                               _x0( m.radiationLength() ),
                                               _lambda( m.interactionLength() ) {}

      MaterialData( const IMaterial& m )  : _name( m.name() ),
                                            _Z( m.Z() ),
                                            _A( m.A() ),
                                            _rho( m.density() ),
                                            _x0( m.radiationLength() ),
                                            _lambda( m.interactionLength() ) {}
      
       MaterialData& operator=(const MaterialData& m){
        if ( this != &m )  {
          _name = m._name ;
          _Z = m._Z ;
          _A = m._A  ;
          _rho = m._rho ;
          _x0 = m._x0 ;
          _lambda = m._lambda ;
        }
        return *this ;
      }

      MaterialData& operator=(const IMaterial& m){
        if ( this != &m )  {
          _name = m.name() ;
          _Z = m.Z() ;
          _A = m.A() ;
          _rho = m.density() ;
          _x0 = m.radiationLength() ;
          _lambda = m.interactionLength() ;
        }
        return *this ;
      }

      MaterialData& operator=(const Geometry::Material& m){
      
        if( m.isValid() ) {

          _name = m.name() ;
          _Z = m.Z() ;
          _A = m.A() ;
          _rho = m.density() ;
          _x0 = m.radLength() ;
          _lambda = m.intLength() ;

        }  else {

          _name= "unknown";
          _Z = -1.  ;
          _A =  0. ;
          _rho = 0. ;
          _x0 = 0. ;
          _lambda = 0. ;
        }

        return *this ;
      }

      bool isValid() const { return ( _Z > 0. ) ; }

      virtual ~MaterialData() {}

      virtual std::string name() const { return _name ; }
      
      virtual double Z() const {  return _Z ; } 
      
      virtual double A() const { return _A ; } 
      
      virtual double density() const {  return _rho ; }
      
      virtual double radiationLength() const { return _x0 ; } 
      
      virtual double interactionLength() const  { return _lambda ; }

    };


  } /* namespace */
} /* namespace */

#endif /* DDRec_Material_H */