GeoPixelIBLFwdServices.cxx

Go to the documentation of this file.

/*
  Copyright (C) 2002-2023 CERN for the benefit of the ATLAS collaboration
*/
 
 
// Build IBL fwd services (wavy shape)
// This is built one time per layer. 
 
#include "GeoPixelIBLFwdServices.h"
 
#include "GeoPixelIBLFwdSvcCADModel.h"
#include "GeoPixelIBLFwdSvcModel1.h"
 
#include "GeoModelKernel/GeoPhysVol.h"
#include "GeoModelKernel/GeoMaterial.h"
 
 
#include <algorithm>
#include <iostream> 
#include <iomanip> 
using std::max;
 
GeoPixelIBLFwdServices::GeoPixelIBLFwdServices(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,
                                               int section)
  : GeoVPixelFactory(ddmgr, mgr, sqliteReader, mapFPV, mapAX),
  m_bFwdComplexGeometry_CAD(false),
  m_bFwdComplexGeometry_Mod1(false),
  m_section(section),
  m_supportPhysA(nullptr),
  m_supportPhysC(nullptr),
  m_xformSupportA(nullptr),
  m_xformSupportC(nullptr)
{
}
 
 
bool GeoPixelIBLFwdServices::isComplexShapeDefined()
{
  // check if the wavy shape option is defined in the DB (based on the definition of the material)
 
  bool bFwd02Defined=false;
  int nbBarrelSvc = m_gmt_mgr->PixelServiceElements("barrel");
  for(int i=0; i<nbBarrelSvc; i++){
    std::string volName = m_gmt_mgr->PixelServiceName("barrel",i);
    if(volName.substr(0,9)=="IBL_Fwd02")bFwd02Defined=true;
  }
  int nbSimpleSvc = m_gmt_mgr->PixelServiceElements("simple");
  for(int i=0; i<nbSimpleSvc; i++){
    std::string volName = m_gmt_mgr->PixelServiceName("simple",i);
    if(volName.substr(0,9)=="IBL_Fwd02")bFwd02Defined=true;
  }
 
  if(m_mat_mgr->hasMaterial("pix::IBL_Fwd02_Cable_Wvy") &&!bFwd02Defined) m_bFwdComplexGeometry_CAD = true;
  if(m_mat_mgr->hasMaterial("pix::IBL_Fwd02_Cable_Wvy_M1")  && !bFwd02Defined) m_bFwdComplexGeometry_Mod1 = true;
 
  m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - CAD tool design / S. model1 : "<<  m_bFwdComplexGeometry_CAD<<"/"<<m_bFwdComplexGeometry_Mod1<<endmsg;
 
 
  return m_bFwdComplexGeometry_CAD||m_bFwdComplexGeometry_Mod1;
 
}
 
 
GeoVPhysVol* GeoPixelIBLFwdServices::Build()
{
 
  m_gmt_mgr->msg(MSG::INFO) <<"GeoPixelIBLFwdService - CAD tool design / S. model1 : "<<  m_bFwdComplexGeometry_CAD<<"/"<<m_bFwdComplexGeometry_Mod1<<endmsg;
  
  if(m_bFwdComplexGeometry_CAD) {
 
    m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - CAD tool design"<<endmsg;
    GeoPixelIBLFwdSvcCADModel fwdSvc = GeoPixelIBLFwdSvcCADModel(m_DDmgr, m_gmt_mgr ,m_sqliteReader, m_mapFPV, m_mapAX, m_section);
    
    fwdSvc.Build();
    m_supportPhysA = fwdSvc.getSupportA();
    m_supportPhysC = fwdSvc.getSupportC();
    m_xformSupportA = fwdSvc.getSupportTrfA();
    m_xformSupportC = fwdSvc.getSupportTrfC();
 
    if(m_supportPhysA==nullptr) {
      m_gmt_mgr->msg(MSG::ERROR) <<"Build IBL fwd services - CAD tool design : the PixelGeoModel package is not meant to be used to build this design"<<endmsg;
      m_gmt_mgr->msg(MSG::ERROR) <<"                                         : check that the GeoTorus object is defined in GeoModelKernel/GeoModel"<<endmsg;
    }
    return nullptr;
  }
 
 
  if(m_bFwdComplexGeometry_Mod1) {
 
    m_gmt_mgr->msg(MSG::INFO) <<"Build IBL fwd services - S. Menke  design"<<endmsg;
    GeoPixelIBLFwdSvcModel1 fwdSvc = GeoPixelIBLFwdSvcModel1(m_DDmgr, m_gmt_mgr, m_sqliteReader, m_mapFPV, m_mapAX,  m_section);
    
    fwdSvc.Build();
    m_supportPhysA = fwdSvc.getSupportA();
    m_supportPhysC = fwdSvc.getSupportC();
    m_xformSupportA = fwdSvc.getSupportTrfA();
    m_xformSupportC = fwdSvc.getSupportTrfC();
 
    return nullptr;
  }
    
  return nullptr;
 
}