#include <ZdcHardcodeGeometryLoader.h>

Public Types
typedef CaloSubdetectorGeometry *	ReturnType

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

ReturnType	load ()

void	setAddRPD (bool flag)

	ZdcHardcodeGeometryLoader (const ZdcTopology &ht)

	~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

bool	m_zdcAddRPD

float	theEMSectiondX

float	theEMSectiondY

float	theEMSectiondZ

float	theHADSectiondX

float	theHADSectiondY

float	theHADSectiondZ

float	theLUMSectiondX

float	theLUMSectiondY

float	theLUMSectiondZ

float	theRPDSectiondX

float	theRPDSectiondY

float	theRPDSectiondZ

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()

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader ( const ZdcTopology & ht )

explicit

Definition at line 10 of file ZdcHardcodeGeometryLoader.cc.

References init().

10 : extTopology(&ht) { init(); }

ZdcHardcodeGeometryLoader::init

void init()

Definition: ZdcHardcodeGeometryLoader.cc:12

ZdcHardcodeGeometryLoader::extTopology

const ZdcTopology * extTopology

Definition: ZdcHardcodeGeometryLoader.h:33

◆ ~ZdcHardcodeGeometryLoader()

ZdcHardcodeGeometryLoader::~ZdcHardcodeGeometryLoader ( )

inline

Definition at line 21 of file ZdcHardcodeGeometryLoader.h.

21 {}

Member Function Documentation

◆ fill()

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

private

Definition at line 38 of file ZdcHardcodeGeometryLoader.cc.

References extTopology, ZdcTopology::firstCell(), relativeConstraints::geom, EcalPhiSymFlatTableProducers_cfi::id, ZdcGeometry::k_NumberOfParametersPerShape, ZdcGeometry::k_NumberOfShapes, ZdcTopology::kSizeForDenseIndexing(), ZdcTopology::lastCell(), makeCell(), hgcalPlots::section, 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.emplace_back(id);
     id = HcalZDCDetId(section, false, ichannel);
     if (extTopology->valid(id))
       zdcIds.emplace_back(id);
   }
   if (geom->cornersMgr() == nullptr)
     geom->allocateCorners(extTopology->kSizeForDenseIndexing());
   if (geom->parMgr() == nullptr)
     geom->allocatePar(ZdcGeometry::k_NumberOfParametersPerShape * ZdcGeometry::k_NumberOfShapes,
                       ZdcGeometry::k_NumberOfParametersPerShape);
 
   edm::LogVerbatim("ZdcGeometry") << "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);
   }
 }

◆ init()

void ZdcHardcodeGeometryLoader::init ( void )

private

Definition at line 12 of file ZdcHardcodeGeometryLoader.cc.

Referenced by ZdcHardcodeGeometryLoader().

12 {}

◆ load() [1/2]

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

Definition at line 14 of file ZdcHardcodeGeometryLoader.cc.

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

                                                                                                {
   ReturnType hg(new ZdcGeometry(extTopology));
   if (subdet == HcalZDCDetId::SubdetectorId) {
     fill(HcalZDCDetId::EM, hg);
     fill(HcalZDCDetId::LUM, hg);
     fill(HcalZDCDetId::HAD, hg);
     if (m_zdcAddRPD)
       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;
 }

◆ load() [2/2]

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load ( )

Definition at line 28 of file ZdcHardcodeGeometryLoader.cc.

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

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

◆ makeCell()

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

private

Definition at line 67 of file ZdcHardcodeGeometryLoader.cc.

References cms::cuda::assert(), funct::cos(), hcalRecHitTable_cff::detId, cuy::dh, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, HcalZDCDetId::EM, relativeConstraints::geom, CaloCellGeometry::getParmPtr(), HcalZDCDetId::HAD, ZdcGeometry::k_NumberOfParametersPerShape, HcalZDCDetId::LUM, HcalZDCDetId::RPD, hgcalPlots::section, x, y, z, ppsFastLocalSimulation_cfi::z0, ecaldqm::zside(), and geometryCSVtoXML::zz.

Referenced by fill().

                                                                                          {
   double zside(detId.zside());
 
   const HcalZDCDetId::Section section(detId.section());
 
   const int channel(detId.channel());
 
   edm::LogVerbatim("ZdcGeometry") << "ZDCGeometry::Cell: " << detId << " Section " << section << " channel " << 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);
 }