DDHGCalNoTaperEndcap Class Reference

Inheritance diagram for DDHGCalNoTaperEndcap:

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

Private Member Functions
int	createQuarter (DDCompactView &cpv, int xQuadrant, int yQuadrant, int startCopyNo)

Private Attributes
std::string	m_childName
std::string	m_idNameSpace
int	m_incrCopyNo
int	m_invert
int	m_n
double	m_rMax
double	m_rMin
double	m_startAngle
int	m_startCopyNo
double	m_tiltAngle
double	m_xyoffset
double	m_zoffset

Detailed Description

Definition at line 21 of file DDHGCalNoTaperEndcap.cc.

Constructor & Destructor Documentation

DDHGCalNoTaperEndcap::DDHGCalNoTaperEndcap

(

void

)

Definition at line 50 of file DDHGCalNoTaperEndcap.cc.

                                           {
  edm::LogVerbatim("HGCalGeom") << "DDHGCalNoTaperEndcap test: Creating an instance";
}

Member Function Documentation

int DDHGCalNoTaperEndcap::createQuarter ( DDCompactView &  cpv, int  xQuadrant, int  yQuadrant, int  startCopyNo  )

private

Definition at line 102 of file DDHGCalNoTaperEndcap.cc.

References funct::abs(), DDcreateRotationMatrix(), DDrot(), SpecificationBuilder_cfi::parent(), DDCompactView::position(), idealTransformation::rotation, mathSSE::sqrt(), AlCaHLTBitMon_QueryRunRegistry::string, and theta().

                                                                                                         {
  int copyNo = startCopyNo;
  double tiltAngle = m_tiltAngle;
  double xphi = xQuadrant * tiltAngle;
  double yphi = yQuadrant * tiltAngle;
  double theta = 90._deg;
  double phiX = 0.0;
  double phiY = theta;
  double phiZ = 3 * theta;
  double offsetZ = m_zoffset;
  double offsetXY = m_xyoffset;

  double offsetX = xQuadrant * 0.5 * offsetXY;
  double offsetY = yQuadrant * 0.5 * offsetXY;

#ifdef EDM_ML_DEBUG
  int rowmax(0), column(0);
#endif
  while (std::abs(offsetX) < m_rMax) {
#ifdef EDM_ML_DEBUG
    column++;
    int row(0);
#endif
    while (std::abs(offsetY) < m_rMax) {
#ifdef EDM_ML_DEBUG
      row++;
#endif
      double limit1 = sqrt((offsetX + 0.5 * xQuadrant * offsetXY) * (offsetX + 0.5 * xQuadrant * offsetXY) +
                           (offsetY + 0.5 * yQuadrant * offsetXY) * (offsetY + 0.5 * yQuadrant * offsetXY));
      double limit2 = sqrt((offsetX - 0.5 * xQuadrant * offsetXY) * (offsetX - 0.5 * xQuadrant * offsetXY) +
                           (offsetY - 0.5 * yQuadrant * offsetXY) * (offsetY - 0.5 * yQuadrant * offsetXY));
      // Make sure we do not add supermodules in rMin area
      if (limit2 > m_rMin && limit1 < m_rMax) {
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << m_childName << " copyNo = " << copyNo << " (" << column << "," << row
                                      << "): offsetX,Y = " << offsetX << "," << offsetY << " limit=" << limit1 << ":"
                                      << limit2 << " rMin, rMax = " << m_rMin << "," << m_rMax;
#endif
        DDRotation rotation;
        std::string rotstr("NULL");

        // Check if we've already created the rotation matrix
        rotstr = "R";
        rotstr += std::to_string(copyNo);
        rotation = DDRotation(DDName(rotstr));
        if (!rotation) {
          rotation = DDrot(DDName(rotstr, m_idNameSpace),
                           std::make_unique<DDRotationMatrix>(
                               *DDcreateRotationMatrix(theta, phiX, theta + yphi, phiY, -yphi, phiZ) *
                               (*DDcreateRotationMatrix(theta + xphi, phiX, 90._deg, 90._deg, xphi, 0.0))));
        }

        DDTranslation tran(offsetX, offsetY, offsetZ);
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << "Module " << copyNo << ": location = " << tran << " Rotation " << rotation;
#endif
        DDName parentName = parent().name();
        cpv.position(DDName(m_childName), parentName, copyNo, tran, rotation);

        copyNo += m_incrCopyNo;
      } else {
#ifdef EDM_ML_DEBUG
        edm::LogVerbatim("HGCalGeom") << " (" << column << "," << row << "): offsetX,Y = " << offsetX << "," << offsetY
                                      << " is out of limit=" << limit1 << ":" << limit2 << " rMin, rMax = " << m_rMin
                                      << "," << m_rMax;
#endif
      }

      yphi += yQuadrant * 2. * tiltAngle;
      offsetY += yQuadrant * offsetXY;
    }
#ifdef EDM_ML_DEBUG
    if (row > rowmax)
      rowmax = row;
#endif
    xphi += xQuadrant * 2. * tiltAngle;
    yphi = yQuadrant * tiltAngle;
    offsetY = yQuadrant * 0.5 * offsetXY;
    offsetX += xQuadrant * offsetXY;
  }
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << rowmax << " rows and " << column << " columns in quadrant " << xQuadrant << ":"
                                << yQuadrant;
#endif
  return copyNo;
}

void DDHGCalNoTaperEndcap::execute ( DDCompactView &  cpv )

override

Definition at line 82 of file DDHGCalNoTaperEndcap.cc.

                                                     {
  int lastCopyNo = m_startCopyNo;
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Create quarter 1:1";
#endif
  lastCopyNo = createQuarter(cpv, 1, 1, lastCopyNo);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Create quarter -1:1";
#endif
  lastCopyNo = createQuarter(cpv, -1, 1, lastCopyNo);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Create quarter -1:-1";
#endif
  lastCopyNo = createQuarter(cpv, -1, -1, lastCopyNo);
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Create quarter 1:-1";
#endif
  createQuarter(cpv, 1, -1, lastCopyNo);
}

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

override

Definition at line 54 of file DDHGCalNoTaperEndcap.cc.

References DDCurrentNamespace::ns(), and SpecificationBuilder_cfi::parent().

                                                                      {
  m_tiltAngle = nArgs["tiltAngle"];
  m_invert = int(nArgs["invert"]);
  m_rMin = int(nArgs["rMin"]);
  m_rMax = int(nArgs["rMax"]);
  m_zoffset = nArgs["zoffset"];
  m_xyoffset = nArgs["xyoffset"];
  m_n = int(nArgs["n"]);
  m_startCopyNo = int(nArgs["startCopyNo"]);
  m_incrCopyNo = int(nArgs["incrCopyNo"]);
  m_childName = sArgs["ChildName"];
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "Tilt Angle " << m_tiltAngle << " Invert " << m_invert << " R " << m_rMin << ":"
                                << m_rMax << " Offset " << m_zoffset << ":" << m_xyoffset << " N " << m_n << " Copy "
                                << m_startCopyNo << ":" << m_incrCopyNo << " Child " << m_childName;
#endif

  m_idNameSpace = DDCurrentNamespace::ns();
#ifdef EDM_ML_DEBUG
  edm::LogVerbatim("HGCalGeom") << "DDHGCalNoTaperEndcap: NameSpace " << m_idNameSpace << "\tParent "
                                << parent().name();
#endif
}

Member Data Documentation

std::string DDHGCalNoTaperEndcap::m_childName

private

Definition at line 46 of file DDHGCalNoTaperEndcap.cc.

std::string DDHGCalNoTaperEndcap::m_idNameSpace

private

Definition at line 47 of file DDHGCalNoTaperEndcap.cc.

int DDHGCalNoTaperEndcap::m_incrCopyNo

private

Definition at line 45 of file DDHGCalNoTaperEndcap.cc.

int DDHGCalNoTaperEndcap::m_invert

private

Definition at line 38 of file DDHGCalNoTaperEndcap.cc.

int DDHGCalNoTaperEndcap::m_n

private

Definition at line 43 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_rMax

private

Definition at line 40 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_rMin

private

Definition at line 39 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_startAngle

private

Definition at line 36 of file DDHGCalNoTaperEndcap.cc.

int DDHGCalNoTaperEndcap::m_startCopyNo

private

Definition at line 44 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_tiltAngle

private

Definition at line 37 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_xyoffset

private

Definition at line 42 of file DDHGCalNoTaperEndcap.cc.

double DDHGCalNoTaperEndcap::m_zoffset

private

Definition at line 41 of file DDHGCalNoTaperEndcap.cc.