CastorHardcodeGeometryLoader Class Reference

#include <CastorHardcodeGeometryLoader.h>

Public Member Functions
	CastorHardcodeGeometryLoader ()
	CastorHardcodeGeometryLoader (const CastorTopology &ht)
virtual std::unique_ptr< CaloSubdetectorGeometry >	load (DetId::Detector det, int subdet)
std::unique_ptr< CaloSubdetectorGeometry >	load ()
virtual	~CastorHardcodeGeometryLoader ()

Private Member Functions
void	fill (HcalCastorDetId::Section section, CaloSubdetectorGeometry *cg)
void	init ()
void	makeCell (const HcalCastorDetId &detId, CaloSubdetectorGeometry *geom) const

Private Attributes
const CastorTopology *	extTopology
float	theEMSectiondX
float	theEMSectiondY
float	theEMSectiondZ
float	theHADSectiondX
float	theHADSectiondY
float	theHADSectiondZ
CastorTopology *	theTopology

Detailed Description

Definition at line 11 of file CastorHardcodeGeometryLoader.h.

Constructor & Destructor Documentation

CastorHardcodeGeometryLoader() [1/2]

CastorHardcodeGeometryLoader::CastorHardcodeGeometryLoader

(

)

Definition at line 12 of file CastorHardcodeGeometryLoader.cc.

References init().

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

CastorHardcodeGeometryLoader() [2/2]

CastorHardcodeGeometryLoader::CastorHardcodeGeometryLoader ( const CastorTopology &  ht )

explicit

Definition at line 17 of file CastorHardcodeGeometryLoader.cc.

References init().

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

~CastorHardcodeGeometryLoader()

virtual CastorHardcodeGeometryLoader::~CastorHardcodeGeometryLoader ( )

inlinevirtual

Definition at line 15 of file CastorHardcodeGeometryLoader.h.

References theTopology.

{ delete theTopology; };

Member Function Documentation

fill()

void CastorHardcodeGeometryLoader::fill ( HcalCastorDetId::Section  section, CaloSubdetectorGeometry *  cg  )

private

Definition at line 47 of file CastorHardcodeGeometryLoader.cc.

References extTopology, CastorTopology::firstCell(), relativeConstraints::geom, l1ctLayer2EG_cff::id, CastorGeometry::k_NumberOfParametersPerShape, CastorGeometry::k_NumberOfShapes, HcalCastorDetId::kNumberSectorsPerEnd, HcalCastorDetId::kSizeForDenseIndexing, CastorTopology::lastCell(), makeCell(), hgcalPlots::section, and CastorTopology::valid().

Referenced by load().

                                                                                                     {
  if (geom->cornersMgr() == nullptr)
    geom->allocateCorners(HcalCastorDetId::kSizeForDenseIndexing);
  if (geom->parMgr() == nullptr)
    geom->allocatePar(CastorGeometry::k_NumberOfShapes * CastorGeometry::k_NumberOfParametersPerShape,
                      CastorGeometry::k_NumberOfParametersPerShape);

  // start by making the new HcalDetIds
  std::vector<HcalCastorDetId> castorIds;

  const int firstCell(extTopology->firstCell(section));
  const int lastCell(extTopology->lastCell(section));

  for (int imodule(firstCell); imodule <= lastCell; ++imodule) {
    for (int isector(1); isector <= HcalCastorDetId::kNumberSectorsPerEnd; ++isector) {
      const HcalCastorDetId id(section, false, isector, imodule);
      if (extTopology->valid(id))
        castorIds.emplace_back(id);
    }
  }
  //   edm::LogInfo("CastorHardcodeGeometry")
  //      << "Number of Castor DetIds made: " << section
  //      << " " << castorIds.size();

  // for each new HcalCastorDetId, make a CaloCellGeometry

  for (std::vector<HcalCastorDetId>::const_iterator castorIdItr(castorIds.begin()); castorIdItr != castorIds.end();
       ++castorIdItr) {
    makeCell(*castorIdItr, geom);
  }
}

init()

void CastorHardcodeGeometryLoader::init ( void  )

private

Definition at line 22 of file CastorHardcodeGeometryLoader.cc.

References dXEMPlate, dXHADPlate, dYEMPlate, dYHADPlate, theEMSectiondX, theEMSectiondY, theEMSectiondZ, theHADSectiondX, theHADSectiondY, and theHADSectiondZ.

Referenced by CastorHardcodeGeometryLoader().

                                        {
  theEMSectiondX = 2. * dXEMPlate;
  theEMSectiondY = 2. * dYEMPlate;
  theEMSectiondZ = 101.0;
  theHADSectiondX = 2. * dXHADPlate;
  theHADSectiondY = 2. * dYHADPlate;
  theHADSectiondZ = 1212.;
}

load() [1/2]

std::unique_ptr< CaloSubdetectorGeometry > CastorHardcodeGeometryLoader::load ( DetId::Detector  det, int  subdet  )

virtual

Definition at line 31 of file CastorHardcodeGeometryLoader.cc.

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

Referenced by CastorHardcodeGeometryEP::produce().

                                                                                                     {
  std::unique_ptr<CaloSubdetectorGeometry> hg(new CastorGeometry(extTopology));
  if (subdet == HcalCastorDetId::SubdetectorId) {
    fill(HcalCastorDetId::EM, hg.get());
    fill(HcalCastorDetId::HAD, hg.get());
  }
  return hg;
}

load() [2/2]

std::unique_ptr< CaloSubdetectorGeometry > CastorHardcodeGeometryLoader::load

(

)

Definition at line 40 of file CastorHardcodeGeometryLoader.cc.

References HcalCastorDetId::EM, extTopology, fill(), and HcalCastorDetId::HAD.

                                                                          {
  std::unique_ptr<CaloSubdetectorGeometry> hg(new CastorGeometry(extTopology));
  fill(HcalCastorDetId::EM, hg.get());
  fill(HcalCastorDetId::HAD, hg.get());
  return hg;
}

makeCell()

void CastorHardcodeGeometryLoader::makeCell ( const HcalCastorDetId &  detId, CaloSubdetectorGeometry *  geom  ) const

private

Definition at line 79 of file CastorHardcodeGeometryLoader.cc.

References timingPdfMaker::can, funct::cos(), hcalRecHitTable_cff::detId, cuy::dh, HGC3DClusterGenMatchSelector_cfi::dR, PVValHelper::dx, PVValHelper::dy, PVValHelper::dz, HcalCastorDetId::EM, relativeConstraints::geom, CaloCellGeometry::getParmPtr(), CastorGeometry::k_NumberOfParametersPerShape, HcalCastorDetId::kNumberSectorsPerEnd, M_PI, callgraph::module, phi, hgcalPlots::section, nano_mu_digi_cff::sector, Validation_hcalonly_cfi::sign, funct::sin(), mathSSE::sqrt(), ecaldqm::zside(), and geometryCSVtoXML::zz.

Referenced by fill().

                                                                                                             {
  const double zside(1.0 * detId.zside());
  const HcalCastorDetId::Section section(detId.section());
  const int module(detId.module());
  const int isect(detId.sector());
  const double sector(1.0 * isect);

  // length units are cm

  const double sign(0 == isect % 2 ? -1 : 1);

  //********* HERE ARE HARDWIRED GEOMETRY NUMBERS ACTUALLY USED *****

  static const double an(atan(1.));  //angle of cant w.r.t. beam
  static const double can(cos(an));
  static const double san(sin(an));
  static const double dxlEM(1.55 / 2.);  //halflength of side near beam
  static const double dxhEM(5.73 / 2.);  //halflength of side away from beam
  static const double dhEM(14.26 / 2.);  //halflength of 2nd longest side
  static const double dR(0.1 + 2. * dhEM * san * dxlEM / (dxhEM - dxlEM));
  static const double dhHAD(19.88 / 2.);  //halflength of 2nd longest side

  static const double dxhHAD(dxhEM * (2. * dhHAD * san + dR) /
                             (2. * dhEM * san + dR));  //halflength of side away from beam
  static const double zm(1439.0);                      //z of start of EM
  static const double dzEM(5.45 / 2);                  // halflength in z of EM
  static const double dzHAD(10.075 / 2);               // halflength in z of HAD

  //*****************************************************************

  const double dxl(sign * dxlEM);  // same for EM and HAD

  const double dxh(sign * (section == HcalCastorDetId::EM ? dxhEM : dxhHAD));
  const double dh(section == HcalCastorDetId::EM ? dhEM : dhHAD);
  const double dz(section == HcalCastorDetId::EM ? dzEM : dzHAD);

  const double delz(dh * can);
  const double dy(dh * san);
  const double dx((dxl + dxh) / 2.);
  const double leg(dR + dy);
  const double len(sqrt(leg * leg + dx * dx));

  static const double dphi(2. * M_PI / (1.0 * HcalCastorDetId::kNumberSectorsPerEnd));

  const double fphi(atan(dx / (dR + dy)));

  const double phi(0 == isect % 2 ? (sector - 1.) * dphi - fphi : sector * dphi - fphi);

  const double sphi(sin(phi));
  const double cphi(cos(phi));

  const double xc(len * cphi);
  const double yc(len * sphi);
  const double zc(
      zside *
      (zm + delz + (module < 3 ? (1. * module - 1.0) * 2. * dzEM : 4. * dzEM + (1. * (module - 2) - 1) * 2. * dzHAD)));

  const GlobalPoint fc(xc, yc, zc);

  std::vector<CCGFloat> zz;
  zz.reserve(CastorGeometry::k_NumberOfParametersPerShape);
  zz.emplace_back(dxl);
  zz.emplace_back(dxh);
  zz.emplace_back(dh);
  zz.emplace_back(dz);
  zz.emplace_back(an);
  zz.emplace_back(dR);

  geom->newCell(fc, fc, fc, CaloCellGeometry::getParmPtr(zz, geom->parMgr(), geom->parVecVec()), detId);
}

Member Data Documentation

extTopology

const CastorTopology* CastorHardcodeGeometryLoader::extTopology

private

Definition at line 26 of file CastorHardcodeGeometryLoader.h.

Referenced by fill(), and load().

theEMSectiondX

float CastorHardcodeGeometryLoader::theEMSectiondX

private

Definition at line 28 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theEMSectiondY

float CastorHardcodeGeometryLoader::theEMSectiondY

private

Definition at line 29 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theEMSectiondZ

float CastorHardcodeGeometryLoader::theEMSectiondZ

private

Definition at line 30 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theHADSectiondX

float CastorHardcodeGeometryLoader::theHADSectiondX

private

Definition at line 31 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theHADSectiondY

float CastorHardcodeGeometryLoader::theHADSectiondY

private

Definition at line 32 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theHADSectiondZ

float CastorHardcodeGeometryLoader::theHADSectiondZ

private

Definition at line 33 of file CastorHardcodeGeometryLoader.h.

Referenced by init().

theTopology

CastorTopology* CastorHardcodeGeometryLoader::theTopology

private

Definition at line 25 of file CastorHardcodeGeometryLoader.h.

Referenced by ~CastorHardcodeGeometryLoader().