DDTOBAxCableAlgo.cc

Go to the documentation of this file.

// File: DDTOBAxCableAlgo.cc
// Description: Equipping the axial cylinder of TOB with cables etc

#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Core/interface/DDCurrentNamespace.h"
#include "DetectorDescription/Core/interface/DDSplit.h"
#include "DetectorDescription/Core/interface/DDTypes.h"
#include "DetectorDescription/Core/interface/DDAlgorithm.h"
#include "DetectorDescription/Core/interface/DDAlgorithmFactory.h"
#include <CLHEP/Units/GlobalPhysicalConstants.h>
#include <CLHEP/Units/SystemOfUnits.h>

#include <string>
#include <vector>

using namespace std;

class DDTOBAxCableAlgo : public DDAlgorithm {
public:
  //Constructor and Destructor
  DDTOBAxCableAlgo();
  ~DDTOBAxCableAlgo() override;

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

  void execute(DDCompactView& cpv) override;

private:
  string idNameSpace;  // Namespace of this and ALL sub-parts

  vector<string> sectorNumber;  // Id. Number of the sectors

  double sectorRin;                 // Inner radius of service sectors
  double sectorRout;                // Outer radius of service sectors
  double sectorDz;                  // Sector half-length
  double sectorDeltaPhi_B;          // Sector B phi width [A=C=0.5*(360/sectors)]
  vector<double> sectorStartPhi;    // Starting phi for the service sectors
  vector<string> sectorMaterial_A;  // Material for the A sectors
  vector<string> sectorMaterial_B;  // Material for the B sectors
  vector<string> sectorMaterial_C;  // Material for the C sectors
};

DDTOBAxCableAlgo::DDTOBAxCableAlgo() : sectorRin(0), sectorRout(0), sectorDeltaPhi_B(0) {
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo info: Creating an instance";
}

DDTOBAxCableAlgo::~DDTOBAxCableAlgo() {}

void DDTOBAxCableAlgo::initialize(const DDNumericArguments& nArgs,
                                  const DDVectorArguments& vArgs,
                                  const DDMapArguments&,
                                  const DDStringArguments& sArgs,
                                  const DDStringVectorArguments& vsArgs) {
  idNameSpace = DDCurrentNamespace::ns();
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: Parent " << parent().name() << " NameSpace " << idNameSpace;

  sectorNumber = vsArgs["SectorNumber"];
  sectorRin = nArgs["SectorRin"];
  sectorRout = nArgs["SectorRout"];
  sectorDz = nArgs["SectorDz"];
  sectorDeltaPhi_B = nArgs["SectorDeltaPhi_B"];
  sectorStartPhi = vArgs["SectorStartPhi"];
  sectorMaterial_A = vsArgs["SectorMaterial_A"];
  sectorMaterial_B = vsArgs["SectorMaterial_B"];
  sectorMaterial_C = vsArgs["SectorMaterial_C"];

  for (int i = 0; i < (int)(sectorNumber.size()); i++)
    LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: sectorNumber[" << i << "] = " << sectorNumber[i];

  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: Axial Service Sectors half-length " << sectorDz << "\tRin "
                      << sectorRin << "\tRout = " << sectorRout << "\tPhi of sectors position:";
  for (int i = 0; i < (int)(sectorNumber.size()); i++)
    LogDebug("TOBGeom") << "\t[" << i << "]\tPhi = " << sectorStartPhi[i];
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: List of materials for the sectors/3 parts";
  //
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 A";
  for (int i = 0; i < (int)(sectorNumber.size()); i++)
    LogDebug("TOBGeom") << "\t[" << i << "]\tsectorMaterial_A = " << sectorMaterial_A[i];
  //
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 B";
  for (int i = 0; i < (int)(sectorNumber.size()); i++)
    LogDebug("TOBGeom") << "\t[" << i << "]\tsectorMaterial_B = " << sectorMaterial_B[i];
  //
  LogDebug("TOBGeom") << "DDTOBAxCableAlgo debug: Sector/3 C";
  for (int i = 0; i < (int)(sectorNumber.size()); i++)
    LogDebug("TOBGeom") << "\t[" << i << "]\tsectorMaterial_C = " << sectorMaterial_C[i];
}

void DDTOBAxCableAlgo::execute(DDCompactView& cpv) {
  LogDebug("TOBGeom") << "==>> Constructing DDTOBAxCableAlgo...";
  DDName tubsName(parent().name());

  // Loop over sectors (sectorNumber vector)
  for (int i = 0; i < (int)(sectorNumber.size()); i++) {
    DDSolid solid;
    string name;
    double dz, rin, rout, startphi, widthphi, deltaphi;

    // Axial Services
    // Each sector is divided in 3 parts from phi[i] to phi[i+1]

    widthphi = ((i + 1 == (int)(sectorStartPhi.size())) ? (sectorStartPhi[0] + CLHEP::twopi) - sectorStartPhi[i]
                                                        : (sectorStartPhi[i + 1] - sectorStartPhi[i]));

    // First Part: A
    name = "TOBAxService_" + sectorNumber[i] + "A";
    dz = sectorDz;
    rin = sectorRin;
    rout = sectorRout;
    startphi = sectorStartPhi[i];
    deltaphi = 0.5 * (widthphi - sectorDeltaPhi_B);

    solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, startphi, deltaphi);

    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " Tubs made of "
                        << sectorMaterial_A[i] << " from " << startphi / CLHEP::deg << " to "
                        << (startphi + deltaphi) / CLHEP::deg << " with Rin " << rin << " Rout " << rout << " ZHalf "
                        << dz;

    DDName sectorMatName(DDSplit(sectorMaterial_A[i]).first, DDSplit(sectorMaterial_A[i]).second);
    DDMaterial sectorMatter(sectorMatName);
    DDLogicalPart sectorLogic(DDName(name, idNameSpace), sectorMatter, solid);

    cpv.position(DDName(name, idNameSpace), tubsName, i + 1, DDTranslation(), DDRotation());
    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " number " << i + 1
                        << " positioned in " << tubsName << " with no translation and no rotation";

    // Second Part: B
    name = "TOBAxService_" + sectorNumber[i] + "B";
    startphi += deltaphi;
    deltaphi = sectorDeltaPhi_B;

    solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, startphi, deltaphi);

    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " Tubs made of "
                        << sectorMaterial_B[i] << " from " << startphi / CLHEP::deg << " to "
                        << (startphi + deltaphi) / CLHEP::deg << " with Rin " << rin << " Rout " << rout << " ZHalf "
                        << dz;

    sectorMatName = DDName(DDSplit(sectorMaterial_B[i]).first, DDSplit(sectorMaterial_B[i]).second);
    sectorMatter = DDMaterial(sectorMatName);
    sectorLogic = DDLogicalPart(DDName(name, idNameSpace), sectorMatter, solid);

    cpv.position(DDName(name, idNameSpace), tubsName, i + 1, DDTranslation(), DDRotation());
    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " number " << i + 1
                        << " positioned in " << tubsName << " with no translation and no rotation";

    // Third Part: C
    name = "TOBAxService_" + sectorNumber[i] + "C";
    startphi += deltaphi;
    deltaphi = 0.5 * (widthphi - sectorDeltaPhi_B);

    solid = DDSolidFactory::tubs(DDName(name, idNameSpace), dz, rin, rout, startphi, deltaphi);

    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " Tubs made of "
                        << sectorMaterial_C[i] << " from " << startphi / CLHEP::deg << " to "
                        << (startphi + deltaphi) / CLHEP::deg << " with Rin " << rin << " Rout " << rout << " ZHalf "
                        << dz;

    sectorMatName = DDName(DDSplit(sectorMaterial_C[i]).first, DDSplit(sectorMaterial_C[i]).second);
    sectorMatter = DDMaterial(sectorMatName);
    sectorLogic = DDLogicalPart(DDName(name, idNameSpace), sectorMatter, solid);

    cpv.position(DDName(name, idNameSpace), tubsName, i + 1, DDTranslation(), DDRotation());
    LogDebug("TOBGeom") << "DDTOBAxCableAlgo test: " << DDName(name, idNameSpace) << " number " << i + 1
                        << " positioned in " << tubsName << " with no translation and no rotation";
  }

  LogDebug("TOBGeom") << "<<== End of DDTOBAxCableAlgo construction ...";
}

DEFINE_EDM_PLUGIN(DDAlgorithmFactory, DDTOBAxCableAlgo, "track:DDTOBAxCableAlgo");