Inheritance diagram for DDHGCalModule:

Public Member Functions
	DDHGCalModule ()

void	execute (DDCompactView &cpv) override

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

	~DDHGCalModule () 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	idName

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	waferW

std::vector< double >	zFront

double	zMinBlock

Detailed Description

Definition at line 31 of file DDHGCalModule.cc.

Constructor & Destructor Documentation

DDHGCalModule::DDHGCalModule ( )

Definition at line 72 of file DDHGCalModule.cc.

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

DDHGCalModule::~DDHGCalModule ( )

override

Definition at line 78 of file DDHGCalModule.cc.

78 {}

Member Function Documentation

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

protected

Definition at line 166 of file DDHGCalModule.cc.

References zMuMuMuonUserData::alpha, angle_units::operators::convertRadToDeg(), filterCSVwithJSON::copy, funct::cos(), DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, mps_fire::i, cuy::ii, dqmdumpme::k, hgcalTopologyTester_cfi::layers, Skims_PA_cff::name, DDBase< N, C >::name(), names, angle_units::piRadians(), DDSolidFactory::polyhedra(), DDCompactView::position(), diffTwoXMLs::r1, photonAnalyzer_cfi::rMax, makeMuonMisalignmentScenario::rot, edm::second(), volumeBasedMagneticField_1103l_cfi::sectors, AlCaHLTBitMon_QueryRunRegistry::string, DDSolidFactory::tubs(), and geometryCSVtoXML::zz.

                                                                                    {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: \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") << "DDHGCalModule: 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 = geant_units::piRadians / 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 * geant_units::piRadians, pgonZ, pgonRin, pgonRout);
         glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: " << 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 * geant_units::piRadians);
         glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: " << solid.name() << " Tubs made of " << matName
                                       << " of dimensions " << rinB << ", " << routF << ", " << 0.5 * thick[ii]
                                       << ", 0.0, 360.0";
         edm::LogVerbatim("HGCalGeom") << "DDHGCalModule test position in: " << glog.name() << " number " << copy;
 #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") << "DDHGCalModule: " << 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]) < 0.00001) {
     } else 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 if (thickTot < layerThick[i]) {
       edm::LogWarning("HGCalGeom") << "Thickness of the partition " << layerThick[i] << " does not match with "
                                    << thickTot << " of the components\n";
     }
   }  // End of loop over blocks
 }

void DDHGCalModule::execute ( DDCompactView & cpv )

override

Definition at line 148 of file DDHGCalModule.cc.

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

                                               {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalModule...";
 #endif
   copies.clear();
   constructLayers(parent(), cpv);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << copies.size() << " different wafer copy numbers";
   int k(0);
   for (std::unordered_set<int>::const_iterator itr = copies.begin(); itr != copies.end(); ++itr, ++k)
     edm::LogVerbatim("HGCalGeom") << "Copy[" << k << "] : " << (*itr);
 #endif
   copies.clear();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "<<== End of DDHGCalModule construction ...";
 #endif
 }

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

override

Definition at line 80 of file DDHGCalModule.cc.

References dbl_to_int(), mps_fire::i, createfilelist::int, hgcalTopologyTester_cfi::layers, names, DDCurrentNamespace::ns(), and volumeBasedMagneticField_1103l_cfi::sectors.

                                                                       {
   wafer = vsArgs["WaferName"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: " << 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"];
   for (unsigned int i = 0; i < materials.size(); ++i) {
     copyNumber.emplace_back(1);
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: " << 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") << "DDHGCalModule: " << 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") << "DDHGCalModule: " << 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"];
   sectors = (int)(nArgs["Sectors"]);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: zStart " << zMinBlock << " rFineCoarse " << rMaxFine
                                 << " wafer width " << waferW << " sectors " << sectors;
 #endif
   slopeB = vArgs["SlopeBottom"];
   slopeT = vArgs["SlopeTop"];
   zFront = vArgs["ZFront"];
   rMaxFront = vArgs["RMaxFront"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: 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") << "DDHGCalModule: NameSpace " << idNameSpace;
 #endif
 }

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

protected

Definition at line 268 of file DDHGCalModule.cc.

References filterCSVwithJSON::copy, distTCMET_cfi::corner, DDBase< N, C >::ddname(), DDSplit(), DEFINE_EDM_PLUGIN, PVValHelper::dx, PVValHelper::dy, createfilelist::int, HGCalParameters::k_CornerSize, Skims_PA_cff::name, nin, EgHLTOffHistBins_cfi::nr, DDCompactView::position(), idealTransformation::rotation, findQualityFiles::rr, mathSSE::sqrt(), funct::tan(), and HGCalGeomTools::waferCorner().

                                                                                                       {
   double dx = 0.5 * waferW;
   double dy = 3.0 * dx * tan(30._deg);
   double rr = 2.0 * dx * tan(30._deg);
   int ncol = (int)(2.0 * rout / waferW) + 1;
   int nrow = (int)(rout / (waferW * tan(30._deg))) + 1;
   int incm(0), inrm(0), kount(0), ntot(0), nin(0), nfine(0), ncoarse(0);
 #ifdef EDM_ML_DEBUG
   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;
         std::pair<int, int> corner = HGCalGeomTools::waferCorner(xpos, ypos, dx, rr, rin, rout, true);
         ++ntot;
         if (corner.first > 0) {
           int copy = inr * 100 + inc;
           if (nc < 0)
             copy += 10000;
           if (nr < 0)
             copy += 100000;
           if (inc > incm)
             incm = inc;
           if (inr > inrm)
             inrm = inr;
           kount++;
           if (copies.count(copy) == 0)
             copies.insert(copy);
           if (corner.first == (int)(HGCalParameters::k_CornerSize)) {
             double rpos = std::sqrt(xpos * xpos + ypos * ypos);
             DDTranslation tran(xpos, ypos, 0.0);
             DDRotation rotation;
             ++nin;
             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 (rpos < rMaxFine)
               ++nfine;
             else
               ++ncoarse;
 #ifdef EDM_ML_DEBUG
             edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: " << name << " number " << copy << " positioned in "
                                           << glog.ddname() << " at " << tran << " with " << rotation;
 #endif
           }
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalModule: # of columns " << incm << " # of rows " << inrm << " and " << nin
                                 << ":" << kount << ":" << ntot << " wafers (" << nfine << ":" << ncoarse << ") for "
                                 << glog.ddname() << " R " << rin << ":" << rout;
 #endif
 }

double DDHGCalModule::rMax ( double z )

protected

Definition at line 249 of file DDHGCalModule.cc.

References dqmdumpme::k, and alignCSCRings::r.

                                    {
   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

std::unordered_set<int> DDHGCalModule::copies

private

Definition at line 69 of file DDHGCalModule.cc.

std::vector<int> DDHGCalModule::copyNumber

private

Definition at line 54 of file DDHGCalModule.cc.

std::string DDHGCalModule::idName

private

Definition at line 67 of file DDHGCalModule.cc.

std::string DDHGCalModule::idNameSpace

private

Definition at line 68 of file DDHGCalModule.cc.

std::vector<int> DDHGCalModule::layers

private

Definition at line 55 of file DDHGCalModule.cc.

std::vector<int> DDHGCalModule::layerSense

private

Definition at line 58 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::layerThick

private

Definition at line 56 of file DDHGCalModule.cc.

std::vector<int> DDHGCalModule::layerType

private

Definition at line 57 of file DDHGCalModule.cc.

std::vector<std::string> DDHGCalModule::materials

private

Definition at line 51 of file DDHGCalModule.cc.

std::vector<std::string> DDHGCalModule::names

private

Definition at line 52 of file DDHGCalModule.cc.

double DDHGCalModule::rMaxFine

private

Definition at line 60 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::rMaxFront

private

Definition at line 66 of file DDHGCalModule.cc.

int DDHGCalModule::sectors

private

Definition at line 62 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::slopeB

private

Definition at line 63 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::slopeT

private

Definition at line 64 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::thick

private

Definition at line 53 of file DDHGCalModule.cc.

std::vector<std::string> DDHGCalModule::wafer

private

Definition at line 50 of file DDHGCalModule.cc.

double DDHGCalModule::waferW

private

Definition at line 61 of file DDHGCalModule.cc.

std::vector<double> DDHGCalModule::zFront

private

Definition at line 65 of file DDHGCalModule.cc.

double DDHGCalModule::zMinBlock

private

Definition at line 59 of file DDHGCalModule.cc.

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation