/data/refman/pasoursint/CMSSW_4_1_8_patch9/src/Geometry/ForwardGeometry/src/CastorGeometry.cc

Go to the documentation of this file.

#include "Geometry/CaloGeometry/interface/CaloGenericDetId.h"
#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/ForwardGeometry/interface/CastorGeometry.h"
#include "Geometry/ForwardGeometry/interface/IdealCastorTrapezoid.h"
#include "CastorGeometryData.h"
#include <algorithm>

CastorGeometry::CastorGeometry() :
   theTopology( new CastorTopology ), 
   lastReqDet_(DetId::Detector(0)),
   lastReqSubdet_(0),
   m_ownsTopology ( true )
{
}

CastorGeometry::CastorGeometry( const CastorTopology* topology ) :
   theTopology(topology), 
   lastReqDet_(DetId::Detector(0)),
   lastReqSubdet_(0),
   m_ownsTopology ( false )
{
}


CastorGeometry::~CastorGeometry() 
{
   if( m_ownsTopology ) delete theTopology ;
}

const std::vector<DetId>& 
CastorGeometry::getValidDetIds( DetId::Detector det,
                                int             subdet ) const 
{
   const std::vector<DetId>& baseIds ( CaloSubdetectorGeometry::getValidDetIds() ) ;
   if( det    == DetId::Detector( 0 ) &&
       subdet == 0                        )
   {
      return baseIds ;
   }
   
   if( lastReqDet_    != det    ||
       lastReqSubdet_ != subdet    ) 
   {
      lastReqDet_     = det    ;
      lastReqSubdet_  = subdet ;
      m_validIds.clear();
      m_validIds.reserve( baseIds.size() ) ;
   }

   if( m_validIds.empty() ) 
   {
      for( unsigned int i ( 0 ) ; i != baseIds.size() ; ++i ) 
      {
         const DetId id ( baseIds[i] );
         if( id.det()      == det    &&
             id.subdetId() == subdet    )
         { 
            m_validIds.push_back( id ) ;
         }
      }
      std::sort(m_validIds.begin(),m_validIds.end());
   }
   return m_validIds;
}

/*  NOTE only draft implementation at the moment
    what about dead volumes?
*/

DetId 
CastorGeometry::getClosestCell(const GlobalPoint& r) const
{
   DetId returnId ( 0 ) ;
   const std::vector<DetId>& detIds ( getValidDetIds() ) ;
   for( std::vector<DetId>::const_iterator it ( detIds.begin() ) ;
        it != detIds.end(); ++it )
   {
      const CaloCellGeometry& cell ( *getGeometry( *it ) ) ;
      if( cell.inside( r ) )
      {
         returnId = *it ;
         break ;
      }
   }
   return returnId ;
}



unsigned int
CastorGeometry::alignmentTransformIndexLocal( const DetId& id )
{
   const CaloGenericDetId gid ( id ) ;

   assert( gid.isCastor() ) ;

   return 0 ;
}

unsigned int
CastorGeometry::alignmentTransformIndexGlobal( const DetId& id )
{
   return (unsigned int)DetId::Calo - 1 ;
}

std::vector<HepGeom::Point3D<double> > 
CastorGeometry::localCorners( const double* pv,
                              unsigned int  i,
                              HepGeom::Point3D<double> &   ref )
{
   return ( calogeom::IdealCastorTrapezoid::localCorners( pv, ref ) ) ;
}

CaloCellGeometry* 
CastorGeometry::newCell( const GlobalPoint& f1 ,
                         const GlobalPoint& f2 ,
                         const GlobalPoint& f3 ,
                         CaloCellGeometry::CornersMgr* mgr,
                         const double*      parm ,
                         const DetId&       detId   ) 
{
   const CaloGenericDetId cgid ( detId ) ;
   
   assert( cgid.isCastor() ) ;

   return ( new calogeom::IdealCastorTrapezoid( f1, mgr, parm ) ) ;
}