Muon::TgcHitClusteringObj Struct Reference

#include <TgcHitClustering.h>

Collaboration diagram for Muon::TgcHitClusteringObj:

Public Types
using	HitList = TgcClusterObj3D::HitList

Public Member Functions
	TgcHitClusteringObj (const TgcIdHelper *tgcIdHelp)
bool	cluster (const HitList &col)
bool	cluster (HitList &filteredHits, std::vector< HitList > &finalClusts)
bool	buildClusters3D ()
void	dump () const
const HitList &	bestEtaCluster () const
const HitList &	bestPhiCluster () const

Public Attributes
const TgcIdHelper *	m_tgcIdHelper {nullptr}
std::vector< HitList >	clustersEta {}
std::vector< HitList >	clustersPhi {}
std::vector< TgcClusterObj3D >	clusters3D {}

Detailed Description

Definition at line 47 of file TgcHitClustering.h.

Member Typedef Documentation

HitList

using Muon::TgcHitClusteringObj::HitList = TgcClusterObj3D::HitList

Definition at line 49 of file TgcHitClustering.h.

Constructor & Destructor Documentation

TgcHitClusteringObj()

Muon::TgcHitClusteringObj::TgcHitClusteringObj ( const TgcIdHelper *  tgcIdHelp )

inline

Definition at line 51 of file TgcHitClustering.h.

                                                        :
        m_tgcIdHelper(tgcIdHelp) {}

Member Function Documentation

bestEtaCluster()

const HitList & Muon::TgcHitClusteringObj::bestEtaCluster

(

)

const

Definition at line 12 of file TgcHitClustering.cxx.

                                                             {      
        return clustersEta.size() ? clustersEta.front() : dummy;
    }

bestPhiCluster()

const HitList & Muon::TgcHitClusteringObj::bestPhiCluster

(

)

const

Definition at line 15 of file TgcHitClustering.cxx.

                                                             {      
        return clustersPhi.size() ? clustersPhi.front() : dummy;
    }

buildClusters3D()

bool Muon::TgcHitClusteringObj::buildClusters3D

(

)

Definition at line 65 of file TgcHitClustering.cxx.

                                            {
    
    if( clustersPhi.empty() || clustersEta.empty() ) return false;
    
    const TgcPrepData* etaHit = bestEtaCluster().front(); 
 
    const MuonGM::TgcReadoutElement* detEl = etaHit->detectorElement();
 
    // now loop over eta and phi clusters and form space points
    for(HitList& eit : clustersEta) {      
 
        const TgcPrepData* firstEta = eit.front();
        const TgcPrepData* lastEta = eit.back();
      
        for(HitList& pit : clustersPhi) {
 
          const TgcPrepData* firstPhi = pit.front();
          const TgcPrepData* lastPhi = pit.back();
 
          TgcClusterObj3D cl3D{eit, pit};
          using Edge = TgcClusterObj3D::Edge;
 
          detEl->spacePointPosition( firstPhi->identify(), firstEta->identify(), cl3D.getEdge(Edge::LowEtaLowPhi) );
          if( lastPhi != firstPhi ) { 
            detEl->spacePointPosition( lastPhi->identify(), firstEta->identify(), cl3D.getEdge(Edge::LowEtaHighPhi)); 
          } else {
            cl3D.getEdge(Edge::LowEtaHighPhi) = cl3D.getEdge(Edge::LowEtaLowPhi);
          }
          if( lastEta != firstEta ) {
            detEl->spacePointPosition( firstPhi->identify(),lastEta->identify(),  cl3D.getEdge(Edge::HighEtaLowPhi) );
            if( lastPhi != firstPhi ) { 
              detEl->spacePointPosition( lastPhi->identify(),lastEta->identify(), cl3D.getEdge(Edge::HighEtaHighPhi) ); 
            } else { 
              cl3D.getEdge(Edge::HighEtaHighPhi) = cl3D.getEdge(Edge::HighEtaLowPhi);
            }
          }else{
            cl3D.getEdge(Edge::HighEtaLowPhi) = cl3D.getEdge(Edge::LowEtaLowPhi);
            cl3D.getEdge(Edge::HighEtaHighPhi) = cl3D.getEdge(Edge::LowEtaHighPhi);
          }
          clusters3D.push_back(std::move(cl3D));
      }
    }
    return true;
  }

cluster() [1/2]

bool Muon::TgcHitClusteringObj::cluster

(

const HitList &

col

)

cluster() [2/2]

bool Muon::TgcHitClusteringObj::cluster	(	HitList &	filteredHits,
		std::vector< HitList > &	finalClusts
	)

Create new clusters if the cluster is more separated by more than one channel or if the gas gaps don't match

Definition at line 18 of file TgcHitClustering.cxx.

                                                                       {
      if (filteredHits.empty()) return false;
      finalClusts.reserve(filteredHits.size());
      std::sort(filteredHits.begin(), filteredHits.end(), 
                [this](const TgcPrepData* h1,const TgcPrepData* h2) {
                  const int gap1{m_tgcIdHelper->gasGap(h1->identify())}, gap2{m_tgcIdHelper->gasGap(h2->identify())};
                  if (gap1 != gap2) return gap1 < gap2;
                  return m_tgcIdHelper->channel(h1->identify()) < m_tgcIdHelper->channel(h2->identify());
                });
      
      for (const TgcPrepData* prd : filteredHits) {
          if (finalClusts.empty() || 
              m_tgcIdHelper->channel(finalClusts.back().back()->identify()) +1 != m_tgcIdHelper->channel(prd->identify()) ||
              m_tgcIdHelper->gasGap(finalClusts.back().back()->identify()) != m_tgcIdHelper->gasGap(prd->identify())) {
              finalClusts.emplace_back();
          }  
          finalClusts.back().push_back(prd);
      }
      std::stable_sort(finalClusts.begin(),finalClusts.end(), 
                       []( const HitList& a, const HitList& b){
                            return a.size() < b.size();
                       });
      return true;
  }

dump()

void Muon::TgcHitClusteringObj::dump

(

)

const

Member Data Documentation

clusters3D

std::vector<TgcClusterObj3D> Muon::TgcHitClusteringObj::clusters3D {}

Definition at line 69 of file TgcHitClustering.h.

clustersEta

std::vector<HitList> Muon::TgcHitClusteringObj::clustersEta {}

Definition at line 67 of file TgcHitClustering.h.

clustersPhi

std::vector<HitList> Muon::TgcHitClusteringObj::clustersPhi {}

Definition at line 68 of file TgcHitClustering.h.

m_tgcIdHelper

const TgcIdHelper* Muon::TgcHitClusteringObj::m_tgcIdHelper {nullptr}

Definition at line 66 of file TgcHitClustering.h.

---

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

• TgcHitClustering.h
• TgcHitClustering.cxx