Skip to content

[EMCAL-1157] Update clusterizer algorithm to use relative time window… #15822

New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

Sign up for GitHub

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Draft
mhemmer-cern wants to merge 1 commit into
base: dev
Choose a base branch
from
+7 −5
Draft

[EMCAL-1157] Update clusterizer algorithm to use relative time window… #15822

Show file tree
Hide file tree
Changes from all commits
Commits
Show all changes
1 commit
Select commit
File filter

Filter by extension

Filter by extension
Conversations
Failed to load comments.
Loading
Jump to
Jump to file
Failed to load files.
Loading
Diff view
Diff view
3 changes: 2 additions & 1 deletion Detectors/EMCAL/reconstruction/include/EMCALReconstruction/Clusterizer.h
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -143,7 +143,8 @@ class Clusterizer
/// \param[in,out] clusterInputs Cells/digits of prototype cluster
/// \param row Row number from neighbor search in recursion step
/// \param column Column number for neighbor search in recursion step
void getClusterFromNeighbours(std::vector<InputwithIndex>& clusterInputs, int row, int column);
/// \param seedTime Timestamp of the seed cell, used as fixed reference for the time cut
void getClusterFromNeighbours(std::vector<InputwithIndex>& clusterInputs, int row, int column, double seedTime);

/// \brief Get row (phi) and column (eta) of a cell/digit, values corresponding to topology
/// \param input Input object (cell/digit)
Expand Down
9 changes: 5 additions & 4 deletions Detectors/EMCAL/reconstruction/src/Clusterizer.cxx
View file
Open in desktop
Original file line number Diff line number Diff line change
Expand Up @@ -45,7 +45,7 @@ void Clusterizer<InputType>::initialize(double timeCut, double timeMin, double t

//____________________________________________________________________________
template <class InputType>
void Clusterizer<InputType>::getClusterFromNeighbours(std::vector<InputwithIndex>& clusterInputs, int row, int column)
void Clusterizer<InputType>::getClusterFromNeighbours(std::vector<InputwithIndex>& clusterInputs, int row, int column, double seedTime)
{
// Recursion 0, add seed cell/digit to cluster
if (!clusterInputs.size()) {
Expand All @@ -71,8 +71,8 @@ void Clusterizer<InputType>::getClusterFromNeighbours(std::vector<InputwithIndex
if (mDoEnergyGradientCut && (mInputMap[row + rowDiffs[dir]][column + colDiffs[dir]].mInput->getEnergy() > mInputMap[row][column].mInput->getEnergy() + mGradientCut)) {
continue;
}
if (not(TMath::Abs(mInputMap[row + rowDiffs[dir]][column + colDiffs[dir]].mInput->getTimeStamp() - mInputMap[row][column].mInput->getTimeStamp()) > mTimeCut)) {
getClusterFromNeighbours(clusterInputs, row + rowDiffs[dir], column + colDiffs[dir]);
if (not(TMath::Abs(mInputMap[row + rowDiffs[dir]][column + colDiffs[dir]].mInput->getTimeStamp() - seedTime) > mTimeCut)) {
getClusterFromNeighbours(clusterInputs, row + rowDiffs[dir], column + colDiffs[dir], seedTime);
// Add the cell/digit to the current cluster -- if we end up here, the selected cluster fulfills the condition
clusterInputs.emplace_back(mInputMap[row + rowDiffs[dir]][column + colDiffs[dir]]);
}
Expand Down Expand Up @@ -182,7 +182,8 @@ void Clusterizer<InputType>::findClusters(const gsl::span<InputType const>& inpu

// Seed is found, form cluster recursively
std::vector<InputwithIndex> clusterInputs;
getClusterFromNeighbours(clusterInputs, row, column);
double seedTime = mInputMap[row][column].mInput->getTimeStamp();
getClusterFromNeighbours(clusterInputs, row, column, seedTime);

// Add cells/digits for current cluster to cell/digit index vector
int inputIndexStart = mInputIndices.size();
Expand Down
Loading