FPGATrackSimSectorBank Class Reference

#include <FPGATrackSimSectorBank.h>

Collaboration diagram for FPGATrackSimSectorBank:

Public Member Functions
	FPGATrackSimSectorBank (std::string const &filepath)
size_t	getNLayers () const
size_t	getNSectors () const
std::vector< module_t > const &	getModules (sector_t sector) const
std::vector< sector_t >	getSectors (unsigned layer, module_t module) const
sector_t	findSector (std::vector< module_t > const &modules) const
sector_t	findSector (std::vector< std::vector< const FPGATrackSimHit * >> const &hits) const
const std::vector< double > &	getQOverPtBins () const
size_t	getNQOverPtBins () const
bool	isAbsQOverPtBinning () const
void	storeQOverPtBinning (const std::vector< double > &qOverPtBins, bool absBinning)

Private Member Functions
void	readHeader (std::ifstream &fin)
void	readSectors (std::ifstream &fin)
void	makeInverseMap ()

Private Attributes
size_t	m_nLayers = 0UL
std::vector< std::vector< module_t > >	m_s2mMap
std::vector< std::multimap< module_t, sector_t > >	m_m2sMap
std::vector< double >	m_qOverPtBins
bool	m_absQOverPtBinning = false

Detailed Description

Definition at line 29 of file FPGATrackSimSectorBank.h.

Constructor & Destructor Documentation

FPGATrackSimSectorBank()

FPGATrackSimSectorBank::FPGATrackSimSectorBank

(

std::string const &

filepath

)

Definition at line 23 of file FPGATrackSimSectorBank.cxx.

{
    // Open the file for reading
  ANA_MSG_INFO("Reading " << filepath);
    std::ifstream fin(filepath);
    if (!fin.is_open())
    {
      ANA_MSG_ERROR("Couldn't open " << filepath);
      throw ("FPGATrackSimSectorBank couldn't open " + filepath);
    }
    
    // Retrieve the number of sectors and layers, resize m_s2mMap
    readHeader(fin);
 
    // Read the modules from each sector
    readSectors(fin);
 
    // Make m2sMap
    makeInverseMap();
}

Member Function Documentation

findSector() [1/2]

sector_t FPGATrackSimSectorBank::findSector

(

std::vector< module_t > const &

modules

)

const

Definition at line 121 of file FPGATrackSimSectorBank.cxx.

{
    assert(modules.size() == m_nLayers);
 
    std::set<sector_t> sectors_good; // List of matching sectors so far.
    for (size_t layer = 0; layer < m_nLayers; layer++)
    {
        std::set<sector_t> sectors_new;
        auto range = m_m2sMap[layer].equal_range(modules[layer]);
 
        for (auto it_m2s = range.first; it_m2s != range.second; it_m2s++)
            if (layer == 0 || sectors_good.count(it_m2s->second))
                sectors_new.insert(it_m2s->second);
 
        if (sectors_new.empty()) return -1;
        sectors_good = sectors_new;
    }
    
    assert(sectors_good.size() == 1);
    return *(sectors_good.begin());
}

findSector() [2/2]

sector_t FPGATrackSimSectorBank::findSector

(

std::vector< std::vector< const FPGATrackSimHit * >> const &

hits

)

const

Definition at line 146 of file FPGATrackSimSectorBank.cxx.

{
    std::vector<module_t> modules(hits.size());
 
    for (size_t i = 0; i < hits.size(); i++)
    {
        if (hits[i].empty())
        {
            modules[i] = MODULE_BADMODULE;
        }
        else if (hits[i].size() == 1)
        {
            modules[i] = hits[i][0]->getIdentifierHash();
    }
        else
        {
            std::unordered_map<unsigned, unsigned> hashCount;
            for (const FPGATrackSimHit* h : hits[i]) hashCount[h->getIdentifierHash()]++;
            modules[i] = std::max_element(hashCount.begin(), hashCount.end())->first;
        }
    }
 
    return findSector(modules);
}

getModules()

std::vector<module_t> const& FPGATrackSimSectorBank::getModules ( sector_t  sector ) const

inline

Definition at line 38 of file FPGATrackSimSectorBank.h.

{ return m_s2mMap[sector]; }

getNLayers()

size_t FPGATrackSimSectorBank::getNLayers ( ) const

inline

Definition at line 35 of file FPGATrackSimSectorBank.h.

{ return m_nLayers; }

getNQOverPtBins()

size_t FPGATrackSimSectorBank::getNQOverPtBins ( ) const

inline

Definition at line 51 of file FPGATrackSimSectorBank.h.

{return m_qOverPtBins.size(); }

getNSectors()

size_t FPGATrackSimSectorBank::getNSectors ( ) const

inline

Definition at line 36 of file FPGATrackSimSectorBank.h.

{ return m_s2mMap.size(); }

getQOverPtBins()

const std::vector<double>& FPGATrackSimSectorBank::getQOverPtBins ( ) const

inline

Definition at line 50 of file FPGATrackSimSectorBank.h.

{ return m_qOverPtBins; }

getSectors()

std::vector< sector_t > FPGATrackSimSectorBank::getSectors	(	unsigned	layer,
		module_t	module
	)		const

Definition at line 110 of file FPGATrackSimSectorBank.cxx.

{
    std::vector<sector_t> out;
    auto range = m_m2sMap[layer].equal_range(module);
    for (auto it_m2s = range.first; it_m2s != range.second; it_m2s++)
        out.push_back(it_m2s->second);
 
    return out;
}

isAbsQOverPtBinning()

bool FPGATrackSimSectorBank::isAbsQOverPtBinning ( ) const

inline

Definition at line 52 of file FPGATrackSimSectorBank.h.

{ return m_absQOverPtBinning; }

makeInverseMap()

void FPGATrackSimSectorBank::makeInverseMap ( )

private

Definition at line 96 of file FPGATrackSimSectorBank.cxx.

{
    m_m2sMap.resize(m_nLayers);
    for (size_t layer = 0; layer < m_nLayers; layer++)
    {
        for (sector_t sector = 0; (size_t)sector < m_s2mMap.size(); sector++)
        {
            module_t moduleID = m_s2mMap[sector][layer];
            m_m2sMap[layer].insert(std::make_pair(moduleID, sector));
        }
    }
}

readHeader()

void FPGATrackSimSectorBank::readHeader ( std::ifstream &  fin )

private

Definition at line 46 of file FPGATrackSimSectorBank.cxx.

{
    bool ok = true;
    std::string line;
    int nSectors;
 
    ok = ok && std::getline(fin, line);
    ANA_MSG_DEBUG(line);
 
    std::istringstream sline(line);
    ok = ok && (sline >> nSectors >> m_nLayers);
 
    if (!ok) throw "Error reading header";
 
    m_s2mMap.resize(nSectors, std::vector<module_t>(m_nLayers));
}

readSectors()

void FPGATrackSimSectorBank::readSectors ( std::ifstream &  fin )

private

Definition at line 64 of file FPGATrackSimSectorBank.cxx.

{
    std::string line;
    bool ok = true;
    int sector = 0; // make sure sectors are numbered in order
    int sector_file, UNUSED, coverage;
 
    while (std::getline(fin, line))
    {
        if (line.empty() || line[0] == '!') continue;
        ANA_MSG_DEBUG(line);
        std::istringstream sline(line);
 
        ok = ok && (sline >> sector_file);
        ok = ok && (sector_file == sector);
 
        for (unsigned i = 0; i < m_nLayers; i++)
            ok = ok && (sline >> m_s2mMap[sector][i]);
 
        ok = ok && (sline >> UNUSED >> coverage); // both these are not stored
 
        if (!ok) break;
        sector++;
    }
 
    if (!ok) ANA_MSG_ERROR("Error reading file on line: " << line);
}

storeQOverPtBinning()

void FPGATrackSimSectorBank::storeQOverPtBinning	(	const std::vector< double > &	qOverPtBins,
		bool	absBinning
	)

Definition at line 173 of file FPGATrackSimSectorBank.cxx.

{
  m_absQOverPtBinning = absBinning;
 
  // Retrieve q/pt binning information. Fallback to no binning if not set.
  // Note that this behavior on whether or not we are using |q/pt|.
  m_qOverPtBins = qOverPtBins;
  if (m_qOverPtBins.size() == 0) {
    ANA_MSG_FATAL("q/pt bin information not set in bank tag!");
  }
}

Member Data Documentation

m_absQOverPtBinning

bool FPGATrackSimSectorBank::m_absQOverPtBinning = false

private

Definition at line 64 of file FPGATrackSimSectorBank.h.

m_m2sMap

std::vector<std::multimap<module_t, sector_t> > FPGATrackSimSectorBank::m_m2sMap

private

Definition at line 60 of file FPGATrackSimSectorBank.h.

m_nLayers

size_t FPGATrackSimSectorBank::m_nLayers = 0UL

private

Definition at line 58 of file FPGATrackSimSectorBank.h.

m_qOverPtBins

std::vector<double> FPGATrackSimSectorBank::m_qOverPtBins

private

Definition at line 63 of file FPGATrackSimSectorBank.h.

m_s2mMap

std::vector<std::vector<module_t> > FPGATrackSimSectorBank::m_s2mMap

private

Definition at line 59 of file FPGATrackSimSectorBank.h.

---

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

• FPGATrackSimSectorBank.h
• FPGATrackSimSectorBank.cxx