#include <FPGATrackSimPlaneMap.h>

Collaboration diagram for FPGATrackSimPlaneMap:

Public Member Functions
	FPGATrackSimPlaneMap (const std::string &filepath, 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	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 {}

Detailed Description

Definition at line 61 of file FPGATrackSimPlaneMap.h.

Constructor & Destructor Documentation

◆ FPGATrackSimPlaneMap()

FPGATrackSimPlaneMap::FPGATrackSimPlaneMap	(	const std::string &	filepath,
		unsigned	region,
		unsigned	stage,
		std::vector< int >	layerOverrides = `std::vector<int>()`
	)

Definition at line 27 of file FPGATrackSimPlaneMap.cxx.

                                                                                                                                      :
     m_map(static_cast<int>(SiliconTech::nTechs),
           vector<vector<LayerSection>>(static_cast<int>(DetectorZone::nZones))
     ),
     m_layerOverrides(std::move(layerOverrides))
 {
     ANA_MSG_INFO("Reading " << filepath);
     ifstream fin(filepath);
     if (!fin.is_open())
     {
         ANA_MSG_FATAL("Couldn't open " << filepath);
         throw ("FPGATrackSimPlaneMap Couldn't open " + filepath);
     }
  
     // 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_INFO("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 64 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);
     ANA_MSG_VERBOSE(line);
     istringstream sline(line);
     ok = ok && (sline >> geoKey);
     m_diskIndex = Remappings::diskIndices(geoKey);
  
     ANA_MSG_INFO("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);
 }

◆ getCoordLayer()

uint32_t FPGATrackSimPlaneMap::getCoordLayer ( uint32_t coord ) const

inline

Definition at line 101 of file FPGATrackSimPlaneMap.h.

101 { return m_coordLayer[coord]; }

◆ getCoordOffset()

uint32_t FPGATrackSimPlaneMap::getCoordOffset ( size_t logiLayer ) const

inline

Definition at line 95 of file FPGATrackSimPlaneMap.h.

95 { return m_coordOffset[logiLayer]; }

◆ getDetType()

SiliconTech FPGATrackSimPlaneMap::getDetType ( size_t logiLayer ) const

inline

Definition at line 82 of file FPGATrackSimPlaneMap.h.

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

◆ getDim()

uint32_t FPGATrackSimPlaneMap::getDim ( size_t logiLayer ) const

inline

Definition at line 81 of file FPGATrackSimPlaneMap.h.

81 { return m_dimension[logiLayer]; }

◆ getLayer()

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

inline

Definition at line 134 of file FPGATrackSimPlaneMap.h.

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

◆ getLayerInfo()

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

inline

Definition at line 107 of file FPGATrackSimPlaneMap.h.

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

◆ getLayerSection() [1/2]

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

inline

Definition at line 131 of file FPGATrackSimPlaneMap.h.

132 { 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 120 of file FPGATrackSimPlaneMap.h.

         {
             try {
                 return m_map.at(static_cast<int>(siTech)).at(static_cast<int>(zone)).at(physLayer);
             } catch ( std::out_of_range const & e)
             {
                 // TMP HACK FIX REMOVE ME
                 return m_map.at(1).at(0).at(0);
             }
         }

◆ getNCoords()

uint32_t FPGATrackSimPlaneMap::getNCoords ( ) const

inline

Definition at line 79 of file FPGATrackSimPlaneMap.h.

79 { return m_nCoords; }

◆ getNDetLayers()

uint32_t FPGATrackSimPlaneMap::getNDetLayers ( ) const

inline

Definition at line 77 of file FPGATrackSimPlaneMap.h.

77 { return m_nDetLayers; }

◆ getNLogiLayers()

uint32_t FPGATrackSimPlaneMap::getNLogiLayers ( ) const

inline

Definition at line 78 of file FPGATrackSimPlaneMap.h.

78 { return m_nLogiLayers; }

◆ getNSections()

uint32_t FPGATrackSimPlaneMap::getNSections ( size_t logiLayer ) const

inline

Definition at line 80 of file FPGATrackSimPlaneMap.h.

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

◆ getOverrideLayer()

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

private

Definition at line 258 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 137 of file FPGATrackSimPlaneMap.h.

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

◆ isEC()

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

inline

Definition at line 111 of file FPGATrackSimPlaneMap.h.

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

◆ isPixel()

bool FPGATrackSimPlaneMap::isPixel ( int pl ) const

inline

Definition at line 110 of file FPGATrackSimPlaneMap.h.

110 { return m_dimension[pl] == 2; }

◆ isSCT()

bool FPGATrackSimPlaneMap::isSCT ( int pl ) const

inline

Definition at line 109 of file FPGATrackSimPlaneMap.h.

109 { return m_dimension[pl] == 1; }

◆ layerName()

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

Definition at line 266 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 237 of file FPGATrackSimPlaneMap.cxx.

 {
     if (hit.isMapped()) return;
     
     // 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?
     m_moduleRelabel->remap(hit);
  
     const LayerSection &pinfo = getLayerSection(hit.getDetType(), hit.getDetectorZone(), hit.getPhysLayer());
     hit.setSection(pinfo.section);
     hit.setLayer(pinfo.layer);
     if (!hit.isMapped()) // failsafe if for some reason someone calls this on a clustered hit again, or something
         hit.setHitType(HitType::mapped);
 }

◆ readLayers()

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

private

Definition at line 157 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 137 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 158 of file FPGATrackSimPlaneMap.h.

◆ m_coordOffset

std::vector<uint32_t> FPGATrackSimPlaneMap::m_coordOffset

private

Definition at line 154 of file FPGATrackSimPlaneMap.h.

◆ m_dimension

std::vector<uint32_t> FPGATrackSimPlaneMap::m_dimension

private

Definition at line 153 of file FPGATrackSimPlaneMap.h.

◆ m_diskIndex

std::vector<uint32_t> FPGATrackSimPlaneMap::m_diskIndex

private

Definition at line 170 of file FPGATrackSimPlaneMap.h.

◆ m_layerInfo

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

private

Definition at line 163 of file FPGATrackSimPlaneMap.h.

◆ m_layerOverrides

std::vector<int> FPGATrackSimPlaneMap::m_layerOverrides

private

Definition at line 166 of file FPGATrackSimPlaneMap.h.

◆ m_map

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

private

Definition at line 161 of file FPGATrackSimPlaneMap.h.

◆ m_moduleRelabel

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

private

Definition at line 173 of file FPGATrackSimPlaneMap.h.

◆ m_nCoords

uint32_t FPGATrackSimPlaneMap::m_nCoords = 0

private

Definition at line 151 of file FPGATrackSimPlaneMap.h.

◆ m_nDetLayers

uint32_t FPGATrackSimPlaneMap::m_nDetLayers = 0

private

Definition at line 150 of file FPGATrackSimPlaneMap.h.

◆ m_nLogiLayers

uint32_t FPGATrackSimPlaneMap::m_nLogiLayers = 0

private

Definition at line 149 of file FPGATrackSimPlaneMap.h.

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

Public Member Functions

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ FPGATrackSimPlaneMap()

Member Function Documentation

◆ allocateMap()

◆ getCoordLayer()

◆ getCoordOffset()

◆ getDetType()

◆ getDim()

◆ getLayer()

◆ getLayerInfo()

◆ getLayerSection() [1/2]

◆ getLayerSection() [2/2]

◆ getNCoords()

◆ getNDetLayers()

◆ getNLogiLayers()

◆ getNSections()

◆ getOverrideLayer()

◆ getSection()

◆ isEC()

◆ isPixel()

◆ isSCT()

◆ layerName()

◆ map()

◆ readLayers()

◆ seek()

Member Data Documentation

◆ m_coordLayer

◆ m_coordOffset

◆ m_dimension

◆ m_diskIndex

◆ m_layerInfo

◆ m_layerOverrides

◆ m_map

◆ m_moduleRelabel

◆ m_nCoords

◆ m_nDetLayers

◆ m_nLogiLayers