Inheritance diagram for DDHGCalTBModule:

Public Member Functions
	DDHGCalTBModule ()

void	execute (DDCompactView &cpv) override

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

	~DDHGCalTBModule () override

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

void	positionSensitive (DDLogicalPart &glog, int type, double rin, double rout, DDCompactView &cpv)

double	rMax (double z)

Private Attributes
double	absorbH_

double	absorbW_

std::unordered_set< int >	copies_

std::vector< int >	copyNumber_

std::vector< std::string >	covers_

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	waferGap_

double	waferW_

std::vector< double >	zFront_

double	zMinBlock_

Detailed Description

Definition at line 24 of file DDHGCalTBModule.cc.

Constructor & Destructor Documentation

DDHGCalTBModule::DDHGCalTBModule ( )

Definition at line 69 of file DDHGCalTBModule.cc.

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

DDHGCalTBModule::~DDHGCalTBModule ( )

override

Definition at line 75 of file DDHGCalTBModule.cc.

75 {}

Member Function Documentation

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

protected

Definition at line 174 of file DDHGCalTBModule.cc.

References DDSolidFactory::box(), filterCSVwithJSON::copy, DDBase< N, C >::ddname(), DDSplit(), dqmdumpme::first, mps_fire::i, cuy::ii, Skims_PA_cff::name, DDBase< N, C >::name(), angle_units::piRadians(), DDCompactView::position(), diffTwoXMLs::r1, photonAnalyzer_cfi::rMax, makeMuonMisalignmentScenario::rot, edm::second(), AlCaHLTBitMon_QueryRunRegistry::string, DDSolidFactory::tubs(), and geometryCSVtoXML::zz.

                                                                                      {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule 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];
     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") << "DDHGCalTBModule: 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) {
         DDSolid solid = DDSolidFactory::box(DDName(name, idNameSpace_), absorbW_, absorbH_, 0.5 * thick_[ii]);
         glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule test: " << solid.name() << " box of dimension " << absorbW_
                                       << ":" << absorbH_ << ":" << 0.5 * thick_[ii];
 #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") << "DDHGCalTBModule: " << solid.name() << " Tubs made of " << matName
                                       << " of dimensions " << rinB << ", " << routF << ", " << 0.5 * thick_[ii]
                                       << ", 0.0, 360.0";
 #endif
         positionSensitive(glog, layerSense_[ly], 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") << "DDHGCalTBModule test: " << 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 DDHGCalTBModule::execute ( DDCompactView & cpv )

override

Definition at line 159 of file DDHGCalTBModule.cc.

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

                                                 {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "==>> Constructing DDHGCalTBModule...";
 #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 DDHGCalTBModule construction ...";
 #endif
 }

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

override

Definition at line 77 of file DDHGCalTBModule.cc.

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

                                                                         {
   wafer_ = vsArgs["WaferName"];
   covers_ = vsArgs["CoverName"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: " << wafer_.size() << " wafers";
   unsigned int i(0);
   for (auto wafer : wafer_) {
     edm::LogVerbatim("HGCalGeom") << "Wafer[" << i << "] " << wafer;
     ++i;
   }
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: " << covers_.size() << " covers";
   i = 0;
   for (auto cover : covers_) {
     edm::LogVerbatim("HGCalGeom") << "Cover[" << i << "] " << cover;
     ++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") << "DDHGCalTBModule: " << 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") << "DDHGCalTBModule: " << 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") << "DDHGCalTBModule: " << 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"];
   absorbW_ = nArgs["absorberW"];
   absorbH_ = nArgs["absorberH"];
   sectors_ = (int)(nArgs["Sectors"]);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: zStart " << zMinBlock_ << " rFineCoarse " << rMaxFine_
                                 << " wafer width " << waferW_ << " gap among wafers " << waferGap_ << " absorber width "
                                 << absorbW_ << " absorber height " << absorbH_ << " sectors " << sectors_;
 #endif
   slopeB_ = vArgs["SlopeBottom"];
   slopeT_ = vArgs["SlopeTop"];
   zFront_ = vArgs["ZFront"];
   rMaxFront_ = vArgs["RMaxFront"];
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: 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") << "DDHGCalTBModule: NameSpace " << idNameSpace_;
 #endif
 }

void DDHGCalTBModule::positionSensitive	(	DDLogicalPart &	glog,
		int	type,
		double	rin,
		double	rout,
		DDCompactView &	cpv
	)

protected

Definition at line 261 of file DDHGCalTBModule.cc.

References filterCSVwithJSON::copy, DDBase< N, C >::ddname(), DDSplit(), DEFINE_EDM_PLUGIN, PVValHelper::dx, PVValHelper::dy, dqmdumpme::first, createfilelist::int, dqmdumpme::k, Skims_PA_cff::name, EgHLTOffHistBins_cfi::nr, DDCompactView::position(), idealTransformation::rotation, findQualityFiles::rr, edm::second(), mathSSE::sqrt(), and funct::tan().

                                                                                                                   {
   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), kount(0);
   double xc[6], yc[6];
 #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;
         xc[0] = xpos + dx;
         yc[0] = ypos - 0.5 * rr;
         xc[1] = xpos + dx;
         yc[1] = ypos + 0.5 * rr;
         xc[2] = xpos;
         yc[2] = ypos + rr;
         xc[3] = xpos - dx;
         yc[3] = ypos + 0.5 * rr;
         xc[4] = xpos + dx;
         yc[4] = ypos - 0.5 * rr;
         xc[5] = xpos;
         yc[5] = ypos - rr;
         bool cornerAll(true);
         for (int k = 0; k < 6; ++k) {
           double rpos = std::sqrt(xc[k] * xc[k] + yc[k] * yc[k]);
           if (rpos < rin || rpos > rout)
             cornerAll = false;
         }
         if (cornerAll) {
           double rpos = std::sqrt(xpos * xpos + ypos * ypos);
           DDTranslation tran(xpos, ypos, 0.0);
           DDRotation rotation;
           int copy = inr * 100 + inc;
           if (nc < 0)
             copy += 10000;
           if (nr < 0)
             copy += 100000;
           DDName name;
           if (type == 1) {
             name = (rpos < rMaxFine_) ? DDName(DDSplit(wafer_[0]).first, DDSplit(wafer_[0]).second)
                                       : DDName(DDSplit(wafer_[1]).first, DDSplit(wafer_[1]).second);
           } else {
             name = DDName(DDSplit(covers_[type - 2]).first, DDSplit(covers_[type - 2]).second);
           }
           cpv.position(name, glog.ddname(), copy, tran, rotation);
           if (inc > incm)
             incm = inc;
           if (inr > inrm)
             inrm = inr;
           kount++;
           if (copies_.count(copy) == 0 && type == 1)
             copies_.insert(copy);
 #ifdef EDM_ML_DEBUG
           edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: " << name << " number " << copy << " positioned in "
                                         << glog.ddname() << " at " << tran << " with " << rotation;
 #endif
         }
       }
     }
   }
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "DDHGCalTBModule: # of columns " << incm << " # of rows " << inrm << " and " << kount
                                 << " wafers for " << glog.ddname();
 #endif
 }

double DDHGCalTBModule::rMax ( double z )

protected

Definition at line 242 of file DDHGCalTBModule.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

double DDHGCalTBModule::absorbH_

private

Definition at line 58 of file DDHGCalTBModule.cc.

double DDHGCalTBModule::absorbW_

private

Definition at line 57 of file DDHGCalTBModule.cc.

std::unordered_set<int> DDHGCalTBModule::copies_

private

Definition at line 66 of file DDHGCalTBModule.cc.

std::vector<int> DDHGCalTBModule::copyNumber_

private

Definition at line 48 of file DDHGCalTBModule.cc.

std::vector<std::string> DDHGCalTBModule::covers_

private

Definition at line 44 of file DDHGCalTBModule.cc.

std::string DDHGCalTBModule::idName_

private

Definition at line 64 of file DDHGCalTBModule.cc.

std::string DDHGCalTBModule::idNameSpace_

private

Definition at line 65 of file DDHGCalTBModule.cc.

std::vector<int> DDHGCalTBModule::layers_

private

Definition at line 49 of file DDHGCalTBModule.cc.

std::vector<int> DDHGCalTBModule::layerSense_

private

Definition at line 52 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::layerThick_

private

Definition at line 50 of file DDHGCalTBModule.cc.

std::vector<int> DDHGCalTBModule::layerType_

private

Definition at line 51 of file DDHGCalTBModule.cc.

std::vector<std::string> DDHGCalTBModule::materials_

private

Definition at line 45 of file DDHGCalTBModule.cc.

std::vector<std::string> DDHGCalTBModule::names_

private

Definition at line 46 of file DDHGCalTBModule.cc.

double DDHGCalTBModule::rMaxFine_

private

Definition at line 54 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::rMaxFront_

private

Definition at line 63 of file DDHGCalTBModule.cc.

int DDHGCalTBModule::sectors_

private

Definition at line 59 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::slopeB_

private

Definition at line 60 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::slopeT_

private

Definition at line 61 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::thick_

private

Definition at line 47 of file DDHGCalTBModule.cc.

std::vector<std::string> DDHGCalTBModule::wafer_

private

Definition at line 43 of file DDHGCalTBModule.cc.

double DDHGCalTBModule::waferGap_

private

Definition at line 56 of file DDHGCalTBModule.cc.

double DDHGCalTBModule::waferW_

private

Definition at line 55 of file DDHGCalTBModule.cc.

std::vector<double> DDHGCalTBModule::zFront_

private

Definition at line 62 of file DDHGCalTBModule.cc.

double DDHGCalTBModule::zMinBlock_

private

Definition at line 53 of file DDHGCalTBModule.cc.

Public Member Functions

Protected Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation