#include <SCT_Barrel.h>

Inheritance diagram for SCT_Barrel:

Collaboration diagram for SCT_Barrel:

Public Member Functions
	SCT_Barrel (const std::string &name, InDetDD::SCT_DetectorManager detectorManager, SCT_GeometryManager geometryManager, SCT_MaterialManager materials, GeoModelIO::ReadGeoModel sqliteReader, std::shared_ptr< std::map< std::string, GeoFullPhysVol * > > mapFPV, std::shared_ptr< std::map< std::string, GeoAlignableTransform * > > mapAX)
virtual GeoVPhysVol *	build (SCT_Identifier id)
int	numLayers () const
double	innerRadius () const
double	outerRadius () const
double	length () const
const std::string &	getName () const

Protected Member Functions
double	epsilon () const

Protected Attributes
const GeoLogVol *	m_logVolume {}
GeoModelIO::ReadGeoModel *	m_sqliteReader {}
std::shared_ptr< std::map< std::string, GeoFullPhysVol * > >	m_mapFPV
std::shared_ptr< std::map< std::string, GeoAlignableTransform * > >	m_mapAX
InDetDD::SCT_DetectorManager *	m_detectorManager {}
SCT_GeometryManager *	m_geometryManager {}
SCT_MaterialManager *	m_materials {}

Private Member Functions
void	getParameters ()
virtual const GeoLogVol *	preBuild ()
void	buildThermalShield (GeoFullPhysVol *parent)
void	buildEMIShield (GeoFullPhysVol *parent)

Private Attributes
int	m_numLayers = 0
double	m_innerRadius = 0.0
double	m_outerRadius = 0.0
double	m_length = 0.0
double	m_thermalShieldEndWallThickness = 0.0
double	m_zClearance = 0.0
bool	m_isOldGeometry = false
std::string	m_name

Static Private Attributes
static const double	s_epsilon = 1.0e-6 * Gaudi::Units::mm

Detailed Description

Definition at line 16 of file SCT_Barrel.h.

Constructor & Destructor Documentation

◆ SCT_Barrel()

SCT_Barrel::SCT_Barrel	(	const std::string &	name,
		InDetDD::SCT_DetectorManager *	detectorManager,
		SCT_GeometryManager *	geometryManager,
		SCT_MaterialManager *	materials,
		GeoModelIO::ReadGeoModel *	sqliteReader,
		std::shared_ptr< std::map< std::string, GeoFullPhysVol * > >	mapFPV,
		std::shared_ptr< std::map< std::string, GeoAlignableTransform * > >	mapAX )

Definition at line 44 of file SCT_Barrel.cxx.

  : SCT_UniqueComponentFactory(name, detectorManager, geometryManager, materials, sqliteReader, std::move(mapFPV), std::move(mapAX))
{
  getParameters();
  if(!m_sqliteReader) {
    m_logVolume = SCT_Barrel::preBuild();
  }
}

Member Function Documentation

◆ build()

GeoVPhysVol * SCT_Barrel::build ( SCT_Identifier id )

virtual

Implements SCT_UniqueComponentFactory.

Definition at line 96 of file SCT_Barrel.cxx.

{
    GeoFullPhysVol * barrel=nullptr;
    if(!m_sqliteReader)
    {
         barrel = new GeoFullPhysVol(m_logVolume);
        
        // Old geometries are no longer supported - give up now if one is requested
        if(m_isOldGeometry) {
            MsgStream log(Athena::getMessageSvc(), "SCT_Barrel");
            log<< MSG::INFO <<"SCT_Barrel Old barrel geometry versions are not supported"<<endmsg;
            return barrel;
        }
        
        // There is only one type of module. So we create it just the once and pass it to the layers.
        SCT_Module module("Module", m_detectorManager, m_geometryManager, m_materials, m_sqliteReader, m_mapFPV, m_mapAX);
        
        // Create the interlinks
        SCT_InterLink interLink("InterLink", m_detectorManager, m_geometryManager, m_materials);
        
        // Calculte the length of the layer cylinder. This is the barrel length less the thermal
        // shield and interlink width.
        // This is only used for 'OldGeometry". In new geometry, layer length is set internally,
        // and is equal to support cylinder length
        double layerLength = m_length - 2*m_thermalShieldEndWallThickness - 2*interLink.length();
        
        // We reduce to allow some alignment clearance
        layerLength -= 2*m_zClearance;
        
        for (int iLayer = 0; iLayer < m_numLayers; iLayer++) {
            
            // Create the layers
            
            layerLength = 0.;
            SCT_Layer layer(std::format("Layer{}",iLayer), iLayer, &module, m_detectorManager, m_geometryManager, m_materials, m_sqliteReader, m_mapFPV, m_mapAX);
            barrel->add(new GeoNameTag(std::format("Layer#{}",iLayer)));
            barrel->add(new GeoIdentifierTag(iLayer)); // Identifier layer= iLayer
            id.setLayerDisk(iLayer);
            GeoAlignableTransform * transform = new GeoAlignableTransform(GeoTrf::Transform3D::Identity());
            barrel->add(transform);
            GeoVPhysVol * layerPV = layer.build(id);
            barrel->add(layerPV);
            // Store alignable transform
            m_detectorManager->addAlignableTransform(2, id.getWaferId(), transform, layerPV);
            layerLength = std::max(layerLength,layer.length());
        }
        
        // Build and place the interlinks
        double interLinkZPos = 0.;
        interLinkZPos = 0.5 * layerLength + m_zClearance + 0.5 * interLink.length();
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(+interLinkZPos)));
        barrel->add(interLink.getVolume());
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(-interLinkZPos)));
        barrel->add(interLink.getVolume());
        
        // Build and place the cooling spiders
        double spiderZPos = 0.;
        SCT_Spider spider("Spider", m_detectorManager, m_geometryManager, m_materials);
        spiderZPos =  interLinkZPos + 0.5*interLink.length() + 0.5*spider.length();
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(+spiderZPos)));
        barrel->add(spider.getVolume());
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(-spiderZPos)));
        barrel->add(spider.getVolume());
        
        // Build and place the thermal shield.
        buildThermalShield(barrel);
        
        // Build and place the EMI shield (inner thermal shield).
        buildEMIShield(barrel);
        
        // Build and place SCT to Pixel attachment
        SCT_PixelAttachment pixelAttachment("AttachmentPixelToSCT", m_detectorManager, m_geometryManager, m_materials);
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(+pixelAttachment.zPosition()))); // +ve z
        barrel->add(pixelAttachment.getVolume());
        barrel->add(new GeoTransform(GeoTrf::TranslateZ3D(-pixelAttachment.zPosition()))); // -ve z
        barrel->add(pixelAttachment.getVolume());
        
        // Extra Material
        InDetDD::ExtraMaterial xMat(m_geometryManager->distortedMatManager());
        xMat.add(barrel, "SCTBarrel");
        
    }
    else {
      // There is only one type of module. So we create it just the once and pass it to the layers.
      SCT_Module module("Module", m_detectorManager, m_geometryManager, nullptr, m_sqliteReader, m_mapFPV, m_mapAX);
 
      for (int iLayer = 0; iLayer < m_numLayers; iLayer++) {
    // Create the layers
    SCT_Layer layer(std::format("Layer{}",iLayer), iLayer, &module, m_detectorManager, m_geometryManager, m_materials, m_sqliteReader, m_mapFPV, m_mapAX);
    id.setLayerDisk(iLayer);
    layer.build(id); //MB to verify
    // Store alignable transform
    m_detectorManager->addAlignableTransform(2, id.getWaferId(), (*m_mapAX)[std::format("Layer#{}",iLayer)], (*m_mapFPV)[std::format("Layer#{}",iLayer)]);
      }
    }
    
    return barrel;
}

◆ buildEMIShield()

void SCT_Barrel::buildEMIShield ( GeoFullPhysVol * parent )

private

Definition at line 255 of file SCT_Barrel.cxx.

{
  
  const SCT_BarrelParameters * parameters = m_geometryManager->barrelParameters();
 
  // Parameters of cylinder 
  double innerRadius       = parameters->emiShieldInnerRadius();
  double deltaR            = parameters->emiShieldDeltaR();
  double outerRadius       = innerRadius + deltaR;
  double length            = 2 * parameters->emiShieldZMax();
  std::string materialName = parameters->emiShieldMaterial();
 
  // Parameters of Electrical Shield Joint
  double jointDeltaR = 0;
  double jointRPhi = 0;
  std::string jointMaterialName;
  if(!m_isOldGeometry) {
    jointDeltaR            = parameters->emiJointDeltaR();
    jointRPhi              = parameters->emiJointRPhi();
    jointMaterialName = parameters->emiJointMaterial();
  }
 
  // Parameters of Pixel Attachment - needed for cut-out
  double pixelAttachmentLength  = parameters->pixelAttachmentDeltaZ();
  double pixelAttachmentZpos    = parameters->pixelAttachmentZMin() + 0.5 * pixelAttachmentLength; 
 
  // Build cylinder (with cut-outs)
  const GeoShape * emiShieldShape  = nullptr;
  const GeoMaterial * material;
  const GeoTube * emiShieldTube  = new GeoTube(innerRadius, outerRadius,  0.5*length);
  if (m_isOldGeometry) {
    emiShieldShape = emiShieldTube;
    material = m_materials->getMaterial(materialName);
  } else {
    const GeoTube* cutOut = new GeoTube(innerRadius, outerRadius, 0.5*pixelAttachmentLength);
    const GeoShape* emiTemp = (GeoShape*)&(emiShieldTube->subtract(*cutOut << GeoTrf::TranslateZ3D(pixelAttachmentZpos)));
    emiShieldShape = (GeoShape*)&emiTemp->subtract(*cutOut << GeoTrf::TranslateZ3D(-pixelAttachmentZpos));
    double emiVolume = emiShieldTube->volume() - 2. * cutOut->volume();
    material = m_materials->getMaterialForVolume(materialName, emiVolume);
  }
  const GeoLogVol  * emiShieldLog = new GeoLogVol("EMIShield", emiShieldShape, material);
  GeoPhysVol * emiShield = new GeoPhysVol(emiShieldLog);
  parent->add(emiShield);
 
  // Build electrical shield joints (with cut-outs)
  if (!m_isOldGeometry) {
    double dphi = jointRPhi / outerRadius;
    const GeoTubs* emiJointTubs = new GeoTubs(outerRadius, outerRadius+jointDeltaR, 0.5*length,
                                              -0.5 * dphi * Gaudi::Units::radian, dphi * Gaudi::Units::radian);
    const GeoTubs* jointCutOut = new GeoTubs(outerRadius, outerRadius+jointDeltaR, 0.5*pixelAttachmentLength,
                                             -0.5 * dphi * Gaudi::Units::radian, dphi * Gaudi::Units::radian);
    const GeoShape* jointTemp = (GeoShape*)&(emiJointTubs->subtract(*jointCutOut << GeoTrf::TranslateZ3D(pixelAttachmentZpos)));
    const GeoShape* emiJointShape = (GeoShape*)&jointTemp->subtract(*jointCutOut << GeoTrf::TranslateZ3D(-pixelAttachmentZpos));
    double jointVolume = emiJointTubs->volume() - 2. * jointCutOut->volume();
    const GeoMaterial * jointMaterial = m_materials->getMaterialForVolume(jointMaterialName, jointVolume);
    const GeoLogVol  * emiJointLog = new GeoLogVol("EMIShieldJoint", emiJointShape, jointMaterial);
    GeoPhysVol * emiJoint = new GeoPhysVol(emiJointLog);
    // Place 3 copies
    for (int i=0; i<3; i++) {
      double angle = (90. + i * 120.) * Gaudi::Units::degree;
      parent->add(new GeoTransform(GeoTrf::RotateZ3D(angle)));
      parent->add(emiJoint);
    }
  }
}

◆ buildThermalShield()

void SCT_Barrel::buildThermalShield ( GeoFullPhysVol * parent )

private

Definition at line 195 of file SCT_Barrel.cxx.

{
 
  // The thermal shield is now in 3 parts:
  // (a) outer cylinder
  // (b) bulkheads
  // (c) end panels
  // The inner cylinder is called 'EMI shield' for backwards compatibility
 
  const SCT_BarrelParameters * parameters = m_geometryManager->barrelParameters();
 
  double cylinderOuterRadius   = parameters->thermalShieldOuterRadius();
  double cylinderInnerRadius   = parameters->thermalShieldInnerRadius();
  double cylinderLength        = parameters->cylinderLength();
  double bulkheadInnerRadius   = parameters->thermalShieldBulkheadInnerRadius();
  double bulkheadOuterRadius   = parameters->thermalShieldBulkheadOuterRadius();
  double bulkheadThickness     = parameters->thermalShieldEndCapCylThickness();
  double endPanelInnerRadius   = parameters->thermalShieldEndPanelInnerRadius();
  double endPanelOuterRadius   = parameters->thermalShieldEndPanelOuterRadius();
  double endPanelThickness     = parameters->thermalShieldEndCapThickness();
  double endPanelZMax          = parameters->thermalShieldEndZMax();
  
  std::string cylinderMaterialName = parameters->thermalShieldMaterialCyl();
  std::string bulkheadMaterialName = parameters->thermalShieldMaterialOuterSect();
  std::string endPanelMaterialName = parameters->thermalShieldMaterialInnerSect();
 
 
  // The outer cylinder part of thermal shield.
  const GeoTube * cylinderShape  = new GeoTube(cylinderInnerRadius, cylinderOuterRadius, 0.5*cylinderLength);
  const GeoMaterial* cylinderMaterial = m_materials->getMaterialForVolume(cylinderMaterialName,cylinderShape->volume());
  const GeoLogVol  * cylinderLog = new GeoLogVol("ThShieldOuterCyl", cylinderShape, cylinderMaterial);
  GeoPhysVol * cylinder = new GeoPhysVol(cylinderLog);
  parent->add(cylinder);
 
  // The bulkheads
  const GeoTube * bulkheadShape  = new GeoTube(bulkheadInnerRadius, bulkheadOuterRadius, 0.5*bulkheadThickness);
  const GeoMaterial* bulkheadMaterial = m_materials->getMaterialForVolume(bulkheadMaterialName,bulkheadShape->volume());
  const GeoLogVol  * bulkheadLog = new GeoLogVol("ThShieldBulkhead", bulkheadShape, bulkheadMaterial);
  GeoPhysVol * bulkhead = new GeoPhysVol(bulkheadLog);
  GeoTransform * bulkheadPosPlus  =  new GeoTransform(GeoTrf::TranslateZ3D(+(endPanelZMax-endPanelThickness-0.5*bulkheadThickness)));
  GeoTransform * bulkheadPosMinus =  new GeoTransform(GeoTrf::TranslateZ3D(-(endPanelZMax-endPanelThickness-0.5*bulkheadThickness)));
  parent->add(bulkheadPosPlus);
  parent->add(bulkhead);
  parent->add(bulkheadPosMinus);
  parent->add(bulkhead);
 
  // The end panels
  const GeoTube * endPanelShape  = new GeoTube(endPanelInnerRadius, endPanelOuterRadius,  0.5*endPanelThickness);
  const GeoMaterial* endPanelMaterial = m_materials->getMaterialForVolume(endPanelMaterialName,endPanelShape->volume());
  const GeoLogVol  * endPanelLog = new GeoLogVol("ThShieldEndPanel", endPanelShape, endPanelMaterial);
  GeoPhysVol * endPanel = new GeoPhysVol(endPanelLog);
  GeoTransform * endPanelPosPlus  =  new GeoTransform(GeoTrf::TranslateZ3D(+(endPanelZMax-0.5*endPanelThickness)));
  GeoTransform * endPanelPosMinus =  new GeoTransform(GeoTrf::TranslateZ3D(-(endPanelZMax-0.5*endPanelThickness)));
  parent->add(endPanelPosPlus);
  parent->add(endPanel);
  parent->add(endPanelPosMinus);
  parent->add(endPanel);
 
}

◆ epsilon()

double SCT_ComponentFactory::epsilon ( ) const

protectedinherited

Definition at line 28 of file SCT_ComponentFactory.cxx.

{
  return s_epsilon;
}

◆ getName()

const std::string & SCT_ComponentFactory::getName ( ) const

inlineinherited

Definition at line 35 of file SCT_ComponentFactory.h.

35{return m_name;}

SCT_ComponentFactory::m_name

std::string m_name

Definition SCT_ComponentFactory.h:46

◆ getParameters()

void SCT_Barrel::getParameters ( )

private

Definition at line 61 of file SCT_Barrel.cxx.

{
    const SCT_BarrelParameters * parameters = m_geometryManager->barrelParameters();
    
    if(!m_sqliteReader){
        m_innerRadius = parameters->barrelInnerRadius();
        m_outerRadius = parameters->barrelOuterRadius();
        m_length =      parameters->barrelLength();
        
        // Used in old geometry
        m_thermalShieldEndWallThickness = parameters->thermalShieldEndCapThickness();
        
        // Clearannce in z between layer and interlink.
        m_zClearance = 1.*Gaudi::Units::mm;
        
        // Layer internal structure and services depend on geometry version
        m_isOldGeometry = parameters->isOldGeometry();
    }
    m_numLayers =   parameters->numLayers();
    // Set numerology
    m_detectorManager->numerology().setNumLayers(m_numLayers);
 
}

◆ innerRadius()

double SCT_Barrel::innerRadius ( ) const

inline

Definition at line 31 of file SCT_Barrel.h.

31{return m_innerRadius;}

◆ length()

double SCT_Barrel::length ( ) const

inline

Definition at line 33 of file SCT_Barrel.h.

33{return m_length;}

◆ numLayers()

int SCT_Barrel::numLayers ( ) const

inline

Definition at line 30 of file SCT_Barrel.h.

30{return m_numLayers;}

◆ outerRadius()

double SCT_Barrel::outerRadius ( ) const

inline

Definition at line 32 of file SCT_Barrel.h.

32{return m_outerRadius;}

◆ preBuild()

const GeoLogVol * SCT_Barrel::preBuild ( )

privatevirtual

Implements SCT_UniqueComponentFactory.

Definition at line 86 of file SCT_Barrel.cxx.

{
  // Create the barrel volume
  // Tube envelope containing the barrel.
  const GeoTube * barrelEnvelopeShape = new GeoTube(m_innerRadius, m_outerRadius, 0.5 * m_length);
  GeoLogVol * barrelLog = new GeoLogVol(getName(), barrelEnvelopeShape, m_materials->gasMaterial());
  return barrelLog;
}

Member Data Documentation

◆ m_detectorManager

InDetDD::SCT_DetectorManager* SCT_ComponentFactory::m_detectorManager {}

protectedinherited

Definition at line 38 of file SCT_ComponentFactory.h.

38{};

◆ m_geometryManager

SCT_GeometryManager* SCT_ComponentFactory::m_geometryManager {}

protectedinherited

Definition at line 39 of file SCT_ComponentFactory.h.

39{};

◆ m_innerRadius

double SCT_Barrel::m_innerRadius = 0.0

private

Definition at line 45 of file SCT_Barrel.h.

◆ m_isOldGeometry

bool SCT_Barrel::m_isOldGeometry = false

private

Definition at line 53 of file SCT_Barrel.h.

◆ m_length

double SCT_Barrel::m_length = 0.0

private

Definition at line 47 of file SCT_Barrel.h.

◆ m_logVolume

const GeoLogVol* SCT_UniqueComponentFactory::m_logVolume {}

protectedinherited

Definition at line 86 of file SCT_ComponentFactory.h.

86{};

◆ m_mapAX

std::shared_ptr<std::map<std::string, GeoAlignableTransform*> > SCT_UniqueComponentFactory::m_mapAX

protectedinherited

Definition at line 92 of file SCT_ComponentFactory.h.

◆ m_mapFPV

std::shared_ptr<std::map<std::string, GeoFullPhysVol*> > SCT_UniqueComponentFactory::m_mapFPV

protectedinherited

Definition at line 91 of file SCT_ComponentFactory.h.

◆ m_materials

SCT_MaterialManager* SCT_ComponentFactory::m_materials {}

protectedinherited

Definition at line 40 of file SCT_ComponentFactory.h.

40{};

◆ m_name

std::string SCT_ComponentFactory::m_name

privateinherited

Definition at line 46 of file SCT_ComponentFactory.h.

◆ m_numLayers

int SCT_Barrel::m_numLayers = 0

private

Definition at line 43 of file SCT_Barrel.h.

◆ m_outerRadius

double SCT_Barrel::m_outerRadius = 0.0

private

Definition at line 46 of file SCT_Barrel.h.

◆ m_sqliteReader

GeoModelIO::ReadGeoModel* SCT_UniqueComponentFactory::m_sqliteReader {}

protectedinherited

Definition at line 87 of file SCT_ComponentFactory.h.

87{};

◆ m_thermalShieldEndWallThickness

double SCT_Barrel::m_thermalShieldEndWallThickness = 0.0

private

Definition at line 49 of file SCT_Barrel.h.

◆ m_zClearance

double SCT_Barrel::m_zClearance = 0.0

private

Definition at line 51 of file SCT_Barrel.h.

◆ s_epsilon

const double SCT_ComponentFactory::s_epsilon = 1.0e-6 * Gaudi::Units::mm

staticprivateinherited

Definition at line 47 of file SCT_ComponentFactory.h.

The documentation for this class was generated from the following files:

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Static Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ SCT_Barrel()

Member Function Documentation

◆ build()

◆ buildEMIShield()

◆ buildThermalShield()

◆ epsilon()

◆ getName()

◆ getParameters()

◆ innerRadius()

◆ length()

◆ numLayers()

◆ outerRadius()

◆ preBuild()

Member Data Documentation

◆ m_detectorManager

◆ m_geometryManager

◆ m_innerRadius

◆ m_isOldGeometry

◆ m_length

◆ m_logVolume

◆ m_mapAX

◆ m_mapFPV

◆ m_materials

◆ m_name

◆ m_numLayers

◆ m_outerRadius

◆ m_sqliteReader

◆ m_thermalShieldEndWallThickness

◆ m_zClearance

◆ s_epsilon