ZdcHardcodeGeometryLoader.cc

Go to the documentation of this file.

#include "Geometry/CaloGeometry/interface/CaloCellGeometry.h"
#include "Geometry/ForwardGeometry/interface/ZdcHardcodeGeometryLoader.h"
#include "Geometry/ForwardGeometry/interface/IdealZDCTrapezoid.h"
#include "Geometry/ForwardGeometry/interface/ZdcGeometry.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <algorithm>

typedef CaloCellGeometry::CCGFloat CCGFloat;

ZdcHardcodeGeometryLoader::ZdcHardcodeGeometryLoader(const ZdcTopology& ht) : extTopology(&ht) { init(); }

void ZdcHardcodeGeometryLoader::init() {}

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load(DetId::Detector /*det*/, int subdet) {
  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;
}

ZdcHardcodeGeometryLoader::ReturnType ZdcHardcodeGeometryLoader::load() {
  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;
}

void ZdcHardcodeGeometryLoader::fill(HcalZDCDetId::Section section, ReturnType geom) {
  // 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);
  }
}

void ZdcHardcodeGeometryLoader::makeCell(const HcalZDCDetId& detId, ReturnType geom) const {
  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);
}