Cluster splitter algorithm to be run on GPUs, using the same cellular automaton-based approach as Topo-Automaton Clustering. More...

#include <TopoAutomatonSplitting.h>

Inheritance diagram for TopoAutomatonSplitting:

Collaboration diagram for TopoAutomatonSplitting:

Public Member Functions
	TopoAutomatonSplitting (const std::string &type, const std::string &name, const IInterface *parent)

virtual StatusCode	initialize () override

virtual StatusCode	initialize_non_CUDA () override
	Initialization that does not invoke CUDA functions. More...

virtual StatusCode	initialize_CUDA () override
	Initialization that invokes CUDA functions. More...

virtual StatusCode	execute (const EventContext &ctx, const CaloRecGPU::ConstantDataHolder &constant_data, CaloRecGPU::EventDataHolder &event_data, void *temporary_buffer) const override

virtual StatusCode	finalize () override

virtual	~TopoAutomatonSplitting ()=default

void	handle (const Incident &incident) override

Protected Member Functions
void	record_times (const size_t event_num, const std::vector< size_t > &times) const

template<class ... Args>
void	record_times (const size_t event_num, const size_t &value) const

template<class ... Args>
void	record_times (const size_t event_num, const size_t &value, Args &&... args) const

void	print_times (const std::string &header, const size_t time_size) const

Protected Attributes
std::shared_mutex	m_timeMutex
	Mutex that is locked when recording times. More...

std::vector< size_t > m_times	ATLAS_THREAD_SAFE
	Vector to hold execution times to be recorded if necessary. More...

std::vector< size_t > m_eventNumbers	ATLAS_THREAD_SAFE
	Vector to hold the event numbers to be recorded if necessary. More...

Gaudi::Property< bool >	m_measureTimes
	If `true`, times are recorded to the file given by `m_timeFileName`. More...

Gaudi::Property< std::string >	m_timeFileName
	File to which times should be saved. More...

Private Member Functions
void	record_times_helper (const size_t) const

template<class Arg >
void	record_times_helper (const size_t index, Arg &&arg) const

template<class ... Args>
void	record_times_helper (size_t index, Args &&... args) const

Private Attributes
Gaudi::Property< std::vector< std::string > >	m_samplingNames {this, "SamplingNames", {}, "Name(s) of Calorimeter Samplings to consider for local maxima"}
	vector of names of the calorimeter samplings to consider for seeds. More...

Gaudi::Property< std::vector< std::string > >	m_secondarySamplingNames {this, "SecondarySamplingNames", {}, "Name(s) of secondary Calorimeter Samplings to consider for local maxima"}
	vector of names of the secondary calorimeter samplings to consider. More...

Gaudi::Property< int >	m_nCells {this, "NumberOfCellsCut", 4, "Local maxima need at least this number of neighbors to become seeds. Must be a non-negative number (obviously)."}
	local maxima need at least this number of neighbors to become seeds More...

Gaudi::Property< float >	m_minEnergy {this, "EnergyCut", 500 * CLHEP::MeV, "Minimal energy for a local max"}
	local maxima need at least this energy content More...

Gaudi::Property< bool >	m_shareBorderCells {this, "ShareBorderCells", false, "Whether or not to share cells at the boundary between two clusters"}
	share cells at the border between two local maxima More...

Gaudi::Property< float >	m_emShowerScale {this, "EMShowerScale", 5 * CLHEP::cm, "Typical EM shower distance for which the energy density should drop to 1/e"}
	typical EM shower scale to use for distance criteria in shared cells More...

Gaudi::Property< bool >	m_absOpt {this, "WeightingOfNegClusters", false, "Should absolute value be used to identify potential seed cells"}
	if set to true, splitter only looks at absolute value of Energy in order to identify potential seed cells More...

Gaudi::Property< bool >	m_treatL1PredictedCellsAsGood {this, "TreatL1PredictedCellsAsGood", true, "Treat bad cells with dead OTX if predicted from L1 as good"}
	if set to true treat cells with a dead OTX which can be predicted by L1 trigger info as good instead of bad cells More...

Gaudi::Property< std::string >	m_neighborOptionString
	type of neighbor relations to use. More...

Gaudi::Property< bool >	m_restrictHECIWandFCalNeighbors
	if set to true limit the neighbors in HEC IW and FCal2&3. More...

Gaudi::Property< bool >	m_restrictPSNeighbors
	if set to true limit the neighbors in presampler Barrel and Endcap. More...

TASplitting::TASOptionsHolder	m_options
	Options for the algorithm, held in a GPU-friendly way. More...

ServiceHandle< IGPUKernelSizeOptimizerSvc >	m_kernelSizeOptimizer { this, "KernelSizeOptimizer", "GPUKernelSizeOptimizerSvc", "CUDA kernel size optimization service." }
	Handle to the CUDA kernel block and grid size optimization service. More...

Detailed Description

Cluster splitter algorithm to be run on GPUs, using the same cellular automaton-based approach as Topo-Automaton Clustering.

Author: Nuno Fernandes nuno..nosp@m.dos..nosp@m.santo.nosp@m.s.fe.nosp@m.rnand.nosp@m.es@c.nosp@m.ern.c.nosp@m.h

Date: 25 August 2022

Definition at line 31 of file TopoAutomatonSplitting.h.

Constructor & Destructor Documentation

◆ TopoAutomatonSplitting()

TopoAutomatonSplitting::TopoAutomatonSplitting	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 24 of file TopoAutomatonSplitting.cxx.

                                                                                                                        :
   base_class(type, name, parent),
   CaloGPUTimed(this)
 {
 }

◆ ~TopoAutomatonSplitting()

virtual TopoAutomatonSplitting::~TopoAutomatonSplitting ( )

virtualdefault

Member Function Documentation

◆ execute()

StatusCode TopoAutomatonSplitting::execute	(	const EventContext &	ctx,
		const CaloRecGPU::ConstantDataHolder &	constant_data,
		CaloRecGPU::EventDataHolder &	event_data,
		void *	temporary_buffer
	)		const

overridevirtual

Definition at line 208 of file TopoAutomatonSplitting.cxx.

 {
  
   using clock_type = boost::chrono::thread_clock;
   auto time_cast = [](const auto & before, const auto & after)
   {
     return boost::chrono::duration_cast<boost::chrono::microseconds>(after - before).count();
   };
  
   static_assert(sizeof(TopoAutomatonSplittingTemporaries) <= sizeof(ClusterMomentsArr), "We store the temporaries in the cluster moments, so the sizes must be compatible!");
   
   const auto start = clock_type::now();
   const auto preprocessing_end = clock_type::now();
   
   fillNeighbours(event_data, constant_data, m_options, *(m_kernelSizeOptimizer.get()), m_measureTimes);
   
   const auto after_neighs = clock_type::now();
   
   findLocalMaxima(event_data, constant_data, m_options, *(m_kernelSizeOptimizer.get()), m_measureTimes);
   
   const auto after_maxima = clock_type::now();
   
   excludeSecondaryMaxima(event_data, constant_data, m_options, *(m_kernelSizeOptimizer.get()), m_measureTimes);
   
   const auto after_secondary_maxima = clock_type::now();
   
   splitClusterGrowing(event_data, constant_data, m_options, *(m_kernelSizeOptimizer.get()), m_measureTimes);
   
   const auto after_growing = clock_type::now();
   
   cellWeightingAndFinalization(event_data, constant_data, m_options, *(m_kernelSizeOptimizer.get()), m_measureTimes);
  
   const auto end = clock_type::now();
   
   if (m_measureTimes)
     {
       record_times(ctx.evt(),
                    time_cast(start, preprocessing_end),
                    time_cast(preprocessing_end, after_neighs),
                    time_cast(after_neighs, after_maxima),
                    time_cast(after_maxima, after_secondary_maxima),
                    time_cast(after_secondary_maxima, after_growing),
                    time_cast(after_growing, end)
                   );
     }
  
   return StatusCode::SUCCESS;
  
  
 }

◆ finalize()

StatusCode TopoAutomatonSplitting::finalize ( )

overridevirtual

Definition at line 261 of file TopoAutomatonSplitting.cxx.

 {
   if (m_measureTimes)
     {
       print_times("Preprocessing Fill_List_of_Intra-Cluster_Neighbours Find_Local_Maxima Find_Secondary_Maxima Splitter_Tag_Propagation Cell_Weighting_And_Finalization", 6);
     }
   return StatusCode::SUCCESS;
 }

◆ handle()

void CaloGPUCUDAInitialization::handle ( const Incident & incident )

inlineoverrideinherited

Definition at line 66 of file CaloGPUCUDAInitialization.h.

   {
     const bool is_multiprocess = (Gaudi::Concurrency::ConcurrencyFlags::numProcs() > 0);
     if (is_multiprocess && incident.type() == AthenaInterprocess::UpdateAfterFork::type())
       {
         if (!this->initialize_CUDA().isSuccess())
         {
           throw GaudiException("Failed to perform the CUDA initialization!",
                                "CaloGPUCUDAInitialization::handle",
                                StatusCode::FAILURE);
         }
       }
   }

◆ initialize()

virtual StatusCode TopoAutomatonSplitting::initialize ( )

inlineoverridevirtual

Reimplemented from CaloGPUCUDAInitialization.

Definition at line 38 of file TopoAutomatonSplitting.h.

   {
     return CaloGPUCUDAInitialization::initialize();
   }

◆ initialize_CUDA()

StatusCode TopoAutomatonSplitting::initialize_CUDA ( )

overridevirtual

Initialization that invokes CUDA functions.

Reimplemented from CaloGPUCUDAInitialization.

Definition at line 200 of file TopoAutomatonSplitting.cxx.

 {
   m_options.sendToGPU();
   register_kernels( *(m_kernelSizeOptimizer.get()) );
   
   return StatusCode::SUCCESS;
 }

◆ initialize_non_CUDA()

StatusCode TopoAutomatonSplitting::initialize_non_CUDA ( )

overridevirtual

Initialization that does not invoke CUDA functions.

Reimplemented from CaloGPUCUDAInitialization.

Definition at line 30 of file TopoAutomatonSplitting.cxx.

 {
   m_options.allocate();
  
   using PackType = decltype(m_options.m_options->valid_sampling_primary);
  
   static_assert(CaloCell_ID::getNumberOfSamplings() <= sizeof(PackType) * CHAR_BIT,  "We are assuming that we have fewer samplings that bits per int!");
  
   auto get_option_from_string = [](const std::string & str, bool & failed)
   {
     failed = false;
     CRGPU_RECURSIVE_MACRO(
             CRGPU_CHEAP_STRING_TO_ENUM( str, CaloCell_ID,
                                         PreSamplerB,
                                         EMB1,
                                         EMB2,
                                         EMB3,
                                         PreSamplerE,
                                         EME1,
                                         EME2,
                                         EME3,
                                         HEC0,
                                         HEC1,
                                         HEC2,
                                         HEC3,
                                         TileBar0,
                                         TileBar1,
                                         TileBar2,
                                         TileGap1,
                                         TileGap2,
                                         TileGap3,
                                         TileExt0,
                                         TileExt1,
                                         TileExt2,
                                         FCAL0,
                                         FCAL1,
                                         FCAL2
                                       )
     )
     else
       {
         failed = true;
         return CaloCell_ID::Unknown;
       }
   };
  
  
   auto process_sampling = [&get_option_from_string](const std::vector<std::string> & sampling_names, std::string & invalid_names, PackType & sampling_option)
   {
     sampling_option = 0;
     for (const std::string & samp_name : sampling_names)
       {
         bool failed = false;
         const PackType sampling = (PackType) get_option_from_string(samp_name, failed);
  
         if (failed)
           {
             if (invalid_names.size() == 0)
               {
                 invalid_names = "'" + samp_name + "'";
               }
             else
               {
                 invalid_names += ", '" + samp_name + "'";
               }
           }
         else
           {
             sampling_option |= ((PackType) 1) << sampling;
           }
       }
   };
  
   std::string invalid_names;
  
   process_sampling(m_samplingNames, invalid_names, m_options.m_options->valid_sampling_primary);
  
   if (invalid_names.size() > 0)
     {
       ATH_MSG_ERROR( "Calorimeter samplings " << invalid_names
                      << " are not a valid Calorimeter sampling name and will be ignored! "
                      << "Valid names are: "
                      << "PreSamplerB, EMB1, EMB2, EMB3, "
                      << "PreSamplerE, EME1, EME2, EME3, "
                      << "HEC0, HEC1, HEC2, HEC3, "
                      << "TileBar0, TileBar1, TileBar2, "
                      << "TileGap1, TileGap2, TileGap3, "
                      << "TileExt0, TileExt1, TileExt2, "
                      << "FCAL0, FCAL1, FCAL2."  );
     }
  
   invalid_names.clear();
  
   process_sampling(m_secondarySamplingNames, invalid_names, m_options.m_options->valid_sampling_secondary);
  
   if (invalid_names.size() > 0)
     {
       ATH_MSG_ERROR( "Calorimeter samplings " << invalid_names
                      << " are not a valid Calorimeter sampling name and will be ignored! "
                      << "Valid names are: "
                      << "PreSamplerB, EMB1, EMB2, EMB3, "
                      << "PreSamplerE, EME1, EME2, EME3, "
                      << "HEC0, HEC1, HEC2, HEC3, "
                      << "TileBar0, TileBar1, TileBar2, "
                      << "TileGap1, TileGap2, TileGap3, "
                      << "TileExt0, TileExt1, TileExt2, "
                      << "FCAL0, FCAL1, FCAL2."  );
     }
  
   auto get_neighbour_option_from_string = [](const std::string & str, bool & failed)
   {
     failed = false;
     CRGPU_RECURSIVE_MACRO(
             CRGPU_CHEAP_STRING_TO_ENUM( str, LArNeighbours,
                                         prevInPhi,
                                         nextInPhi,
                                         prevInEta,
                                         nextInEta,
                                         faces2D,
                                         corners2D,
                                         all2D,
                                         prevInSamp,
                                         nextInSamp,
                                         upAndDown,
                                         prevSubDet,
                                         nextSubDet,
                                         all3D,
                                         corners3D,
                                         all3DwithCorners,
                                         prevSuperCalo,
                                         nextSuperCalo,
                                         super3D
                                       )
     )
     //I know Topological Clustering only supports a subset of those,
     //but this is supposed to be a general data exporting tool...
     else
       {
         failed = true;
         return LArNeighbours::super3D;
       }
   };
  
   bool neigh_failed = false;
   m_options.m_options->neighbour_options = (unsigned int) get_neighbour_option_from_string(m_neighborOptionString, neigh_failed);
  
   if (neigh_failed)
     {
       ATH_MSG_ERROR("Invalid Neighbour Option: " << m_neighborOptionString);
     }
  
   //We must repeat this printing part because ATH_MSG_ERROR
   //is a macro that apparently calls a this->msg(...) function.
   //Of course it won't work within a lambda...
  
   m_options.m_options->min_num_cells = m_nCells;
   m_options.m_options->min_maximum_energy = m_minEnergy;
   m_options.m_options->EM_shower_scale = m_emShowerScale;
   m_options.m_options->share_border_cells = m_shareBorderCells;
   m_options.m_options->use_absolute_energy = m_absOpt;
   m_options.m_options->treat_L1_predicted_as_good = m_treatL1PredictedCellsAsGood;
  
   m_options.m_options->limit_HECIW_and_FCal_neighs = m_restrictHECIWandFCalNeighbors;
   m_options.m_options->limit_PS_neighs = m_restrictPSNeighbors;
   
   ATH_CHECK( m_kernelSizeOptimizer.retrieve() );
   
   return StatusCode::SUCCESS;
 }

◆ print_times()

void CaloGPUTimed::print_times	(	const std::string &	header,
		const size_t	time_size
	)		const

inlineprotectedinherited

Definition at line 143 of file CaloGPUTimed.h.

   {
     std::shared_lock<std::shared_mutex> lock(m_timeMutex);
     
     if (m_timeFileName.size() == 0)
       {
         return;
       }
     
     std::vector<size_t> indices(m_eventNumbers.size());
     
     std::iota(indices.begin(), indices.end(), 0);
     std::sort(indices.begin(), indices.end(), [&](size_t a, size_t b)
     {
       return m_eventNumbers[a] < m_eventNumbers[b];
     }
              );
     std::ofstream out(m_timeFileName);
  
     out << "Event_Number Total " << header << "\n";
  
     for (const size_t idx : indices)
       {
         out << m_eventNumbers[idx] << " ";
  
         size_t total = 0;
  
         for (size_t i = 0; i < time_size; ++i)
           {
             total += m_times[idx * time_size + i];
           }
  
         out << total << " ";
  
         for (size_t i = 0; i < time_size; ++i)
           {
             out << m_times[idx * time_size + i] << " ";
           }
         out << "\n";
       }
  
     out << std::endl;
  
     out.close();
   }

◆ record_times() [1/3]

template<class ... Args>

void CaloGPUTimed::record_times	(	const size_t	event_num,
		const size_t &	value
	)		const

inlineprotectedinherited

Definition at line 105 of file CaloGPUTimed.h.

   {
     const size_t time_size = 1;
  
     size_t old_size;
  
     {
       std::unique_lock<std::shared_mutex> lock(m_timeMutex);
       old_size = m_times.size();
       m_times.resize(old_size + time_size);
       m_eventNumbers.push_back(event_num);
     }
     {
       std::shared_lock<std::shared_mutex> lock(m_timeMutex);
       record_times_helper(old_size, value);
     }
   }

◆ record_times() [2/3]

template<class ... Args>

void CaloGPUTimed::record_times	(	const size_t	event_num,
		const size_t &	value,
		Args &&...	args
	)		const

inlineprotectedinherited

Definition at line 124 of file CaloGPUTimed.h.

   {
     const size_t time_size = sizeof...(args) + 1;
  
     size_t old_size;
  
     {
       std::unique_lock<std::shared_mutex> lock(m_timeMutex);
       old_size = m_times.size();
       m_times.resize(old_size + time_size);
       m_eventNumbers.push_back(event_num);
     }
     {
       std::shared_lock<std::shared_mutex> lock(m_timeMutex);
       record_times_helper(old_size, value, std::forward<Args>(args)...);
     }
  
   }

◆ record_times() [3/3]

void CaloGPUTimed::record_times	(	const size_t	event_num,
		const std::vector< size_t > &	times
	)		const

inlineprotectedinherited

Definition at line 86 of file CaloGPUTimed.h.

   {
     size_t old_size;
     {
       std::unique_lock<std::shared_mutex> lock(m_timeMutex);
       old_size = m_times.size();
       m_times.resize(old_size + times.size());
       m_eventNumbers.push_back(event_num);
     }
     {
       std::shared_lock<std::shared_mutex> lock(m_timeMutex);
       for (size_t i = 0; i < times.size(); ++i)
         {
           m_times[old_size + i] = times[i];
         }
     }
   }

◆ record_times_helper() [1/3]

template<class Arg >

void CaloGPUTimed::record_times_helper	(	const size_t	index,
		Arg &&	arg
	)		const

inlineprivateinherited

Definition at line 70 of file CaloGPUTimed.h.

   {
     // coverity[missing_lock]
     m_times[index] = std::forward<Arg>(arg);
     
     //This is called within a function that holds the lock itself.
   }

◆ record_times_helper() [2/3]

void CaloGPUTimed::record_times_helper ( const size_t ) const

inlineprivateinherited

Definition at line 64 of file CaloGPUTimed.h.

   {
     //Do nothing
   }

◆ record_times_helper() [3/3]

template<class ... Args>

void CaloGPUTimed::record_times_helper	(	size_t	index,
		Args &&...	args
	)		const

inlineprivateinherited

Definition at line 79 of file CaloGPUTimed.h.

   {
     (record_times_helper(index++, std::forward<Args>(args)), ...);
   }

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

std::vector<size_t> m_times CaloGPUTimed::ATLAS_THREAD_SAFE

mutableprotectedinherited

Vector to hold execution times to be recorded if necessary.

Definition at line 35 of file CaloGPUTimed.h.

◆ ATLAS_THREAD_SAFE [2/2]

std::vector<size_t> m_eventNumbers CaloGPUTimed::ATLAS_THREAD_SAFE

mutableprotectedinherited

Vector to hold the event numbers to be recorded if necessary.

Definition at line 40 of file CaloGPUTimed.h.

◆ m_absOpt

Gaudi::Property<bool> TopoAutomatonSplitting::m_absOpt {this, "WeightingOfNegClusters", false, "Should absolute value be used to identify potential seed cells"}

private

if set to true, splitter only looks at absolute value of Energy in order to identify potential seed cells

Definition at line 134 of file TopoAutomatonSplitting.h.

◆ m_emShowerScale

Gaudi::Property<float> TopoAutomatonSplitting::m_emShowerScale {this, "EMShowerScale", 5 * CLHEP::cm, "Typical EM shower distance for which the energy density should drop to 1/e"}

private

typical EM shower scale to use for distance criteria in shared cells

a shared cell is included in both clusters neighboring the cell with weights depending on the cluster energies and the distance of the shared cell to the cluster centroids. The distance is measured in units of this property to roughly describe the exponential slope of the energy density distribution for em showers. The exact choice of this property is not critical but should roughly match the Moliere radius in the LArEM since here the sharing of cells has the biggest use case.

Definition at line 129 of file TopoAutomatonSplitting.h.

◆ m_kernelSizeOptimizer

ServiceHandle<IGPUKernelSizeOptimizerSvc> TopoAutomatonSplitting::m_kernelSizeOptimizer { this, "KernelSizeOptimizer", "GPUKernelSizeOptimizerSvc", "CUDA kernel size optimization service." }

private

Handle to the CUDA kernel block and grid size optimization service.

Definition at line 198 of file TopoAutomatonSplitting.h.

◆ m_measureTimes

Gaudi::Property<bool> CaloGPUTimed::m_measureTimes

protectedinherited

If true, times are recorded to the file given by m_timeFileName.

Defaults to false.

Definition at line 46 of file CaloGPUTimed.h.

◆ m_minEnergy

Gaudi::Property<float> TopoAutomatonSplitting::m_minEnergy {this, "EnergyCut", 500 * CLHEP::MeV, "Minimal energy for a local max"}

private

local maxima need at least this energy content

potential seed cells have to pass this cut on the energy content.

Definition at line 93 of file TopoAutomatonSplitting.h.

◆ m_nCells

Gaudi::Property<int> TopoAutomatonSplitting::m_nCells {this, "NumberOfCellsCut", 4, "Local maxima need at least this number of neighbors to become seeds. Must be a non-negative number (obviously)."}

private

local maxima need at least this number of neighbors to become seeds

each cell above the energy cut having at least this many neighbors in the parent cluster and only neighbors with smaller energy seed a split cluster.

Definition at line 86 of file TopoAutomatonSplitting.h.

◆ m_neighborOptionString

Gaudi::Property<std::string> TopoAutomatonSplitting::m_neighborOptionString

private

Initial value:

{this, "NeighborOption", "super3D",

"Neighbor option to be used for cell neighborhood relations"}

type of neighbor relations to use.

The CaloIdentifier package defines different types of neighbors for the calorimeter cells. Currently supported neighbor relations for topological clustering are:

"all2D" for all cells in the same layer (sampling or module) of one calorimeter subsystem. Note that endcap and barrel will be unconnected in this case even for the LAREM.

"all3D" for all cells in the same calorimeter. This means all the "all2D" neighbors for each cell plus the cells in adjacent samplings overlapping at least partially in \(\eta\) and \(\phi\) with the cell. Note that endcap and barrel will be connected in this case for the LAREM.

"super3D" for all cells. This means all the "all3D" neighbors for each cell plus the cells in adjacent samplings from other subsystems overlapping at least partially in \(\eta\) and \(\phi\) with the cell. All calorimeters are connected in this case.

The default setting is "super3D".

Definition at line 165 of file TopoAutomatonSplitting.h.

◆ m_options

TASplitting::TASOptionsHolder TopoAutomatonSplitting::m_options

private

Options for the algorithm, held in a GPU-friendly way.

Definition at line 195 of file TopoAutomatonSplitting.h.

◆ m_restrictHECIWandFCalNeighbors

Gaudi::Property<bool> TopoAutomatonSplitting::m_restrictHECIWandFCalNeighbors

private

Initial value:

{this, "RestrictHECIWandFCalNeighbors",

false, "Limit the neighbors in HEC IW and FCal2&3"}

if set to true limit the neighbors in HEC IW and FCal2&3.

The cells in HEC IW and FCal2&3 get very large in terms of eta and phi. Since this might pose problems on certain jet algorithms one might need to avoid expansion in eta and phi for those cells. If this property is set to true the 2d neighbors of these cells are not used - only the next sampling neighbors are probed.

Definition at line 177 of file TopoAutomatonSplitting.h.

◆ m_restrictPSNeighbors

Gaudi::Property<bool> TopoAutomatonSplitting::m_restrictPSNeighbors

private

Initial value:

{this, "RestrictPSNeighbors",

false, "Limit the neighbors in presampler Barrel and Endcap"}

if set to true limit the neighbors in presampler Barrel and Endcap.

The presampler cells add a lot of PileUp in the Hilum samples. With this option set to true the presampler cells do not expand the cluster in the presampler layer. Only the next sampling is used as valid neighbor source.

Definition at line 187 of file TopoAutomatonSplitting.h.

◆ m_samplingNames

Gaudi::Property<std::vector<std::string> > TopoAutomatonSplitting::m_samplingNames {this, "SamplingNames", {}, "Name(s) of Calorimeter Samplings to consider for local maxima"}

private

vector of names of the calorimeter samplings to consider for seeds.

The default is to use all calorimeter samplings. Excluding a sampling from this vector prevents the definition of a seed cell in this sampling. Cells in those samplings are still used and incorporated in the topo clusters (both on the neighbor and the cell level) they can therefore even expand a cluster but not seed one ...

Definition at line 68 of file TopoAutomatonSplitting.h.

◆ m_secondarySamplingNames

Gaudi::Property<std::vector<std::string> > TopoAutomatonSplitting::m_secondarySamplingNames {this, "SecondarySamplingNames", {}, "Name(s) of secondary Calorimeter Samplings to consider for local maxima"}

private

vector of names of the secondary calorimeter samplings to consider.

Samplings in this list will be considered for local maxima only if no local max in the primary list is overlapping. By default this list is empty

Definition at line 77 of file TopoAutomatonSplitting.h.

◆ m_shareBorderCells

Gaudi::Property<bool> TopoAutomatonSplitting::m_shareBorderCells {this, "ShareBorderCells", false, "Whether or not to share cells at the boundary between two clusters"}

private

share cells at the border between two local maxima

this property needs to be set to true in order to treat cells which would be included in 2 clusters (for more then 2 the 2 with the largest E for the current seed cells are used) as shared cells. Shared cells are first excluded from the clustering and then clustered after all normal cells are clustered. The shared clusters are added to the 2 clusters they neighbor with the weights \(w_1 = E_1/(E_1+r E_2)\) and \(w_2 = 1-w_1\), where \(E_{1,2}\) are the current energies of the 2 neighboring clusters without the shared cells and \(r=\exp(d_1-d_2)\) is the ratio of the expected dependencies on the distances \(d_i\) (in units of a typical em shower scale) of each shared cell to the cluster centers. If the property is set to false the border cells are included in the normal clustering and the cluster with the largest E for the current seed cells gets the current border cell.

Definition at line 114 of file TopoAutomatonSplitting.h.

◆ m_timeFileName

Gaudi::Property<std::string> CaloGPUTimed::m_timeFileName

protectedinherited

File to which times should be saved.

Definition at line 50 of file CaloGPUTimed.h.

◆ m_timeMutex

std::shared_mutex CaloGPUTimed::m_timeMutex

mutableprotectedinherited

Mutex that is locked when recording times.

Definition at line 32 of file CaloGPUTimed.h.

◆ m_treatL1PredictedCellsAsGood

Gaudi::Property<bool> TopoAutomatonSplitting::m_treatL1PredictedCellsAsGood {this, "TreatL1PredictedCellsAsGood", true, "Treat bad cells with dead OTX if predicted from L1 as good"}

private

if set to true treat cells with a dead OTX which can be predicted by L1 trigger info as good instead of bad cells

Definition at line 139 of file TopoAutomatonSplitting.h.

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

Public Member Functions

Protected Member Functions

Protected Attributes

Private Member Functions

Private Attributes

Detailed Description

Constructor & Destructor Documentation

◆ TopoAutomatonSplitting()

◆ ~TopoAutomatonSplitting()

Member Function Documentation

◆ execute()

◆ finalize()

◆ handle()

◆ initialize()

◆ initialize_CUDA()

◆ initialize_non_CUDA()

◆ print_times()

◆ record_times() [1/3]

◆ record_times() [2/3]

◆ record_times() [3/3]

◆ record_times_helper() [1/3]

◆ record_times_helper() [2/3]

◆ record_times_helper() [3/3]

Member Data Documentation

◆ ATLAS_THREAD_SAFE [1/2]

◆ ATLAS_THREAD_SAFE [2/2]

◆ m_absOpt

◆ m_emShowerScale

◆ m_kernelSizeOptimizer

◆ m_measureTimes

◆ m_minEnergy

◆ m_nCells

◆ m_neighborOptionString

◆ m_options

◆ m_restrictHECIWandFCalNeighbors

◆ m_restrictPSNeighbors

◆ m_samplingNames

◆ m_secondarySamplingNames

◆ m_shareBorderCells

◆ m_timeFileName

◆ m_timeMutex

◆ m_treatL1PredictedCellsAsGood