DetectorDescription/Parser/src/DDDividedCons.cc - master - githubmirror/cmssw

DetectorDescription/Parser/src/DDDividedCons.cc
218 строк8 KB
Cms Build
Clang-Format
11 май 2019, 12:09
11 май 2019, 12:0902cea70
#include "DetectorDescription/Parser/src/DDDividedCons.h"
#include "DetectorDescription/Core/interface/DDRotationMatrix.h"
#include "DetectorDescription/Core/interface/DDAxes.h"
#include "DetectorDescription/Core/interface/DDLogicalPart.h"
#include "DetectorDescription/Core/interface/DDMaterial.h"
#include "DetectorDescription/Core/interface/DDName.h"
#include "DetectorDescription/Core/interface/DDSolid.h"
#include "DetectorDescription/Core/interface/DDTransform.h"
#include "DataFormats/Math/interface/GeantUnits.h"
#include "DetectorDescription/Parser/src/DDDividedGeometryObject.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"

#include <iostream>
#include <string>
#include <utility>

class DDCompactView;

using namespace geant_units::operators;

DDDividedConsRho::DDDividedConsRho(const DDDivision& div, DDCompactView* cpv)
    : DDDividedGeometryObject::DDDividedGeometryObject(div, cpv) {
  checkParametersValidity();
  setType("DivisionConsRho");
  DDCons msol = (DDCons)(div_.parent().solid());

  if (msol.rInPlusZ() == 0.) {
    std::cout << "WARNING - DDDividedConsRho, rOutMinusZ = 0. " << std::endl
              << "          Width is calculated as that of rOutMinusZ !" << std::endl;
  }

  if (divisionType_ == DivWIDTH) {
    compNDiv_ = calculateNDiv(msol.rOutMinusZ() - msol.rInMinusZ(), div_.width(), div_.offset());
  } else if (divisionType_ == DivNDIV) {
    DDCons msol = (DDCons)(div_.parent().solid());
    compWidth_ = calculateWidth(msol.rOutMinusZ() - msol.rInMinusZ(), div_.nReplicas(), div_.offset());
  }
}

double DDDividedConsRho::getMaxParameter(void) const {
  DDCons msol = (DDCons)(div_.parent().solid());
  return msol.rOutMinusZ() - msol.rInMinusZ();
}

DDRotation DDDividedConsRho::makeDDRotation(const int copyNo) const {
  DDRotation myddrot;  // sets to identity.
  return myddrot;
}

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

DDLogicalPart DDDividedConsRho::makeDDLogicalPart(const int copyNo) const {
  DDName solname(div_.parent().ddname().name() + "_DIVCHILD" + std::to_string(copyNo), div_.parent().ddname().ns());
  DDSolid ddcons(solname);
  DDMaterial usemat(div_.parent().material());
  DDCons msol = (DDCons)(div_.parent().solid());

  double pRMin1 = msol.rInMinusZ() + div_.offset() + compWidth_ * copyNo;
  double pRMax1 = msol.rInMinusZ() + div_.offset() + compWidth_ * (copyNo + 1);

  //width at Z Plus
  //- double compWidth_Plus =
  //   compWidth_ * ( msol.rOutPlusZ()/ msol.rInPlusZ())
  //-         / ( msol.rOutMinusZ() - msol.rInMinusZ());
  double compWidth_Plus = calculateWidth(msol.rOutPlusZ() - msol.rInPlusZ(), compNDiv_, div_.offset());
  double pRMin2 = msol.rInPlusZ() + div_.offset() + compWidth_Plus * copyNo;
  double pRMax2 = msol.rInPlusZ() + div_.offset() + compWidth_Plus * (copyNo + 1);
  double pDz = msol.zhalf();

  double pSPhi = msol.phiFrom();
  double pDPhi = msol.deltaPhi();

  ddcons = DDSolidFactory::cons(DDName(solname), pDz, pRMin1, pRMax1, pRMin2, pRMax2, pSPhi, pDPhi);

  DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);
  return ddlp;
}

DDDividedConsPhi::DDDividedConsPhi(const DDDivision& div, DDCompactView* cpv)
    : DDDividedGeometryObject::DDDividedGeometryObject(div, cpv) {
  checkParametersValidity();
  setType("DivisionConsPhi");
  DDCons msol = (DDCons)(div_.parent().solid());

  if (divisionType_ == DivWIDTH) {
    DDCons msol = (DDCons)(div_.parent().solid());
    //If you divide a tube of 360 degrees the offset displaces the starting angle, but you still fill the 360 degrees
    if (msol.deltaPhi() == 360._deg) {
      compNDiv_ = calculateNDiv(msol.deltaPhi(), div_.width(), 0.);
    } else {
      compNDiv_ = calculateNDiv(msol.deltaPhi(), div_.width(), div_.offset());
    }
  } else if (divisionType_ == DivNDIV) {
    DDCons msol = (DDCons)(div_.parent().solid());
    if (msol.deltaPhi() == 360._deg) {
      compWidth_ = calculateWidth(msol.deltaPhi(), div_.nReplicas(), 0.);
    } else {
      compWidth_ = calculateWidth(msol.deltaPhi(), div_.nReplicas(), div_.offset());
    }
  }
}

double DDDividedConsPhi::getMaxParameter(void) const {
  DDCons msol = (DDCons)(div_.parent().solid());
  return msol.deltaPhi();
}

DDRotation DDDividedConsPhi::makeDDRotation(const int copyNo) const {
  DDRotation myddrot;  // sets to identity.
  double posi = (copyNo - 1) * compWidth_;
  // how to name the rotation??
  // i hate this crap :-)
  DDName ddrotname(div_.parent().ddname().name() + "_DIVCHILD_ROT" + std::to_string(copyNo),
                   div_.parent().ddname().ns());
  myddrot = DDrot(ddrotname, changeRotMatrix(posi));

  return myddrot;
}

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

DDLogicalPart DDDividedConsPhi::makeDDLogicalPart(const int copyNo) const {
  DDName solname(div_.parent().ddname().name() + "_DIVCHILD", div_.parent().ddname().ns());
  DDSolid ddcons(solname);
  DDMaterial usemat(div_.parent().material());
  DDCons msol = (DDCons)(div_.parent().solid());

  if (!ddcons.isDefined().second) {
    double pRMin1 = msol.rInMinusZ();
    double pRMax1 = msol.rOutMinusZ();
    double pRMin2 = msol.rInPlusZ();
    double pRMax2 = msol.rOutPlusZ();
    double pDz = msol.zhalf();

    //- already rotated  double pSPhi = div_.offset() + copyNo*compWidth_;
    double pSPhi = div_.offset() + msol.phiFrom();
    double pDPhi = compWidth_;
    ddcons = DDSolidFactory::cons(DDName(solname), pDz, pRMin1, pRMax1, pRMin2, pRMax2, pSPhi, pDPhi);
  }

  DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);

  return ddlp;
}

DDDividedConsZ::DDDividedConsZ(const DDDivision& div, DDCompactView* cpv)
    : DDDividedGeometryObject::DDDividedGeometryObject(div, cpv) {
  checkParametersValidity();

  DDCons msol = (DDCons)(div_.parent().solid());
  setType("DivisionConsZ");

  if (divisionType_ == DivWIDTH) {
    DDCons msol = (DDCons)(div_.parent().solid());
    compNDiv_ = calculateNDiv(2 * msol.zhalf(), div_.width(), div_.offset());
  } else if (divisionType_ == DivNDIV) {
    DDCons msol = (DDCons)(div_.parent().solid());
    compWidth_ = calculateWidth(2 * msol.zhalf(), div_.nReplicas(), div_.offset());
  }
}

double DDDividedConsZ::getMaxParameter(void) const {
  DDCons msol = (DDCons)(div_.parent().solid());
  return 2 * msol.zhalf();
}

DDRotation DDDividedConsZ::makeDDRotation(const int copyNo) const {
  DDRotation myddrot;  // sets to identity.
  return myddrot;
}

DDTranslation DDDividedConsZ::makeDDTranslation(const int copyNo) const {
  DDTranslation translation;

  DDCons motherCons = (DDCons)(div_.parent().solid());
  double posi = -motherCons.zhalf() + div_.offset() + compWidth_ / 2 + copyNo * compWidth_;
  translation.SetZ(posi);

  return translation;
}

DDLogicalPart DDDividedConsZ::makeDDLogicalPart(const int copyNo) const {
  DDName solname(div_.parent().ddname().name() + "_DIVCHILD" + std::to_string(copyNo), div_.parent().ddname().ns());
  DDSolid ddcons(solname);
  DDMaterial usemat(div_.parent().material());
  DDCons msol = (DDCons)(div_.parent().solid());

  double mHalfLength = msol.zhalf();
  double aRInner = (msol.rInPlusZ() - msol.rInMinusZ()) / (2 * mHalfLength);
  double bRInner = (msol.rInPlusZ() + msol.rInMinusZ()) / 2;
  double aROuter = (msol.rOutPlusZ() - msol.rOutMinusZ()) / (2 * mHalfLength);
  double bROuter = (msol.rOutPlusZ() + msol.rOutMinusZ()) / 2;
  double xMinusZ = -mHalfLength + div_.offset() + compWidth_ * copyNo;
  double xPlusZ = -mHalfLength + div_.offset() + compWidth_ * (copyNo + 1);

  double pDz = compWidth_ / 2.;
  double pSPhi = msol.phiFrom();
  double pDPhi = msol.deltaPhi();

  ddcons = DDSolidFactory::cons(DDName(solname),
                                pDz,
                                aRInner * xMinusZ + bRInner,
                                aROuter * xMinusZ + bROuter,
                                aRInner * xPlusZ + bRInner,
                                aROuter * xPlusZ + bROuter,
                                pSPhi,
                                pDPhi);

  DDLogicalPart ddlp = DDLogicalPart(solname, usemat, ddcons);

  return ddlp;
}