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:
      static constexpr const double epsilon = 1e-4;
      
    public:
 
      MaterialManager(Volume world);
 
#if defined(G__ROOT)
      MaterialManager() = default ;
#else
      MaterialManager() = delete ;
#endif
 
      ~MaterialManager();
 
      class ScanData  {
      public:
        const MaterialVec&  materials;
        const PlacementVec& places;
      };
      
      const MaterialVec& materialsBetween(const Vector3D& p0,
                                          const Vector3D& p1,
                                          double eps = MaterialManager::epsilon);
      const ScanData     entriesBetween(const Vector3D& p0,
                                        const Vector3D& p1,
                                        double eps = MaterialManager::epsilon);
 
      const PlacementVec& placementsBetween(const Vector3D& p0, const Vector3D& p1 , double eps = MaterialManager::epsilon );
      
      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