SurfaceManager.cpp

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
//
//==========================================================================
#include "DDRec/SurfaceManager.h"
 
#include "DDRec/SurfaceHelper.h"
#include "DD4hep/VolumeManager.h"
#include "DD4hep/Detector.h"
 
#include <sstream>
 
namespace dd4hep {
  
  using namespace detail ;
 
  namespace rec {
    
 
    SurfaceManager::SurfaceManager(const Detector& theDetector){
 
      // have to make sure the volume manager is populated once in order to have
      // the volumeIDs attached to the DetElements
 
      VolumeManager::getVolumeManager(theDetector);
 
      initialize(theDetector) ;
    }
    
    SurfaceManager::~SurfaceManager(){
      // nothing to do
    }
    
    
    const SurfaceMap* SurfaceManager::map( const std::string name ) const {
 
      SurfaceMapsMap::const_iterator it = _map.find( name ) ;
 
      if( it != _map.end() ){
 
        return & it->second ;
      }
 
      return 0 ;
    }
 
    void SurfaceManager::initialize(const Detector& description) {
      
      const std::vector<std::string>& types = description.detectorTypes() ;
 
      for(unsigned i=0,N=types.size();i<N;++i){
 
        const std::vector<DetElement>& dets = description.detectors( types[i] ) ;  
 
        for(unsigned j=0,M=dets.size();j<M;++j){
 
          std::string name = dets[j].name() ;
 
          SurfaceHelper surfH( dets[j] ) ;
      
          const SurfaceList& detSL = surfH.surfaceList() ;
  
          // add an empty map for this detector in case there are no surfaces attached 
          _map.emplace(name , SurfaceMap());
 
          for( SurfaceList::const_iterator it = detSL.begin() ; it != detSL.end() ; ++it ){
            ISurface* surf =  *it ;
        
            // enter surface into map for this detector
            _map[ name ].emplace(surf->id(), surf );
 
            // enter surface into map for detector type
            _map[ types[i] ].emplace(surf->id(), surf );
 
            // enter surface into world map 
            _map[ "world" ].emplace(surf->id(), surf );
 
          }
        }
      }
 
    }
 
    std::string SurfaceManager::toString() const {
      
      std::stringstream sstr ;
       
      sstr << "--------  SurfaceManager contains the following maps : --------- " << std::endl ;
 
      for( SurfaceMapsMap::const_iterator mi = _map.begin() ; mi != _map.end() ; ++mi ) {
    
        sstr << "  key: " <<  mi->first << " \t number of surfaces : " << mi->second.size() << std::endl ; 
      }
      sstr << "---------------------------------------------------------------- " << std::endl ;
 
      return sstr.str() ;
    }
 
 
  } // namespace
}// namespace