DetectorDesign.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2021 CERN for the benefit of the ATLAS collaboration
*/
 
#include "ReadoutGeometryBase/DetectorDesign.h"
#include "ReadoutGeometryBase/SiIntersect.h"
 
namespace InDetDD {
 
// Constructor with parameters:
DetectorDesign::DetectorDesign(double thickness,
                               bool phiSymmetric,
                               bool etaSymmetric,
                               bool depthSymmetric,
                               InDetDD::CarrierType carrierType,
                               int readoutSide):
    m_etaAxis(zAxis),
    m_phiAxis(yAxis),
    m_depthAxis(xAxis),
    m_thickness(thickness),
    m_carrierType(carrierType),
    m_phiSymmetric(phiSymmetric),
    m_etaSymmetric(etaSymmetric),
    m_depthSymmetric(depthSymmetric),
    m_readoutSidePosDepth(readoutSide > 0) {
}
 
// Constructor with axis parameters:
DetectorDesign::DetectorDesign(double thickness,
                               bool phiSymmetric,
                               bool etaSymmetric,
                               bool depthSymmetric,
                               InDetDD::CarrierType carrierType,
                               int readoutSide,
                               Axis stripDirection,
                               Axis thicknessDirection):
    m_etaAxis(stripDirection),
    m_depthAxis(thicknessDirection),
    m_thickness(thickness),
    m_carrierType(carrierType),
    m_phiSymmetric(phiSymmetric),
    m_etaSymmetric(etaSymmetric),
    m_depthSymmetric(depthSymmetric),
    m_readoutSidePosDepth(readoutSide > 0) {
    if (stripDirection == thicknessDirection) {
        throw std::runtime_error(
                  "ERROR: DetectorDesign called with phi and thickness directions equal");
    }
    // phiAxis is "the other one"
    m_phiAxis = static_cast<Axis> ((xAxis + yAxis + zAxis) - (stripDirection + thicknessDirection));
 
    
 
}
 
DetectorDesign::~DetectorDesign()= default;
 
 
 
Amg::Vector3D DetectorDesign::sensorCenter() const {
    return Amg::Vector3D(0., 0., 0.);
}
 
//For future: Should update to use appropriate message streams in place of cout
void DetectorDesign::setSymmetry(bool phiSymmetric, bool etaSymmetric,
                                 bool depthSymmetric) {
    // Flags can be changed from true to false but not false to true.
    if (m_phiSymmetric) {
        m_phiSymmetric = phiSymmetric;
    }
    else if (phiSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xPhi axis direction ignored."
                  << std::endl;
    }
 
    if (m_etaSymmetric) {
        m_etaSymmetric = etaSymmetric;
    }
    else if (etaSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xEta axis direction ignored."
                  << std::endl;
    }
 
    if (m_depthSymmetric) {
        m_depthSymmetric = depthSymmetric;
    }
    else if (depthSymmetric) {
        std::cout <<
            "SiDetectorDesign: WARNING! Attempt to allow swapping of xDepth axis direction ignored."
                  << std::endl;
    }
}
 
DetectorShape DetectorDesign::shape() const {
    // Default is Box.
    return InDetDD::Box;
}
 
DetectorType DetectorDesign::type() const {
    // Default is Undefined.
    return InDetDD::Undefined;
}
 
SiIntersect DetectorDesign::inDetector(const SiLocalPosition &localPosition,
                                       double phiTol, double etaTol) const {
    double etaDist = 0;
    double phiDist = 0;
 
    distanceToDetectorEdge(localPosition, etaDist, phiDist);
 
    SiIntersect state;
 
    if (phiDist < -phiTol || etaDist < -etaTol) {
        state.setOut();
        return state;
    }
 
    if (phiDist > phiTol && etaDist > etaTol) {
        state.setIn();
        return state;
    }
 
    // Near boundary.
    state.setNearBoundary();
    return state;
}
 
  const Amg::Transform3D DetectorDesign::moduleShift() const {
    return Amg::Transform3D::Identity();
  }
 
} // namespace InDetDD