Gep::EratioMaker Class Reference

#include <EratioMaker.h>

Collaboration diagram for Gep::EratioMaker:

Public Member Functions
	EratioMaker (const GepCellMap &caloCellsMap, unsigned int etaWindowHalfSize=8, unsigned int phiWindowHalfSize=1)
const EratioObj	makeEratio (const ROOT::Math::PtEtaPhiEVector obj) const

Private Member Functions
std::vector< std::vector< GepCaloCell > >	makeWindow (const ROOT::Math::PtEtaPhiEVector &obj) const
std::pair< double, double >	findLocalMaxima (const std::vector< std::vector< GepCaloCell > > &window) const
double	computeEratio (double E1, double E2) const

Private Attributes
const GepCellMap &	m_caloCellsMap
const unsigned int	m_etaWindowHalfSize
const unsigned int	m_phiWindowHalfSize

Detailed Description

Definition at line 27 of file EratioMaker.h.

Constructor & Destructor Documentation

EratioMaker()

Gep::EratioMaker::EratioMaker	(	const GepCellMap &	caloCellsMap,
		unsigned int	etaWindowHalfSize = 8,
		unsigned int	phiWindowHalfSize = 1 )

Definition at line 15 of file EratioMaker.cxx.

    : m_caloCellsMap{caloCellsMap},
      m_etaWindowHalfSize{etaWindowHalfSize},
      m_phiWindowHalfSize{phiWindowHalfSize}
{
 
}

Member Function Documentation

computeEratio()

double Gep::EratioMaker::computeEratio ( double E1, double E2 ) const

private

Definition at line 158 of file EratioMaker.cxx.

{
    return (E1 + E2 > 0) ? (E1 - E2) / (E1 + E2) : -999;
}

findLocalMaxima()

std::pair< double, double > Gep::EratioMaker::findLocalMaxima ( const std::vector< std::vector< GepCaloCell > > & window ) const

private

Definition at line 103 of file EratioMaker.cxx.

{
 
    struct LocalMax {
        double energy;
        int i;
        int j;
    };
 
    std::vector<LocalMax> localMaxima;
 
    int rows = window.size();
    int cols = window[0].size();
 
    for (int i = 0; i < rows; ++i) {
        for (int j = 0; j < cols; ++j) {
            double current = window[i][j].e;
            bool isMax = true;
 
            for (int di = -1; di <= 1; ++di) {
                for (int dj = -1; dj <= 1; ++dj) {
                    if (di == 0 && dj == 0) continue;
                    int ni = i + di;
                    int nj = j + dj;
                    if (ni >= 0 && ni < rows && nj >= 0 && nj < cols) {
                        if (window[ni][nj].e >= current) {
                            isMax = false;
                            break;
                        }
                    }
                }
                if (!isMax) break;
            }
 
            if (isMax) localMaxima.push_back({current, i, j });
        }
    }
 
    if (localMaxima.empty()) return {0., 0.};
 
    // Sort local maximas in descendent fashion
    std::sort(localMaxima.begin(), localMaxima.end(),
                [](const LocalMax& a, const LocalMax& b){ return a.energy > b.energy; });
 
    double E1 = localMaxima[0].energy;
    double E2 = (localMaxima.size() > 1) ? localMaxima[1].energy : 0.0;
 
    return {E1, E2};
 
}

makeEratio()

const EratioObj Gep::EratioMaker::makeEratio

(

const ROOT::Math::PtEtaPhiEVector

obj

)

const

Definition at line 167 of file EratioMaker.cxx.

{
    auto window = makeWindow(seed);
    auto [E1, E2] = findLocalMaxima(window);
    double eratio = computeEratio(E1, E2);
 
    EratioObj obj;
    obj.seedEta = seed.eta();
    obj.seedPhi = seed.phi();
    obj.E1 = E1;
    obj.E2 = E2;
    obj.Eratio = eratio;
 
    return obj;
}

makeWindow()

std::vector< std::vector< GepCaloCell > > Gep::EratioMaker::makeWindow ( const ROOT::Math::PtEtaPhiEVector & obj ) const

private

Definition at line 31 of file EratioMaker.cxx.

{
 
    // Window dimensions set based on configured half sizes
    unsigned int etaWindowSize = 2 * m_etaWindowHalfSize + 1;
    unsigned int phiWindowSize = 2 * m_phiWindowHalfSize + 1;  
 
    std::vector<std::vector<GepCaloCell>> window(etaWindowSize, std::vector<GepCaloCell>(phiWindowSize));
 
    auto cellMap = m_caloCellsMap.getCellMap();
    
    // 1. Find cell that contains the seed 
    const GepCaloCell* centerCell = nullptr;
    unsigned int centerCellId = 0;
    bool exactMatchFound = false;
    
    for (auto it = cellMap->begin(); it != cellMap->end(); ++it) {
        const auto& cell = it->second;
        if (cell.sampling != 1) continue; // Only ECAL layer 1 
        
        // Verify if seed coordinates fall inside the limits of some cell
        if (cell.etaMin < seed.eta() && seed.eta() < cell.etaMax &&
            cell.phiMin < seed.phi() && seed.phi() < cell.phiMax) {
            centerCell = &cell;
            centerCellId = cell.id;
            exactMatchFound = true;
            break;
        }
    }
    
    // If not exact match, return empty window
    if (!exactMatchFound || !centerCell) {
        return window;
    }
 
    // 2. Construct 17x3 (rows x cols) window based on cell ID relations starting from seed
    // Rows: -8 to +8 from central cell in eta (17 total)
    // Columns: -1, 0, +1 from central cell in phi (3 total)
    
    // Construct each window row
    for (int etaOffset = -8; etaOffset <= 8; etaOffset++) {
        unsigned int windowRow = etaOffset + 8; // Convertir a índice 0..16
        
        // Get base ID for this row (central cell of the phi column)
        unsigned int baseId = static_cast<unsigned int>(static_cast<int>(centerCellId) + (etaOffset * 512));
        
        // Loop phi columns (-1, 0, +1)
        for (int phiOffset = -1; phiOffset <= 1; phiOffset++) {
            int windowCol = phiOffset + 1; 
            
            // Compute desired cell ID 
            unsigned int cellId = static_cast<unsigned int>(static_cast<int>(baseId) + phiOffset * 2);
            
            // Look for that cell within the map
            auto it = cellMap->find(cellId);
            if (it != cellMap->end()) {
                const auto& cell = it->second;
                if (cell.sampling == 1) {
                    window[windowRow][windowCol] = cell;
                }
            }
        }
    }
 
    return window;
}

Member Data Documentation

m_caloCellsMap

const GepCellMap& Gep::EratioMaker::m_caloCellsMap

private

Definition at line 45 of file EratioMaker.h.

m_etaWindowHalfSize

const unsigned int Gep::EratioMaker::m_etaWindowHalfSize

private

Definition at line 46 of file EratioMaker.h.

m_phiWindowHalfSize

const unsigned int Gep::EratioMaker::m_phiWindowHalfSize

private

Definition at line 47 of file EratioMaker.h.

---

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

• EratioMaker.h
• EratioMaker.cxx