#include <GeoPixelServices.h>

Inheritance diagram for GeoPixelServices:

Collaboration diagram for GeoPixelServices:

Public Member Functions
	GeoPixelServices (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, InDetDD::Zone *envelopeZone=0)
	GeoPixelServices (const GeoPixelServices &)=delete
GeoPixelServices &	operator= (const GeoPixelServices &)=delete
	~GeoPixelServices ()
virtual GeoVPhysVol *	Build () override
void	initialize (const std::string &)
void	initializeOld (const std::string &)
InDetDD::VolumeBuilder *	getBuilder ()
InDetDD::VolumeBuilder *	getServMatBuilder ()

Protected Attributes
PixelGeometryManager *	m_gmt_mgr
InDetMaterialManager *	m_mat_mgr
InDetDD::PixelDetectorManager *	m_DDmgr
GeoModelIO::ReadGeoModel *	m_sqliteReader
std::shared_ptr< std::map< std::string, GeoFullPhysVol * > >	m_mapFPV
std::shared_ptr< std::map< std::string, GeoAlignableTransform * > >	m_mapAX
const double	m_epsilon

Private Attributes
std::vector< const InDetDD::ServiceVolume * >	m_services
InDetDD::VolumeBuilder *	m_pixServBuilder
InDetDD::VolumeBuilder *	m_servMatBuilder
std::vector< double >	m_layerShift

Detailed Description

Definition at line 17 of file GeoPixelServices.h.

Constructor & Destructor Documentation

◆ GeoPixelServices() [1/2]

GeoPixelServices::GeoPixelServices	(	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,
		InDetDD::Zone *	envelopeZone = 0 )

Definition at line 114 of file GeoPixelServices.cxx.

  : GeoVPixelFactory(ddmgr, mgr, sqliteReader, std::move(mapFPV), std::move(mapAX)),
    m_pixServBuilder(nullptr),
    m_servMatBuilder(nullptr),
    m_layerShift(0)
{
  
  // Get dimensions of barrel and endcap
  double barrelRmin = m_gmt_mgr->PixelBarrelRMin();
  double barrelRmax = m_gmt_mgr->PixelBarrelRMax();
  double barrelZmin = -m_gmt_mgr->PixelBarrelHalfLength();
  double barrelZmax = m_gmt_mgr->PixelBarrelHalfLength();
  double endcapRmin = m_gmt_mgr->PixelEndcapRMin();
  double endcapRmax = m_gmt_mgr->PixelEndcapRMax();
  double endcapZmin = m_gmt_mgr->PixelEndcapZMin();
  double endcapZmax = m_gmt_mgr->PixelEndcapZMax();
 
  double pixelRmin, pixelRmax, pixelZmax;
  if (!pixZone) {
    pixelRmin = m_gmt_mgr->PixelRMin();
    pixelRmax = m_gmt_mgr->PixelRMax();
    pixelZmax = m_gmt_mgr->PixelHalfLength();
    pixZone = new InDetDD::TubeZone("Pixel",-pixelZmax,pixelZmax,pixelRmin,pixelRmax);
  }
  else {
    if (m_gmt_mgr->PixelSimpleEnvelope()) {
      pixelRmin = m_gmt_mgr->PixelRMin();
      pixelRmax = m_gmt_mgr->PixelRMax();
      pixelZmax = m_gmt_mgr->PixelHalfLength();
    }
    else {
      pixelZmax = m_gmt_mgr->PixelEnvelopeZ(0);
      pixelRmin = m_gmt_mgr->PixelEnvelopeRMin(0);
      pixelRmax = m_gmt_mgr->PixelEnvelopeRMax(0);
    }
  }  
  
  InDetDD::UnboundedZone topZone("Mother");
  pixZone->add(new InDetDD::TubeZone("C",-endcapZmax,-endcapZmin,endcapRmin,endcapRmax, true));
  pixZone->add(new InDetDD::TubeZone("B",barrelZmin,barrelZmax,barrelRmin,barrelRmax));
  pixZone->add(new InDetDD::TubeZone("A",endcapZmin,endcapZmax,endcapRmin,endcapRmax));
  topZone.add(pixZone);
 
  // Collect the layer shifts / IBL 2mm shift
  m_layerShift.clear();
  for(int ii = 0; ii < m_gmt_mgr->PixelBarrelNLayer(); ii++){
    m_gmt_mgr->SetCurrentLD(ii);
    m_layerShift.push_back(m_gmt_mgr->PixelLayerGlobalShift());
  }
 
 
  // We process all tables. 
  bool barrelPresent   = m_gmt_mgr->partPresent("Barrel");
  bool endcapAPresent  = m_gmt_mgr->partPresent("EndcapA");
  bool endcapCPresent  = m_gmt_mgr->partPresent("EndcapC");
 
  if (barrelPresent) initialize("barrel");
  if (endcapAPresent || endcapCPresent) initialize("endcap");
  initialize("simple");
  initialize("envelope");
 
  // all services are added at a time
  m_pixServBuilder = new InDetDD::VolumeBuilder(topZone, m_services);
  m_pixServBuilder->setMaterialManager(m_mat_mgr);
  if (m_gmt_mgr->athenaComps()->serviceBuilderTool()) {
    const std::vector<const InDetDD::ServiceVolume *> & services = m_gmt_mgr->athenaComps()->serviceBuilderTool()->getServices();
    m_servMatBuilder = new InDetDD::VolumeBuilder(topZone, services);
    m_servMatBuilder->setMaterialManager(m_mat_mgr);
  }
}

◆ GeoPixelServices() [2/2]

GeoPixelServices::GeoPixelServices ( const GeoPixelServices & )

delete

◆ ~GeoPixelServices()

GeoPixelServices::~GeoPixelServices ( )

Definition at line 190 of file GeoPixelServices.cxx.

{
  // delete contents of m_multiParam
  for (unsigned int i = 0; i <  m_services.size(); ++i) {
    delete  m_services[i];
  }
  delete m_pixServBuilder;
  delete m_servMatBuilder;
}

Member Function Documentation

◆ Build()

GeoVPhysVol * GeoPixelServices::Build ( )

overridevirtual

Implements GeoVPixelFactory.

Definition at line 356 of file GeoPixelServices.cxx.

{
  return nullptr;
}

◆ getBuilder()

InDetDD::VolumeBuilder * GeoPixelServices::getBuilder ( )

inline

Definition at line 31 of file GeoPixelServices.h.

31{return m_pixServBuilder;}

◆ getServMatBuilder()

InDetDD::VolumeBuilder * GeoPixelServices::getServMatBuilder ( )

inline

Definition at line 32 of file GeoPixelServices.h.

32{return m_servMatBuilder;}

◆ initialize()

void GeoPixelServices::initialize ( const std::string & a )

Definition at line 201 of file GeoPixelServices.cxx.

{
  // For "simple" services use ServiceVolumeMaker to read in the table. The barrel and endcap can also be read in 
  // with ServiceVolumeMaker but there is still use of parameters specified by radiation length which is not handled by ServiceVolumeMaker
  // so we use the old version. 
 
  if (a!="simple") {
    return initializeOld(a);
  }
 
  InDetDD::ServiceVolumeSchema schema;
  std::string label;
  if (a=="simple") {
    schema.setSimpleSchema();
    label = "Svc";
  } else if (a=="barrel") { // Not used yet. initializeOld() is used instead
    schema.setPixelSchema();
    label = "Brl";
  } else if (a=="endcap") { // Not used yet. initializeOld() is used instead
    schema.setPixelSchema();
    label = "EC";
  } else {
    m_gmt_mgr->msg(MSG::ERROR) << "Unrecognized service table type: " << a << endmsg;
    return;
  } 
 
  IRDBRecordset_ptr table = m_gmt_mgr->getPixelServiceRecordset(a);
 
  InDetDD::ServiceVolumeMaker volMaker(label, std::move(table), schema, m_gmt_mgr->athenaComps());
  for (unsigned int i = 0; i < volMaker.numElements(); ++i) {
    m_services.push_back(volMaker.make(i));
  }
}

◆ initializeOld()

void GeoPixelServices::initializeOld ( const std::string & a )

Definition at line 237 of file GeoPixelServices.cxx.

{
  //
  // Loop over the service elements:
  //
  int numServices =  m_gmt_mgr->PixelServiceElements(a);
  for(int ii = 0; ii < numServices; ii++) {
    
 
    // Will return <0 if the element doesn't belong to current inside/outside zone
    // This check is removed as it now checks accroding to the dimensions.
    //if (m_gmt_mgr->PixelServiceFrameNum(a, ii) < 0) continue;
    //
    //
    // Retrieve/calculate the parameters for the volume.
    //
    InDetDD::ServiceVolume param;
    param.setMaterial(m_gmt_mgr->PixelServiceMaterial(a, ii));
    param.setRmin(m_gmt_mgr->PixelServiceRMin(a, ii));
    param.setRmax(m_gmt_mgr->PixelServiceRMax(a, ii));
    param.setZmin(m_gmt_mgr->PixelServiceZMin(a, ii));
    param.setZmax(m_gmt_mgr->PixelServiceZMax(a, ii));
    param.setZsymm(m_gmt_mgr->PixelServiceZsymm(a, ii));
    param.setVolName(m_gmt_mgr->PixelServiceName(a, ii));
 
    double zShift=0.;           // the famous IBL Z shift
    int iShiftIndex = m_gmt_mgr->PixelServiceShift(a, ii);
    // FIXME: The magic number 100 here should be explained...
    // cppcheck-suppress negativeContainerIndex
    if(iShiftIndex>0) zShift=m_layerShift[iShiftIndex-100];
    param.setZShift(zShift);
 
    // Service envelopes
    int envNum=m_gmt_mgr->PixelServiceEnvelopeNum(a, ii);
    param.setEnvelopeNum(envNum);
    int envParent=m_gmt_mgr->PixelServiceParentEnvelopeNum(a, ii);
    param.setParentEnvelope(envParent);
    int volId = m_gmt_mgr->PixelServiceFrameNum(a, ii);
 
    bool needsRotation = false;
 
    // For TUBE there is no need to read the rest 
    std::string shapeType = m_gmt_mgr->PixelServiceShape(a, ii);
    if (!shapeType.empty() && shapeType != "TUBE") {
      double rmin2 = m_gmt_mgr->PixelServiceRMin2(a, ii);
      double rmax2 = m_gmt_mgr->PixelServiceRMax2(a, ii);
      int repeat = m_gmt_mgr->PixelServiceRepeat(a, ii);
      double phiLoc =  m_gmt_mgr->PixelServicePhiLoc(a, ii);
      double phiWidth =  m_gmt_mgr->PixelServiceWidth(a, ii);
      
      // Can be in degree or Gaudi::Units::mm. Usually it is Gaudi::Units::deg expect for BOX, TRAP and ROD shape
      // Geometry manager makes no assumptions about units. So we must interpret here.
      if (shapeType == "BOX" || shapeType == "ROD" || shapeType == "ROD2" || shapeType == "TRAP") {
    phiWidth *= Gaudi::Units::mm;
      } else {
    phiWidth *= Gaudi::Units::degree;
      }
      
      if (phiWidth == 0) phiWidth = 2*Gaudi::Units::pi;
      if (rmin2 <= 0) rmin2 = param.rmin(); 
      if (rmax2 <= 0) rmax2 = param.rmax(); 
      if (repeat == 0) repeat = 1;
      
      if (shapeType == "PGON"  || shapeType == "PGON2" || 
      shapeType == "CONE"  || shapeType == "CONS" || 
      shapeType == "PGON3" || shapeType == "PGON31" || shapeType == "PGON4") {
    if ((rmin2 != param.rmin()) || (rmax2 != param.rmax())) {
      needsRotation = true;
    }
      }
      
      int sides = 0;
      int nCopies = 1;
      if (shapeType == "PGON"  || shapeType == "PGON2" ||
      shapeType == "PGON3" || shapeType == "PGON31" || shapeType == "PGON4") {
    sides = repeat;
      } else {
    nCopies = repeat;
      }
      
      // Force nCopies to 1 for TUBE and CONE 
      if (shapeType.empty() || shapeType == "TUBE" || shapeType == "CONE") {
    nCopies = 1;
      }
      
      param.setShapeType(shapeType);
      param.setRmin2(rmin2);
      param.setRmax2(rmax2);
      param.setPhiLoc(phiLoc);
      param.setPhiWidth(phiWidth);
      param.setSides(sides);
      param.setNCopies(nCopies);
    }
 
    param.setNeedsRotation(needsRotation);
    
    //
    // If zin is 0... (within 10^-5) this is a volume symmetric around
    // the origin
    //
    if(std::abs(param.zmin()) < 0.000001) {
      param.setZmin(-param.zmax());
      param.setZsymm(false);
    }   
 
    std::string label;
    if(envNum>0) label="Env";
    else if (a == "barrel") label = "Brl";
    else if (a == "endcap") label = "EC";
    else if (a == "simple") label = "Svc";
    param.setLabel(label, volId); 
 
    // Fill m_multiVolumes;
    m_services.push_back(new InDetDD::ServiceVolume(param));
    
  }
}

◆ operator=()

GeoPixelServices & GeoPixelServices::operator= ( const GeoPixelServices & )

delete

Member Data Documentation

◆ m_DDmgr

InDetDD::PixelDetectorManager* GeoVPixelFactory::m_DDmgr

protectedinherited

Definition at line 45 of file GeoVPixelFactory.h.

◆ m_epsilon

const double GeoVPixelFactory::m_epsilon

protectedinherited

Definition at line 49 of file GeoVPixelFactory.h.

◆ m_gmt_mgr

PixelGeometryManager* GeoVPixelFactory::m_gmt_mgr

protectedinherited

Definition at line 43 of file GeoVPixelFactory.h.

◆ m_layerShift

std::vector<double> GeoPixelServices::m_layerShift

private

Definition at line 39 of file GeoPixelServices.h.

◆ m_mapAX

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

protectedinherited

Definition at line 48 of file GeoVPixelFactory.h.

◆ m_mapFPV

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

protectedinherited

Definition at line 47 of file GeoVPixelFactory.h.

◆ m_mat_mgr

InDetMaterialManager* GeoVPixelFactory::m_mat_mgr

protectedinherited

Definition at line 44 of file GeoVPixelFactory.h.

◆ m_pixServBuilder

InDetDD::VolumeBuilder* GeoPixelServices::m_pixServBuilder

private

Definition at line 37 of file GeoPixelServices.h.

◆ m_services

std::vector<const InDetDD::ServiceVolume *> GeoPixelServices::m_services

private

Definition at line 36 of file GeoPixelServices.h.

◆ m_servMatBuilder

InDetDD::VolumeBuilder* GeoPixelServices::m_servMatBuilder

private

Definition at line 38 of file GeoPixelServices.h.

◆ m_sqliteReader

GeoModelIO::ReadGeoModel* GeoVPixelFactory::m_sqliteReader

protectedinherited

Definition at line 46 of file GeoVPixelFactory.h.

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

Public Member Functions

Protected Attributes

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ GeoPixelServices() [1/2]

◆ GeoPixelServices() [2/2]

◆ ~GeoPixelServices()

Member Function Documentation

◆ Build()

◆ getBuilder()

◆ getServMatBuilder()

◆ initialize()

◆ initializeOld()

◆ operator=()

Member Data Documentation

◆ m_DDmgr

◆ m_epsilon

◆ m_gmt_mgr

◆ m_layerShift

◆ m_mapAX

◆ m_mapFPV

◆ m_mat_mgr

◆ m_pixServBuilder

◆ m_services

◆ m_servMatBuilder

◆ m_sqliteReader