In This Package:

TransportSvcIntersections.h

Go to the documentation of this file.

// $Id: TransportSvcIntersections.h,v 1.5 2007/09/13 14:19:16 wouter Exp $ 
// ============================================================================
// CVS tag $Name: v2r3 $
// ============================================================================
#ifndef     __DETDESC_TRANSPORTSVC_TRANSPORTSVCINTERSECTIONS_H__ 
#define     __DETDESC_TRANSPORTSVC_TRANSPORTSVCINTERSECTIONS_H__  1

#include "TransportSvc.h" 

// ============================================================================
// ============================================================================

// ============================================================================
// ============================================================================
unsigned long TransportSvc::intersections
( const Gaudi::XYZPoint&   point               , 
  const Gaudi::XYZVector&  vect                , 
  const ISolid::Tick&      tickMin             , 
  const ISolid::Tick&      tickMax             , 
  ILVolume::Intersections& intersept           , 
  double                   threshold           , 
  IGeometryInfo*           alternativeGeometry , 
  IGeometryInfo*           guessGeometry       ) 
{

  try { 
    
    intersept.clear();
    
    if( tickMin >= tickMax && vect.mag2() <= 0 ) { return 0;}
    
    // Set the top level geometry, because this is what is in the cache
    IGeometryInfo* topGeometry = alternativeGeometry ? alternativeGeometry : standardGeometry() ;
 
    Gaudi::XYZPoint point1( point + vect * tickMin ) ; 
    Gaudi::XYZPoint point2( point + vect * tickMax ) ; 
    if( point1              == m_prevPoint1          && 
        point2              == m_prevPoint2          &&
        threshold           == m_previousThreshold   &&
        topGeometry         == m_previousTopGeometry && 
        guessGeometry       == m_previousGuess         ) 
      {
        intersept.reserve( m_localIntersections.size() ); 
        intersept.insert(intersept.begin(),
                         m_localIntersections.begin() , 
                         m_localIntersections.end()); 
        return intersept.size();  
      }
    
    IGeometryInfo* giLocal = 0 ; 
    IGeometryInfo* gi      = 
      // try the guess geometry
     ( 0 != guessGeometry       && 
        0 != ( giLocal = findLocalGI( point1 , 
                                      point2 , 
                                      guessGeometry ,
                                      alternativeGeometry ) ) ) ? 
      guessGeometry       : 
      // try the previous geometry (only if top-geometry matches)
      ( previousGeometry() && topGeometry == m_previousTopGeometry &&
        0 != ( giLocal = findLocalGI( point1 , 
                                      point2 , 
                                      previousGeometry() ,
                                      topGeometry ) ) ) ? 
      // just take the top geometry
      previousGeometry() :
      ( 0 != ( giLocal = findLocalGI( point1 , 
                                      point2 , 
                                      topGeometry,
                                      topGeometry ) ) ) ? 
      
      topGeometry : 0 ;
    
    if( 0 == giLocal )
      throw TransportSvcException( "TransportSvc::(1) unable to locate points",
                                   StatusCode::FAILURE );
    if( 0 == gi )
      throw TransportSvcException( "TransportSvc::(2) unable to locate points",
                                   StatusCode::FAILURE );

    setPreviousGeometry( giLocal ); 
    
    const ILVolume* lv = giLocal->lvolume();   
    lv->intersectLine
      ( giLocal->toLocalMatrix() * point , 
        giLocal->toLocalMatrix() * vect  , 
        intersept                 , 
        tickMin                   , 
        tickMax                   , 
        threshold                 ); 
    
    m_prevPoint1           = point1              ;
    m_prevPoint2           = point2              ; 
    m_previousThreshold    = threshold           ;
    m_previousGuess        = guessGeometry       ;
    m_previousTopGeometry  = topGeometry ;
    
    m_localIntersections.clear();
    m_localIntersections.reserve( intersept.size()); 
    m_localIntersections.insert(m_localIntersections.begin(),
                                intersept.begin() , 
                                intersept.end());
  }
  catch ( const GaudiException& Exception ) 
    {
      m_prevPoint1           = Gaudi::XYZPoint() ;
      m_prevPoint2           = Gaudi::XYZPoint() ; 
      m_previousThreshold    = -1000.0      ; 
      m_previousGuess        = 0            ; 
      m_previousTopGeometry  = 0            ;
      m_localIntersections.clear()          ; 
      
      std::string message
        ("TransportSvc::intersection(...), exception caught; Params: ");
      {
        std::ostringstream ost;
        ost << "Point=" << point 
            << ",Vect=" << vect
            << ",Tick=" << tickMin << "/" << tickMax  
            << "Thrsh=" << threshold; 
        if( 0 != alternativeGeometry ) { ost << "A.Geo!=NULL" ; }  
        if( 0 != guessGeometry ) { ost << "G.Geo!=NULL" ; }  
        message += ost.str();             
        Assert( false , message , Exception );
      }
    }
  catch( ... ) 
    {
      m_prevPoint1           = Gaudi::XYZPoint() ;
      m_prevPoint2           = Gaudi::XYZPoint() ; 
      m_previousThreshold    = -1000.0      ; 
      m_previousGuess        = 0            ; 
      m_previousTopGeometry  = 0            ;
      m_localIntersections.clear()          ; 
      
      std::string message
        ("TransportSvc::intersection(...), unknown exception caught; Params: ");
      {
        std::ostringstream ost;
        ost << "Point=" << point 
            << ",Vect=" << vect
            << ",Tick=" << tickMin << "/" << tickMax  
            << "Thrsh=" << threshold; 
        if( 0 != alternativeGeometry ) { ost << "A.Geo!=NULL" ; }  
        if( 0 != guessGeometry ) { ost << "G.Geo!=NULL" ; }  
        message += ost.str();             
        Assert( false , message );
      }
    }

  return intersept.size() ;
};

// ============================================================================
// The End 
// ============================================================================
#endif  //  __DETDESC_TRANSPORTSVC_TRANSPORTSVCINTERSECTIONS_H__ 
// ============================================================================

---