FPGATrackSimPlaneMap Class Reference

#include <FPGATrackSimPlaneMap.h>

Collaboration diagram for FPGATrackSimPlaneMap:

Public Member Functions
	FPGATrackSimPlaneMap (std::ifstream &fin, unsigned region, unsigned stage, std::vector< int > layerOverrides=std::vector< int >())
uint32_t	getNDetLayers () const
uint32_t	getNLogiLayers () const
uint32_t	getNCoords () const
uint32_t	getNSections (size_t logiLayer) const
uint32_t	getDim (size_t logiLayer) const
SiliconTech	getDetType (size_t logiLayer) const
uint32_t	getCoordOffset (size_t logiLayer) const
uint32_t	getCoordLayer (uint32_t coord) const
const LayerInfo &	getLayerInfo (uint32_t layer, uint32_t section) const
bool	isSCT (int pl) const
bool	isPixel (int pl) const
bool	isEC (uint32_t layer, uint32_t section) const
std::string	layerName (uint32_t layer, uint32_t section) const
const LayerSection &	getLayerSection (SiliconTech siTech, DetectorZone zone, uint32_t physLayer) const
const LayerSection &	getLayerSection (LayerInfo const &mi) const
int	getLayer (SiliconTech siTech, DetectorZone zone, uint32_t physLayer) const
unsigned	getSection (SiliconTech siTech, DetectorZone zone, uint32_t physLayer) const
void	doRemap (FPGATrackSimHit &hit) const
void	map (FPGATrackSimHit &hit) const

Private Member Functions
void	allocateMap (std::ifstream &fin, uint32_t stage)
void	seek (std::ifstream &fin, unsigned region)
void	readLayers (std::ifstream &fin, uint32_t stage)
int	getOverrideLayer (SiliconTech si, DetectorZone dz, int physLayer)

Private Attributes
uint32_t	m_nLogiLayers = 0
uint32_t	m_nDetLayers = 0
uint32_t	m_nCoords = 0
std::vector< uint32_t >	m_dimension
std::vector< uint32_t >	m_coordOffset
std::vector< uint32_t >	m_coordLayer
std::vector< std::vector< std::vector< LayerSection > > >	m_map
std::vector< std::vector< LayerInfo > >	m_layerInfo
std::vector< int >	m_layerOverrides
std::vector< uint32_t >	m_diskIndex
std::unique_ptr< FPGATrackSimModuleRelabel >	m_moduleRelabel {}
const LayerSection	m_inValidLayerSection {}

Detailed Description

Definition at line 61 of file FPGATrackSimPlaneMap.h.

Constructor & Destructor Documentation

FPGATrackSimPlaneMap()

FPGATrackSimPlaneMap::FPGATrackSimPlaneMap	(	std::ifstream &	fin,
		unsigned	region,
		unsigned	stage,
		std::vector< int >	layerOverrides = std::vector<int>() )

Definition at line 26 of file FPGATrackSimPlaneMap.cxx.

                                                                                                                           :
    m_map(static_cast<int>(SiliconTech::nTechs),
          vector<vector<LayerSection>>(static_cast<int>(DetectorZone::nZones))
    ),
    m_layerOverrides(std::move(layerOverrides))
{
    // Reads the header of the file to resize all the vector members
    allocateMap(fin, stage);
 
    // Seek to the correct region
    seek(fin, region);
 
    // Reads the rest of the file to populate all the member vectors
    readLayers(fin, stage);
 
    // Postprocessing on coordinate indices
    for (uint32_t l = 0; l < m_nLogiLayers; l++)
    {
        m_coordOffset[l] = m_nCoords;
        m_nCoords += m_dimension[l];
        for (uint32_t i = 0; i < m_dimension[l]; i++)
            m_coordLayer.push_back(l);
    }
 
    ANA_MSG_DEBUG("Using " << m_nLogiLayers << " logical layers and " << m_nCoords << " coordinates");
}

Member Function Documentation

allocateMap()

void FPGATrackSimPlaneMap::allocateMap ( std::ifstream & fin, uint32_t stage )

private

Definition at line 53 of file FPGATrackSimPlaneMap.cxx.

{
    // Initialize read variables
    vector<int> layerCounts((int)SiliconTech::nTechs * static_cast<int>(DetectorZone::nZones)); // pixel_barrel, pixel_EC, SCT_barrel, SCT_EC
    std::string line, silicon, detReg, layerKey, geoKey;
    bool ok = true;
    // Read Geometry Version
    ok = ok && getline(fin, line);
    while (line.empty() || line[0] == '!')
    {
        ok = ok && getline(fin, line);
    }
    ANA_MSG_VERBOSE(line);
    istringstream sline(line);
    ok = ok && (sline >> geoKey);
    m_diskIndex = Remappings::diskIndices(geoKey);
 
    ANA_MSG_DEBUG("Allocating map for geometry " << geoKey <<" diskIndex size="<<m_diskIndex.size());
    m_moduleRelabel = std::make_unique<FPGATrackSimModuleRelabel>(geoKey, false);
 
 
 
    // Read number of logical layers
    if (stage == 1)
    {
        ok = ok && getline(fin, line);
        ANA_MSG_VERBOSE(line);
        istringstream sline(line);
        ok = ok && (sline >> m_nLogiLayers >> layerKey);
        ok = ok && (layerKey == "logical_s1");
        ok = ok && getline(fin, line); // skip the 2nd stage line
    }
    else if (stage == 2)
    {
        ok = ok && getline(fin, line); // skip the 1st stage line
        ok = ok && getline(fin, line);
        ANA_MSG_VERBOSE(line);
        istringstream sline(line);
        ok = ok && (sline >> m_nLogiLayers >> layerKey);
        ok = ok && (layerKey == "logical_s2");
    }
    else ANA_MSG_ERROR("Bad stage " << stage << ". Stage should be 1 or 2");
 
    if (!m_layerOverrides.empty()) m_nLogiLayers = m_layerOverrides.size();
 
    // Read number of physical layers
    int nHeaderLines = 6;
    for (int i = 0; i < nHeaderLines; i++)
    {
        if (!ok) break;
        ok = ok && getline(fin, line);
        ANA_MSG_VERBOSE(line);
 
        istringstream sline(line);
        ok = ok && (sline >> layerCounts[i] >> silicon >> detReg);
 
        ok = ok && ( (i < (nHeaderLines/2) && silicon == "pixel") || (i >= (nHeaderLines/2) && silicon == "SCT") );
        ok = ok && ( (i % (nHeaderLines/2) == 0 && detReg == "barrel") || (i % (nHeaderLines/2) != 0 && (detReg == "endcap+" || detReg == "endcap-")) );
    }
    
    if (!ok) ANA_MSG_FATAL("Error reading layer counts");
    m_nDetLayers = std::accumulate(layerCounts.begin(), layerCounts.end(), 0);
 
    // Resize objects
    m_map[(int)SiliconTech::pixel][(int)DetectorZone::barrel].resize(layerCounts[0]);
    m_map[(int)SiliconTech::pixel][(int)DetectorZone::posEndcap].resize(layerCounts[1]);
    m_map[(int)SiliconTech::pixel][(int)DetectorZone::negEndcap].resize(layerCounts[2]);
    m_map[(int)SiliconTech::strip][(int)DetectorZone::barrel].resize(layerCounts[3]);
    m_map[(int)SiliconTech::strip][(int)DetectorZone::posEndcap].resize(layerCounts[4]);
    m_map[(int)SiliconTech::strip][(int)DetectorZone::negEndcap].resize(layerCounts[5]);
 
    m_dimension.resize(m_nLogiLayers);
    m_layerInfo.resize(m_nLogiLayers);
    m_coordOffset.resize(m_nLogiLayers);
}

doRemap()

void FPGATrackSimPlaneMap::doRemap

(

FPGATrackSimHit &

hit

)

const

Definition at line 227 of file FPGATrackSimPlaneMap.cxx.

{
   if (!hit.isRemapped()) {
       m_moduleRelabel->remap(hit);
    }
}

getCoordLayer()

uint32_t FPGATrackSimPlaneMap::getCoordLayer ( uint32_t coord ) const

inline

Definition at line 98 of file FPGATrackSimPlaneMap.h.

{ return m_coordLayer[coord]; }

getCoordOffset()

uint32_t FPGATrackSimPlaneMap::getCoordOffset ( size_t logiLayer ) const

inline

Definition at line 92 of file FPGATrackSimPlaneMap.h.

{ return m_coordOffset[logiLayer]; }

getDetType()

SiliconTech FPGATrackSimPlaneMap::getDetType ( size_t logiLayer ) const

inline

Definition at line 79 of file FPGATrackSimPlaneMap.h.

{ return m_dimension[logiLayer] == 2 ? SiliconTech::pixel : SiliconTech::strip; }

getDim()

uint32_t FPGATrackSimPlaneMap::getDim ( size_t logiLayer ) const

inline

Definition at line 78 of file FPGATrackSimPlaneMap.h.

{ return m_dimension[logiLayer]; }

getLayer()

int FPGATrackSimPlaneMap::getLayer ( SiliconTech siTech, DetectorZone zone, uint32_t physLayer ) const

inline

Definition at line 132 of file FPGATrackSimPlaneMap.h.

            { return getLayerSection(siTech, zone, physLayer).layer; }

getLayerInfo()

const LayerInfo & FPGATrackSimPlaneMap::getLayerInfo ( uint32_t layer, uint32_t section ) const

inline

Definition at line 104 of file FPGATrackSimPlaneMap.h.

{ return m_layerInfo.at(layer).at(section); }

getLayerSection() [1/2]

const LayerSection & FPGATrackSimPlaneMap::getLayerSection ( LayerInfo const & mi ) const

inline

Definition at line 129 of file FPGATrackSimPlaneMap.h.

            { return getLayerSection(mi.siTech, mi.zone, mi.physLayer); }

getLayerSection() [2/2]

const LayerSection & FPGATrackSimPlaneMap::getLayerSection ( SiliconTech siTech, DetectorZone zone, uint32_t physLayer ) const

inline

Definition at line 117 of file FPGATrackSimPlaneMap.h.

        {
            if ( static_cast<size_t>(siTech) < m_map.size() ) {
                if ( static_cast<size_t>(zone) < m_map.at(static_cast<size_t>(siTech)).size() ) {
                    if ( physLayer < m_map.at(static_cast<size_t>(siTech)).at(static_cast<size_t>(zone)).size() ) {
                        return m_map[static_cast<size_t>(siTech)][static_cast<size_t>(zone)][physLayer];
                    }
                }
            }
            return m_inValidLayerSection;
        }

getNCoords()

uint32_t FPGATrackSimPlaneMap::getNCoords ( ) const

inline

Definition at line 76 of file FPGATrackSimPlaneMap.h.

{ return m_nCoords; }

getNDetLayers()

uint32_t FPGATrackSimPlaneMap::getNDetLayers ( ) const

inline

Definition at line 74 of file FPGATrackSimPlaneMap.h.

{ return m_nDetLayers; }

getNLogiLayers()

uint32_t FPGATrackSimPlaneMap::getNLogiLayers ( ) const

inline

Definition at line 75 of file FPGATrackSimPlaneMap.h.

{ return m_nLogiLayers; }

getNSections()

uint32_t FPGATrackSimPlaneMap::getNSections ( size_t logiLayer ) const

inline

Definition at line 77 of file FPGATrackSimPlaneMap.h.

{ return m_layerInfo[logiLayer].size(); }

getOverrideLayer()

int FPGATrackSimPlaneMap::getOverrideLayer ( SiliconTech si, DetectorZone dz, int physLayer )

private

Definition at line 259 of file FPGATrackSimPlaneMap.cxx.

{
    int code = static_cast<int>(si) * 1000 + static_cast<int>(dz) * 100 + physLayer;
    for (unsigned i = 0; i < m_layerOverrides.size(); i++)
        if (code == m_layerOverrides[i]) return i;
    return -1;
}

getSection()

unsigned FPGATrackSimPlaneMap::getSection ( SiliconTech siTech, DetectorZone zone, uint32_t physLayer ) const

inline

Definition at line 135 of file FPGATrackSimPlaneMap.h.

            { return getLayerSection(siTech, zone, physLayer).section; }

isEC()

bool FPGATrackSimPlaneMap::isEC ( uint32_t layer, uint32_t section ) const

inline

Definition at line 108 of file FPGATrackSimPlaneMap.h.

{ return getLayerInfo(layer, section).zone != DetectorZone::barrel; }

isPixel()

bool FPGATrackSimPlaneMap::isPixel ( int pl ) const

inline

Definition at line 107 of file FPGATrackSimPlaneMap.h.

{ return m_dimension[pl] == 2; }

isSCT()

bool FPGATrackSimPlaneMap::isSCT ( int pl ) const

inline

Definition at line 106 of file FPGATrackSimPlaneMap.h.

{ return m_dimension[pl] == 1; }

layerName()

std::string FPGATrackSimPlaneMap::layerName	(	uint32_t	layer,
		uint32_t	section ) const

Definition at line 267 of file FPGATrackSimPlaneMap.cxx.

{
    std::string out;
    out += (isPixel(layer) ? "PIX" : "SCT");
    out += (isEC(layer, section) ? "_EC_" : "_B_");
    out += std::to_string(getLayerInfo(layer, section).physLayer);
    return out;
}

map()

void FPGATrackSimPlaneMap::map

(

FPGATrackSimHit &

hit

)

const

Definition at line 234 of file FPGATrackSimPlaneMap.cxx.

{
    // re-assign layers in the pixel endcap to be each individual disk
    // technically this returns a success/fail but I'm not sure we need it?
    // This should only happen if the hit is not already mapped.
    const LayerSection &pinfo = getLayerSection(hit.getDetType(), hit.getDetectorZone(), hit.getPhysLayer());
    hit.setSection(pinfo.section);
    hit.setLayer(pinfo.layer);
    if (!hit.isMapped() && hit.getHitType() == HitType::unmapped) { // failsafe if for some reason someone calls this on a clustered hit again, or something
        hit.setHitType(HitType::mapped);
    }    
    // Special case if this is a spacepoint (now possible because slicing engine happens after SP formation).
    if (hit.getHitType() == HitType::spacepoint) {
        const LayerSection &pinfo_sp = getLayerSection(hit.getPairedDetType(), hit.getPairedDetZone(), hit.getPairedPhysLayer());
        hit.setPairedSection(pinfo_sp.section);
        hit.setPairedLayer(pinfo_sp.layer);
    }
}

readLayers()

void FPGATrackSimPlaneMap::readLayers ( std::ifstream & fin, uint32_t stage )

private

Definition at line 150 of file FPGATrackSimPlaneMap.cxx.

{
    // Initialize read variables
    int BEC{}, physLayer{}, physDisk{}, logiLayer{}, logiLayer1{}, logiLayer2{}, stereo{}, nDim{};
    string line, silicon, planeKey1, planeKey2, stereoKey;
    bool ok = true;
 
    uint32_t linesRead = 0;
 
    // Each line is a detector layer
    while (getline(fin, line))
    {
        if (line.empty() || line[0] == '!') continue;
        ANA_MSG_VERBOSE(line);
        istringstream sline(line);
 
        ok = ok && (sline >> silicon >> BEC >> physDisk >> physLayer);
 
        // Detector
        DetectorZone zone  = DetectorZone::undefined;
        if (BEC == 0) zone = DetectorZone::barrel;
        else if (BEC == 1) zone = DetectorZone::posEndcap;
        else if (BEC == 2) zone = DetectorZone::negEndcap;
        else ANA_MSG_FATAL("Bad Detector zone" << BEC);
 
        // Case on silicon type
        int sil = 0;
        SiliconTech siTech = SiliconTech::pixel;
        if (silicon == "pixel")
        {
            sil = static_cast<int>(SiliconTech::pixel);
            nDim = 2;
            stereo = -1;
            ok = ok && (sline >> planeKey1 >> logiLayer1 >> planeKey2 >> logiLayer2);
            ok = ok && (planeKey1 == "plane1" && planeKey2 == "plane2");
        }
        else if (silicon == "SCT")
        {
            sil = static_cast<int>(SiliconTech::strip);
            siTech = SiliconTech::strip;
 
            nDim = 1;
            ok = ok && (sline >> stereoKey >> stereo >> planeKey1 >> logiLayer1 >> planeKey2 >> logiLayer2);
            ok = ok && (planeKey1 == "plane1" && planeKey2 == "plane2" && stereoKey == "stereo");
        }
        else ANA_MSG_FATAL("Bad silicon type " << silicon);
 
 
        if (!ok) break;
 
        // Logical layer
        if (!m_layerOverrides.empty()) logiLayer = getOverrideLayer(siTech, zone, physLayer);
        else if (stage == 1) logiLayer = logiLayer1;
        else if (stage == 2) logiLayer = logiLayer2;
 
        // Check number of layers
        linesRead++;
        if (logiLayer >= (int)m_nLogiLayers) ANA_MSG_FATAL("Logical layer " << logiLayer << " exceeds expected number " << m_nLogiLayers);
 
        // Store to map
        m_map.at(sil).at(BEC).at(physLayer).layer = logiLayer;
        if (logiLayer >= 0)
        {
            m_dimension[logiLayer] = nDim;
            m_map[sil][BEC][physLayer].section = m_layerInfo[logiLayer].size(); // i.e. index into m_layerInfo[logiLayer] entry below
            m_layerInfo[logiLayer].push_back({ siTech, zone, physLayer, physDisk, stereo});
        }
        if (m_nDetLayers == linesRead) break;
    }
 
    if (!ok) ANA_MSG_FATAL("Error reading file");
}

seek()

void FPGATrackSimPlaneMap::seek ( std::ifstream & fin, unsigned region )

private

Definition at line 130 of file FPGATrackSimPlaneMap.cxx.

{
    std::string line, key;
    unsigned region_read = -1;
 
    while (getline(fin, line))
    {
        if (line.empty() || line[0] == '!') continue;
        ANA_MSG_VERBOSE(line);
        istringstream sline(line);
 
        bool ok = static_cast<bool>(sline >> key >> region_read);
        if (ok && key == "region" && region_read == region) break;
    }
 
    ANA_MSG_DEBUG("seek() ended at:" << line);
}

Member Data Documentation

m_coordLayer

std::vector<uint32_t> FPGATrackSimPlaneMap::m_coordLayer

private

Definition at line 156 of file FPGATrackSimPlaneMap.h.

m_coordOffset

std::vector<uint32_t> FPGATrackSimPlaneMap::m_coordOffset

private

Definition at line 152 of file FPGATrackSimPlaneMap.h.

m_dimension

std::vector<uint32_t> FPGATrackSimPlaneMap::m_dimension

private

Definition at line 151 of file FPGATrackSimPlaneMap.h.

m_diskIndex

std::vector<uint32_t> FPGATrackSimPlaneMap::m_diskIndex

private

Definition at line 168 of file FPGATrackSimPlaneMap.h.

m_inValidLayerSection

const LayerSection FPGATrackSimPlaneMap::m_inValidLayerSection {}

private

Definition at line 174 of file FPGATrackSimPlaneMap.h.

{};

m_layerInfo

std::vector<std::vector<LayerInfo> > FPGATrackSimPlaneMap::m_layerInfo

private

Definition at line 161 of file FPGATrackSimPlaneMap.h.

m_layerOverrides

std::vector<int> FPGATrackSimPlaneMap::m_layerOverrides

private

Definition at line 164 of file FPGATrackSimPlaneMap.h.

m_map

std::vector<std::vector<std::vector<LayerSection> > > FPGATrackSimPlaneMap::m_map

private

Definition at line 159 of file FPGATrackSimPlaneMap.h.

m_moduleRelabel

std::unique_ptr<FPGATrackSimModuleRelabel> FPGATrackSimPlaneMap::m_moduleRelabel {}

private

Definition at line 171 of file FPGATrackSimPlaneMap.h.

{};

m_nCoords

uint32_t FPGATrackSimPlaneMap::m_nCoords = 0

private

Definition at line 149 of file FPGATrackSimPlaneMap.h.

m_nDetLayers

uint32_t FPGATrackSimPlaneMap::m_nDetLayers = 0

private

Definition at line 148 of file FPGATrackSimPlaneMap.h.

m_nLogiLayers

uint32_t FPGATrackSimPlaneMap::m_nLogiLayers = 0

private

Definition at line 147 of file FPGATrackSimPlaneMap.h.

---

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

• FPGATrackSimPlaneMap.h
• FPGATrackSimPlaneMap.cxx