#include <ZdcHardcodeGeometryLoader.h>

Public Types
typedef CaloSubdetectorGeometry *	ReturnType

Public Member Functions
virtual ReturnType	load (DetId::Detector det, int subdet)

ReturnType	load ()

	ZdcHardcodeGeometryLoader ()

	ZdcHardcodeGeometryLoader (const ZdcTopology &ht)

virtual	~ZdcHardcodeGeometryLoader ()

Private Member Functions
void	fill (HcalZDCDetId::Section section, CaloSubdetectorGeometry *cg)

void	init ()

void	makeCell (const HcalZDCDetId &detId, CaloSubdetectorGeometry *geom) const

Private Attributes
const ZdcTopology *	extTopology

float	theEMSectiondX

float	theEMSectiondY

float	theEMSectiondZ

float	theHADSectiondX

float	theHADSectiondY

float	theHADSectiondZ

float	theLUMSectiondX

float	theLUMSectiondY

float	theLUMSectiondZ

ZdcTopology *	theTopology

Detailed Description

Author: E. Garcia - UIC

Definition at line 17 of file ZdcHardcodeGeometryLoader.h.

Member Typedef Documentation

typedef CaloSubdetectorGeometry* ZdcHardcodeGeometryLoader::ReturnType

Definition at line 20 of file ZdcHardcodeGeometryLoader.h.

Constructor & Destructor Documentation

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader ( )

Definition at line 10 of file ZdcHardcodeGeometryLoader.cc.

References init().

                                                      :
    theTopology ( new ZdcTopology ) ,
    extTopology ( theTopology )
 {
   init();
 }

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader ( const ZdcTopology & ht )

explicit

Definition at line 17 of file ZdcHardcodeGeometryLoader.cc.

References init().

                                                                           : 
    theTopology( 0   ) ,
    extTopology( &ht ) 
 {
   init();
 }

virtual ZdcHardcodeGeometryLoader::~ZdcHardcodeGeometryLoader ( )

inlinevirtual

Definition at line 24 of file ZdcHardcodeGeometryLoader.h.

References theTopology.

24 { delete theTopology ; }

ZdcHardcodeGeometryLoader::theTopology

ZdcTopology * theTopology

Definition: ZdcHardcodeGeometryLoader.h:35

Member Function Documentation

void ZdcHardcodeGeometryLoader::fill	(	HcalZDCDetId::Section	section,
		CaloSubdetectorGeometry *	cg
	)

private

Definition at line 51 of file ZdcHardcodeGeometryLoader.cc.

References CaloSubdetectorGeometry::allocateCorners(), CaloSubdetectorGeometry::allocatePar(), CaloSubdetectorGeometry::cornersMgr(), extTopology, ZdcTopology::firstCell(), ZdcGeometry::k_NumberOfParametersPerShape, ZdcGeometry::k_NumberOfShapes, HcalZDCDetId::kSizeForDenseIndexing, ZdcTopology::lastCell(), makeCell(), CaloSubdetectorGeometry::parMgr(), and ZdcTopology::valid().

Referenced by load().

 {
   // start by making the new HcalDetIds
   std::vector<HcalZDCDetId> zdcIds;
   HcalZDCDetId id;
   int firstCell = extTopology->firstCell(section);
   int lastCell = extTopology->lastCell(section);
   for(int ichannel = firstCell; ichannel <= lastCell; ++ichannel) {
     id = HcalZDCDetId(section, true, ichannel);
     if(extTopology->valid(id)) zdcIds.push_back(id);
     id = HcalZDCDetId(section, false, ichannel);
     if(extTopology->valid(id)) zdcIds.push_back(id);
    }
   if( geom->cornersMgr() == 0 ) geom->allocateCorners( HcalZDCDetId::kSizeForDenseIndexing ) ;
   if( geom->parMgr()     == 0 ) geom->allocatePar( 
      ZdcGeometry::k_NumberOfParametersPerShape*ZdcGeometry::k_NumberOfShapes,
      ZdcGeometry::k_NumberOfParametersPerShape ) ;
 
   edm::LogInfo("ZdcHardcodeGeometry") << "Number of ZDC DetIds made: " << section << " " << zdcIds.size();
  
   // for each new HcalZdcDetId, make a CaloCellGeometry
 
  for(std::vector<HcalZDCDetId>::const_iterator zdcIdItr = zdcIds.begin();
      zdcIdItr != zdcIds.end(); ++zdcIdItr)
    {
       makeCell( *zdcIdItr, geom ) ;
    }
 }

void ZdcHardcodeGeometryLoader::init ( void )

private

Definition at line 24 of file ZdcHardcodeGeometryLoader.cc.

Referenced by ZdcHardcodeGeometryLoader().

25 {

26 }

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load	(	DetId::Detector	det,
		int	subdet
	)

virtual

Definition at line 29 of file ZdcHardcodeGeometryLoader.cc.

References HcalZDCDetId::EM, extTopology, fill(), HcalZDCDetId::HAD, HcalZDCDetId::LUM, and HcalZDCDetId::SubdetectorId.

Referenced by ZdcHardcodeGeometryEP::produce().

 {
    ReturnType hg(new ZdcGeometry( extTopology ) );
    if(subdet == HcalZDCDetId::SubdetectorId)
    {
       fill(HcalZDCDetId::EM  ,hg );
       fill(HcalZDCDetId::LUM ,hg );
       fill(HcalZDCDetId::HAD ,hg );
    }
    return hg;
 }

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load ( )

Definition at line 42 of file ZdcHardcodeGeometryLoader.cc.

References HcalZDCDetId::EM, extTopology, fill(), HcalZDCDetId::HAD, and HcalZDCDetId::LUM.

 {
    ReturnType hg(new ZdcGeometry( extTopology ) );
    fill(HcalZDCDetId::EM  ,hg );
    fill(HcalZDCDetId::LUM ,hg );
    fill(HcalZDCDetId::HAD ,hg );
    return hg;
 }

void ZdcHardcodeGeometryLoader::makeCell	(	const HcalZDCDetId &	detId,
		CaloSubdetectorGeometry *	geom
	)		const

private

Definition at line 82 of file ZdcHardcodeGeometryLoader.cc.

References HcalZDCDetId::channel(), funct::cos(), cuy::dh, HcalZDCDetId::EM, benchmark_cfg::fc, CaloCellGeometry::getParmPtr(), HcalZDCDetId::HAD, ZdcGeometry::k_NumberOfParametersPerShape, HcalZDCDetId::LUM, CaloSubdetectorGeometry::newCell(), CaloSubdetectorGeometry::parMgr(), CaloSubdetectorGeometry::parVecVec(), HcalZDCDetId::section(), x, detailsBasic3DVector::y, detailsBasic3DVector::z, and HcalZDCDetId::zside().

Referenced by fill().

 {
    double zside ( detId.zside() ) ;
 
    const HcalZDCDetId::Section section ( detId.section() ) ;
 
    const int channel ( detId.channel() ) ;
 
 //********* Here are all the hardcoded numbers you need to know, in **cm**
 //********* Most are from the zdc.xml and zdclum.xml files ******
 
    static const double x0 ( 0 ) ; // these 3 are for the "mother" volume
    static const double y0 ( 0 ) ;
    static const double z0 ( 14000 ) ;
 
    static const double angEM  ( 0 ) ; // the angles of front face wrt vertical
    static const double angLUM ( 0 ) ;
    static const double angHAD ( atan( 1. ) ) ; // this is 45 deg
 
    // these dimensions are **half**-sizes
 
    static const double dxHAD ( 4.8 ) ;
    static const double dxEM  ( dxHAD/5. ) ;
    static const double dxLUM ( 4.8 ) ; // to be updated when known
 
    static const double dhEM  ( 6.25 ) ;
    static const double dhLUM ( 6.25 ) ; // to be updated when known
    static const double dhHAD ( 6.25 ) ;
 
    static const double dzEM  ( 33.*0.15 ) ;
    static const double dzLUM ( 23.5 ) ; // to be updated when known
    static const double dzHAD ( 0.82*6./cos(angHAD) ) ;
 
    // these are not half-dimensions, they are offsets from nominal
    // for the center-of-front-face points
 
    static const double xOffEM  ( -4.*dxEM ) ; 
    static const double xOffLUM ( 0 ) ; 
    static const double xOffHAD ( 0 ) ; 
 
    static const double yOffEM  ( 0 ) ; 
    static const double yOffLUM ( 0 ) ; 
    static const double yOffHAD ( 0 ) ; 
 
    static const double zOffEM  ( -49.85  - 0.15 ) ; 
    static const double zOffLUM ( -39.555        ) ; 
    static const double zOffHAD ( -29.00         ) ; 
 
 //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
 
    double dx, dh, dz, x, y, z, an ;
 
    if( section==HcalZDCDetId::EM )
    {
       dx = dxEM ;
       dh = dhEM ;
       dz = dzEM ;
       an = angEM ;
       x  = zside*( x0 + xOffEM + ( channel - 1.0 )*dxEM*2. ) ;
       y  = y0 + yOffEM ;
       z  = zside*( z0 + zOffEM ) ;
    }
    else
    {
       if( section==HcalZDCDetId::LUM )
       {
      dx = dxLUM ;
      dh = dhLUM ;
      dz = dzLUM ;
      an = angLUM ;
      x  = zside*( x0 + xOffLUM ) ;
      y  = y0 + yOffLUM ;
      z  = zside*( z0 + zOffLUM + ( channel - 1.0 )*dzLUM*2. ) ;
       }
       else
       {
      assert( section==HcalZDCDetId::HAD ) ;
      dx = dxHAD ;
      dh = dhHAD ;
      dz = dzHAD ;
      an = angHAD ;
      x  = zside*( x0 + xOffHAD ) ;
      y  = y0 + yOffHAD ;
      z  = zside*( z0 + zOffHAD + ( channel - 1.0 )*dzHAD*2. ) ;
       }
    }
 
    const GlobalPoint fc ( x, y, z );
 
    const double dy ( dh*cos( an ) ) ;
 
    std::vector<CCGFloat> zz ;
    zz.reserve( ZdcGeometry::k_NumberOfParametersPerShape ) ;
    zz.push_back( an ) ;
    zz.push_back( dx ) ;
    zz.push_back( dy ) ;
    zz.push_back( dz ) ;
 
    geom->newCell( fc, fc, fc, 
           CaloCellGeometry::getParmPtr( zz, 
                         geom->parMgr(), 
                         geom->parVecVec() ),
           detId ) ;
 }