GeoPixelStaveRing.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
// Build detailed stave support : face plate + carbon foam + cable flex + cooling pipe + end blocks
// This is built one time per layer. 
 
#include "GeoPixelStaveRing.h"
 
#include "GeoPixelSiCrystal.h"
#include "GeoPixelModule.h"
 
#include "GeoModelKernel/GeoBox.h"
#include "GeoModelKernel/GeoTube.h"
#include "GeoModelKernel/GeoSimplePolygonBrep.h"
#include "GeoModelKernel/GeoLogVol.h"
#include "GeoModelKernel/GeoPhysVol.h"
#include "GeoModelKernel/GeoMaterial.h"
#include "GeoModelKernel/GeoNameTag.h"
#include "GaudiKernel/SystemOfUnits.h"
 
#include <algorithm>
using std::max;
using namespace std;
 
GeoPixelStaveRing::GeoPixelStaveRing(InDetDD::PixelDetectorManager* ddmgr,
                                     PixelGeometryManager* mgr,
                     GeoModelIO::ReadGeoModel* sqliteReader,
                                     std::shared_ptr<std::map<std::string, GeoFullPhysVol*>> mapFPV,
                                     std::shared_ptr<std::map<std::string, GeoAlignableTransform*>> mapAX)
  : GeoVPixelFactory (ddmgr, mgr, sqliteReader, mapFPV, mapAX),
    m_physVol (nullptr),
    m_zPosition (0),
    m_innerRadius (0),
    m_outerRadius (0),
    m_ringPosition ("AC"),
    m_ringName ("staveRing")
{
}
 
GeoVPhysVol* GeoPixelStaveRing::SetParametersAndBuild(const std::string& ringName,
                                                      const std::string& ringPos)
{
  m_ringName=ringName;
  m_ringPosition=ringPos;
  m_physVol = Build();
 
  return m_physVol;
}
 
GeoVPhysVol* GeoPixelStaveRing::Build(){
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build detailed stave ring support : "<<m_ringName<<"  "<<m_ringPosition<<endmsg;
 
  double safety = 0.001*Gaudi::Units::mm; 
  bool isBLayer = false;
  if(m_gmt_mgr->GetLD() == 0) isBLayer = true;
  GeoPixelSiCrystal theSensor(m_DDmgr, m_gmt_mgr, m_sqliteReader, m_mapFPV, m_mapAX, isBLayer);
  GeoPixelModule pm(m_DDmgr, m_gmt_mgr, m_sqliteReader,  m_mapFPV, m_mapAX, theSensor);
  if(m_sqliteReader) return nullptr;
  // Ladder geometry
 
  double endBlockFixingPoint= m_gmt_mgr->IBLStaveMechanicalStaveEndBlockFixPoint();
  double totalStaveLength = m_gmt_mgr->IBLStaveLength();
 
  if(m_ringPosition=="AC")
    {
      // Stave ring geometry
      double ringWidth = m_gmt_mgr->IBLSupportRingWidth();
      m_innerRadius = m_gmt_mgr->IBLSupportRingInnerRadius();
      m_outerRadius = m_gmt_mgr->IBLSupportRingOuterRadius();
      double fixingPoint = m_gmt_mgr->IBLSupportMechanicalStaveRingFixPoint();
 
      // SES (dec 2012 - IBL pcache ): values hard coded as the one defined in ATLAS-IBL-02-01-00
      //   to get rid off the overlap between staverings ans ladder with ATLAS-IBL-02-00-00 values
      if(m_innerRadius<29.56) m_innerRadius=29.56;     // instead of 29.2
      if(m_outerRadius>31.85) m_outerRadius=31.85;     // instead of 32.
 
      double endBlockPosition = totalStaveLength*0.5+endBlockFixingPoint;
      double ringPosition = -ringWidth*0.5+fixingPoint;
      m_zPosition = endBlockPosition-ringPosition;
 
      double IPTouterRadius = m_gmt_mgr->IBLServiceGetMaxRadialPosition("IPT","simple",m_zPosition,m_zPosition)+safety;
      m_gmt_mgr->msg(MSG::DEBUG)<<"IBL stave ring "<<m_zPosition<<" "<<m_innerRadius<<"  "<<IPTouterRadius<<endmsg;
      if(IPTouterRadius>m_innerRadius) m_innerRadius=IPTouterRadius;
 
      // create log and phys volumes
      m_gmt_mgr->msg(MSG::DEBUG)<<"-> IBL stave ring "<<m_zPosition<<" "<<m_innerRadius<<"  "<<m_outerRadius<<endmsg;
      const GeoTube* ring_tube = new GeoTube(m_innerRadius,m_outerRadius,ringWidth*0.5);
      const GeoMaterial* ring_material_weight = m_mat_mgr->getMaterialForVolume("pix::StaveRing_IBLwght",ring_tube->volume());
      m_gmt_mgr->msg(MSG::DEBUG)<<"IBL stave ring weighted material : "<<(ring_material_weight==nullptr)<<endmsg;
 
      GeoLogVol * logVol = nullptr;
      if(ring_material_weight)
    logVol = new GeoLogVol(m_ringName,ring_tube,ring_material_weight);
      else {
    const GeoMaterial* ring_material = m_mat_mgr->getMaterial("pix::StaveRing_IBL");
    logVol = new GeoLogVol(m_ringName,ring_tube,ring_material);
      }
      GeoPhysVol * logVolPV = new GeoPhysVol(logVol);
  m_gmt_mgr->msg(MSG::INFO) <<"Build detailed stave ring support - logVol : "<<logVol->getName()<<endmsg;
 
      return logVolPV;
    }
 
  // Stave ring geometry
  double ringWidth = m_gmt_mgr->IBLSupportMidRingWidth();
  m_innerRadius = m_gmt_mgr->IBLSupportMidRingInnerRadius();
  m_outerRadius = m_gmt_mgr->IBLSupportMidRingOuterRadius();
  
  m_zPosition = 0.0;
  
  // create log and phys volumes
  const GeoTube* ring_tube = new GeoTube(m_innerRadius,m_outerRadius,ringWidth*0.5);
  const GeoMaterial* ring_material_weight = m_mat_mgr->getMaterialForVolume("pix::StaveRingMid_IBLwght",ring_tube->volume());
  GeoLogVol *logVol=nullptr;
  if(ring_material_weight)
    logVol = new GeoLogVol(m_ringName,ring_tube,ring_material_weight);
  else {
    const GeoMaterial* ring_material = m_mat_mgr->getMaterial("pix::MiddleRing_IBL");
    logVol = new GeoLogVol(m_ringName,ring_tube,ring_material);
  }
 
  GeoPhysVol * logVolPV = new GeoPhysVol(logVol); 
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build detailed stave ring support - logVol : "<<logVol->getName()<<endmsg;
 
  return logVolPV;
 
}