MaterialManager.h

Go to the documentation of this file.

#ifndef DDRec_MaterialManager_H_
#define DDRec_MaterialManager_H_

#include "DD4hep/Objects.h"
#include "DDSurfaces/Vector3D.h"
#include "DDRec/Material.h"
#include "DD4hep/DD4hepUnits.h"

#include <vector>


class TGeoManager ;

namespace DD4hep {
  namespace DDRec {

    //  export Geometry::Material to this namespace ;
    using  Geometry::Material ;

    typedef std::vector< std::pair< Material, double > > MaterialVec ;
    
    class MaterialManager {

    public:

      MaterialManager();
      
      ~MaterialManager();
      
      const MaterialVec& materialsBetween(const DDSurfaces::Vector3D& p0, const DDSurfaces::Vector3D& p1 , double epsilon=1e-4 ) ;

      const Material& materialAt(const DDSurfaces::Vector3D& pos );


      MaterialData createAveragedMaterial( const MaterialVec& materials ) ;

    protected :

      //cached materials
      MaterialVec  _mV ;
      Material _m ;

      // cached last points
      DDSurfaces::Vector3D _p0 , _p1, _pos ;

      TGeoManager* _tgeoMgr ;
    };

   
    inline std::ostream& operator<<( std::ostream& os , const Material& m ) {
      os << "  " << m.name() << " Z: " << m.Z() << " A: " << m.A() << " density: " << m.density() 
         << " radiationLength: "  <<  m.radLength() 
         << " interactionLength: " << m.intLength() ;
      return os ;
    }
    

    inline std::ostream& operator<<( std::ostream& os , const MaterialVec& m ) {

      for( unsigned i=0,n=m.size() ; i<n ; ++i ) {
        os << "  material: " << m[i].first << " thickness: " << m[i].second << std::endl ; 
      }
      return os ;
    }

  } /* namespace DDRec */
} /* namespace DD4hep */

#endif // DDRec_MaterialManager_H_