DDTBH4Algo.cc

Go to the documentation of this file.

// File: DDTBH4Algo.cc
// Description: Position inside the mother according to (eta,phi)

#include <cmath>
#include <algorithm>

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/DDCurrentNamespace.h"
#include "DetectorDescription/Core/interface/DDSplit.h"
#include "Geometry/EcalTestBeam/plugins/DDTBH4Algo.h"
#include <CLHEP/Units/SystemOfUnits.h>

using CLHEP::deg;

DDTBH4Algo::DDTBH4Algo()
    : m_idNameSpace(""),
      m_BLZBeg(0),
      m_BLZEnd(0),
      m_BLZPiv(0),
      m_BLRadius(0),
      m_VacName(""),
      m_VacMat(""),
      m_vecVacZBeg(),
      m_vecVacZEnd(),
      m_WinName(""),
      m_vecWinMat(),
      m_vecWinZBeg(),
      m_vecWinThick(),
      m_TrgMat(""),
      m_HoleMat(""),
      m_TrgVetoHoleRadius(0),
      m_vecTrgName(),
      m_vecTrgSide(),
      m_vecTrgThick(),
      m_vecTrgPhi(),
      m_vecTrgXOff(),
      m_vecTrgYOff(),
      m_vecTrgZPiv(),
      m_FibFibName(""),
      m_FibCladName(""),
      m_FibFibMat(""),
      m_FibCladMat(""),
      m_FibSide(0),
      m_FibCladThick(0),
      m_FibLength(0),
      m_vecFibPhi(),
      m_vecFibXOff(),
      m_vecFibYOff(),
      m_vecFibZPiv()

{
  edm::LogVerbatim("EcalGeom") << "creating an instance if DDTBH4Algo";
  LogDebug("EcalGeom") << "DDTBH4Algo test: Creating an instance";
}

DDTBH4Algo::~DDTBH4Algo() {}

DDRotation DDTBH4Algo::myrot(const std::string& s, const CLHEP::HepRotation& r) const {
  return DDrot(
      ddname(idNameSpace() + ":" + s),
      std::make_unique<DDRotationMatrix>(r.xx(), r.xy(), r.xz(), r.yx(), r.yy(), r.yz(), r.zx(), r.zy(), r.zz()));
}

DDMaterial DDTBH4Algo::ddmat(const std::string& s) const { return DDMaterial(ddname(s)); }

DDName DDTBH4Algo::ddname(const std::string& s) const {
  const std::pair<std::string, std::string> temp(DDSplit(s));
  return DDName(temp.first, temp.second);
}

void DDTBH4Algo::initialize(const DDNumericArguments& nArgs,
                            const DDVectorArguments& vArgs,
                            const DDMapArguments& mArgs,
                            const DDStringArguments& sArgs,
                            const DDStringVectorArguments& vsArgs) {
  m_idNameSpace = DDCurrentNamespace::ns();
  m_BLZBeg = nArgs["BLZBeg"];
  m_BLZEnd = nArgs["BLZEnd"];
  m_BLZPiv = nArgs["BLZPiv"];
  m_BLRadius = nArgs["BLRadius"];
  m_VacName = sArgs["VacName"];
  m_VacMat = sArgs["VacMat"];
  m_vecVacZBeg = vArgs["VacZBeg"];
  m_vecVacZEnd = vArgs["VacZEnd"];

  m_WinName = sArgs["WinName"];
  m_vecWinMat = vsArgs["WinMat"];
  m_vecWinZBeg = vArgs["WinZBeg"];
  m_vecWinThick = vArgs["WinThick"];

  m_TrgMat = sArgs["TrgMat"];
  m_HoleMat = sArgs["HoleMat"];
  m_TrgVetoHoleRadius = nArgs["TrgVetoHoleRadius"];
  m_vecTrgName = vsArgs["TrgName"];
  m_vecTrgSide = vArgs["TrgSide"];
  m_vecTrgThick = vArgs["TrgThick"];
  m_vecTrgPhi = vArgs["TrgPhi"];
  m_vecTrgXOff = vArgs["TrgXOff"];
  m_vecTrgYOff = vArgs["TrgYOff"];
  m_vecTrgZPiv = vArgs["TrgZPiv"];

  m_FibFibName = sArgs["FibFibName"];
  m_FibCladName = sArgs["FibCladName"];
  m_FibFibMat = sArgs["FibFibMat"];
  m_FibCladMat = sArgs["FibCladMat"];
  m_FibSide = nArgs["FibSide"];
  m_FibCladThick = nArgs["FibCladThick"];
  m_FibLength = nArgs["FibLength"];
  m_vecFibPhi = vArgs["FibPhi"];
  m_vecFibXOff = vArgs["FibXOff"];
  m_vecFibYOff = vArgs["FibYOff"];
  m_vecFibZPiv = vArgs["FibZPiv"];
}

void DDTBH4Algo::execute(DDCompactView& cpv) {
  const unsigned int copyOne(1);

  const double halfZbl((blZEnd() - blZBeg()) / 2.);
  for (unsigned int i(0); i != vecVacZBeg().size(); ++i) {
    DDName vacNameNm(ddname(vacName() + std::to_string(i + 1)));
    const double halfZvac((vecVacZEnd()[i] - vecVacZBeg()[i]) / 2.);
    DDSolid vTubeSolid(DDSolidFactory::tubs(vacNameNm, halfZvac, 0, blRadius(), 0 * deg, 360 * deg));
    const DDLogicalPart vacLog(vacNameNm, vacMat(), vTubeSolid);

    cpv.position(vacLog,
                 parent().name(),
                 1 + i,
                 DDTranslation(0, 0, -halfZbl + halfZvac + vecVacZBeg()[i] - blZBeg()),
                 DDRotation());
  }

  for (unsigned int i(0); i != vecWinZBeg().size(); ++i) {
    DDName wName(ddname(winName() + std::to_string(i + 1)));
    DDSolid wTubeSolid(DDSolidFactory::tubs(wName, vecWinThick()[i] / 2., 0, blRadius(), 0 * deg, 360 * deg));
    const DDLogicalPart wLog(wName, ddmat(vecWinMat()[i]), wTubeSolid);

    const double off(0 < vecWinZBeg()[i] ? vecWinZBeg()[i] : fabs(vecWinZBeg()[i]) - vecWinThick()[i]);

    cpv.position(wLog,
                 parent().name(),
                 1 + i,
                 DDTranslation(0, 0, -halfZbl + vecWinThick()[i] / 2. + off - blZBeg()),
                 DDRotation());
  }

  for (unsigned int i(0); i != vecTrgName().size(); ++i) {
    DDName tName(ddname(vecTrgName()[i]));
    DDSolid tSolid(DDSolidFactory::box(tName, vecTrgSide()[i] / 2., vecTrgSide()[i] / 2., vecTrgThick()[i] / 2.));
    const DDLogicalPart tLog(tName, trgMat(), tSolid);

    if (tName.name() == "VETO") {
      DDName vName(ddname(tName.name() + "Hole"));
      DDSolid vTubeSolid(
          DDSolidFactory::tubs(vName, vecTrgThick()[i] / 2., 0, trgVetoHoleRadius(), 0 * deg, 360 * deg));
      const DDLogicalPart vLog(vName, holeMat(), vTubeSolid);

      cpv.position(vLog, tName, copyOne, DDTranslation(0, 0, 0), DDRotation());
    }

    cpv.position(tLog,
                 parent().name(),
                 copyOne,
                 DDTranslation(vecTrgXOff()[i], vecTrgYOff()[i], vecTrgZPiv()[i] - halfZbl + blZPiv() - blZBeg()),
                 myrot(tName.name() + "Rot", CLHEP::HepRotationZ(vecTrgPhi()[i])));
  }

  DDName pName(fibCladName());
  const double planeWidth(32.5 * fibSide() + 33.5 * fibCladThick());
  const double planeThick(2 * fibSide() + 3 * fibCladThick());
  DDSolid pSolid(DDSolidFactory::box(pName, planeWidth / 2., fibLength() / 2., planeThick / 2.));
  const DDLogicalPart pLog(pName, fibCladMat(), pSolid);

  DDSolid fSolid(DDSolidFactory::box(fibFibName(), fibSide() / 2., fibLength() / 2., fibSide() / 2.));

  const DDLogicalPart fLog(fibFibName(), fibFibMat(), fSolid);

  for (unsigned int j(0); j != 32; ++j) {
    const double xoff(planeWidth / 2. - (1 + j) * fibCladThick() - (1 + j) * fibSide());
    const double zoff(-planeThick / 2 + fibCladThick() + fibSide() / 2.);
    cpv.position(fLog, pName, 1 + j, DDTranslation(xoff, 0, zoff), DDRotation());

    cpv.position(fLog, pName, 33 + j, DDTranslation(xoff + (fibCladThick() + fibSide()) / 2., 0, -zoff), DDRotation());
  }
  for (unsigned int i(0); i != vecFibZPiv().size(); ++i) {
    cpv.position(
        pLog,
        parent().name(),
        1 + i,
        DDTranslation(
            vecFibXOff()[i] - 0.5 * fibSide(), vecFibYOff()[i], vecFibZPiv()[i] - halfZbl + blZPiv() - blZBeg()),
        myrot(pName.name() + "Rot" + std::to_string(i), CLHEP::HepRotationZ(vecFibPhi()[i])));
  }
}