mkfit::LayerOfHits Class Reference

#include <HitStructures.h>

Classes
struct	HitInfo

Public Member Functions
void	beginRegistrationOfHits (const HitVec &hitv)
void	endRegistrationOfHits (bool build_original_to_internal_map)
int	getHitIndexFromOriginal (int i) const
int	getOriginalHitIndex (int i) const
float	hit_phi (int i) const
float	hit_q (int i) const
const Hit *	hitArray () const
bool	is_barrel () const
bool	is_endcap () const
bool	is_pix_lyr () const
bool	is_pixb_lyr () const
bool	is_pixe_lyr () const
bool	is_stereo () const
bool	is_tec_lyr () const
bool	is_tib_lyr () const
bool	is_tid_lyr () const
bool	is_tob_lyr () const
bool	is_within_r_limits (float r) const
WSR_Result	is_within_r_sensitive_region (float r, float dr) const
bool	is_within_z_limits (float z) const
WSR_Result	is_within_z_sensitive_region (float z, float dz) const
int	layer_id () const
const LayerInfo *	layer_info () const
	LayerOfHits ()=default
int	nHits () const
bool	phi_bin_dead (int qi, int pi) const
const vecvecPhiBinDead_t &	phi_bin_deads () const
PhiBinInfo_t	phi_bin_info (int qi, int pi) const
const vecvecPhiBinInfo_t &	phi_bin_infos () const
int	phiBin (float phi) const
int	phiBinChecked (float phi) const
int	phiBinFine (float phi) const
int	phiMaskApply (int in) const
void	printBins ()
int	qBin (float q) const
int	qBinChecked (float q) const
const Hit &	refHit (int i) const
void	registerHit (int idx)
void	reset ()
void	setupLayer (const LayerInfo &li)
void	suckInDeads (const DeadVec &deadv)
void	suckInHits (const HitVec &hitv)
const vecPhiBinInfo_t &	vecPhiBinInfo (float q) const
	~LayerOfHits ()

Private Member Functions
void	empty_phi_bins (int q_bin, int phi_bin_1, int phi_bin_2, uint16_t hit_count)
void	empty_phi_bins_dead (int q_bin, int phi_bin_1, int phi_bin_2)
void	empty_q_bins (int q_bin_1, int q_bin_2, uint16_t hit_count)
void	empty_q_bins_dead (int q_bin_1, int q_bin_2)
void	setup_bins (float qmin, float qmax, float dq)

Private Attributes
const HitVec *	m_ext_hits
std::vector< int >	m_ext_idcs
float	m_fq
std::vector< HitInfo >	m_hit_infos
std::vector< float >	m_hit_phis
std::vector< float >	m_hit_qs
unsigned int *	m_hit_ranks = nullptr
bool	m_is_barrel
const LayerInfo *	m_layer_info = nullptr
int	m_max_ext_idx
int	m_min_ext_idx
int	m_n_hits = 0
int	m_nq = 0
vecvecPhiBinDead_t	m_phi_bin_deads
vecvecPhiBinInfo_t	m_phi_bin_infos
float	m_qmax
float	m_qmin
std::vector< uint32_t >	m_qphifines

Static Private Attributes
static constexpr float	m_fphi = Config::m_nphi / Const::TwoPI
static constexpr float	m_fphi_fine = 1024 / Const::TwoPI
static constexpr int	m_phi_bits = 8
static constexpr int	m_phi_bits_fine = 10
static constexpr int	m_phi_bits_shift = m_phi_bits_fine - m_phi_bits
static constexpr int	m_phi_fine_xmask = ~((1 << m_phi_bits_shift) - 1)
static constexpr int	m_phi_mask = 0xff
static constexpr int	m_phi_mask_fine = 0x3ff

Detailed Description

Definition at line 45 of file HitStructures.h.

Constructor & Destructor Documentation

mkfit::LayerOfHits::LayerOfHits ( )

default

mkfit::LayerOfHits::~LayerOfHits

(

)

Definition at line 11 of file HitStructures.cc.

References m_hit_ranks.

                            {
#ifdef COPY_SORTED_HITS
    free_hits();
#endif
    operator delete[](m_hit_ranks);
  }

Member Function Documentation

void mkfit::LayerOfHits::beginRegistrationOfHits

(

const HitVec &

hitv

)

Definition at line 166 of file HitStructures.cc.

References cms::cuda::assert(), m_ext_hits, m_ext_idcs, m_hit_infos, m_max_ext_idx, m_min_ext_idx, m_n_hits, m_nq, m_qphifines, SiStripPI::max, and SiStripPI::min.

                                                              {
    assert(m_nq > 0 && "setupLayer() was not called.");

    m_ext_hits = &hitv;

    m_n_hits = 0;
    m_hit_infos.clear();
    m_qphifines.clear();
    m_ext_idcs.clear();
    m_min_ext_idx = std::numeric_limits<int>::max();
    m_max_ext_idx = std::numeric_limits<int>::min();
  }

void mkfit::LayerOfHits::empty_phi_bins ( int  q_bin, int  phi_bin_1, int  phi_bin_2, uint16_t  hit_count  )

inlineprivate

Definition at line 156 of file HitStructures.h.

References m_phi_bin_infos.

Referenced by empty_q_bins().

                                                                                     {
      for (int pb = phi_bin_1; pb < phi_bin_2; ++pb) {
        m_phi_bin_infos[q_bin][pb] = {hit_count, hit_count};
      }
    }

void mkfit::LayerOfHits::empty_phi_bins_dead ( int  q_bin, int  phi_bin_1, int  phi_bin_2  )

inlineprivate

Definition at line 168 of file HitStructures.h.

References m_phi_bin_deads.

Referenced by empty_q_bins_dead().

                                                                      {
      for (int pb = phi_bin_1; pb < phi_bin_2; ++pb) {
        m_phi_bin_deads[q_bin][pb] = false;
      }
    }

void mkfit::LayerOfHits::empty_q_bins ( int  q_bin_1, int  q_bin_2, uint16_t  hit_count  )

inlineprivate

Definition at line 162 of file HitStructures.h.

References empty_phi_bins(), and mkfit::Config::m_nphi.

Referenced by endRegistrationOfHits(), and suckInHits().

                                                                    {
      for (int qb = q_bin_1; qb < q_bin_2; ++qb) {
        empty_phi_bins(qb, 0, Config::m_nphi, hit_count);
      }
    }

void mkfit::LayerOfHits::empty_q_bins_dead ( int  q_bin_1, int  q_bin_2  )

inlineprivate

Definition at line 174 of file HitStructures.h.

References empty_phi_bins_dead(), and mkfit::Config::m_nphi.

Referenced by suckInDeads().

                                                     {
      for (int qb = q_bin_1; qb < q_bin_2; ++qb) {
        empty_phi_bins_dead(qb, 0, Config::m_nphi);
      }
    }

void mkfit::LayerOfHits::endRegistrationOfHits

(

bool

build_original_to_internal_map

)

Definition at line 195 of file HitStructures.cc.

References empty_q_bins(), mps_fire::i, dqmiolumiharvest::j, isotrackApplyRegressor::k, m_ext_idcs, m_hit_infos, m_hit_phis, m_hit_qs, m_hit_ranks, m_max_ext_idx, m_min_ext_idx, m_n_hits, m_nq, m_phi_bin_infos, m_phi_bits_shift, m_phi_fine_xmask, m_phi_mask_fine, m_qphifines, gpuVertexFinder::printf(), RADIX_UNSIGNED, RadixSort::RelinquishRanks(), RadixSort::Sort(), and mkfit::Config::usePhiQArrays.

                                                                             {
    m_n_hits = m_ext_idcs.size();
    if (m_n_hits == 0)
      return;

    // radix
    operator delete[](m_hit_ranks);
    {
      RadixSort sort;
      sort.Sort(&m_qphifines[0], m_n_hits, RADIX_UNSIGNED);
      m_hit_ranks = sort.RelinquishRanks();
    }

    // copy q/phi

#ifdef COPY_SORTED_HITS
    if (m_capacity < m_n_hits) {
      free_hits();
      alloc_hits(m_n_hits);
    }
#endif

    if (Config::usePhiQArrays) {
      m_hit_phis.resize(m_n_hits);
      m_hit_qs.resize(m_n_hits);
    }

    int curr_qphi = -1;
    empty_q_bins(0, m_nq, 0);

    for (int i = 0; i < m_n_hits; ++i) {
      int j = m_hit_ranks[i];  // index in intermediate
      int k = m_ext_idcs[j];   // index in external hit_vec

#ifdef COPY_SORTED_HITS
      memcpy(&m_hits[i], &hitv[k], sizeof(Hit));
#endif

      if (Config::usePhiQArrays) {
        m_hit_phis[i] = m_hit_infos[j].phi;
        m_hit_qs[i] = m_hit_infos[j].q;
      }

      // Combined q-phi bin with fine part masked off
      const int jqphi = m_qphifines[j] & m_phi_fine_xmask;

      const int phi_bin = (jqphi & m_phi_mask_fine) >> m_phi_bits_shift;
      const int q_bin = jqphi >> 16;

      // Fill the bin info
      if (jqphi != curr_qphi) {
        m_phi_bin_infos[q_bin][phi_bin] = {i, i};
        curr_qphi = jqphi;
      }

      m_phi_bin_infos[q_bin][phi_bin].second++;

      // m_hit_ranks[i] will never be used again - use it to point to external/original index.
      m_hit_ranks[i] = k;
    }

    if (build_original_to_internal_map) {
      if (m_max_ext_idx - m_min_ext_idx + 1 > 8 * m_n_hits) {
        // If this happens we might:
        // a) Use external indices for everything. -- *** We are now. ***
        // b) Build these maps for seeding layers only.
        // c) Have a flag in hit-on-track that tells us if the hit index has been remapped,
        //    essentially, if it is a seed hit. This might be smart anyway.
        //    One could use index < -256 or something similar.

        printf(
            "LayerOfHits::endRegistrationOfHits() original_to_internal index map vector is largish: m_n_hits=%d, "
            "map_vector_size=%d\n",
            m_n_hits,
            m_max_ext_idx - m_min_ext_idx + 1);
      }

      m_ext_idcs.resize(m_max_ext_idx - m_min_ext_idx + 1);
      for (int i = 0; i < m_n_hits; ++i) {
        m_ext_idcs[m_hit_ranks[i] - m_min_ext_idx] = i;
      }
    }

    // We can release m_hit_infos and m_qphifines -- and realloc on next BeginInput.
    // m_qphifines could still be used as pre-selection in selectHitIndices().
  }

int mkfit::LayerOfHits::getHitIndexFromOriginal ( int  i ) const

inline

Definition at line 95 of file HitStructures.h.

References m_ext_idcs, and m_min_ext_idx.

{ return m_ext_idcs[i - m_min_ext_idx]; }

int mkfit::LayerOfHits::getOriginalHitIndex ( int  i ) const

inline

Definition at line 97 of file HitStructures.h.

References mps_fire::i, and m_hit_ranks.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return m_hit_ranks[i]; }

float mkfit::LayerOfHits::hit_phi ( int  i ) const

inline

Definition at line 92 of file HitStructures.h.

References mps_fire::i, and m_hit_phis.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return m_hit_phis[i]; }

float mkfit::LayerOfHits::hit_q ( int  i ) const

inline

Definition at line 91 of file HitStructures.h.

References mps_fire::i, and m_hit_qs.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return m_hit_qs[i]; }

const Hit* mkfit::LayerOfHits::hitArray ( ) const

inline

Definition at line 104 of file HitStructures.h.

References m_ext_hits.

Referenced by mkfit::MkFinder::addBestHit(), mkfit::MkFinder::findCandidates(), and mkfit::MkFinder::findCandidatesCloneEngine().

{ return m_ext_hits->data(); }

bool mkfit::LayerOfHits::is_barrel ( ) const

inline

Definition at line 118 of file HitStructures.h.

References m_is_barrel.

Referenced by mkfit::MkFinder::bkFitFitTracks(), mkfit::MkFinder::findCandidates(), mkfit::MkFinder::findCandidatesCloneEngine(), printBins(), and mkfit::MkFinder::selectHitIndices().

{ return m_is_barrel; }

bool mkfit::LayerOfHits::is_endcap ( ) const

inline

Definition at line 119 of file HitStructures.h.

References m_is_barrel.

{ return !m_is_barrel; }

bool mkfit::LayerOfHits::is_pix_lyr ( ) const

inline

Definition at line 137 of file HitStructures.h.

References mkfit::LayerInfo::is_pix_lyr(), and m_layer_info.

Referenced by mkfit::MkFinder::findCandidates(), and mkfit::MkFinder::findCandidatesCloneEngine().

{ return m_layer_info->is_pix_lyr(); }

bool mkfit::LayerOfHits::is_pixb_lyr ( ) const

inline

Definition at line 135 of file HitStructures.h.

References mkfit::LayerInfo::is_pixb_lyr(), and m_layer_info.

{ return m_layer_info->is_pixb_lyr(); }

bool mkfit::LayerOfHits::is_pixe_lyr ( ) const

inline

Definition at line 136 of file HitStructures.h.

References mkfit::LayerInfo::is_pixe_lyr(), and m_layer_info.

{ return m_layer_info->is_pixe_lyr(); }

bool mkfit::LayerOfHits::is_stereo ( ) const

inline

Definition at line 132 of file HitStructures.h.

References mkfit::LayerInfo::is_stereo(), and m_layer_info.

Referenced by mkfit::MkFinder::bkFitFitTracks().

{ return m_layer_info->is_stereo(); }

bool mkfit::LayerOfHits::is_tec_lyr ( ) const

inline

Definition at line 141 of file HitStructures.h.

References mkfit::LayerInfo::is_tec_lyr(), and m_layer_info.

{ return m_layer_info->is_tec_lyr(); }

bool mkfit::LayerOfHits::is_tib_lyr ( ) const

inline

Definition at line 138 of file HitStructures.h.

References mkfit::LayerInfo::is_tib_lyr(), and m_layer_info.

{ return m_layer_info->is_tib_lyr(); }

bool mkfit::LayerOfHits::is_tid_lyr ( ) const

inline

Definition at line 140 of file HitStructures.h.

References mkfit::LayerInfo::is_tid_lyr(), and m_layer_info.

{ return m_layer_info->is_tid_lyr(); }

bool mkfit::LayerOfHits::is_tob_lyr ( ) const

inline

Definition at line 139 of file HitStructures.h.

References mkfit::LayerInfo::is_tob_lyr(), and m_layer_info.

{ return m_layer_info->is_tob_lyr(); }

bool mkfit::LayerOfHits::is_within_r_limits ( float  r ) const

inline

Definition at line 122 of file HitStructures.h.

References mkfit::LayerInfo::is_within_r_limits(), and m_layer_info.

{ return m_layer_info->is_within_r_limits(r); }

WSR_Result mkfit::LayerOfHits::is_within_r_sensitive_region ( float  r, float  dr  ) const

inline

Definition at line 128 of file HitStructures.h.

References mkfit::LayerInfo::is_within_r_sensitive_region(), and m_layer_info.

Referenced by mkfit::MkFinder::selectHitIndices().

                                                                     {
      return m_layer_info->is_within_r_sensitive_region(r, dr);
    }

bool mkfit::LayerOfHits::is_within_z_limits ( float  z ) const

inline

Definition at line 121 of file HitStructures.h.

References mkfit::LayerInfo::is_within_z_limits(), and m_layer_info.

{ return m_layer_info->is_within_z_limits(z); }

WSR_Result mkfit::LayerOfHits::is_within_z_sensitive_region ( float  z, float  dz  ) const

inline

Definition at line 124 of file HitStructures.h.

References mkfit::LayerInfo::is_within_z_sensitive_region(), and m_layer_info.

Referenced by mkfit::MkFinder::selectHitIndices().

                                                                     {
      return m_layer_info->is_within_z_sensitive_region(z, dz);
    }

int mkfit::LayerOfHits::layer_id ( ) const

inline

Definition at line 116 of file HitStructures.h.

References mkfit::LayerInfo::layer_id(), and m_layer_info.

Referenced by mkfit::MkFinder::addBestHit(), mkfit::MkFinder::findCandidates(), mkfit::MkFinder::findCandidatesCloneEngine(), and mkfit::MkFinder::selectHitIndices().

{ return m_layer_info->layer_id(); }

const LayerInfo* mkfit::LayerOfHits::layer_info ( ) const

inline

Definition at line 115 of file HitStructures.h.

References m_layer_info.

Referenced by mkfit::MkFinder::bkFitFitTracks(), mkfit::MkFinder::bkFitFitTracksBH(), and mkfit::MkFinder::selectHitIndices().

{ return m_layer_info; }

int mkfit::LayerOfHits::nHits ( void  ) const

inline

Definition at line 68 of file HitStructures.h.

References m_n_hits.

Referenced by mkfit::MkBuilder::find_tracks_in_layers(), mkfit::MkBuilder::findTracksBestHit(), and mkfit::MkBuilder::findTracksStandard().

{ return m_n_hits; }

bool mkfit::LayerOfHits::phi_bin_dead ( int  qi, int  pi  ) const

inline

Definition at line 89 of file HitStructures.h.

References m_phi_bin_deads, and pi.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return m_phi_bin_deads[qi][pi]; }

const vecvecPhiBinDead_t& mkfit::LayerOfHits::phi_bin_deads ( ) const

inline

Definition at line 87 of file HitStructures.h.

References m_phi_bin_deads.

{ return m_phi_bin_deads; }

PhiBinInfo_t mkfit::LayerOfHits::phi_bin_info ( int  qi, int  pi  ) const

inline

Definition at line 88 of file HitStructures.h.

References m_phi_bin_infos, and pi.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return m_phi_bin_infos[qi][pi]; }

const vecvecPhiBinInfo_t& mkfit::LayerOfHits::phi_bin_infos ( ) const

inline

Definition at line 86 of file HitStructures.h.

References m_phi_bin_infos.

{ return m_phi_bin_infos; }

int mkfit::LayerOfHits::phiBin ( float  phi ) const

inline

Definition at line 78 of file HitStructures.h.

References m_phi_bits_shift, and phiBinFine().

Referenced by phiBinChecked(), mkfit::MkFinder::selectHitIndices(), and suckInDeads().

{ return phiBinFine(phi) >> m_phi_bits_shift; }

int mkfit::LayerOfHits::phiBinChecked ( float  phi ) const

inline

Definition at line 80 of file HitStructures.h.

References m_phi_mask, and phiBin().

{ return phiBin(phi) & m_phi_mask; }

int mkfit::LayerOfHits::phiBinFine ( float  phi ) const

inline

Definition at line 77 of file HitStructures.h.

References m_fphi_fine, and mkfit::Const::PI.

Referenced by phiBin(), registerHit(), and suckInHits().

{ return std::floor(m_fphi_fine * (phi + Const::PI)); }

int mkfit::LayerOfHits::phiMaskApply ( int  in ) const

inline

Definition at line 82 of file HitStructures.h.

References m_phi_mask.

Referenced by mkfit::MkFinder::selectHitIndices().

{ return in & m_phi_mask; }

void mkfit::LayerOfHits::printBins

(

)

Definition at line 355 of file HitStructures.cc.

References first, is_barrel(), mkfit::Config::m_nphi, m_nq, m_phi_bin_infos, gpuVertexFinder::printf(), and edm::second().

                              {
    for (int qb = 0; qb < m_nq; ++qb) {
      printf("%c bin %d\n", is_barrel() ? 'Z' : 'R', qb);
      for (int pb = 0; pb < Config::m_nphi; ++pb) {
        if (pb % 8 == 0)
          printf(" Phi %4d: ", pb);
        printf("%5d,%4d   %s",
               m_phi_bin_infos[qb][pb].first,
               m_phi_bin_infos[qb][pb].second,
               ((pb + 1) % 8 == 0) ? "\n" : "");
      }
    }
  }

int mkfit::LayerOfHits::qBin ( float  q ) const

inline

Definition at line 72 of file HitStructures.h.

References m_fq, and m_qmin.

Referenced by qBinChecked(), and vecPhiBinInfo().

{ return (q - m_qmin) * m_fq; }

int mkfit::LayerOfHits::qBinChecked ( float  q ) const

inline

Definition at line 74 of file HitStructures.h.

References m_nq, and qBin().

Referenced by registerHit(), mkfit::MkFinder::selectHitIndices(), suckInDeads(), and suckInHits().

{ return std::clamp(qBin(q), 0, m_nq - 1); }

const Hit& mkfit::LayerOfHits::refHit ( int  i ) const

inline

Definition at line 103 of file HitStructures.h.

References mps_fire::i, and m_ext_hits.

Referenced by mkfit::MkFinder::addBestHit(), mkfit::MkFinder::bkFitFitTracks(), mkfit::MkFinder::bkFitFitTracksBH(), mkfit::MkFinder::findCandidates(), mkfit::MkFinder::findCandidatesCloneEngine(), mkfit::findSeedsByRoadSearch(), mkfit::MkBuilder::findTracksStandard(), mkfit::MkFinder::selectHitIndices(), and mkfit::MkFinder::updateWithLastHit().

{ return (*m_ext_hits)[i]; }

void mkfit::LayerOfHits::registerHit

(

int

idx

)

Definition at line 179 of file HitStructures.cc.

References h, m_ext_idcs, m_hit_infos, m_is_barrel, m_max_ext_idx, m_min_ext_idx, m_qphifines, SiStripPI::max, SiStripPI::min, mkfit::Hit::phi(), mkfit::LayerOfHits::HitInfo::phi, phiBinFine(), mkfit::LayerOfHits::HitInfo::q, qBinChecked(), mkfit::Hit::r(), mkfit::Config::usePhiQArrays, and mkfit::Hit::z().

                                       {
    const Hit &h = (*m_ext_hits)[idx];

    m_ext_idcs.push_back(idx);
    m_min_ext_idx = std::min(m_min_ext_idx, idx);
    m_max_ext_idx = std::max(m_max_ext_idx, idx);

    HitInfo hi = {h.phi(), m_is_barrel ? h.z() : h.r()};

    m_qphifines.push_back(phiBinFine(hi.phi) + (qBinChecked(hi.q) << 16));

    if (Config::usePhiQArrays) {
      m_hit_infos.emplace_back(hi);
    }
  }

void mkfit::LayerOfHits::reset ( void  )

inline

Definition at line 55 of file HitStructures.h.

{}

void mkfit::LayerOfHits::setup_bins ( float  qmin, float  qmax, float  dq  )

private

Definition at line 30 of file HitStructures.cc.

References reco::ceil(), m_fq, m_nq, m_phi_bin_deads, m_phi_bin_infos, m_qmax, and m_qmin.

Referenced by setupLayer().

                                                               {
    // Define layer with min/max and number of bins along q.

    if (dq < 0) {
      m_nq = (int)-dq;
      m_qmin = qmin;
      m_qmax = qmax;
    } else {
      float extent = qmax - qmin;
      m_nq = std::ceil(extent / dq);
      float extra = 0.5f * (m_nq * dq - extent);
      m_qmin = qmin - extra;
      m_qmax = qmax + extra;
    }
    m_fq = m_nq / (qmax - qmin);  // used in e.g. qbin = (q_hit - m_qmin) * m_fq;

    m_phi_bin_infos.resize(m_nq);
    m_phi_bin_deads.resize(m_nq);
  }

void mkfit::LayerOfHits::setupLayer

(

const LayerInfo &

li

)

Definition at line 50 of file HitStructures.cc.

References cms::cuda::assert(), mkfit::LayerInfo::is_barrel(), m_is_barrel, m_layer_info, mkfit::LayerInfo::q_bin(), mkfit::LayerInfo::rin(), mkfit::LayerInfo::rout(), setup_bins(), mkfit::LayerInfo::zmax(), and mkfit::LayerInfo::zmin().

                                                  {
    // Note, LayerInfo::q_bin( ==>  > 0 - bin width, < 0 - number of bins

    assert(m_layer_info == nullptr && "setupLayer() already called.");

    m_layer_info = &li;

    m_is_barrel = m_layer_info->is_barrel();

    if (m_is_barrel)
      setup_bins(li.zmin(), li.zmax(), li.q_bin());
    else
      setup_bins(li.rin(), li.rout(), li.q_bin());
  }

void mkfit::LayerOfHits::suckInDeads

(

const DeadVec &

deadv

)

Definition at line 139 of file HitStructures.cc.

References cms::cuda::assert(), ztail::d, empty_q_bins_dead(), mkfit::Config::m_nphi, m_nq, m_phi_bin_deads, phiBin(), and qBinChecked().

                                                    {
    assert(m_nq > 0 && "setupLayer() was not called.");

    empty_q_bins_dead(0, m_nq);

    for (const auto &d : deadv) {
      int q_bin_1 = qBinChecked(d.q1);
      int q_bin_2 = qBinChecked(d.q2) + 1;
      int phi_bin_1 = phiBin(d.phi1);
      int phi_bin_2 = phiBin(d.phi2) + 1;
      for (int q_bin = q_bin_1; q_bin < q_bin_2; q_bin++) {
        if (phi_bin_1 > phi_bin_2) {
          for (int pb = phi_bin_1; pb < Config::m_nphi; pb++) {
            m_phi_bin_deads[q_bin][pb] = true;
          }
          for (int pb = 0; pb < phi_bin_2; pb++) {
            m_phi_bin_deads[q_bin][pb] = true;
          }
        } else {
          for (int pb = phi_bin_1; pb < phi_bin_2; pb++) {
            m_phi_bin_deads[q_bin][pb] = true;
          }
        }
      }
    }
  }

void mkfit::LayerOfHits::suckInHits

(

const HitVec &

hitv

)

Definition at line 67 of file HitStructures.cc.

References cms::cuda::assert(), empty_q_bins(), h, mps_fire::i, dqmiolumiharvest::j, m_ext_hits, m_hit_infos, m_hit_phis, m_hit_qs, m_hit_ranks, m_is_barrel, m_n_hits, m_nq, m_phi_bin_infos, m_phi_bits_shift, m_phi_fine_xmask, m_phi_mask_fine, m_qphifines, mkfit::Hit::phi(), mkfit::LayerOfHits::HitInfo::phi, phiBinFine(), mkfit::LayerOfHits::HitInfo::q, qBinChecked(), mkfit::Hit::r(), RADIX_UNSIGNED, RadixSort::RelinquishRanks(), RadixSort::Sort(), mkfit::Config::usePhiQArrays, and mkfit::Hit::z().

                                                 {
    assert(m_nq > 0 && "setupLayer() was not called.");

    m_n_hits = hitv.size();
    m_ext_hits = &hitv;

#ifdef COPY_SORTED_HITS
    if (m_capacity < m_n_hits) {
      free_hits();
      alloc_hits(m_n_hits);
    }
#endif

    if (Config::usePhiQArrays) {
      m_hit_phis.resize(m_n_hits);
      m_hit_qs.resize(m_n_hits);
      m_hit_infos.resize(m_n_hits);
    }
    m_qphifines.resize(m_n_hits);

    for (int i = 0; i < m_n_hits; ++i) {
      const Hit &h = hitv[i];

      HitInfo hi = {h.phi(), m_is_barrel ? h.z() : h.r()};

      m_qphifines[i] = phiBinFine(hi.phi) + (qBinChecked(hi.q) << 16);

      if (Config::usePhiQArrays) {
        m_hit_infos[i] = hi;
      }
    }

    operator delete[](m_hit_ranks);
    {
      RadixSort sort;
      sort.Sort(&m_qphifines[0], m_n_hits, RADIX_UNSIGNED);
      m_hit_ranks = sort.RelinquishRanks();
    }

    int curr_qphi = -1;
    empty_q_bins(0, m_nq, 0);

    for (int i = 0; i < m_n_hits; ++i) {
      int j = m_hit_ranks[i];

#ifdef COPY_SORTED_HITS
      memcpy(&m_hits[i], &hitv[j], sizeof(Hit));
#endif

      if (Config::usePhiQArrays) {
        m_hit_phis[i] = m_hit_infos[j].phi;
        m_hit_qs[i] = m_hit_infos[j].q;
      }

      // Combined q-phi bin with fine part masked off
      const int jqphi = m_qphifines[j] & m_phi_fine_xmask;

      const int phi_bin = (jqphi & m_phi_mask_fine) >> m_phi_bits_shift;
      const int q_bin = jqphi >> 16;

      // Fill the bin info
      if (jqphi != curr_qphi) {
        m_phi_bin_infos[q_bin][phi_bin] = {i, i};
        curr_qphi = jqphi;
      }

      m_phi_bin_infos[q_bin][phi_bin].second++;
    }
  }

const vecPhiBinInfo_t& mkfit::LayerOfHits::vecPhiBinInfo ( float  q ) const

inline

Definition at line 84 of file HitStructures.h.

References m_phi_bin_infos, and qBin().

{ return m_phi_bin_infos[qBin(q)]; }

Member Data Documentation

const HitVec* mkfit::LayerOfHits::m_ext_hits

private

Definition at line 187 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), hitArray(), refHit(), and suckInHits().

std::vector<int> mkfit::LayerOfHits::m_ext_idcs

private

Definition at line 190 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), getHitIndexFromOriginal(), and registerHit().

constexpr float mkfit::LayerOfHits::m_fphi = Config::m_nphi / Const::TwoPI

staticprivate

Definition at line 145 of file HitStructures.h.

constexpr float mkfit::LayerOfHits::m_fphi_fine = 1024 / Const::TwoPI

staticprivate

Definition at line 148 of file HitStructures.h.

Referenced by phiBinFine().

float mkfit::LayerOfHits::m_fq

private

Definition at line 204 of file HitStructures.h.

Referenced by qBin(), and setup_bins().

std::vector<HitInfo> mkfit::LayerOfHits::m_hit_infos

private

Definition at line 213 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), registerHit(), and suckInHits().

std::vector<float> mkfit::LayerOfHits::m_hit_phis

private

Definition at line 199 of file HitStructures.h.

Referenced by endRegistrationOfHits(), hit_phi(), and suckInHits().

std::vector<float> mkfit::LayerOfHits::m_hit_qs

private

Definition at line 200 of file HitStructures.h.

Referenced by endRegistrationOfHits(), hit_q(), and suckInHits().

unsigned int* mkfit::LayerOfHits::m_hit_ranks = nullptr

private

Definition at line 189 of file HitStructures.h.

Referenced by endRegistrationOfHits(), getOriginalHitIndex(), suckInHits(), and ~LayerOfHits().

bool mkfit::LayerOfHits::m_is_barrel

private

Definition at line 206 of file HitStructures.h.

Referenced by is_barrel(), is_endcap(), registerHit(), setupLayer(), and suckInHits().

const LayerInfo* mkfit::LayerOfHits::m_layer_info = nullptr

private

Definition at line 203 of file HitStructures.h.

Referenced by is_pix_lyr(), is_pixb_lyr(), is_pixe_lyr(), is_stereo(), is_tec_lyr(), is_tib_lyr(), is_tid_lyr(), is_tob_lyr(), is_within_r_limits(), is_within_r_sensitive_region(), is_within_z_limits(), is_within_z_sensitive_region(), layer_id(), layer_info(), and setupLayer().

int mkfit::LayerOfHits::m_max_ext_idx

private

Definition at line 191 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), and registerHit().

int mkfit::LayerOfHits::m_min_ext_idx

private

Definition at line 191 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), getHitIndexFromOriginal(), and registerHit().

int mkfit::LayerOfHits::m_n_hits = 0

private

Definition at line 192 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), nHits(), and suckInHits().

int mkfit::LayerOfHits::m_nq = 0

private

Definition at line 205 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), printBins(), qBinChecked(), setup_bins(), suckInDeads(), and suckInHits().

vecvecPhiBinDead_t mkfit::LayerOfHits::m_phi_bin_deads

private

Definition at line 196 of file HitStructures.h.

Referenced by empty_phi_bins_dead(), phi_bin_dead(), phi_bin_deads(), setup_bins(), and suckInDeads().

vecvecPhiBinInfo_t mkfit::LayerOfHits::m_phi_bin_infos

private

Definition at line 195 of file HitStructures.h.

Referenced by empty_phi_bins(), endRegistrationOfHits(), phi_bin_info(), phi_bin_infos(), printBins(), setup_bins(), suckInHits(), and vecPhiBinInfo().

constexpr int mkfit::LayerOfHits::m_phi_bits = 8

staticprivate

Definition at line 147 of file HitStructures.h.

constexpr int mkfit::LayerOfHits::m_phi_bits_fine = 10

staticprivate

Definition at line 150 of file HitStructures.h.

constexpr int mkfit::LayerOfHits::m_phi_bits_shift = m_phi_bits_fine - m_phi_bits

staticprivate

Definition at line 151 of file HitStructures.h.

Referenced by endRegistrationOfHits(), phiBin(), and suckInHits().

constexpr int mkfit::LayerOfHits::m_phi_fine_xmask = ~((1 << m_phi_bits_shift) - 1)

staticprivate

Definition at line 152 of file HitStructures.h.

Referenced by endRegistrationOfHits(), and suckInHits().

constexpr int mkfit::LayerOfHits::m_phi_mask = 0xff

staticprivate

Definition at line 146 of file HitStructures.h.

Referenced by phiBinChecked(), and phiMaskApply().

constexpr int mkfit::LayerOfHits::m_phi_mask_fine = 0x3ff

staticprivate

Definition at line 149 of file HitStructures.h.

Referenced by endRegistrationOfHits(), and suckInHits().

float mkfit::LayerOfHits::m_qmax

private

Definition at line 204 of file HitStructures.h.

Referenced by setup_bins().

float mkfit::LayerOfHits::m_qmin

private

Definition at line 204 of file HitStructures.h.

Referenced by qBin(), and setup_bins().

std::vector<uint32_t> mkfit::LayerOfHits::m_qphifines

private

Definition at line 214 of file HitStructures.h.

Referenced by beginRegistrationOfHits(), endRegistrationOfHits(), registerHit(), and suckInHits().