#include <L1TopoLUT.h>

Collaboration diagram for LVL1MUCTPIPHASE1::L1TopoLUT:

Classes
struct	L1TopoLUTKey

struct	L1TopoLUTKeyHasher

Public Member Functions
bool	initializeLUT (const std::string &barrelFileName, const std::string &ecfFileName, const std::string &side0LUTFileName, const std::string &side1LUTFileName)

L1TopoCoordinates	getCoordinates (const unsigned short &side, const unsigned short &subsystem, const unsigned short &sectorID, const unsigned short &roi) const

unsigned short	getBarrelROI (unsigned short side, unsigned short sector, unsigned short ieta, unsigned short iphi) const

std::vector< std::string >	getErrors () const

float	getCompactedValue_eta (unsigned short subsystem, unsigned short side, unsigned short sectorID, unsigned short roi)

float	getCompactedValue_phi (unsigned short subsystem, unsigned short side, unsigned short sectorID, unsigned short roi)

Protected Member Functions
void	fillFromPtree (const boost::property_tree::ptree &node, std::map< unsigned short, std::vector< float >> &theLut) const

void	fillFromPtree (const boost::property_tree::ptree &node, std::map< unsigned short, std::map< unsigned short, std::pair< unsigned short, unsigned short >>> &theLut) const

bool	initializeCompactedLUTs (const std::string &side0LUTFileName, const std::string &side1LUTFileName)

bool	initializeLUT (const std::string &inFileName, const bool &isBarrel)

bool	initializeJSON (const std::string &inFileName, bool side)

bool	initializeJSONForSubsystem (pt::ptree &root, const std::string &nodeName, bool side, unsigned short subsystem)

Protected Attributes
std::map< unsigned short, std::vector< float > >	m_barrel_eta_lookup0

std::map< unsigned short, std::vector< float > >	m_barrel_eta_lookup1

std::map< unsigned short, std::vector< float > >	m_barrel_phi_lookup0

std::map< unsigned short, std::vector< float > >	m_barrel_phi_lookup1

std::map< unsigned short, std::vector< float > >	m_endcap_eta_lookup0

std::map< unsigned short, std::vector< float > >	m_endcap_eta_lookup1

std::map< unsigned short, std::vector< float > >	m_endcap_phi_lookup0

std::map< unsigned short, std::vector< float > >	m_endcap_phi_lookup1

std::map< unsigned short, std::vector< float > >	m_forward_eta_lookup0

std::map< unsigned short, std::vector< float > >	m_forward_eta_lookup1

std::map< unsigned short, std::vector< float > >	m_forward_phi_lookup0

std::map< unsigned short, std::vector< float > >	m_forward_phi_lookup1

std::map< unsigned short, std::map< unsigned short, std::pair< unsigned short, unsigned short > > >	m_barrel_encoding0

std::map< unsigned short, std::map< unsigned short, std::pair< unsigned short, unsigned short > > >	m_barrel_encoding1

std::map< unsigned short, unsigned short >	m_barrel_reverse_encoding0

std::map< unsigned short, unsigned short >	m_barrel_reverse_encoding1

std::unordered_map< L1TopoLUTKey, L1TopoCoordinates, L1TopoLUTKeyHasher >	m_encoding

std::vector< std::string >	m_errors

Detailed Description

Definition at line 40 of file L1TopoLUT.h.

Member Function Documentation

◆ fillFromPtree() [1/2]

void LVL1MUCTPIPHASE1::L1TopoLUT::fillFromPtree	(	const boost::property_tree::ptree &	node,
		std::map< unsigned short, std::map< unsigned short, std::pair< unsigned short, unsigned short >>> &	theLut
	)		const

protected

Definition at line 36 of file L1TopoLUT.cxx.

   {
     theLut.clear();
     for (const boost::property_tree::ptree::value_type& outer_elem : node){
       
       std::map<unsigned short,std::pair<unsigned short, unsigned short>> buff;
       for (const boost::property_tree::ptree::value_type& inner_elem : outer_elem.second) {
         std::vector<unsigned short> codeBuff;
         for (const boost::property_tree::ptree::value_type& code_elem : inner_elem.second) {
           codeBuff.push_back( code_elem.second.get_value<unsigned short>() );
         }
         if ( codeBuff.size() != 2 ) {
           //stricter error checking/handling comes later
           continue;
         }
         buff[std::stoi(inner_elem.first)] = { codeBuff[0], codeBuff[1] } ; 
       }
       //cannot get around stoi as outer layer may be non-continuous (some indices/keys are not present)
       theLut[std::stoi(outer_elem.first)] = buff;
     }
     return;
   }

◆ fillFromPtree() [2/2]

void LVL1MUCTPIPHASE1::L1TopoLUT::fillFromPtree	(	const boost::property_tree::ptree &	node,
		std::map< unsigned short, std::vector< float >> &	theLut
	)		const

protected

Definition at line 13 of file L1TopoLUT.cxx.

   {
     theLut.clear();
     size_t outerIndex = 0;
     for (const boost::property_tree::ptree::value_type& outer_elem : node){
       
       std::vector<float> buff(outer_elem.second.size()); 
       size_t index = 0;
       for (const boost::property_tree::ptree::value_type& inner_elem : outer_elem.second) {
         buff[index] = inner_elem.second.get_value<float>(); 
         ++index; //count index manually to avoid some stoi
       }
       //cannot get around stoi as outer layer may be non-continuous (some indices/keys are not present)
       if (outer_elem.first.length() > 0) { //string type key
         theLut[std::stoi(outer_elem.first)] = buff;
       } else {
         theLut[outerIndex] = buff;
         ++outerIndex; //ok to only count in this case, mixed keys/indices at same depth are not valid 
       }
     }
     return;
   }

◆ getBarrelROI()

unsigned short LVL1MUCTPIPHASE1::L1TopoLUT::getBarrelROI	(	unsigned short	side,
		unsigned short	sector,
		unsigned short	ieta,
		unsigned short	iphi
	)		const

Definition at line 232 of file L1TopoLUT.cxx.

                                                                                                                                    {
     unsigned short etaPhiSubWord = ((ieta&0xF) << 9) | ((sector&0x1F) << 3) | (iphi&0x7);
     unsigned short retval = 0;
     if (side) {
       retval = m_barrel_reverse_encoding0.at(etaPhiSubWord);
     } else {
       retval = m_barrel_reverse_encoding1.at(etaPhiSubWord);
     }
     return retval;
   }

◆ getCompactedValue_eta()

float LVL1MUCTPIPHASE1::L1TopoLUT::getCompactedValue_eta	(	unsigned short	subsystem,
		unsigned short	side,
		unsigned short	sectorID,
		unsigned short	roi
	)

Definition at line 244 of file L1TopoLUT.cxx.

                                                                                                                                    {
     //note: scheme used for "side" argument here is
     // 0 -> C side
     // 1 -> A side
     float retval = std::numeric_limits<float>::quiet_NaN();
     try {
       if (subsystem == 0) { //barrel
         //barrel is very irregular, so need to retrieve indices for MUCTPI->L1Topo encoding scheme first
         unsigned short lutCode = (side ? m_barrel_encoding0 : m_barrel_encoding1)[sectorID][roi].first;
         retval = (side ? m_barrel_eta_lookup0 : m_barrel_eta_lookup1)[sectorID][lutCode];
       } else if (subsystem == 1) { //endcap
         retval = (side ? m_endcap_eta_lookup0 : m_endcap_eta_lookup1)[roi >> 2][sectorID & 0x1];
       } else if (subsystem == 2) { //forward
         retval = (side ? m_forward_eta_lookup0 : m_forward_eta_lookup1)[sectorID & 0x1][roi >> 2];
       } else {
         std::stringstream ess;
         ess << "[TrigT1MuctpiPhase1/L1TopoLUT] Unknown muon subsystem value: " << subsystem;
         throw std::out_of_range(ess.str());
       }
     } catch (const std::exception& e) {
       std::stringstream ess;
       ess << "[TrigT1MuctpiPhase1/L1TopoLUT] Failed to look up eta value for subsystem=" << subsystem
           << ", side="<< (side?"A":"C")
           << ", sectorID="<< sectorID
           << ", roi="<< roi
           << "\n  (caught exception: " << e.what() <<")";
       throw std::out_of_range(ess.str());
     }
     return retval;
   }

◆ getCompactedValue_phi()

float LVL1MUCTPIPHASE1::L1TopoLUT::getCompactedValue_phi	(	unsigned short	subsystem,
		unsigned short	side,
		unsigned short	sectorID,
		unsigned short	roi
	)

Definition at line 275 of file L1TopoLUT.cxx.

                                                                                                                                    {
     //note: scheme used for "side" argument here is
     // 0 -> C side
     // 1 -> A side
     float retval = std::numeric_limits<float>::quiet_NaN();
     try {
       if (subsystem == 0) { //barrel
         //barrel is very irregular, so need to retrieve indices for MUCTPI->L1Topo encoding scheme first
         unsigned short lutCode = (side ? m_barrel_encoding0 : m_barrel_encoding1)[sectorID][roi].second;
         retval = (side ? m_barrel_phi_lookup0 : m_barrel_phi_lookup1)[sectorID][lutCode];
       } else if (subsystem == 1) { //endcap
         retval = (side ? m_endcap_phi_lookup0 : m_endcap_phi_lookup1)[sectorID][roi & 0x3];
       } else if (subsystem == 2) { //forward
         retval = (side ? m_forward_phi_lookup0 : m_forward_phi_lookup1)[roi & 0x3][sectorID];
       } else {
         std::stringstream ess;
         ess << "[TrigT1MuctpiPhase1/L1TopoLUT] Unknown muon subsystem value: " << subsystem;
         throw std::out_of_range(ess.str());
       }
     } catch (const std::exception& e) {
       std::stringstream ess;
       ess << "[TrigT1MuctpiPhase1/L1TopoLUT] Failed to look up phi value for subsystem=" << subsystem
           << ", side="<< (side?"A":"C")
           << ", sectorID="<< sectorID
           << ", roi="<< roi
           << "\n  (caught exception: " << e.what() <<")";
       throw std::out_of_range(ess.str());
     }
     return retval;
   }

◆ getCoordinates()

L1TopoCoordinates LVL1MUCTPIPHASE1::L1TopoLUT::getCoordinates	(	const unsigned short &	side,
		const unsigned short &	subsystem,
		const unsigned short &	sectorID,
		const unsigned short &	roi
	)		const

Definition at line 218 of file L1TopoLUT.cxx.

                                                            {
     L1TopoLUTKey key = {side, subsystem, sectorID, roi};
     auto itr = m_encoding.find(key);
     if (itr == m_encoding.end())
     {
       L1TopoCoordinates null;
       return null;
     }
     return itr->second;
   }

◆ getErrors()

std::vector<std::string> LVL1MUCTPIPHASE1::L1TopoLUT::getErrors ( ) const

inline

Definition at line 55 of file L1TopoLUT.h.

55 {return m_errors;}

◆ initializeCompactedLUTs()

bool LVL1MUCTPIPHASE1::L1TopoLUT::initializeCompactedLUTs	(	const std::string &	side0LUTFileName,
		const std::string &	side1LUTFileName
	)

protected

Definition at line 60 of file L1TopoLUT.cxx.

                                                          {    
     bool success = true;
     boost::property_tree::ptree root0;
     boost::property_tree::read_json(side0LUTFileName.c_str(), root0);
     boost::property_tree::ptree root1;
     boost::property_tree::read_json(side1LUTFileName.c_str(), root1);
  
     //Barrel encoding
     fillFromPtree(root0.get_child("encode_lookup_barrel") , m_barrel_encoding0); 
     fillFromPtree(root1.get_child("encode_lookup_barrel") , m_barrel_encoding1);
     
     //Barrel, A-side
     fillFromPtree(root0.get_child("lookup_barrel_eta") , m_barrel_eta_lookup0);
     fillFromPtree(root0.get_child("lookup_barrel_phi") , m_barrel_phi_lookup0);
     //Barrel, C-side
     fillFromPtree(root1.get_child("lookup_barrel_eta") , m_barrel_eta_lookup1);
     fillFromPtree(root1.get_child("lookup_barrel_phi") , m_barrel_phi_lookup1);
     //Endcap, A-side
     fillFromPtree(root0.get_child("lookup_endcap_eta") , m_endcap_eta_lookup0);
     fillFromPtree(root0.get_child("lookup_endcap_phi") , m_endcap_phi_lookup0);
     //Endcap, C-side
     fillFromPtree(root1.get_child("lookup_endcap_eta") , m_endcap_eta_lookup1);
     fillFromPtree(root1.get_child("lookup_endcap_phi") , m_endcap_phi_lookup1);
     //Forward, A-side
     fillFromPtree(root0.get_child("lookup_forward_eta") , m_forward_eta_lookup0);
     fillFromPtree(root0.get_child("lookup_forward_phi") , m_forward_phi_lookup0);
     //Forward, C-side
     fillFromPtree(root1.get_child("lookup_forward_eta") , m_forward_eta_lookup1);
     fillFromPtree(root1.get_child("lookup_forward_phi") , m_forward_phi_lookup1);
     
     
     if (m_barrel_encoding1.size() && m_barrel_encoding1.size() &&
         m_barrel_phi_lookup1.size() &&  m_barrel_phi_lookup0.size() && 
         m_endcap_phi_lookup1.size() &&  m_endcap_phi_lookup0.size() &&
         m_forward_phi_lookup1.size() &&  m_forward_phi_lookup0.size() ) {
       return success;
     } 
     return false;
   }

◆ initializeJSON()

bool LVL1MUCTPIPHASE1::L1TopoLUT::initializeJSON	(	const std::string &	inFileName,
		bool	side
	)

protected

◆ initializeJSONForSubsystem()

bool LVL1MUCTPIPHASE1::L1TopoLUT::initializeJSONForSubsystem	(	pt::ptree &	root,
		const std::string &	nodeName,
		bool	side,
		unsigned short	subsystem
	)

protected

◆ initializeLUT() [1/2]

bool LVL1MUCTPIPHASE1::L1TopoLUT::initializeLUT	(	const std::string &	barrelFileName,
		const std::string &	ecfFileName,
		const std::string &	side0LUTFileName,
		const std::string &	side1LUTFileName
	)

Definition at line 101 of file L1TopoLUT.cxx.

                                                    {
     m_encoding.clear();
     m_errors.clear();
     
     bool success = true;
  
     // read the compacted LUTs containing the simplified eta/phi values
     // used for encoding/decoding for/at L1Topo
     success = success && initializeCompactedLUTs(side0LUTFileName, side1LUTFileName);
  
     //read remaining info from the json files 
     //(the input labels are swapped between A-side and C-side)
     //success = success && initializeJSON(side1LUTFileName, 0);
     //success = success && initializeJSON(side0LUTFileName, 1);
     
     //read the barrel file
     success = success && initializeLUT(barrelFileName, true);
  
     //read the endcap+forward file
     success = success && initializeLUT(ecfFileName, false);
  
     return success;
   }

◆ initializeLUT() [2/2]

bool LVL1MUCTPIPHASE1::L1TopoLUT::initializeLUT	(	const std::string &	inFileName,
		const bool &	isBarrel
	)

protected

Definition at line 128 of file L1TopoLUT.cxx.

                                                                                  {
     //map between L1TopoLUTKey -> map between eta/phi -> index
     std::unordered_map<L1TopoLUTKey, std::map<float, unsigned short>, L1TopoLUTKeyHasher> sector_eta_indices, sector_phi_indices;
  
     std::ifstream inFile(inFileName.c_str());
     if (!inFile) return false;
     while (!inFile.eof() && inFile.good()) {
       unsigned short side;
       unsigned short subsystem=0;
       unsigned short sectorID;
       unsigned short roi;
       double eta_min, eta_max;
       double phi_min, phi_max;
  
       inFile >> side;
       if (inFile.eof()) { break; }
  
       if (!isBarrel) {
           inFile >> subsystem;      
           //in the EC+FW file, EC=0 and FW=1, so increment to distinguish from barrel
           subsystem++;
           side = !side; //EC+FW file uses inverted side numbers w.r.t. Barrel...
       } 
       //at this point the "side" convention is
       // 0 -> C side
       // 1 -> A side
       // the LUT json file naming, however, if inverse (follows EC+FW file...)
       
       inFile >> sectorID;
       inFile >> roi;
       inFile >> eta_min;
       inFile >> eta_max;
       inFile >> phi_min;
       inFile >> phi_max;
  
       double eta = getCompactedValue_eta(subsystem, side, sectorID, roi);
       double phi = getCompactedValue_phi(subsystem, side, sectorID, roi);
       
       unsigned short ieta = 0;
       unsigned short iphi = 0;
       if (subsystem == 0) { //barrel
         //note: interpretation of values 0/1 for side flipped on purpose here (different conventions)!
         ieta = (side ? m_barrel_encoding0 : m_barrel_encoding1)[sectorID][roi].first;
         iphi = (side ? m_barrel_encoding0 : m_barrel_encoding1)[sectorID][roi].second;
       } else if (subsystem == 1) { //endcap
         ieta  = (roi>>2); //eta part of roi
         ieta |= 1 << 6; //detector code
         iphi  = roi & 0x3; //phi part of roi
         iphi |= (sectorID & 0x3F) << 2; 
       } else if (subsystem == 2) { //forward
         ieta  = (roi>>2); //eta part of roi
         ieta |= 1 << 5; //detector code
         iphi  = (roi & 0x3) << 1; //phi part of roi
         iphi |= (sectorID & 0x1F) << 3; 
       } 
       //truncate to 8bit values each:
       ieta &= 0xFF;
       iphi &= 0xFF;
      
       if (subsystem == 0) {
         //populate reverse map for barrel
         unsigned short etaPhiSubWord = ((ieta&0xF) << 9) | ((sectorID&0x1F) << 3) | (iphi&0x7);
         if (side) {
           m_barrel_reverse_encoding0[etaPhiSubWord] = roi;
         } else {
           m_barrel_reverse_encoding1[etaPhiSubWord] = roi;
         }
       }
      
       //reminder in case of debugging needs:
       //compacted values should be close to these averages
       //(albeit some are different enough to cause mismatches in Topo!)
       //float avgEta = (eta_min+eta_max)/2.;
       //float avgPhi = (phi_min+phi_max)/2.;
      
      
       L1TopoLUTKey key = {side, subsystem, sectorID, roi};
       L1TopoCoordinates value = {eta, phi, eta_min, eta_max, phi_min, phi_max, ieta, iphi};
       
       if (m_encoding.find(key) != m_encoding.end())
       {
           m_errors.push_back("Duplicate key found in L1TopoLUT: "+key.info());
           return false;
       }
       m_encoding[key] = value;
     }
     return true;
   }

Member Data Documentation

◆ m_barrel_encoding0

std::map<unsigned short,std::map<unsigned short, std::pair<unsigned short, unsigned short> > > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_encoding0

protected

Definition at line 129 of file L1TopoLUT.h.

◆ m_barrel_encoding1

std::map<unsigned short,std::map<unsigned short, std::pair<unsigned short, unsigned short> > > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_encoding1

protected

Definition at line 130 of file L1TopoLUT.h.

◆ m_barrel_eta_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_eta_lookup0

protected

Definition at line 116 of file L1TopoLUT.h.

◆ m_barrel_eta_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_eta_lookup1

protected

Definition at line 117 of file L1TopoLUT.h.

◆ m_barrel_phi_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_phi_lookup0

protected

Definition at line 118 of file L1TopoLUT.h.

◆ m_barrel_phi_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_phi_lookup1

protected

Definition at line 119 of file L1TopoLUT.h.

◆ m_barrel_reverse_encoding0

std::map<unsigned short, unsigned short> LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_reverse_encoding0

protected

Definition at line 133 of file L1TopoLUT.h.

◆ m_barrel_reverse_encoding1

std::map<unsigned short, unsigned short> LVL1MUCTPIPHASE1::L1TopoLUT::m_barrel_reverse_encoding1

protected

Definition at line 134 of file L1TopoLUT.h.

◆ m_encoding

std::unordered_map<L1TopoLUTKey, L1TopoCoordinates, L1TopoLUTKeyHasher> LVL1MUCTPIPHASE1::L1TopoLUT::m_encoding

protected

Definition at line 136 of file L1TopoLUT.h.

◆ m_endcap_eta_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_endcap_eta_lookup0

protected

Definition at line 120 of file L1TopoLUT.h.

◆ m_endcap_eta_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_endcap_eta_lookup1

protected

Definition at line 121 of file L1TopoLUT.h.

◆ m_endcap_phi_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_endcap_phi_lookup0

protected

Definition at line 122 of file L1TopoLUT.h.

◆ m_endcap_phi_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_endcap_phi_lookup1

protected

Definition at line 123 of file L1TopoLUT.h.

◆ m_errors

std::vector<std::string> LVL1MUCTPIPHASE1::L1TopoLUT::m_errors

protected

Definition at line 137 of file L1TopoLUT.h.

◆ m_forward_eta_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_forward_eta_lookup0

protected

Definition at line 124 of file L1TopoLUT.h.

◆ m_forward_eta_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_forward_eta_lookup1

protected

Definition at line 125 of file L1TopoLUT.h.

◆ m_forward_phi_lookup0

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_forward_phi_lookup0

protected

Definition at line 126 of file L1TopoLUT.h.

◆ m_forward_phi_lookup1

std::map<unsigned short,std::vector<float> > LVL1MUCTPIPHASE1::L1TopoLUT::m_forward_phi_lookup1

protected

Definition at line 127 of file L1TopoLUT.h.

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

Classes

Public Member Functions

Protected Member Functions

Protected Attributes

Detailed Description

Member Function Documentation

◆ fillFromPtree() [1/2]

◆ fillFromPtree() [2/2]

◆ getBarrelROI()

◆ getCompactedValue_eta()

◆ getCompactedValue_phi()

◆ getCoordinates()

◆ getErrors()

◆ initializeCompactedLUTs()

◆ initializeJSON()

◆ initializeJSONForSubsystem()

◆ initializeLUT() [1/2]

◆ initializeLUT() [2/2]

Member Data Documentation

◆ m_barrel_encoding0

◆ m_barrel_encoding1

◆ m_barrel_eta_lookup0

◆ m_barrel_eta_lookup1

◆ m_barrel_phi_lookup0

◆ m_barrel_phi_lookup1

◆ m_barrel_reverse_encoding0

◆ m_barrel_reverse_encoding1

◆ m_encoding

◆ m_endcap_eta_lookup0

◆ m_endcap_eta_lookup1

◆ m_endcap_phi_lookup0

◆ m_endcap_phi_lookup1

◆ m_errors

◆ m_forward_eta_lookup0

◆ m_forward_eta_lookup1

◆ m_forward_phi_lookup0

◆ m_forward_phi_lookup1