AsicBackgroundDetector Class Reference

Marks hits as background based on the result of a filter. More...

#include <AsicBackgroundDetector.h>

Inheritance diagram for AsicBackgroundDetector:

Collaboration diagram for AsicBackgroundDetector:

Public Types
using	IOTypes = std::tuple< AIOTypes... >
	Types that should be served to apply on invokation.
using	IOVectors = std::tuple< std::vector< AIOTypes >... >
	Vector types that should be served to apply on invokation.

Public Member Functions
	AsicBackgroundDetector ()=default
	Default constructor.
void	initialize () final
	Access database here:
void	beginRun () final
	Reload channel map if needed.
std::string	getDescription () final
	Short description of the findlet.
virtual void	exposeParameters (ModuleParamList *moduleParamList, const std::string &prefix) final
	Expose the parameters to a module.
virtual void	apply (std::vector< CDCWireHit > &wireHits) final
	Main algorithm marking hit as background.
virtual void	apply (ToVector< AIOTypes > &... ioVectors)=0
	Main function executing the algorithm.
void	beginEvent () override
	Receive and dispatch signal for the start of a new event.
void	endRun () override
	Receive and dispatch signal for the end of the run.
void	terminate () override
	Receive and dispatch Signal for termination of the event processing.

Protected Types
using	ToVector = typename ToVectorImpl< T >::Type
	Short hand for ToRangeImpl.

Protected Member Functions
void	addProcessingSignalListener (ProcessingSignalListener *psl)
	Register a processing signal listener to be notified.
int	getNProcessingSignalListener ()
	Get the number of currently registered listeners.

Private Types
using	Super = Findlet< CDCWireHit & >
	Type of the base class.

Private Member Functions
void	applyAsicFilter (std::vector< CDCWireHit * > &wireHits)
	Algorithm marking hit as background for each CDC ASIC. More...

Private Attributes
std::unique_ptr< DBArray< CDCChannelMap > >	m_channelMapFromDB
	Channel map retrieved from DB.
std::map< int, std::pair< int, int > >	m_map
	map from ewire to board/channel ID
size_t	m_minimal_hit_number {4}
	min. number of hits in ASIC for background check
double	m_deviation_from_median {10.}
	distance from median TDC, to be considered as bg.
size_t	m_nsignal_max {4}
	max. number of signal-like hits in ASIC for background check
size_t	m_max_asic_error_messages {1}
	max. number of logged error messages for number of hits per ASIC check
std::vector< ProcessingSignalListener * >	m_subordinaryProcessingSignalListeners
	References to subordinary signal processing listener contained in this findlet.
bool	m_initialized = false
	Flag to keep track whether initialization happend before.
bool	m_terminated = false
	Flag to keep track whether termination happend before.
std::string	m_initializedAs
	Name of the type during initialisation.

Detailed Description

Marks hits as background based on the result of a filter.

Definition at line 26 of file AsicBackgroundDetector.h.

Member Function Documentation

applyAsicFilter()

void applyAsicFilter ( std::vector< CDCWireHit * > &  wireHits )

private

Algorithm marking hit as background for each CDC ASIC.

Extra information:

This is abnormal situation, detected for few runs, related to CDC unpacker. Hits are to be marked as background.

Definition at line 93 of file AsicBackgroundDetector.cc.

{
 
  if (wireHits.size() < m_minimal_hit_number) {
    return;
  };
 
  if (wireHits.size() > 8) {
    auto eWire = wireHits[0]->getWireID().getEWire();
    auto board = m_map[eWire].first;
 
    if (m_max_asic_error_messages > 0) {
      B2ERROR("Number of hits per asic should not exceed 8, observe too many hits."
              << LogVar("nHits", wireHits.size())
              << LogVar("Board ID", board));
      m_max_asic_error_messages -= 1;
    } else {
      B2WARNING("Number of hits per asic should not exceed 8, observe too many hits."
                << LogVar("nHits", wireHits.size())
                << LogVar("Board ID", board));
    }
 
    for (auto& hit : wireHits) {
      (*hit)->setBackgroundFlag();
      (*hit)->setTakenFlag();
    }
    return;
  }
 
  // compute median time:
  vector<short> times;
  for (auto& hit : wireHits) {
    short tdc = hit->getHit()->getTDCCount();
    times.push_back(tdc);
  }
  sort(times.begin(), times.end());
  int mid = times.size() / 2;
  double median = times.size() % 2 == 0 ? (times[mid] + times[mid - 1]) / 2 : times[mid];
 
  size_t nbg = 0;
  int adcOnMedian = 0;
  int adcOffMedian = 0;
  for (auto& hit : wireHits) {
    int adc = hit->getHit()->getADCCount();
    if (fabs(hit->getHit()->getTDCCount() - median) < m_deviation_from_median) {
      nbg++;
      if (adc > adcOnMedian) adcOnMedian = adc;
    } else {
      if (adc > adcOffMedian) adcOffMedian = adc;
    }
  }
 
  if ((nbg < times.size())     // at least 1 hit with different TDC ("signal")
      && (nbg > 1) // more than one candidate hits
      && (nbg > times.size() - m_nsignal_max) // a few background hits
      && (adcOnMedian < adcOffMedian) // triggered by large ADC "signal"
     ) {
 
    // mark hits too close to the median time as background:
    for (auto& hit : wireHits) {
      if (fabs(hit->getHit()->getTDCCount() - median) < m_deviation_from_median) {
        (*hit)->setBackgroundFlag();
        (*hit)->setTakenFlag();
      }
    }
  }
}

---

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

• tracking/trackFindingCDC/findlets/minimal/include/AsicBackgroundDetector.h
• tracking/trackFindingCDC/findlets/minimal/src/AsicBackgroundDetector.cc