#include <ZdcHardcodeGeometryLoader.h>

Public Types
typedef CaloSubdetectorGeometry *	ReturnType

Public Member Functions
ReturnType	load ()

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

	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

float	theRPDSectiondX

float	theRPDSectiondY

float	theRPDSectiondZ

ZdcTopology *	theTopology

Detailed Description

Author: E. Garcia - UIC

Definition at line 16 of file ZdcHardcodeGeometryLoader.h.

Member Typedef Documentation

◆ ReturnType

typedef CaloSubdetectorGeometry* ZdcHardcodeGeometryLoader::ReturnType

Definition at line 18 of file ZdcHardcodeGeometryLoader.h.

Constructor & Destructor Documentation

◆ ZdcHardcodeGeometryLoader() [1/2]

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader ( )

Definition at line 10 of file ZdcHardcodeGeometryLoader.cc.

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

References init().

◆ ZdcHardcodeGeometryLoader() [2/2]

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader ( const ZdcTopology & ht )

explicit

Definition at line 14 of file ZdcHardcodeGeometryLoader.cc.

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

References init().

◆ ~ZdcHardcodeGeometryLoader()

virtual ZdcHardcodeGeometryLoader::~ZdcHardcodeGeometryLoader ( )

inlinevirtual

Definition at line 22 of file ZdcHardcodeGeometryLoader.h.

22 { delete theTopology; }

References theTopology.

Member Function Documentation

◆ fill()

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

private

Definition at line 40 of file ZdcHardcodeGeometryLoader.cc.

                                                                                  {
   // 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.emplace_back(id);
     id = HcalZDCDetId(section, false, ichannel);
     if (extTopology->valid(id))
       zdcIds.emplace_back(id);
   }
   if (geom->cornersMgr() == nullptr)
     geom->allocateCorners(HcalZDCDetId::kSizeForDenseIndexing);
   if (geom->parMgr() == nullptr)
     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);
   }
 }

References extTopology, ZdcTopology::firstCell(), relativeConstraints::geom, triggerObjects_cff::id, ZdcGeometry::k_NumberOfParametersPerShape, ZdcGeometry::k_NumberOfShapes, HcalZDCDetId::kSizeForDenseIndexing, ZdcTopology::lastCell(), makeCell(), hgcalPlots::section, and ZdcTopology::valid().

Referenced by load().

◆ init()

void ZdcHardcodeGeometryLoader::init ( void )

private

Definition at line 18 of file ZdcHardcodeGeometryLoader.cc.

18 {}

Referenced by ZdcHardcodeGeometryLoader().

◆ load() [1/2]

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load ( )

Definition at line 31 of file ZdcHardcodeGeometryLoader.cc.

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

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

◆ load() [2/2]

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

virtual

Definition at line 20 of file ZdcHardcodeGeometryLoader.cc.

                                                                                                {
   ReturnType hg(new ZdcGeometry(extTopology));
   if (subdet == HcalZDCDetId::SubdetectorId) {
     fill(HcalZDCDetId::EM, hg);
     fill(HcalZDCDetId::LUM, hg);
     fill(HcalZDCDetId::HAD, hg);
     //      fill(HcalZDCDetId::RPD ,hg ); // line will be needed in next version of code, but is commented out for now in order to ensure run1 compatability.
   }
   return hg;
 }

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

Referenced by ZdcHardcodeGeometryEP::produce().

◆ makeCell()

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

private

Definition at line 69 of file ZdcHardcodeGeometryLoader.cc.

                                                                                          {
   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 angRPD(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 dxRPD(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 dhRPD(6.25);  // to be updated when known
  
   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));
   static const double dzRPD(23.5);  // to be updated when known
  
   // 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 xOffRPD(0);
  
   static const double yOffEM(0);
   static const double yOffLUM(0);
   static const double yOffHAD(0);
   static const double yOffRPD(0);
  
   static const double zOffEM(-49.85 - 0.15);
   static const double zOffLUM(-39.555);
   static const double zOffHAD(-29.00);
   static const double zOffRPD(-39.555);
  
   //++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
  
   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 if (section == HcalZDCDetId::RPD) {
     dx = dxRPD;
     dh = dhRPD;
     dz = dzRPD;
     an = angRPD;
     x = zside * (x0 + xOffRPD);
     y = y0 + yOffRPD;
     z = zside * (z0 + zOffRPD + (channel - 1.0) * dzRPD * 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.emplace_back(an);
   zz.emplace_back(dx);
   zz.emplace_back(dy);
   zz.emplace_back(dz);
  
   geom->newCell(fc, fc, fc, CaloCellGeometry::getParmPtr(zz, geom->parMgr(), geom->parVecVec()), detId);
 }