HcalParametersFromDD Class Reference

#include <HcalParametersFromDD.h>

Public Member Functions
bool	build (const DDCompactView *, HcalParameters &)
bool	build (const cms::DDCompactView &, HcalParameters &)
	HcalParametersFromDD ()=default

Private Member Functions
bool	build (HcalParameters &)
void	rescale (std::vector< double > &, const double s)

Detailed Description

Definition at line 9 of file HcalParametersFromDD.h.

Constructor & Destructor Documentation

HcalParametersFromDD()

HcalParametersFromDD::HcalParametersFromDD ( )

default

Member Function Documentation

build() [1/3]

bool HcalParametersFromDD::build	(	const DDCompactView *	cpv,
		HcalParameters &	php
	)

Definition at line 56 of file HcalParametersFromDD.cc.

References newFWLiteAna::build, dbl_to_int(), HcalParameters::dzVcal, HcalParameters::etagroup, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaRange, HcalParameters::etaTable, Exception, DDFilteredView::firstChild(), relativeConstraints::geom, HcalParameters::gparHF, HcalParameters::HBGains, HcalParameters::HBShift, HcalParameters::HEGains, HcalParameters::HEShift, HcalParameters::HFGains, HcalParameters::HFShift, mps_fire::i, HcalGeomParameters::k_ScaleFromDDDToG4, HcalParameters::LayerItem::layer, HcalParameters::Layer0Wt, HcalParameters::LayerItem::layerGroup, HcalParameters::layerGroupEtaRec, HcalParameters::layerGroupEtaSim, HcalParameters::maxDepth, DDFilteredView::mergedSpecifics(), HcalParameters::modHB, HcalParameters::modHE, HcalParameters::noff, convertSQLiteXML::ok, HcalParameters::phibin, HcalParameters::phigroup, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHO, HcalParameters::rTable, AlCaHLTBitMon_QueryRunRegistry::string, pfDeepBoostedJetPreprocessParams_cfi::sv, HcalParameters::topologyMode, findQualityFiles::v, and DDCompactView::vector().

Referenced by HcalParametersDBBuilder::beginRun().

                                                                              {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalParametersFromDD::build(const DDCompactView*, HcalParameters&) is called";
#endif
  //Special parameters at simulation level
  std::string attribute = "OnlyForHcalSimNumbering";
  DDSpecificsHasNamedValueFilter filter1{attribute};
  DDFilteredView fv1(*cpv, filter1);
  bool ok = fv1.firstChild();

  const int nEtaMax = 100;

  if (ok) {
    std::unique_ptr<HcalGeomParameters> geom = std::make_unique<HcalGeomParameters>();
    geom->loadGeometry(fv1, php);
    php.modHB = geom->getModHalfHBHE(0);
    php.modHE = geom->getModHalfHBHE(1);
    php.dzVcal = geom->getConstDzHF();
    geom->getConstRHO(php.rHO);

    php.phioff = cpv->vector("phioff");
    php.etaTable = cpv->vector("etaTable");
    php.rTable = cpv->vector("rTable");
    rescale(php.rTable, HcalGeomParameters::k_ScaleFromDDDToG4);
    php.phibin = cpv->vector("phibin");
    php.phitable = cpv->vector("phitable");
    for (unsigned int i = 1; i <= nEtaMax; ++i) {
      std::stringstream sstm;
      sstm << "layerGroupSimEta" << i;
      std::string tempName = sstm.str();
      auto const& v = cpv->vector(tempName);
      if (!v.empty()) {
        HcalParameters::LayerItem layerGroupEta;
        layerGroupEta.layer = i;
        layerGroupEta.layerGroup = dbl_to_int(v);
        php.layerGroupEtaSim.emplace_back(layerGroupEta);
      }
    }
    php.etaMin = dbl_to_int(cpv->vector("etaMin"));
    php.etaMax = dbl_to_int(cpv->vector("etaMax"));
    php.etaRange = cpv->vector("etaRange");
    php.gparHF = cpv->vector("gparHF");
    rescale(php.gparHF, HcalGeomParameters::k_ScaleFromDDDToG4);
    php.noff = dbl_to_int(cpv->vector("noff"));
    php.Layer0Wt = cpv->vector("Layer0Wt");
    php.HBGains = cpv->vector("HBGains");
    php.HBShift = dbl_to_int(cpv->vector("HBShift"));
    php.HEGains = cpv->vector("HEGains");
    php.HEShift = dbl_to_int(cpv->vector("HEShift"));
    php.HFGains = cpv->vector("HFGains");
    php.HFShift = dbl_to_int(cpv->vector("HFShift"));
    php.maxDepth = dbl_to_int(cpv->vector("MaxDepth"));
  } else {
    throw cms::Exception("HcalParametersFromDD") << "Not found " << attribute.c_str() << " but needed.";
  }
  //Special parameters at reconstruction level
  attribute = "OnlyForHcalRecNumbering";
  DDSpecificsHasNamedValueFilter filter2{attribute};
  DDFilteredView fv2(*cpv, filter2);
  ok = fv2.firstChild();
  if (ok) {
    DDsvalues_type sv(fv2.mergedSpecifics());
    int topoMode = getTopologyMode("TopologyMode", sv, true);
    int trigMode = getTopologyMode("TriggerMode", sv, false);
    php.topologyMode = ((trigMode & 0xFF) << 8) | (topoMode & 0xFF);
    php.etagroup = dbl_to_int(cpv->vector("etagroup"));
    php.phigroup = dbl_to_int(cpv->vector("phigroup"));
    for (unsigned int i = 1; i <= nEtaMax; ++i) {
      std::stringstream sstm;
      sstm << "layerGroupRecEta" << i;
      std::string tempName = sstm.str();
      auto const& v = cpv->vector(tempName);
      if (!v.empty()) {
        HcalParameters::LayerItem layerGroupEta;
        layerGroupEta.layer = i;
        layerGroupEta.layerGroup = dbl_to_int(v);
        php.layerGroupEtaRec.emplace_back(layerGroupEta);
      }
    }
  } else {
    throw cms::Exception("HcalParametersFromDD") << "Not found " << attribute.c_str() << " but needed.";
  }

  return build(php);
}

build() [2/3]

bool HcalParametersFromDD::build	(	const cms::DDCompactView &	cpv,
		HcalParameters &	php
	)

Definition at line 142 of file HcalParametersFromDD.cc.

References newFWLiteAna::build, HcalParameters::dzVcal, HcalParameters::etagroup, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaRange, HcalParameters::etaTable, ALCARECOTkAlBeamHalo_cff::filter, relativeConstraints::geom, cms::DDCompactView::get(), cms::DDFilteredView::get(), cms::DDCompactView::getVector(), HcalParameters::gparHF, HcalParameters::HBGains, HcalParameters::HBShift, HcalParameters::HEGains, HcalParameters::HEShift, HcalParameters::HFGains, HcalParameters::HFShift, mps_fire::i, HcalGeomParameters::k_ScaleFromDD4hepToG4, HcalParameters::LayerItem::layer, HcalParameters::Layer0Wt, HcalParameters::LayerItem::layerGroup, HcalParameters::layerGroupEtaRec, HcalParameters::layerGroupEtaSim, HcalParameters::maxDepth, HcalParameters::modHB, HcalParameters::modHE, HcalParameters::noff, HcalParameters::phibin, HcalParameters::phigroup, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rHO, HcalParameters::rTable, AlCaHLTBitMon_QueryRunRegistry::string, pfDeepBoostedJetPreprocessParams_cfi::sv, to_string(), and HcalParameters::topologyMode.

                                                                                 {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HCalGeom") << "HcalParametersFromDD::build(const cms::DDCompactView*, HcalParameters&) is called";
#endif
  const int nEtaMax = 100;

  std::unique_ptr<HcalGeomParameters> geom = std::make_unique<HcalGeomParameters>();
  geom->loadGeometry(cpv, php);
  php.modHB = geom->getModHalfHBHE(0);
  php.modHE = geom->getModHalfHBHE(1);
  php.dzVcal = geom->getConstDzHF();
  geom->getConstRHO(php.rHO);

  php.phioff = cpv.get<std::vector<double>>("phioff");
  php.etaTable = cpv.get<std::vector<double>>("etaTable");
  php.rTable = cpv.get<std::vector<double>>("rTable");
  php.phibin = cpv.get<std::vector<double>>("phibin");
  php.phitable = cpv.get<std::vector<double>>("phitable");
  php.etaMin = cpv.getVector<int>("etaMin");
  php.etaMax = cpv.getVector<int>("etaMax");
  php.etaRange = cpv.get<std::vector<double>>("etaRange");
  php.gparHF = cpv.get<std::vector<double>>("gparHF");
  php.noff = cpv.getVector<int>("noff");
  php.Layer0Wt = cpv.get<std::vector<double>>("Layer0Wt");
  php.HBGains = cpv.get<std::vector<double>>("HBGains");
  php.HBShift = cpv.getVector<int>("HBShift");
  php.HEGains = cpv.get<std::vector<double>>("HEGains");
  php.HEShift = cpv.getVector<int>("HEShift");
  php.HFGains = cpv.get<std::vector<double>>("HFGains");
  php.HFShift = cpv.getVector<int>("HFShift");
  php.maxDepth = cpv.getVector<int>("MaxDepth");

  rescale(php.rTable, HcalGeomParameters::k_ScaleFromDD4hepToG4);
  rescale(php.gparHF, HcalGeomParameters::k_ScaleFromDD4hepToG4);
  for (unsigned int i = 1; i <= nEtaMax; ++i) {
    HcalParameters::LayerItem layerGroupEta;
    layerGroupEta.layer = i;
    layerGroupEta.layerGroup = cpv.getVector<int>(std::string("layerGroupSimEta") + std::to_string(i));
    if (!layerGroupEta.layerGroup.empty()) {
      php.layerGroupEtaSim.emplace_back(layerGroupEta);
    }
  }

  //Special parameters at reconstruction level
  const cms::DDFilter filter("OnlyForHcalRecNumbering", "HCAL");
  cms::DDFilteredView fv(cpv, filter);

  std::vector<std::string> tempS = fv.get<std::vector<std::string>>("hcal", "TopologyMode");
  std::string sv = (!tempS.empty()) ? tempS[0] : "HcalTopologyMode::SLHC";
  int topoMode = getTopologyMode(sv, true);
  tempS = fv.get<std::vector<std::string>>("hcal", "TriggerMode");
  sv = (!tempS.empty()) ? tempS[0] : "HcalTopologyMode::TriggerMode_2021";
  int trigMode = getTopologyMode(sv, false);
  php.topologyMode = ((trigMode & 0xFF) << 8) | (topoMode & 0xFF);

  php.etagroup = cpv.getVector<int>("etagroup");
  php.phigroup = cpv.getVector<int>("phigroup");

  for (unsigned int i = 1; i <= nEtaMax; ++i) {
    HcalParameters::LayerItem layerGroupEta;
    layerGroupEta.layer = i;
    layerGroupEta.layerGroup = cpv.getVector<int>(std::string("layerGroupRecEta") + std::to_string(i));
    if (!layerGroupEta.layerGroup.empty()) {
      php.layerGroupEtaRec.emplace_back(layerGroupEta);
    }
  }

  return build(php);
}

build() [3/3]

bool HcalParametersFromDD::build ( HcalParameters &  php )

private

Definition at line 212 of file HcalParametersFromDD.cc.

References angle_units::operators::convertMmToCm(), angle_units::operators::convertRadToDeg(), TauDecayModes::dec, HcalParameters::etagroup, HcalParameters::etaMax, HcalParameters::etaMin, HcalParameters::etaRange, HcalParameters::etaTable, HcalParameters::etaTableHF, HcalParameters::gparHF, HcalParameters::HBGains, HcalParameters::HBShift, HcalParameters::HEGains, HcalParameters::HEShift, HcalParameters::HFGains, HcalParameters::HFShift, mps_fire::i, ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::it, isotrackApplyRegressor::k, HcalParameters::Layer0Wt, HcalParameters::layerGroupEtaRec, HcalParameters::layerGroupEtaSim, CrabHelper::log, HcalParameters::maxDepth, HcalParameters::modHB, HcalParameters::modHE, HcalParameters::noff, HcalParameters::phibin, HcalParameters::phigroup, HcalParameters::phioff, HcalParameters::phitable, HcalParameters::rTable, funct::tan(), and HcalParameters::topologyMode.

                                                    {
  php.etaMin[0] = 1;
  if (php.etaMax[1] >= php.etaMin[1])
    php.etaMax[1] = static_cast<int>(php.etaTable.size()) - 1;
  php.etaMax[2] = php.etaMin[2] + static_cast<int>(php.rTable.size()) - 2;

  for (unsigned int i = 0; i < php.rTable.size(); ++i) {
    unsigned int k = php.rTable.size() - i - 1;
    php.etaTableHF.emplace_back(-log(tan(0.5 * atan(php.rTable[k] / php.gparHF[4]))));
  }

#ifdef EDM_ML_DEBUG
  int i(0);
  std::stringstream ss0;
  ss0 << "HcalParametersFromDD: MaxDepth[" << php.maxDepth.size() << "]: ";
  for (const auto& it : php.maxDepth)
    ss0 << it << ", ";
  edm::LogVerbatim("HCalGeom") << ss0.str();
  std::stringstream ss1;
  ss1 << "HcalParametersFromDD: ModHB [" << php.modHB.size() << "]: ";
  for (const auto& it : php.modHB)
    ss1 << it << ", ";
  edm::LogVerbatim("HCalGeom") << ss1.str();
  std::stringstream ss2;
  ss2 << "HcalParametersFromDD: ModHE [" << php.modHE.size() << "]: ";
  for (const auto& it : php.modHE)
    ss2 << it << ", ";
  edm::LogVerbatim("HCalGeom") << ss2.str();
  std::stringstream ss3;
  ss3 << "HcalParametersFromDD: " << php.phioff.size() << " phioff values:";
  std::vector<double>::const_iterator it;
  for (it = php.phioff.begin(), i = 0; it != php.phioff.end(); ++it)
    ss3 << " [" << ++i << "] = " << convertRadToDeg(*it);
  edm::LogVerbatim("HCalGeom") << ss3.str();
  std::stringstream ss4;
  ss4 << "HcalParametersFromDD: " << php.etaTable.size() << " entries for etaTable:";
  for (it = php.etaTable.begin(), i = 0; it != php.etaTable.end(); ++it)
    ss4 << " [" << ++i << "] = " << (*it);
  edm::LogVerbatim("HCalGeom") << ss4.str();
  std::stringstream ss5;
  ss5 << "HcalParametersFromDD: " << php.rTable.size() << " entries for rTable:";
  for (it = php.rTable.begin(), i = 0; it != php.rTable.end(); ++it)
    ss5 << " [" << ++i << "] = " << convertMmToCm(*it);
  edm::LogVerbatim("HCalGeom") << ss5.str();
  std::stringstream ss6;
  ss6 << "HcalParametersFromDD: " << php.phibin.size() << " phibin values:";
  for (it = php.phibin.begin(), i = 0; it != php.phibin.end(); ++it)
    ss6 << " [" << ++i << "] = " << convertRadToDeg(*it);
  edm::LogVerbatim("HCalGeom") << ss6.str();
  std::stringstream ss7;
  ss7 << "HcalParametersFromDD: " << php.phitable.size() << " phitable values:";
  for (it = php.phitable.begin(), i = 0; it != php.phitable.end(); ++it)
    ss7 << " [" << ++i << "] = " << convertRadToDeg(*it);
  edm::LogVerbatim("HCalGeom") << ss7.str();
  edm::LogVerbatim("HCalGeom") << "HcalParametersFromDD: " << php.layerGroupEtaSim.size() << " layerGroupEtaSim blocks"
                               << std::endl;
  std::vector<int>::const_iterator kt;
  for (unsigned int k = 0; k < php.layerGroupEtaSim.size(); ++k) {
    std::stringstream ss8;
    ss8 << "layerGroupEtaSim[" << k << "] Layer " << php.layerGroupEtaSim[k].layer;
    for (kt = php.layerGroupEtaSim[k].layerGroup.begin(), i = 0; kt != php.layerGroupEtaSim[k].layerGroup.end(); ++kt)
      ss8 << " " << ++i << ":" << (*kt);
    edm::LogVerbatim("HCalGeom") << ss8.str();
  }
  std::stringstream ss8;
  ss8 << "HcalParametersFromDD: " << php.etaMin.size() << " etaMin values:";
  for (kt = php.etaMin.begin(), i = 0; kt != php.etaMin.end(); ++kt)
    ss8 << " [" << ++i << "] = " << (*kt);
  edm::LogVerbatim("HCalGeom") << ss8.str();
  std::stringstream ss9;
  ss9 << "HcalParametersFromDD: " << php.etaMax.size() << " etaMax values:";
  for (kt = php.etaMax.begin(), i = 0; kt != php.etaMax.end(); ++kt)
    ss9 << " [" << ++i << "] = " << (*kt);
  edm::LogVerbatim("HCalGeom") << ss9.str();
  std::stringstream ss10;
  ss10 << "HcalParametersFromDD: " << php.etaRange.size() << " etaRange values:";
  for (it = php.etaRange.begin(), i = 0; it != php.etaRange.end(); ++it)
    ss10 << " [" << ++i << "] = " << (*it);
  edm::LogVerbatim("HCalGeom") << ss10.str();
  std::stringstream ss11;
  ss11 << "HcalParametersFromDD: " << php.gparHF.size() << " gparHF values:";
  for (it = php.gparHF.begin(), i = 0; it != php.gparHF.end(); ++it)
    ss11 << " [" << ++i << "] = " << convertMmToCm(*it);
  edm::LogVerbatim("HCalGeom") << ss11.str();
  std::stringstream ss12;
  ss12 << "HcalParametersFromDD: " << php.noff.size() << " noff values:";
  for (kt = php.noff.begin(), i = 0; kt != php.noff.end(); ++kt)
    ss12 << " [" << ++i << "] = " << (*kt);
  edm::LogVerbatim("HCalGeom") << ss12.str();
  std::stringstream ss13;
  ss13 << "HcalParametersFromDD: " << php.Layer0Wt.size() << " Layer0Wt values:";
  for (it = php.Layer0Wt.begin(), i = 0; it != php.Layer0Wt.end(); ++it)
    ss13 << " [" << ++i << "] = " << (*it);
  edm::LogVerbatim("HCalGeom") << ss13.str();
  std::stringstream ss14;
  ss14 << "HcalParametersFromDD: " << php.HBGains.size() << " Shift/Gains values for HB:";
  for (unsigned k = 0; k < php.HBGains.size(); ++k)
    ss14 << " [" << k << "] = " << php.HBShift[k] << ":" << php.HBGains[k];
  edm::LogVerbatim("HCalGeom") << ss14.str();
  std::stringstream ss15;
  ss15 << "HcalParametersFromDD: " << php.HEGains.size() << " Shift/Gains values for HE:";
  for (unsigned k = 0; k < php.HEGains.size(); ++k)
    ss15 << " [" << k << "] = " << php.HEShift[k] << ":" << php.HEGains[k];
  edm::LogVerbatim("HCalGeom") << ss15.str();
  std::stringstream ss16;
  ss16 << "HcalParametersFromDD: " << php.HFGains.size() << " Shift/Gains values for HF:";
  for (unsigned k = 0; k < php.HFGains.size(); ++k)
    ss16 << " [" << k << "] = " << php.HFShift[k] << ":" << php.HFGains[k];
  edm::LogVerbatim("HCalGeom") << ss16.str();
  std::stringstream ss17;
  ss17 << "HcalParametersFromDD: " << php.etagroup.size() << " etagroup values:";
  for (kt = php.etagroup.begin(), i = 0; kt != php.etagroup.end(); ++kt)
    ss17 << " [" << ++i << "] = " << (*kt);
  edm::LogVerbatim("HCalGeom") << ss17.str();
  std::stringstream ss18;
  ss18 << "HcalParametersFromDD: " << php.phigroup.size() << " phigroup values:";
  for (kt = php.phigroup.begin(), i = 0; kt != php.phigroup.end(); ++kt)
    ss18 << " [" << ++i << "] = " << (*kt);
  edm::LogVerbatim("HCalGeom") << ss18.str();
  edm::LogVerbatim("HCalGeom") << "HcalParametersFromDD: " << php.layerGroupEtaRec.size() << " layerGroupEtaRec blocks";
  for (unsigned int k = 0; k < php.layerGroupEtaRec.size(); ++k) {
    std::stringstream ss19;
    ss19 << "layerGroupEtaRec[" << k << "] Layer " << php.layerGroupEtaRec[k].layer;
    for (kt = php.layerGroupEtaRec[k].layerGroup.begin(), i = 0; kt != php.layerGroupEtaRec[k].layerGroup.end(); ++kt)
      ss19 << " " << ++i << ":" << (*kt);
    edm::LogVerbatim("HCalGeom") << ss19.str();
  }
  edm::LogVerbatim("HCalGeom") << "HcalParametersFromDD: (topology|trigger)Mode " << std::hex << php.topologyMode
                               << std::dec;
#endif

  return true;
}

rescale()

void HcalParametersFromDD::rescale ( std::vector< double > &  v, const double  s  )

private

Definition at line 346 of file HcalParametersFromDD.cc.

References create_idmaps::n, alignCSCRings::s, and findQualityFiles::v.

                                                                       {
  std::for_each(v.begin(), v.end(), [s](double& n) { n *= s; });
}