InDet::TotPixelClusterSplitter Class Reference

#include <TotPixelClusterSplitter.h>

Inheritance diagram for InDet::TotPixelClusterSplitter:

Collaboration diagram for InDet::TotPixelClusterSplitter:

Public Member Functions
	TotPixelClusterSplitter (const std::string &type, const std::string &name, const IInterface *parent)
	Constructor. More...
	~TotPixelClusterSplitter ()=default
	Destructor. More...
virtual StatusCode	initialize () override
	AthAlgTool interface methods. More...
virtual StatusCode	finalize () override
virtual std::vector< InDet::PixelClusterParts >	splitCluster (const InDet::PixelCluster &OrigCluster) const override
	take one, give zero or many More...
virtual std::vector< InDet::PixelClusterParts >	splitCluster (const InDet::PixelCluster &OrigCluster, const InDet::PixelClusterSplitProb &spo) const override
	take one, give zero or many - with split probability object More...
void	setMinPixels (unsigned int minPix)
	Set the lower and upper bounds for the number of pixels in clusters to be considered. More...
void	setMaxPixels (unsigned int maxPix)

Private Types
enum	SplitType { PhiSplit = 0, EtaSplit = 1, NoSplit = 2 }

Static Private Member Functions
static int	pixelType (const int PhiIdx, const int EtaIdx)
	Determine a pixel's type. More...

Private Attributes
ServiceHandle< InDetDD::IPixelReadoutManager >	m_pixelReadout {this, "PixelReadoutManager", "PixelReadoutManager", "Pixel readout manager" }
SG::ReadCondHandleKey< PixelChargeCalibCondData >	m_chargeDataKey {this, "PixelChargeCalibCondData", "PixelChargeCalibCondData", "Pixel charge calibration data"}
unsigned int	m_minPixels
	Minimum number of pixels in cluster to consider splitting. More...
unsigned int	m_maxPixels
	Maximum size of cluster for which splitting will be attempted. More...

Detailed Description

Author

Kurt Barry, Tim Nelsen, Andreas Salzburger

Definition at line 33 of file TotPixelClusterSplitter.h.

Member Enumeration Documentation

SplitType

enum InDet::TotPixelClusterSplitter::SplitType

private

Enumerator
PhiSplit
EtaSplit
NoSplit

Definition at line 79 of file TotPixelClusterSplitter.h.

{ PhiSplit = 0, EtaSplit = 1, NoSplit = 2 };

Constructor & Destructor Documentation

TotPixelClusterSplitter()

InDet::TotPixelClusterSplitter::TotPixelClusterSplitter	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Constructor.

Definition at line 17 of file TotPixelClusterSplitter.cxx.

                                                                                                                               :
  base_class(type, name, parent),
  m_minPixels(3),
  m_maxPixels(25)
{
}

~TotPixelClusterSplitter()

InDet::TotPixelClusterSplitter::~TotPixelClusterSplitter ( )

default

Destructor.

Member Function Documentation

finalize()

StatusCode InDet::TotPixelClusterSplitter::finalize ( )

overridevirtual

Definition at line 32 of file TotPixelClusterSplitter.cxx.

                                                  {
  return StatusCode::SUCCESS;
}

initialize()

StatusCode InDet::TotPixelClusterSplitter::initialize ( )

overridevirtual

AthAlgTool interface methods.

Definition at line 24 of file TotPixelClusterSplitter.cxx.

                                                    {
 
  ATH_CHECK(m_pixelReadout.retrieve());
  ATH_CHECK(m_chargeDataKey.initialize());
 
  return StatusCode::SUCCESS;
}

pixelType()

int InDet::TotPixelClusterSplitter::pixelType ( const int  PhiIdx, const int  EtaIdx  )

staticprivate

Determine a pixel's type.

Even numbers have normal eta-pitch, odd numbers have long eta pitch.

Arguments: phiIdx : sensor phi index, in the range [0, 327] etaIdx : sensor eta index, in the range [0, 17]

Return Values: -1 : one or both indicies are out-of-range 0 : short normal pixel 1 : long normal pixel 2 : short ganged pixel 3 : long ganged pixel 4 : short inter-ganged pixel 5 : long inter-ganged pixel

Definition at line 234 of file TotPixelClusterSplitter.cxx.

{
  if (EtaIdx%18 != 0 && EtaIdx%18 != 17)
  {
    if (PhiIdx==152 || PhiIdx==154 || PhiIdx==156 || PhiIdx==158 ||
        PhiIdx==169 || PhiIdx==171 || PhiIdx==173 || PhiIdx==175)
      return 4; // short inter-ganged pixel
    else if (PhiIdx==153 || PhiIdx==155 || PhiIdx==157 || PhiIdx==159 ||
             (PhiIdx > 159 && PhiIdx < 168) ||
             PhiIdx==168 || PhiIdx==170 || PhiIdx==172 || PhiIdx==174)
      return 2; // short ganged pixel
    else
      return 0; // short normal pixel
  }
  else
  {
    if (PhiIdx==152 || PhiIdx==154 || PhiIdx==156 || PhiIdx==158 ||
        PhiIdx==169 || PhiIdx==171 || PhiIdx==173 || PhiIdx==175)
      return 5; // long inter-ganged pixel
    else if (PhiIdx==153 || PhiIdx==155 || PhiIdx==157 || PhiIdx==159 ||
             (PhiIdx > 159 && PhiIdx < 168) ||
             PhiIdx==168 || PhiIdx==170 || PhiIdx==172 || PhiIdx==174)
      return 3; // long ganged pixel
    else
      return 1; // long normal pixel
  }
}

setMaxPixels()

void InDet::TotPixelClusterSplitter::setMaxPixels ( unsigned int  maxPix )

inline

Definition at line 101 of file TotPixelClusterSplitter.h.

  {
    m_maxPixels = maxPix;
    return;
  }

setMinPixels()

void InDet::TotPixelClusterSplitter::setMinPixels ( unsigned int  minPix )

inline

Set the lower and upper bounds for the number of pixels in clusters to be considered.

Definition at line 95 of file TotPixelClusterSplitter.h.

  {
    m_minPixels = minPix;
    return;
  }

splitCluster() [1/2]

std::vector< InDet::PixelClusterParts > InDet::TotPixelClusterSplitter::splitCluster ( const InDet::PixelCluster &  OrigCluster ) const

overridevirtual

take one, give zero or many

Definition at line 37 of file TotPixelClusterSplitter.cxx.

                                                                                       {
 
  std::vector<InDet::PixelClusterParts> Parts;
 
  const std::vector<Identifier> & Rdos = OrigCluster.rdoList();
  const unsigned int NumPixels = static_cast<unsigned int>(Rdos.size());
  const InDetDD::SiDetectorElement * Element = OrigCluster.detectorElement();
 
  // Check cluster-size bounds.
  if (NumPixels < m_minPixels || NumPixels > m_maxPixels) return Parts;
 
  std::vector<InDetDD::SiCellId> CellIds(NumPixels); //NumPixels entries of default constructed SiCellId
  
  auto getCellId=[Element](const Identifier & i){ return Element->cellIdFromIdentifier(i);};
  std::transform(Rdos.begin(),Rdos.end(), CellIds.begin(), getCellId);
 
  SG::ReadCondHandle<PixelChargeCalibCondData> calibDataHandle(m_chargeDataKey);
  const PixelChargeCalibCondData *calibData = *calibDataHandle;
 
  // Detect special pixels and exclude them if necessary.
  // Veto removed on 28-7-2011 by K. Barry - residual-level 
  // studies show RMS improvement when ganged-pixel containing 
  // clusters are split.
  
 
  // Determine maximum and minimum phi and eta indices.
  auto comparePhi=[](const InDetDD::SiCellId & a, const InDetDD::SiCellId & b){ return (a.phiIndex() < b.phiIndex());};
  auto compareEta=[](const InDetDD::SiCellId & a, const InDetDD::SiCellId & b){ return (a.etaIndex() < b.etaIndex());};
  const auto [pMinPhi,pMaxPhi] = std::minmax_element(CellIds.begin(), CellIds.end(), comparePhi);
  const auto [pMinEta,pMaxEta] = std::minmax_element(CellIds.begin(), CellIds.end(), compareEta);
  const auto & MinPhiIdx = pMinPhi->phiIndex();
  const auto & MaxPhiIdx = pMaxPhi->phiIndex();
  const auto & MinEtaIdx = pMinEta->etaIndex();
  const auto & MaxEtaIdx = pMaxEta->etaIndex();
  const int PhiWidth = MaxPhiIdx - MinPhiIdx + 1;
  const int EtaWidth = MaxEtaIdx - MinEtaIdx + 1;
  
  // In the future, might want to add a check here to see if the charge list is empty...
  const std::vector<float> & Charges = OrigCluster.chargeList();
 
  std::vector<float> PhiChargeHist(PhiWidth, 0.);
  unsigned int NumPhiMinima = 0;
  std::vector<unsigned int> PhiMinimaLoc(PhiWidth/2);
  std::vector<float> PhiMinimaQuality(PhiWidth/2);
  if (PhiWidth > 2)
  {
    // Fill up the phi-charge histogram. The values were initialized to zero by the vector constructor.
    for (unsigned int i = 0; i < NumPixels; i++)
      PhiChargeHist[CellIds[i].phiIndex() - MinPhiIdx] += Charges[i];
 
    // Minima search
    for (int i = 1; i < PhiWidth - 1; i++)
    {
      if (PhiChargeHist[i-1] <= PhiChargeHist[i] || PhiChargeHist[i] >= PhiChargeHist[i+1]) continue;
 
      PhiMinimaLoc[NumPhiMinima] = static_cast<unsigned int>(i);
      PhiMinimaQuality[NumPhiMinima++] = (PhiChargeHist[i+1] + PhiChargeHist[i-1] - 2.*PhiChargeHist[i]) /
                                      (PhiChargeHist[i+1] + PhiChargeHist[i-1]);
    }
  }
    
  std::vector<float> EtaChargeHist(EtaWidth, 0.);
  unsigned int NumEtaMinima = 0;
  std::vector<unsigned int> EtaMinimaLoc(EtaWidth/2);
  std::vector<float> EtaMinimaQuality(EtaWidth/2);
  const float LongPixScale = 2./3.;
  if (EtaWidth > 2)
  {
    float Scale;
    // Fill up the eta-charge histogram. The values were initialized to zero by the vector constructor.
    for (unsigned int i = 0; i < NumPixels; i++)
    {
      const int PixType = pixelType(CellIds[i].phiIndex(), CellIds[i].etaIndex());
      if ( (PixType % 2) == 1 )
        Scale = LongPixScale;
      else
        Scale = 1.;
      EtaChargeHist[CellIds[i].etaIndex() - MinEtaIdx] += Scale*Charges[i];
    }
 
    // Minima search
    for (int i = 1; i < EtaWidth - 1; i++)
    {
      if (EtaChargeHist[i-1] <= EtaChargeHist[i] || EtaChargeHist[i] >= EtaChargeHist[i+1]) continue;
 
      EtaMinimaLoc[NumEtaMinima] = static_cast<unsigned int>(i);
      EtaMinimaQuality[NumEtaMinima++] = (EtaChargeHist[i+1] + EtaChargeHist[i-1] - 2.*EtaChargeHist[i]) /
                                         (EtaChargeHist[i+1] + EtaChargeHist[i-1]);
    }
  }
 
  // Select the "best" minimum for splitting
  unsigned int bestMin = 0;
  float bestQual = -1.;
  SplitType Type = NoSplit;
  for (unsigned int i = 0; i < NumPhiMinima; i++)
  {
    if (PhiMinimaQuality[i] <= bestQual) continue;
    Type = PhiSplit;
    bestMin = PhiMinimaLoc[i];
    bestQual = PhiMinimaQuality[i];
  }
  for (unsigned int i = 0; i < NumEtaMinima; i++)
  {
    if (EtaMinimaQuality[i] <= bestQual) continue;
    Type = EtaSplit;
    bestMin = EtaMinimaLoc[i];
    bestQual = EtaMinimaQuality[i];
  }
  const int SplitIndex = static_cast<int>(bestMin);
 
  // Now it's time to do the split!
  if (Type == NoSplit)
  {
    return Parts;
  }
  if (Type != PhiSplit && Type != EtaSplit)
  {
    // Should never happen...major bug someplace if it does
    ATH_MSG_ERROR("<InDet::TotPixelClusterSplitter::splitCluster> : Unrecognized SplitType!");
    return Parts;
  }
 
  int Idx;
  int (InDetDD::SiCellId:: * const IndexFunc) () const = (Type == PhiSplit) ? &InDetDD::SiCellId::phiIndex
                                                                            : &InDetDD::SiCellId::etaIndex;
  const int LowIdx = (Type == PhiSplit) ? MinPhiIdx : MinEtaIdx;
 
  // Variables for PixelClusterParts
  std::vector<Identifier> SplitRdos[2];
  std::vector<int> Totgroups[2];
  std::vector<int> Lvl1groups[2];
 
  const std::vector<int>& OrigTots = OrigCluster.totList();
  const int Lvl1a = OrigCluster.LVL1A();
  
  const AtlasDetectorID* aid = Element->getIdHelper();
  if (aid==nullptr)
  {
    ATH_MSG_ERROR("Could not get ATLASDetectorID");
    return Parts;
  }
  
  if (aid->helper() != AtlasDetectorID::HelperType::Pixel){
    ATH_MSG_ERROR("Could not get PixelID pointer");
    return Parts;
  } 
  const PixelID* pixelIDp=static_cast<const PixelID*>(aid);
  const PixelID& pixelID = *pixelIDp;
 
  
  for (unsigned int i = 0; i < NumPixels; i++)
  {
    Idx = (CellIds[i].*IndexFunc)() - LowIdx;
    if (Idx < SplitIndex)
    {
      SplitRdos[0].push_back(Rdos[i]);
      Totgroups[0].push_back(OrigTots[i]);
      Lvl1groups[0].push_back(Lvl1a);
    }
    else if (Idx > SplitIndex)
    {
      SplitRdos[1].push_back(Rdos[i]);
      Totgroups[1].push_back(OrigTots[i]);
      Lvl1groups[1].push_back(Lvl1a);
    }
    else // only == remains
    {
      for (int j = 0; j < 2; j++)
      {
 
        Identifier pixid = Rdos[i];
        Identifier moduleID = pixelID.wafer_id(pixid);
        IdentifierHash moduleHash = pixelID.wafer_hash(moduleID); // wafer hash
        unsigned int FE = m_pixelReadout->getFE(pixid, moduleID);
        InDetDD::PixelDiodeType type = m_pixelReadout->getDiodeType(pixid);
 
        SplitRdos[j].push_back(Rdos[i]);
        Totgroups[j].push_back(calibData->getToT(type, moduleHash, FE, Charges[i] / 2.0));
        Lvl1groups[j].push_back(Lvl1a);
      }
    }
  }
 
  Parts.emplace_back(SplitRdos[0], Totgroups[0], Lvl1groups[0]);
  Parts.emplace_back(SplitRdos[1], Totgroups[1], Lvl1groups[1]);
    
  return Parts;
}

splitCluster() [2/2]

std::vector< InDet::PixelClusterParts > InDet::TotPixelClusterSplitter::splitCluster ( const InDet::PixelCluster &  OrigCluster, const InDet::PixelClusterSplitProb &  spo  ) const

overridevirtual

take one, give zero or many - with split probability object

Definition at line 227 of file TotPixelClusterSplitter.cxx.

{
  return splitCluster(OrigCluster);
}

Member Data Documentation

m_chargeDataKey

SG::ReadCondHandleKey<PixelChargeCalibCondData> InDet::TotPixelClusterSplitter::m_chargeDataKey {this, "PixelChargeCalibCondData", "PixelChargeCalibCondData", "Pixel charge calibration data"}

private

Definition at line 84 of file TotPixelClusterSplitter.h.

m_maxPixels

unsigned int InDet::TotPixelClusterSplitter::m_maxPixels

private

Maximum size of cluster for which splitting will be attempted.

Definition at line 91 of file TotPixelClusterSplitter.h.

m_minPixels

unsigned int InDet::TotPixelClusterSplitter::m_minPixels

private

Minimum number of pixels in cluster to consider splitting.

Definition at line 88 of file TotPixelClusterSplitter.h.

m_pixelReadout

ServiceHandle<InDetDD::IPixelReadoutManager> InDet::TotPixelClusterSplitter::m_pixelReadout {this, "PixelReadoutManager", "PixelReadoutManager", "Pixel readout manager" }

private

Definition at line 81 of file TotPixelClusterSplitter.h.

---

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

• TotPixelClusterSplitter.h
• TotPixelClusterSplitter.cxx