#include <GeoPixelSiCrystal.h>

Inheritance diagram for GeoPixelSiCrystal:

Collaboration diagram for GeoPixelSiCrystal:

Public Member Functions
	GeoPixelSiCrystal (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, bool isBLayer, bool isModule3D=false, bool even_odd_phi_design=false)
virtual GeoVPhysVol *	Build () override
Identifier	getID ()
bool	GetModule3DFlag ()

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
Identifier	m_id
unsigned int	m_nPhiDesigns =1
std::array< const InDetDD::SiDetectorDesign *, 2 >	m_design {nullptr,nullptr}
bool	m_isBLayer = false
bool	m_isModule3D = false

Detailed Description

Definition at line 19 of file GeoPixelSiCrystal.h.

Constructor & Destructor Documentation

◆ GeoPixelSiCrystal()

GeoPixelSiCrystal::GeoPixelSiCrystal	(	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,
		bool	isBLayer,
		bool	isModule3D = false,
		bool	even_odd_phi_design = false )

Definition at line 44 of file GeoPixelSiCrystal.cxx.

  : GeoVPixelFactory (ddmgr, mgr, sqliteReader, std::move(mapFPV), std::move(mapAX))
{
  // 
  //Builds the design for this crystal
  m_isBLayer = isBLayer;
  m_isModule3D = isModule3D;
 
  // Dimensions
  const double thickness = m_gmt_mgr->PixelBoardThickness(m_isModule3D);
  const double length = m_gmt_mgr->PixelBoardLength(m_isModule3D);
  const double width = m_gmt_mgr->PixelBoardWidth(m_isModule3D);
 
  int circuitsPerPhi = m_gmt_mgr->DesignCircuitsPhi(m_isModule3D); // Warning col/row naming opposite to chip
  int circuitsPerEta = m_gmt_mgr->DesignCircuitsEta(m_isModule3D); // Warning col/row naming opposite to chip
  int cellRowPerCirc = m_gmt_mgr->DesignCellRowsPerCircuit(m_isModule3D);
  int cellColPerCirc = m_gmt_mgr->DesignCellColumnsPerCircuit(m_isModule3D);
  int rowsPerCircuit = m_gmt_mgr->DesignDiodeRowsPerCircuit(m_isModule3D);
  int columnsPerCircuit = m_gmt_mgr->DesignDiodeColumnsPerCircuit(m_isModule3D);
  int readoutSide = m_gmt_mgr->DesignReadoutSide(m_isModule3D);
 
  double pitchEtaLongEnd =  m_gmt_mgr->DesignPitchZLongEnd(m_isModule3D);
  double pitchEtaLong =  m_gmt_mgr->DesignPitchZLong(m_isModule3D);
  double pitchPhi = m_gmt_mgr->DesignPitchRP(m_isModule3D);
  double pitchEta = m_gmt_mgr->DesignPitchZ(m_isModule3D);
 
  auto  readoutTechnology = getPixelReadoutTechnology(rowsPerCircuit,columnsPerCircuit );
  auto circuitsPerPhi_corr = circuitsPerPhi;
  auto rowsPerCircuit_corr = rowsPerCircuit;
  if (readoutTechnology==InDetDD::PixelReadoutTechnology::FEI3 && circuitsPerPhi==1 && rowsPerCircuit==328) {
     // FEI3 has 2x8 circuits and 160 + 4 rows per circuit not 1x8 circuits with 328 rows
     // without this correction it is not possible to use the attribute associated to sub-matrices
     // to assign FE numbers to sub-matrices.
     circuitsPerPhi_corr*=2;
     rowsPerCircuit_corr/=2;
  }
 
  // optionally create different module design for even and odd phi modules, because
  // the front-ends are numbered differently
  m_nPhiDesigns=even_odd_phi_design ? 2 : 1;
  for (unsigned int design_i=0; design_i<m_nPhiDesigns ; ++design_i) {
  // for endcap modules the design for even-phi indices uses a "mirrored" front-end numbering scheme.
  // instead of from top left  0..7 to top right, then bottom right to bottom left 8..15 , the front
  // ends are numbered from bottom left to top left i.e. bottom and top are swapped.
  InDetDD::detail::FENumbering fe_numbering = (even_odd_phi_design && design_i==0 ? InDetDD::detail::FENumbering::kMirror : InDetDD::detail::FENumbering::kRegular);
 
  constexpr auto kNDirections = InDetDD::detail::kNDirections;
  constexpr auto kNPixelLocations = InDetDD::detail::kNPixelLocations;
  PixelDiodeTree diode_tree = InDetDD::detail::makePixelDiodeTree(m_gmt_mgr,
                                                    readoutTechnology,
                                                    std::array<int,kNDirections>{circuitsPerPhi_corr,circuitsPerEta},    // [0]=phi/row, [1]=eta/column
                                                    std::array<int,kNDirections>{rowsPerCircuit_corr,columnsPerCircuit}, // [0]=phi/row, [1]=eta/column
                                                    std::array<std::array<double,kNDirections>,kNPixelLocations>{        // regular/central,longEnd/outer,long/inner
                                                       std::array<double,kNDirections>{pitchPhi,pitchEta},
                                                       std::array<double,kNDirections>{0.,pitchEtaLongEnd},
                                                       std::array<double,kNDirections>{0.,pitchEtaLong}},
                                                    fe_numbering);
 
  std::unique_ptr<PixelModuleDesign> p_barrelDesign2 = std::make_unique<PixelModuleDesign>(thickness,
                                 circuitsPerPhi,
                                 circuitsPerEta,
                                 cellColPerCirc,
                                 cellRowPerCirc,
                                 columnsPerCircuit,
                                 rowsPerCircuit,
                                 std::move(diode_tree),
                                 InDetDD::electrons,
                                 readoutSide,
                                 false,              /* 3D */
                                 InDetDD::Undefined, /* detector type */
                                 readoutTechnology);
 
  // Multiple connections (ganged pixels)
  if (m_gmt_mgr->NumberOfEmptyRows() > 0) {
    int minRow = m_gmt_mgr->EmptyRows(0);
    int maxRow = minRow;
    for (int iConnect = 0; iConnect < m_gmt_mgr->NumberOfEmptyRows(); iConnect++){
      minRow = std::min(minRow, m_gmt_mgr->EmptyRows(iConnect));
      minRow = std::min(minRow, m_gmt_mgr->EmptyRowConnections(iConnect));
      maxRow = std::max(maxRow, m_gmt_mgr->EmptyRows(iConnect));
      maxRow = std::max(maxRow, m_gmt_mgr->EmptyRowConnections(iConnect));
    }
  
    std::vector <int> connections(maxRow-minRow+1);
    
    // We fill them all with a one to one correspondence first.
    for (unsigned int iRow = 0; iRow < connections.size(); iRow++){
      connections[iRow] = iRow +  minRow;
    }
    
    // Now make the connections.
    for (int iConnect = 0; iConnect < m_gmt_mgr->NumberOfEmptyRows(); iConnect++){
      connections[m_gmt_mgr->EmptyRows(iConnect)-minRow] = m_gmt_mgr->EmptyRowConnections(iConnect);
    }
   
    p_barrelDesign2->addMultipleRowConnection(minRow, connections);
 
  } else {
    // No empty rows.
    if (!m_gmt_mgr->ibl()) {
      m_gmt_mgr->msg(MSG::WARNING) << "GeoPixelSiCrystal: No ganged pixels" << endmsg;
    } else {
      if (m_gmt_mgr->msgLvl(MSG::DEBUG)) m_gmt_mgr->msg(MSG::DEBUG) << "GeoPixelSiCrystal: No ganged pixels" << endmsg;
    }
  }
  
  // Check that the active area seems reasonable
 
  if ( (m_gmt_mgr->DesignRPActiveArea(m_isModule3D) > width) ||
       (m_gmt_mgr->DesignZActiveArea(m_isModule3D) >  length) || 
       (width - m_gmt_mgr->DesignRPActiveArea(m_isModule3D) > 4 * Gaudi::Units::mm) || 
       (length - m_gmt_mgr->DesignZActiveArea(m_isModule3D) > 4 * Gaudi::Units::mm) ) { 
    m_gmt_mgr->msg(MSG::WARNING) << "GeoPixelSiCrystal: Active area not consistent with sensor size. Sensor: " 
                   << width/Gaudi::Units::mm << " x " << length/Gaudi::Units::mm << ", Active: " 
                   << m_gmt_mgr->DesignRPActiveArea(m_isModule3D)/Gaudi::Units::mm << " x " << m_gmt_mgr->DesignZActiveArea(m_isModule3D)/Gaudi::Units::mm 
                   << endmsg;
  } else {
    if (m_gmt_mgr->msgLvl(MSG::DEBUG)) m_gmt_mgr->msg(MSG::DEBUG) 
      << "GeoPixelSiCrystal: Sensor: "  
      << width/Gaudi::Units::mm << " x " << length/Gaudi::Units::mm << ", Active: " 
      << m_gmt_mgr->DesignRPActiveArea(m_isModule3D)/Gaudi::Units::mm << " x " << m_gmt_mgr->DesignZActiveArea(m_isModule3D)/Gaudi::Units::mm 
      << endmsg;             
  }
 
 
  m_design.at(design_i) = m_DDmgr->addDesign(std::move(p_barrelDesign2));
  }
  
  
}

Member Function Documentation

◆ Build()

GeoVPhysVol * GeoPixelSiCrystal::Build ( )

overridevirtual

Implements GeoVPixelFactory.

Definition at line 179 of file GeoPixelSiCrystal.cxx.

                                      {
 
  GeoFullPhysVol* siPhys{nullptr};
  int brl_ec= m_gmt_mgr->isEndcap() ? 2*m_gmt_mgr->GetSide() : 0;
 
  if(m_sqliteReader) {
    std::string siName ="SiPhys_" + std::to_string(brl_ec) +"_"+ std::to_string(m_gmt_mgr->GetLD()) +"_"+ std::to_string(m_gmt_mgr->Phi()) +"_"+ std::to_string(m_gmt_mgr->Eta());
    siPhys = (*m_mapFPV)[siName];
  }
  else {
    //(sar) code moved from c'tor..
    const double thickness = m_gmt_mgr->PixelBoardThickness(m_isModule3D);
    const double length = m_gmt_mgr->PixelBoardLength(m_isModule3D);
    const double width = m_gmt_mgr->PixelBoardWidth(m_isModule3D);
    //
    std::string matName = m_gmt_mgr->getMaterialName("Sensor");
    const GeoMaterial* siMat = m_mat_mgr->getMaterial(matName);
    const GeoBox* siBox = new GeoBox(thickness*0.5,width*0.5,length*0.5);
    std::string logname{"siLog"};
    // There is not a strong need to give the blayer a different name but leave it for now. 
    if(m_isBLayer) logname = "siBLayLog";
    auto *logVolume = new GeoLogVol(logname,siBox,siMat);
    //(sar) ...to here
 
 
    siPhys = new GeoFullPhysVol(logVolume);
  }
  // Build the Identifier for the silicon:
  //
  const PixelID * idHelper = m_gmt_mgr->getIdHelper();
  m_id = idHelper->wafer_id(brl_ec,m_gmt_mgr->GetLD(),m_gmt_mgr->Phi(),m_gmt_mgr->Eta());
  SiDetectorElement * element = new SiDetectorElement(m_id, m_design[ m_gmt_mgr->Phi() % m_nPhiDesigns] , siPhys, m_gmt_mgr->commonItems());
  
  // add the element to the manager
  m_DDmgr->addDetectorElement(element);
  return siPhys;
}

◆ getID()

Identifier GeoPixelSiCrystal::getID ( )

inline

Definition at line 28 of file GeoPixelSiCrystal.h.

28{return m_id;}

◆ GetModule3DFlag()

bool GeoPixelSiCrystal::GetModule3DFlag ( )

inline

Definition at line 30 of file GeoPixelSiCrystal.h.

30{ return m_isModule3D; };

Member Data Documentation

◆ m_DDmgr

InDetDD::PixelDetectorManager* GeoVPixelFactory::m_DDmgr

protectedinherited

Definition at line 45 of file GeoVPixelFactory.h.

◆ m_design

std::array<const InDetDD::SiDetectorDesign*,2> GeoPixelSiCrystal::m_design {nullptr,nullptr}

private

Definition at line 39 of file GeoPixelSiCrystal.h.

39{nullptr,nullptr};

◆ 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_id

Identifier GeoPixelSiCrystal::m_id

private

Definition at line 33 of file GeoPixelSiCrystal.h.

◆ m_isBLayer

bool GeoPixelSiCrystal::m_isBLayer = false

private

Definition at line 40 of file GeoPixelSiCrystal.h.

◆ m_isModule3D

bool GeoPixelSiCrystal::m_isModule3D = false

private

Definition at line 41 of file GeoPixelSiCrystal.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_nPhiDesigns

unsigned int GeoPixelSiCrystal::m_nPhiDesigns =1

private

Definition at line 38 of file GeoPixelSiCrystal.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

◆ GeoPixelSiCrystal()

Member Function Documentation

◆ Build()

◆ getID()

◆ GetModule3DFlag()

Member Data Documentation

◆ m_DDmgr

◆ m_design

◆ m_epsilon

◆ m_gmt_mgr

◆ m_id

◆ m_isBLayer

◆ m_isModule3D

◆ m_mapAX

◆ m_mapFPV

◆ m_mat_mgr

◆ m_nPhiDesigns

◆ m_sqliteReader