MaterialManager.h

Go to the documentation of this file.

//==========================================================================
//  AIDA Detector description implementation 
//--------------------------------------------------------------------------
// Copyright (C) Organisation europeenne pour la Recherche nucleaire (CERN)
// All rights reserved.
//
// For the licensing terms see $DD4hepINSTALL/LICENSE.
// For the list of contributors see $DD4hepINSTALL/doc/CREDITS.
//
// Author     : F.Gaede
//
//==========================================================================
#ifndef DDREC_MATERIALMANAGER_H
#define DDREC_MATERIALMANAGER_H

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

#include <vector>


class TGeoManager ;

namespace dd4hep {
  namespace rec {

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

    public:

      MaterialManager(Volume world);

#if defined(G__ROOT)
      MaterialManager() = default ;
#else
      MaterialManager() = delete ;
#endif

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

      const PlacementVec& placementsBetween(const Vector3D& p0, const Vector3D& p1 , double epsilon=1e-4 );
      
      const Material& materialAt(const Vector3D& pos );

      PlacedVolume placementAt(const Vector3D& pos );

      MaterialData createAveragedMaterial( const MaterialVec& materials ) ;

    protected :
      MaterialVec  _mV ;
      Material     _m ;
      PlacedVolume _pv;
      PlacementVec _placeV;
      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 rec */
} /* namespace dd4hep */



#endif // DDREC_MATERIALMANAGER_H