DDDividedGeometryObject.cc

Go to the documentation of this file.

#include "DetectorDescription/Parser/src/DDDividedGeometryObject.h"
#include "DetectorDescription/Core/interface/DDCompactView.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/Utilities/interface/Exception.h"
#include "Math/GenVector/RotationZ.h"

#include <iostream>
#include <utility>

DDDividedGeometryObject::DDDividedGeometryObject(const DDDivision& div, DDCompactView* cpv)
    : div_(div),
      ftype_(),
      compNDiv_(div.nReplicas()),
      compWidth_(div.width()),
      divisionType_(DivNDIVandWIDTH),
      theVoluFirstCopyNo_(1),
      cpv_(cpv) {
  if (div_.nReplicas() == 0 || div_.width() < tolerance()) {
    if (div_.width() < tolerance())
      divisionType_ = DivNDIV;
    else
      divisionType_ = DivWIDTH;
  }
}

std::unique_ptr<DDRotationMatrix> DDDividedGeometryObject::changeRotMatrix(double rotZ) const {
  return std::make_unique<DDRotationMatrix>(ROOT::Math::RotationZ(rotZ));
}

int DDDividedGeometryObject::calculateNDiv(double motherDim, double width, double offset) const {
  return int((motherDim - offset) / width);
}

double DDDividedGeometryObject::calculateWidth(double motherDim, int nDiv, double offset) const {
  return (motherDim - offset) / nDiv;
}

void DDDividedGeometryObject::checkParametersValidity(void) {
  double maxPar = getMaxParameter();
  checkOffset(maxPar);
  checkNDivAndWidth(maxPar);
  if (!div_.parent().isDefined().second) {
    std::string s = "DDDividedGeometryObject::checkParametersValidity() :";
    s += "\n ERROR - the LogicalPart of the parent must be ";
    s += "\n         defined before a division can occur.";
    s += "\n         Parent= " + div_.parent().toString();
    throw cms::Exception("DDException") << s;
  }
}

void DDDividedGeometryObject::checkOffset(double maxPar) {
  if (div_.offset() >= maxPar) {
    std::string s = "DDDividedGeometryObject::checkOffset() IllegalConstruct";
    s += "\nERROR - DDDividedGeometryObject::checkOffset()";
    s += "\n        failed.";
    s += "  Too big an offset.";
    throw cms::Exception("DDException") << s;
  }
}

void DDDividedGeometryObject::checkNDivAndWidth(double maxPar) {
  if ((divisionType_ == DivNDIVandWIDTH) && (div_.offset() + compWidth_ * compNDiv_ - maxPar > tolerance())) {
    std::string s = "ERROR - DDDividedGeometryObject::checkNDivAndWidth()";
    s += "\n        Division of LogicalPart " + div_.parent().name().name();
    s += " has too big an offset.";

    std::cout << compWidth_ << std::endl;
    throw cms::Exception("DDException") << s;
  }
}

const double DDDividedGeometryObject::tolerance(void) {
  // this can come from some global tolerance if you want.
  static const double tol = 1.0 / 1000.00;
  return tol;
}

void DDDividedGeometryObject::setType(const std::string& s) { ftype_ = s; }

const std::string& DDDividedGeometryObject::getType(void) const { return ftype_; }

void DDDividedGeometryObject::execute(void) {
  for (int i = theVoluFirstCopyNo_; i < compNDiv_ + theVoluFirstCopyNo_; ++i) {
    cpv_->position(makeDDLogicalPart(i), div_.parent(), i, makeDDTranslation(i), makeDDRotation(i), &div_);
  }
}

double DDDividedGeometryObject::getMaxParameter(void) const { return 0.0; }

DDRotation DDDividedGeometryObject::makeDDRotation(const int copyNo) const { return DDRotation(); }

DDTranslation DDDividedGeometryObject::makeDDTranslation(const int copyNo) const { return DDTranslation(); }

DDLogicalPart DDDividedGeometryObject::makeDDLogicalPart(const int copyNo) const {
  // just return the parent... this is USELESS
  return div_.parent();
}