/data/refman/pasoursint/CMSSW_4_4_5_patch3/src/Geometry/ForwardGeometry/src/IdealZDCTrapezoid.cc

Go to the documentation of this file.

#include "Geometry/ForwardGeometry/interface/IdealZDCTrapezoid.h"
#include <math.h>

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

IdealZDCTrapezoid::IdealZDCTrapezoid() :
   CaloCellGeometry()
{
}

IdealZDCTrapezoid::IdealZDCTrapezoid( const IdealZDCTrapezoid& idzt ) 
{
   *this = idzt ;
}

IdealZDCTrapezoid& 
IdealZDCTrapezoid::operator=( const IdealZDCTrapezoid& idzt ) 
{
   if( &idzt != this ) CaloCellGeometry::operator=( idzt ) ;
   return *this ;
}

IdealZDCTrapezoid::IdealZDCTrapezoid( const GlobalPoint& faceCenter,
                                      const CornersMgr*  mgr       ,
                                      const CCGFloat*    parm        ) :  
   CaloCellGeometry ( faceCenter, mgr, parm )  
{
}
         
IdealZDCTrapezoid::~IdealZDCTrapezoid() {}

const CCGFloat 
IdealZDCTrapezoid::an() const { return param()[0] ; }

const CCGFloat 
IdealZDCTrapezoid::dx() const 
{
   return param()[1] ; 
}

const CCGFloat 
IdealZDCTrapezoid::dy() const 
{
   return param()[2] ; 
}

const CCGFloat 
IdealZDCTrapezoid::dz() const 
{
   return param()[3] ; 
}

const CCGFloat 
IdealZDCTrapezoid::ta() const 
{
   return tan( an() ) ; 
}

const CCGFloat 
IdealZDCTrapezoid::dt() const 
{
   return dy()*ta() ; 
}

void 
IdealZDCTrapezoid::vocalCorners( Pt3DVec&        vec ,
                                 const CCGFloat* pv  ,
                                 Pt3D&           ref  ) const 
{
   localCorners( vec, pv, ref ) ; 
}

void
IdealZDCTrapezoid::localCorners( Pt3DVec&        lc  ,
                                 const CCGFloat* pv  ,
                                 Pt3D&           ref   )
{
   assert( 8 == lc.size() ) ;
   assert( 0 != pv ) ;

   const CCGFloat an ( pv[0] ) ;
   const CCGFloat dx ( pv[1] ) ;
   const CCGFloat dy ( pv[2] ) ;
   const CCGFloat dz ( pv[3] ) ;
   const CCGFloat ta ( tan( an ) ) ;
   const CCGFloat dt ( dy*ta ) ;

   lc[0] = Pt3D ( -dx, -dy, +dz+dt ) ;
   lc[1] = Pt3D ( -dx, +dy, +dz-dt ) ;
   lc[2] = Pt3D ( +dx, +dy, +dz-dt ) ;
   lc[3] = Pt3D ( +dx, -dy, +dz+dt ) ;
   lc[4] = Pt3D ( -dx, -dy, -dz+dt ) ;
   lc[5] = Pt3D ( -dx, +dy, -dz-dt ) ;
   lc[6] = Pt3D ( +dx, +dy, -dz-dt ) ;
   lc[7] = Pt3D ( +dx, -dy, -dz+dt ) ;

   ref   = 0.25*( lc[0] + lc[1] + lc[2] + lc[3] ) ;
}

const CaloCellGeometry::CornersVec& 
IdealZDCTrapezoid::getCorners() const 
{
   const CornersVec& co ( CaloCellGeometry::getCorners() ) ;
   if( co.uninitialized() ) 
   {
      CaloCellGeometry::CornersVec& corners ( setCorners() ) ;
      const GlobalPoint& p ( getPosition() ) ;
      const CCGFloat zsign ( 0 < p.z() ? 1. : -1. ) ;
      const Pt3D  gf ( p.x(), p.y(), p.z() ) ;

      Pt3D  lf ;
      Pt3DVec lc ( 8, Pt3D(0,0,0) ) ;
      localCorners( lc, param(), lf ) ;
      const Pt3D  lb ( lf.x() , lf.y() , lf.z() - 2.*dz() ) ;
      const Pt3D  ls ( lf.x() - dx(), lf.y(), lf.z() ) ;
      
      const Pt3D  gb ( gf.x() , gf.y() , gf.z() + 2.*zsign*dz() ) ;

      const Pt3D  gs ( gf.x() - zsign*dx(),
                       gf.y() ,
                       gf.z()         ) ;

      const HepGeom::Transform3D tr ( lf, lb, ls,
                                      gf, gb, gs ) ;

      for( unsigned int i ( 0 ) ; i != 8 ; ++i )
      {
         const Pt3D  gl ( tr*lc[i] ) ;
         corners[i] = GlobalPoint( gl.x(), gl.y(), gl.z() ) ;
      }
   }
   return co ;
}

std::ostream& operator<<( std::ostream& s, const IdealZDCTrapezoid& cell ) 
{
   s << "Center: " <<  cell.getPosition() << std::endl ;
   s << "TiltAngle = " << cell.an()*180./M_PI << " deg, dx = " 
     << cell.dx() 
     << ", dy = " << cell.dy() << ", dz = " << cell.dz() << std::endl ;
   return s;
}