TouchableToDetectorElementFast.h

Go to the documentation of this file.

#ifndef TOUCHABLETODETECTORELEMENTFAST_H
#define TOUCHABLETODETECTORELEMENTFAST_H


#include "G4DataHelpers/ITouchableToDetectorElement.h"
#include "GaudiAlg/GaudiTool.h"
#include "DetDesc/ILVolume.h"
#include <list>

class TouchableToDetectorElementFast : public GaudiTool, virtual ITouchableToDetectorElement
{
  public:
  TouchableToDetectorElementFast(const std::string& type,
                             const std::string& name,
                            const IInterface* parent);
  virtual ~TouchableToDetectorElementFast() {};

  virtual StatusCode GetBestDetectorElement(
                          const G4TouchableHistory* inHistory,  // The current particle location
                          const std::vector<std::string>& inPath,// Name(s) of specific detElements, or paths to be searched
                          const IDetectorElement* &outElement,  // output: The best element (may be zero!) 
                          int& outCompatibility );              // output: the goodness of the match (lower is better, <=0 means no match.)

  virtual StatusCode G4VolumeToDetDesc( const G4VPhysicalVolume* inVol,     
                                        const IPVolume* &outVol             
                                       );

 virtual StatusCode ClearCache();

 private:
   
   // Things to store in the cache:
   struct Relation {
     Relation() : mLogVol(0), mPhysVol(0), mRpath(0) {};
     bool isNull() const { return mLogVol==0; };
     const ILVolume*          mLogVol; // The supporting volume
     const IPVolume*          mPhysVol;
     ILVolume::ReplicaPath    mRpath;  // Empty if it is not under the mWorldElement.
   };

   // The Cache:
   const IDetectorElement* mWorldElement;   // From recent calls. Changing this causes cache flushing.
   std::string             mWorldElementName; // The name of above.
   typedef std::map<const G4VPhysicalVolume*,Relation> G4ToRelationMap_t;
   typedef std::map<const IDetectorElement* ,Relation> DetElemToRelationMap_t;
   typedef std::list<const IDetectorElement*>          ElementList_t;
   G4ToRelationMap_t      mG4ToRelationMap;
   DetElemToRelationMap_t mDetElemToRelationMap;
   std::vector<std::string>             mLastSearchPaths;
   ElementList_t                        mElementList;
   
   StatusCode GetRelation(const G4VPhysicalVolume* inVol, Relation* &outRelation );
   StatusCode GetRelation(const IDetectorElement* inElem, Relation* &outRelation );

  int        Compatability(const ILVolume::ReplicaPath& inPlace, const ILVolume::ReplicaPath& inContainer);
  
  template < class T  >
  StatusCode FindObjectsInDirectory(const std::string& dirname, std::list<const T*>& outList);
  
};


#endif /* TOUCHABLETODETECTORELEMENTFAST_H */