#include <DDAHcalModuleAlgo.h>

Inheritance diagram for DDAHcalModuleAlgo:

Public Member Functions
	DDAHcalModuleAlgo ()

void	execute (DDCompactView &cpv) override

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

	~DDAHcalModuleAlgo () override

Protected Member Functions
void	constructLayers (const DDLogicalPart &, DDCompactView &cpv)

void	positionSensitive (DDLogicalPart &glog, DDCompactView &cpv)

Private Attributes
std::unordered_set< int >	copies

std::vector< int >	copyNumber

std::vector< double >	heights

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

std::vector< double >	thick

std::string	tile

std::vector< int >	tileN

std::vector< double >	tileStep

std::vector< double >	widths

double	zMinBlock

Detailed Description

Definition at line 11 of file DDAHcalModuleAlgo.h.

Constructor & Destructor Documentation

DDAHcalModuleAlgo::DDAHcalModuleAlgo ( )

Definition at line 20 of file DDAHcalModuleAlgo.cc.

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

DDAHcalModuleAlgo::~DDAHcalModuleAlgo ( )

override

Definition at line 26 of file DDAHcalModuleAlgo.cc.

26 {}

Member Function Documentation

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

protected

Definition at line 113 of file DDAHcalModuleAlgo.cc.

References DDSolidFactory::box(), popcon2dropbox::copy(), copyNumber, DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, heights, mps_fire::i, idNameSpace, cuy::ii, layers, layerSense, layerThick, layerType, materials, dataset::name, DDBase< N, C >::name(), names, DDCompactView::position(), positionSensitive(), diffTwoXMLs::r1, makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, thick, widths, zMinBlock, and geometryCSVtoXML::zz.

Referenced by execute().

                                                             {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDAHcalModuleAlgo test: \t\tInside Layers";
 #endif
   double zi(zMinBlock);
   int laymin(0);
   for (unsigned int i = 0; i < layers.size(); i++) {
     double zo = zi + layerThick[i];
     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];
       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")
     << "DDAHcalModuleAlgo test: Layer " << ly << ":" << ii << " Front " 
     << zi << " Back " << zo << " superlayer thickness " << layerThick[i];
 #endif
       DDName matName(DDSplit(materials[ii]).first,
                      DDSplit(materials[ii]).second);
       DDMaterial matter(matName);
       DDLogicalPart glog;
       if (layerSense[ly] == 0) {
         DDSolid solid =
             DDSolidFactory::box(DDName(name, idNameSpace), 0.5 * widths[0],
                                 0.5 * heights[0], 0.5 * thick[ii]);
         glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom")
       << "DDAHcalModuleAlgo: " << solid.name() << " Box made of " 
       << matName << " of dimensions " << 0.5 * widths[0] << ", " 
       << 0.5 * heights[0] << ", " << 0.5 * thick[ii];
 #endif
       } else {
         DDSolid solid =
             DDSolidFactory::box(DDName(name, idNameSpace), 0.5 * widths[1],
                                 0.5 * heights[1], 0.5 * thick[ii]);
         glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom")
       << "DDAHcalModuleAlgo: " << solid.name() << " Box made of " 
       << matName << " of dimensions " << 0.5 * widths[1] << ", "
       << 0.5 * heights[1] << ", " << 0.5 * thick[ii];
 #endif
         positionSensitive(glog, cpv);
       }
       DDTranslation r1(0, 0, zz);
       DDRotation rot;
       cpv.position(glog, module, copy, r1, rot);
       ++copyNumber[ii];
 #ifdef EDM_ML_DEBUG
       edm::LogVerbatim("HGCalGeom")
     << "DDAHcalModuleAlgo: " << 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 DDAHcalModuleAlgo::execute ( DDCompactView & cpv )

override

Definition at line 103 of file DDAHcalModuleAlgo.cc.

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

                                                   {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDAHcalModuleAlgo...";
 #endif
   constructLayers(parent(), cpv);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") <<"<<== End of DDAHcalModuleAlgo construction";
 #endif
 }

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

override

Definition at line 28 of file DDAHcalModuleAlgo.cc.

References copyNumber, dbl_to_int(), heights, mps_fire::i, idNameSpace, layers, layerSense, layerThick, layerType, materials, names, DDCurrentNamespace::ns(), thick, tile, tileN, tileStep, widths, and zMinBlock.

                                                                           {
   tile = sArgs["TileName"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDAHcalModuleAlgo: Tile " << tile;
 #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") 
     << "DDAHcalModuleAlgo: " << 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") 
     << "DDAHcalModuleAlgo: " << 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") 
     << "DDAHcalModuleAlgo: " << 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
   widths = vArgs["Widths"];
   heights = vArgs["Heights"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") 
     << "DDAHcalModuleAlgo: " << widths.size() << " sizes for width "
     << "and height:";
   for (unsigned int i = 0; i < widths.size(); ++i)
     edm::LogVerbatim("HGCalGeom") 
       << " [" << i << "] " << widths[i] << ":" << heights[i];
 #endif
   tileN    = dbl_to_int(vArgs["TileN"]);
   tileStep = vArgs["TileStep"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") 
     << "DDAHcalModuleAlgo: " << tileN.size() << " tile positioning parameters";
   for (unsigned int i = 0; i < tileN.size(); ++i)
     edm::LogVerbatim("HGCalGeom") 
       << " [" << i << "] " << tileN[i] << ":" << tileStep[i];
 #endif
   zMinBlock = nArgs["zMinBlock"];
   idNameSpace = DDCurrentNamespace::ns();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom")
     << "DDAHcalModuleAlgo: NameSpace " << idNameSpace;
 #endif
 }

void DDAHcalModuleAlgo::positionSensitive	(	DDLogicalPart &	glog,
		DDCompactView &	cpv
	)

protected

Definition at line 193 of file DDAHcalModuleAlgo.cc.

References popcon2dropbox::copy(), DDBase< N, C >::ddname(), DDSplit(), plotBeamSpotDB::first, dataset::name, create_public_lumi_plots::ncol, DDCompactView::position(), idealTransformation::rotation, edm::second(), tile, tileN, and tileStep.

Referenced by constructLayers().

                                                               {
   int ncol = tileN[0] / 2;
   int nrow = tileN[1] / 2;
 #ifdef EDM_ML_DEBUG
   int kount(0);
   edm::LogVerbatim("HGCalGeom") 
     << glog.ddname() << " Row " << nrow << " Column " << ncol;
 #endif
   for (int nr = -nrow; nr <= nrow; ++nr) {
     int inr = (nr >= 0) ? nr : -nr;
     double ypos =
         (nr >= 0) ? (inr - 0.5) * tileStep[1] : -(inr - 0.5) * tileStep[1];
     for (int nc = -ncol; nc <= ncol; ++nc) {
       int inc = (nc >= 0) ? nc : -nc;
       double xpos =
           (nc >= 0) ? (inc - 0.5) * tileStep[0] : -(inc - 0.5) * tileStep[0];
       if (nr != 0 && nc != 0) {
         DDTranslation tran(xpos, ypos, 0.0);
         DDRotation rotation;
         int copy = inr * 10 + inc;
         if (nc < 0) copy += 100;
         if (nr < 0) copy += 1000;
         DDName name = DDName(DDSplit(tile).first, DDSplit(tile).second);
         cpv.position(name, glog.ddname(), copy, tran, rotation);
 #ifdef EDM_ML_DEBUG
         kount++;
         edm::LogVerbatim("HGCalGeom") 
       << "DDAHcalModuleAlgo: " << name << " number " << copy
       << " positioned in " << glog.ddname() << " at " << tran
       << " with " << rotation;
 #endif
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") 
     << "DDAHcalModuleAlgo: " << kount << " tiles for " << glog.ddname();
 #endif
 }