DDHGCalModuleAlgo Class Reference

Inheritance diagram for DDHGCalModuleAlgo:

Public Member Functions
	DDHGCalModuleAlgo ()
void	execute (DDCompactView &cpv) override
void	initialize (const DDNumericArguments &nArgs, const DDVectorArguments &vArgs, const DDMapArguments &mArgs, const DDStringArguments &sArgs, const DDStringVectorArguments &vsArgs) override

Protected Member Functions
void	constructLayers (const DDLogicalPart &, DDCompactView &cpv)
void	positionSensitive (DDLogicalPart &glog, double rin, double rout, DDCompactView &cpv)
double	rMax (double z)

Private Attributes
std::unordered_set< int >	copies_
std::vector< int >	copyNumber_
std::string	idNameSpace_
std::vector< int >	layers_
std::vector< int >	layerSense_
std::vector< double >	layerThick_
std::vector< int >	layerType_
std::vector< std::string >	materials_
std::vector< std::string >	names_
double	rMaxFine_
std::vector< double >	rMaxFront_
int	sectors_
std::vector< double >	slopeB_
std::vector< double >	slopeT_
std::vector< double >	thick_
std::vector< std::string >	wafer_
double	waferGap_
double	waferW_
std::vector< double >	zFront_
double	zMinBlock_

Static Private Attributes
static constexpr double	tol_ = 0.00001

Detailed Description

Definition at line 32 of file DDHGCalModuleAlgo.cc.

Constructor & Destructor Documentation

DDHGCalModuleAlgo()

DDHGCalModuleAlgo::DDHGCalModuleAlgo

(

)

Definition at line 74 of file DDHGCalModuleAlgo.cc.

                                     {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: Creating an instance";
#endif
}

Member Function Documentation

constructLayers()

void DDHGCalModuleAlgo::constructLayers ( const DDLogicalPart &  module, DDCompactView &  cpv  )

protected

Definition at line 163 of file DDHGCalModuleAlgo.cc.

References isotrackApplyRegressor::alpha, angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, funct::cos(), DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, mps_fire::i, cuy::ii, isotrackApplyRegressor::k, DDBase< N, C >::name(), mergeVDriftHistosByStation::name, DDSolidFactory::polyhedra(), DDCompactView::position(), photonAnalyzer_cfi::rMax, makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, to_string(), DDSolidFactory::tubs(), and geometryCSVtoXML::zz.

                                                                                       {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: \t\tInside Layers";
#endif
  double zi(zMinBlock_);
  int laymin(0);
  const double tol(0.01);
  for (unsigned int i = 0; i < layers_.size(); i++) {
    double zo = zi + layerThick_[i];
    double routF = rMax(zi);
    int laymax = laymin + layers_[i];
    double zz = zi;
    double thickTot(0);
    for (int ly = laymin; ly < laymax; ++ly) {
      int ii = layerType_[ly];
      int copy = copyNumber_[ii];
      double rinB = (layerSense_[ly] == 0) ? (zo * slopeB_[0]) : (zo * slopeB_[1]);
      zz += (0.5 * thick_[ii]);
      thickTot += thick_[ii];

      std::string name = "HGCal" + names_[ii] + std::to_string(copy);
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: Layer " << ly << ":" << ii << " Front " << zi << ", "
                                    << routF << " Back " << zo << ", " << rinB << " superlayer thickness "
                                    << layerThick_[i];
#endif
      DDName matName(DDSplit(materials_[ii]).first, DDSplit(materials_[ii]).second);
      DDMaterial matter(matName);
      DDLogicalPart glog;
      if (layerSense_[ly] == 0) {
        double alpha = 1._pi / sectors_;
        double rmax = routF * cos(alpha) - tol;
        std::vector<double> pgonZ, pgonRin, pgonRout;
        pgonZ.emplace_back(-0.5 * thick_[ii]);
        pgonZ.emplace_back(0.5 * thick_[ii]);
        pgonRin.emplace_back(rinB);
        pgonRin.emplace_back(rinB);
        pgonRout.emplace_back(rmax);
        pgonRout.emplace_back(rmax);
        DDSolid solid =
            DDSolidFactory::polyhedra(DDName(name, idNameSpace_), sectors_, -alpha, 2._pi, pgonZ, pgonRin, pgonRout);
        glog = DDLogicalPart(solid.ddname(), matter, solid);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << solid.name() << " polyhedra of " << sectors_
                                      << " sectors covering " << convertRadToDeg(-alpha) << ":"
                                      << (360.0 + convertRadToDeg(-alpha)) << " with " << pgonZ.size() << " sections";
        for (unsigned int k = 0; k < pgonZ.size(); ++k)
          edm::LogVerbatim("HGCalGeom") << "[" << k << "] z " << pgonZ[k] << " R " << pgonRin[k] << ":" << pgonRout[k];
#endif
      } else {
        DDSolid solid = DDSolidFactory::tubs(DDName(name, idNameSpace_), 0.5 * thick_[ii], rinB, routF, 0.0, 2._pi);
        glog = DDLogicalPart(solid.ddname(), matter, solid);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << solid.name() << " Tubs made of " << matName
                                      << " of dimensions " << rinB << ", " << routF << ", " << 0.5 * thick_[ii]
                                      << ", 0.0, 360.0";
#endif
        positionSensitive(glog, rinB, routF, cpv);
      }
      DDTranslation r1(0, 0, zz);
      DDRotation rot;
      cpv.position(glog, module, copy, r1, rot);
      ++copyNumber_[ii];
#ifdef EDM_ML_DEBUG
      edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << glog.name() << " number " << copy << " positioned in "
                                    << module.name() << " at " << r1 << " with " << rot;
#endif
      zz += (0.5 * thick_[ii]);
    }  // End of loop over layers in a block
    zi = zo;
    laymin = laymax;
    if (fabs(thickTot - layerThick_[i]) > tol_) {
      if (thickTot > layerThick_[i]) {
        edm::LogError("HGCalGeom") << "Thickness of the partition " << layerThick_[i] << " is smaller than thickness "
                                   << thickTot << " of all its components **** ERROR ****\n";
      } else {
        edm::LogWarning("HGCalGeom") << "Thickness of the partition " << layerThick_[i] << " does not match with "
                                     << thickTot << " of the components\n";
      }
    }
  }  // End of loop over blocks
}

execute()

void DDHGCalModuleAlgo::execute ( DDCompactView &  cpv )

override

Definition at line 148 of file DDHGCalModuleAlgo.cc.

References DDHGCalGeom::constructLayers(), and class-composition::parent.

                                                  {
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalModuleAlgo...";
#endif
  copies_.clear();
  constructLayers(parent(), cpv);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << copies_.size() << " different wafer copy numbers";
#endif
  copies_.clear();
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalModuleAlgo construction";
#endif
}

initialize()

void DDHGCalModuleAlgo::initialize ( const DDNumericArguments &  nArgs, const DDVectorArguments &  vArgs, const DDMapArguments &  mArgs, const DDStringArguments &  sArgs, const DDStringVectorArguments &  vsArgs  )

override

Definition at line 80 of file DDHGCalModuleAlgo.cc.

References dbl_to_int(), mps_fire::i, createfilelist::int, and DDCurrentNamespace::ns().

                                                                          {
  wafer_ = vsArgs["WaferName"];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << wafer_.size() << " wafers";
  for (unsigned int i = 0; i < wafer_.size(); ++i)
    edm::LogVerbatim("HGCalGeom") << "Wafer[" << i << "] " << wafer_[i];
#endif
  materials_ = vsArgs["MaterialNames"];
  names_ = vsArgs["VolumeNames"];
  thick_ = vArgs["Thickness"];
  copyNumber_.resize(materials_.size(), 1);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << materials_.size() << " types of volumes";
  for (unsigned int i = 0; i < names_.size(); ++i)
    edm::LogVerbatim("HGCalGeom") << "Volume [" << i << "] " << names_[i] << " of thickness " << thick_[i]
                                  << " filled with " << materials_[i] << " first copy number " << copyNumber_[i];
#endif
  layers_ = dbl_to_int(vArgs["Layers"]);
  layerThick_ = vArgs["LayerThick"];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << layers_.size() << " blocks";
  for (unsigned int i = 0; i < layers_.size(); ++i)
    edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] of thickness " << layerThick_[i] << " with " << layers_[i]
                                  << " layers";
#endif
  layerType_ = dbl_to_int(vArgs["LayerType"]);
  layerSense_ = dbl_to_int(vArgs["LayerSense"]);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << layerType_.size() << " layers";
  for (unsigned int i = 0; i < layerType_.size(); ++i)
    edm::LogVerbatim("HGCalGeom") << "Layer [" << i << "] with material type " << layerType_[i] << " sensitive class "
                                  << layerSense_[i];
#endif
  zMinBlock_ = nArgs["zMinBlock"];
  rMaxFine_ = nArgs["rMaxFine"];
  waferW_ = nArgs["waferW"];
  waferGap_ = nArgs["waferGap"];
  sectors_ = (int)(nArgs["Sectors"]);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: zStart " << zMinBlock_ << " rFineCoarse " << rMaxFine_
                                << " wafer width " << waferW_ << " gap among wafers " << waferGap_ << " sectors "
                                << sectors_;
#endif
  slopeB_ = vArgs["SlopeBottom"];
  slopeT_ = vArgs["SlopeTop"];
  zFront_ = vArgs["ZFront"];
  rMaxFront_ = vArgs["RMaxFront"];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: Bottom slopes " << slopeB_[0] << ":" << slopeB_[1] << " and "
                                << slopeT_.size() << " slopes for top";
  for (unsigned int i = 0; i < slopeT_.size(); ++i)
    edm::LogVerbatim("HGCalGeom") << "Block [" << i << "] Zmin " << zFront_[i] << " Rmax " << rMaxFront_[i] << " Slope "
                                  << slopeT_[i];
#endif
  idNameSpace_ = DDCurrentNamespace::ns();
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: NameSpace " << idNameSpace_;
#endif
}

positionSensitive()

void DDHGCalModuleAlgo::positionSensitive ( DDLogicalPart &  glog, double  rin, double  rout, DDCompactView &  cpv  )

protected

Definition at line 265 of file DDHGCalModuleAlgo.cc.

References filterCSVwithJSON::copy, TCMET_cfi::corner, DDBase< N, C >::ddname(), DDSplit(), PVValHelper::dx, PVValHelper::dy, createfilelist::int, HGCalParameters::k_CornerSize, mergeVDriftHistosByStation::name, EgHLTOffHistBins_cfi::nr, HGCalTypes::packTypeUV(), DDCompactView::position(), idealTransformation::rotation, findQualityFiles::rr, mathSSE::sqrt(), funct::tan(), HGCalGeomTools::waferCorner(), and ALPAKA_ACCELERATOR_NAMESPACE::vertexFinder::ww.

                                                                                                          {
  double ww = (waferW_ + waferGap_);
  double dx = 0.5 * ww;
  double dy = 3.0 * dx * tan(30._deg);
  double rr = 2.0 * dx * tan(30._deg);
  int ncol = (int)(2.0 * rout / ww) + 1;
  int nrow = (int)(rout / (ww * tan(30._deg))) + 1;
  int incm(0), inrm(0);
#ifdef EDM_ML_DEBUG
  int kount(0);
  edm::LogVerbatim("HGCalGeom") << glog.ddname() << " rout " << rout << " Row " << nrow << " Column " << ncol;
#endif
  for (int nr = -nrow; nr <= nrow; ++nr) {
    int inr = (nr >= 0) ? nr : -nr;
    for (int nc = -ncol; nc <= ncol; ++nc) {
      int inc = (nc >= 0) ? nc : -nc;
      if (inr % 2 == inc % 2) {
        double xpos = nc * dx;
        double ypos = nr * dy;
        auto const& corner = HGCalGeomTools::waferCorner(xpos, ypos, dx, rr, rin, rout, true);
        if (corner.first == (int)(HGCalParameters::k_CornerSize)) {
          double rpos = std::sqrt(xpos * xpos + ypos * ypos);
          DDTranslation tran(xpos, ypos, 0.0);
          DDRotation rotation;
          int copy = HGCalTypes::packTypeUV(0, nc, nr);
          DDName name = (rpos < rMaxFine_) ? DDName(DDSplit(wafer_[0]).first, DDSplit(wafer_[0]).second)
                                           : DDName(DDSplit(wafer_[1]).first, DDSplit(wafer_[1]).second);
          cpv.position(name, glog.ddname(), copy, tran, rotation);
          if (inc > incm)
            incm = inc;
          if (inr > inrm)
            inrm = inr;
#ifdef EDM_ML_DEBUG
          kount++;
#endif
          if (copies_.count(copy) == 0)
            copies_.insert(copy);
#ifdef EDM_ML_DEBUG
          edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: " << name << " number " << copy << " positioned in "
                                        << glog.ddname() << " at " << tran << " with " << rotation;
#endif
        }
      }
    }
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalModuleAlgo: # of columns " << incm << " # of rows " << inrm << " and "
                                << kount << " wafers for " << glog.ddname();
#endif
}

rMax()

double DDHGCalModuleAlgo::rMax ( double  z )

protected

Definition at line 246 of file DDHGCalModuleAlgo.cc.

References isotrackApplyRegressor::k.

                                       {
  double r(0);
#ifdef EDM_ML_DEBUG
  unsigned int ik(0);
#endif
  for (unsigned int k = 0; k < slopeT_.size(); ++k) {
    if (z < zFront_[k])
      break;
    r = rMaxFront_[k] + (z - zFront_[k]) * slopeT_[k];
#ifdef EDM_ML_DEBUG
    ik = k;
#endif
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "rMax : " << z << ":" << ik << ":" << r;
#endif
  return r;
}

Member Data Documentation

copies_

std::unordered_set<int> DDHGCalModuleAlgo::copies_

private

Definition at line 71 of file DDHGCalModuleAlgo.cc.

copyNumber_

std::vector<int> DDHGCalModuleAlgo::copyNumber_

private

Definition at line 56 of file DDHGCalModuleAlgo.cc.

idNameSpace_

std::string DDHGCalModuleAlgo::idNameSpace_

private

Definition at line 70 of file DDHGCalModuleAlgo.cc.

layers_

std::vector<int> DDHGCalModuleAlgo::layers_

private

Definition at line 57 of file DDHGCalModuleAlgo.cc.

layerSense_

std::vector<int> DDHGCalModuleAlgo::layerSense_

private

Definition at line 60 of file DDHGCalModuleAlgo.cc.

layerThick_

std::vector<double> DDHGCalModuleAlgo::layerThick_

private

Definition at line 58 of file DDHGCalModuleAlgo.cc.

layerType_

std::vector<int> DDHGCalModuleAlgo::layerType_

private

Definition at line 59 of file DDHGCalModuleAlgo.cc.

materials_

std::vector<std::string> DDHGCalModuleAlgo::materials_

private

Definition at line 53 of file DDHGCalModuleAlgo.cc.

names_

std::vector<std::string> DDHGCalModuleAlgo::names_

private

Definition at line 54 of file DDHGCalModuleAlgo.cc.

rMaxFine_

double DDHGCalModuleAlgo::rMaxFine_

private

Definition at line 62 of file DDHGCalModuleAlgo.cc.

rMaxFront_

std::vector<double> DDHGCalModuleAlgo::rMaxFront_

private

Definition at line 69 of file DDHGCalModuleAlgo.cc.

sectors_

int DDHGCalModuleAlgo::sectors_

private

Definition at line 65 of file DDHGCalModuleAlgo.cc.

slopeB_

std::vector<double> DDHGCalModuleAlgo::slopeB_

private

Definition at line 66 of file DDHGCalModuleAlgo.cc.

slopeT_

std::vector<double> DDHGCalModuleAlgo::slopeT_

private

Definition at line 67 of file DDHGCalModuleAlgo.cc.

thick_

std::vector<double> DDHGCalModuleAlgo::thick_

private

Definition at line 55 of file DDHGCalModuleAlgo.cc.

tol_

constexpr double DDHGCalModuleAlgo::tol_ = 0.00001

staticprivate

Definition at line 50 of file DDHGCalModuleAlgo.cc.

wafer_

std::vector<std::string> DDHGCalModuleAlgo::wafer_

private

Definition at line 52 of file DDHGCalModuleAlgo.cc.

waferGap_

double DDHGCalModuleAlgo::waferGap_

private

Definition at line 64 of file DDHGCalModuleAlgo.cc.

waferW_

double DDHGCalModuleAlgo::waferW_

private

Definition at line 63 of file DDHGCalModuleAlgo.cc.

zFront_

std::vector<double> DDHGCalModuleAlgo::zFront_

private

Definition at line 68 of file DDHGCalModuleAlgo.cc.

zMinBlock_

double DDHGCalModuleAlgo::zMinBlock_

private

Definition at line 61 of file DDHGCalModuleAlgo.cc.