#include <DDHGCalWafer8.h>

Inheritance diagram for DDHGCalWafer8:

Public Member Functions
	DDHGCalWafer8 ()

void	execute (DDCompactView &cpv) override

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

	~DDHGCalWafer8 () override

Private Attributes
std::vector< std::string >	cellNames_

int	cellType_

std::string	material_

double	mouseBite_

std::string	nameSpace_

int	nCells_

double	waferSepar_

double	waferSize_

double	waferT_

Detailed Description

Definition at line 9 of file DDHGCalWafer8.h.

Constructor & Destructor Documentation

DDHGCalWafer8::DDHGCalWafer8 ( )

Definition at line 12 of file DDHGCalWafer8.cc.

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

DDHGCalWafer8::~DDHGCalWafer8 ( )

override

Definition at line 18 of file DDHGCalWafer8.cc.

18 {}

Member Function Documentation

void DDHGCalWafer8::execute ( DDCompactView & cpv )

override

Definition at line 47 of file DDHGCalWafer8.cc.

References cellNames_, cellType_, popcon2dropbox::copy(), DDBase< N, C >::ddname(), DDSplit(), DDSolidFactory::extrudedpolygon(), plotBeamSpotDB::first, gen::k, material_, DDBase< N, C >::name(), nCells_, class-composition::parent, DDCompactView::position(), dttmaxenums::R, alignCSCRings::r, makeMuonMisalignmentScenario::rot, Scenarios_cff::scale, edm::second(), mathSSE::sqrt(), findQualityFiles::v, waferSepar_, waferSize_, waferT_, and zw().

                                               {
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom") << "==>> Executing DDHGCalWafer8...";
 #endif
 
   static const double sqrt3 = std::sqrt(3.0);
   double rM = 0.5 * (waferSize_ + waferSepar_);
   double RM2 = rM / sqrt3;
   double R = waferSize_ / (3.0 * nCells_);
   double r = 0.5 * R * sqrt3;
 
   // First the full cell
   std::vector<double> xM = {rM, 0, -rM, -rM, 0, rM};
   std::vector<double> yM = {RM2, 2 * RM2, RM2, -RM2, -2 * RM2, -RM2};
   std::vector<double> zw = {-0.5 * waferT_, 0.5 * waferT_};
   std::vector<double> zx(2, 0), zy(2, 0), scale(2, 1.0);
   DDName parentName = parent().name();
   DDSolid solid =
       DDSolidFactory::extrudedpolygon(parentName, xM, yM, zw, zx, zy, scale);
   DDName matName(DDSplit(material_).first, DDSplit(material_).second);
   DDMaterial matter(matName);
   DDLogicalPart glog = DDLogicalPart(solid.ddname(), matter, solid);
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom")
       << "DDHGCalWafer8: " << solid.name() << " extruded polygon made of "
       << matName << " z|x|y|s (0) " << zw[0] << ":" << zx[0] << ":" << zy[0]
       << ":" << scale[0] << " z|x|y|s (1) " << zw[1] << ":" << zx[1] << ":"
       << zy[1] << ":" << scale[1] << " and " << xM.size() << " edges";
   for (unsigned int k = 0; k < xM.size(); ++k)
     edm::LogVerbatim("HGCalGeom") << "[" << k << "] " << xM[k] << ":" << yM[k];
 #endif
 
   DDRotation rot;
   for (int u = 0; u < 2 * nCells_; ++u) {
     for (int v = 0; v < 2 * nCells_; ++v) {
       if (((v - u) < nCells_) && (u - v) <= nCells_) {
         int n2 = nCells_ / 2;
         double yp = (u - 0.5 * v - n2) * 2 * r;
         double xp = (1.5 * (v - nCells_) + 1.0) * R;
         int cell(0);
         if ((u == 0) && (v == 0))
           cell = 7;
         else if ((u == 0) && (v == nCells_ - 1))
           cell = 8;
         else if ((u == nCells_) && (v == 2 * nCells_ - 1))
           cell = 9;
         else if ((u == 2 * nCells_ - 1) && (v == 2 * nCells_ - 1))
           cell = 10;
         else if ((u == 2 * nCells_ - 1) && (v == nCells_ - 1))
           cell = 11;
         else if ((u == nCells_) && (v == 0))
           cell = 12;
         else if (u == 0)
           cell = 1;
         else if ((v - u) == (nCells_ - 1))
           cell = 4;
         else if (v == (2 * nCells_ - 1))
           cell = 2;
         else if (u == (2 * nCells_ - 1))
           cell = 5;
         else if ((u - v) == nCells_)
           cell = 3;
         else if (v == 0)
           cell = 6;
         DDTranslation tran(xp, yp, 0);
         int copy = (cellType_ * 100 + v) * 100 + u;
         cpv.position(DDName(cellNames_[cell]), glog, copy, tran, rot);
 #ifdef EDM_ML_DEBUG
         edm::LogVerbatim("HGCalGeom")
             << "DDHGCalWafer8: " << cellNames_[cell] << " number " << copy
             << " position in " << glog.name() << " at " << tran << " with "
             << rot;
 #endif
       }
     }
   }
 }

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

override

Definition at line 20 of file DDHGCalWafer8.cc.

References cellNames_, cellType_, createfilelist::int, gen::k, material_, mouseBite_, nameSpace_, nCells_, DDCurrentNamespace::ns(), class-composition::parent, waferSepar_, waferSize_, and waferT_.

                                                                       {
   waferSize_ = nArgs["WaferSize"];
   waferT_ = nArgs["WaferThick"];
   waferSepar_ = nArgs["SensorSeparation"];
   mouseBite_ = nArgs["MouseBite"];
   nCells_ = (int)(nArgs["NCells"]);
   cellType_ = (int)(nArgs["CellType"]);
   material_ = sArgs["Material"];
   cellNames_ = vsArgs["CellNames"];
   nameSpace_ = DDCurrentNamespace::ns();
 #ifdef EDM_ML_DEBUG
   edm::LogVerbatim("HGCalGeom")
       << "DDHGCalWafer8: Wafer 2r " << waferSize_ << " T " << waferT_
       << " Half Separation " << waferSepar_ << " Mouse bite radius "
       << mouseBite_ << " Cells/Wafer " << nCells_ << " Cell Type " << cellType_
       << " Material " << material_ << " Names " << parent().name()
       << " NameSpace " << nameSpace_ << " # of cells " << cellNames_.size();
   for (unsigned int k = 0; k < cellNames_.size(); ++k)
     edm::LogVerbatim("HGCalGeom")
         << "DDHGCalWafer8: Cell[" << k << "] " << cellNames_[k];
 #endif
 }