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>

typedef CaloCellGeometry::CCGFloat CCGFloat ;
typedef CaloCellGeometry::Pt3D     Pt3D     ;
typedef CaloCellGeometry::Pt3DVec  Pt3DVec  ;

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

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


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

DetId 
CastorGeometry::getClosestCell(const GlobalPoint& r) const
{
  DetId returnId ( 0 ) ;
  const std::vector<DetId>& detIds ( getValidDetIds() ) ;
  for(auto detId : detIds)  {
    auto cell = getGeometry(detId);
    if( nullptr != cell &&
    cell->inside( r ) ) {
      returnId = detId ;
      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 ;
}

void
CastorGeometry::localCorners( Pt3DVec&        lc ,
                  const CCGFloat* pv ,
                  unsigned int  /*i*/,
                  Pt3D&           ref )
{
   IdealCastorTrapezoid::localCorners( lc, pv, ref ) ;
}

void
CastorGeometry::newCell( const GlobalPoint& f1 ,
             const GlobalPoint& /*f2*/ ,
             const GlobalPoint& /*f3*/ ,
             const CCGFloat*    parm ,
             const DetId&       detId   ) 
{
   const CaloGenericDetId cgid ( detId ) ;
   
   assert( cgid.isCastor() ) ;

   const unsigned int di ( cgid.denseIndex() ) ;

   m_cellVec[ di ] = IdealCastorTrapezoid( f1, cornersMgr(), parm ) ;
   addValidID( detId ) ;
}

const CaloCellGeometry* CastorGeometry::getGeometryRawPtr(uint32_t index) const {
  // Modify the RawPtr class
  const CaloCellGeometry* cell(&m_cellVec[index]);
  return (m_cellVec.size() < index ||
      nullptr == cell->param() ? nullptr : cell);
}