ARICHHitRateCounter Class Reference

Class for monitoring beam background hit rates of ARICH. More...

#include <ARICHHitRateCounter.h>

Inheritance diagram for ARICHHitRateCounter:

Classes
struct	TreeStruct
	tree structure More...

Public Member Functions
	ARICHHitRateCounter ()
	Constructor.
virtual void	initialize (TTree *tree) override
	Class initializer: set branch addresses and other staf.
virtual void	clear () override
	Clear time-stamp buffer to prepare for 'accumulate'.
virtual void	accumulate (unsigned timeStamp) override
	Accumulate hits.
virtual void	normalize (unsigned timeStamp) override
	Normalize accumulated hits (e.g.

Private Member Functions
void	setActiveHapds ()
	Sets number of active hapds in each segment.

Private Attributes
TreeStruct	m_rates
	tree variables
std::map< unsigned, TreeStruct >	m_buffer
	average rates in time stamps
StoreArray< ARICHHit >	m_hits
	collection of digits
DBObjPtr< ARICHChannelMask >	m_channelMask
	channel mask
DBObjPtr< ARICHModulesInfo >	m_modulesInfo
	HAPD modules info.
double	m_activeHapds [18] = {0}
	number of active HAPDS in each segment
double	m_activeTotal = 0
	total number of active HAPDS
int	m_segmentMap [420] = {0}
	mapping from module ID to segments

Detailed Description

Class for monitoring beam background hit rates of ARICH.

Definition at line 31 of file ARICHHitRateCounter.h.

Constructor & Destructor Documentation

ARICHHitRateCounter()

ARICHHitRateCounter ( )

inline

Constructor.

Definition at line 60 of file ARICHHitRateCounter.h.

      {}

Member Function Documentation

accumulate()

void accumulate ( unsigned  timeStamp )

overridevirtual

Accumulate hits.

Parameters

timeStamp

time stamp

Implements HitRateBase.

Definition at line 53 of file ARICHHitRateCounter.cc.

    {
      // check if data are available
      if (not m_hits.isValid()) return;
 
      // get buffer element
      auto& rates = m_buffer[timeStamp];
 
      // increment event counter
      rates.numEvents++;
 
      // count and weight hits accoring to channel efficiecny
      for (const auto& hit : m_hits) {
        if (hit.getModule() < 1 || hit.getModule() > 420) continue;
        auto effi = m_modulesInfo->getChannelQE(hit.getModule(), hit.getChannel());
        float wt = std::min(1.0 / effi, 100.);
        rates.segmentRates[m_segmentMap[hit.getModule() - 1]] += wt;
        rates.averageRate += wt;
      }
 
      // set flag to true to indicate the rates are valid
      rates.valid = true;
 
    }

clear()

void clear ( )

overridevirtual

Clear time-stamp buffer to prepare for 'accumulate'.

Implements HitRateBase.

Definition at line 48 of file ARICHHitRateCounter.cc.

    {
      m_buffer.clear();
    }

initialize()

void initialize ( TTree *  tree )

overridevirtual

Class initializer: set branch addresses and other staf.

Parameters

tree	a valid TTree pointer

Implements HitRateBase.

Definition at line 24 of file ARICHHitRateCounter.cc.

    {
      // register collection(s) as optional, your detector might be excluded in DAQ
      m_hits.isOptional();
 
      // set branch address
      tree->Branch("arich", &m_rates, "segmentRates[18]/F:averageRate/F:numEvents/I:valid/O");
 
      // make map of modules to 18 segments
      // (first hapd ring is 6 segments (by sector), for the rest, each segment merges 3 hapd rings (again by sector))
      const int nModInRing[7] = {0, 42, 90, 144, 204, 270, 342};
      int iRing = 0;
      for (int i = 0; i < 420; i++) {
        if (i == nModInRing[iRing + 1]) iRing++;
        int segment = (i - nModInRing[iRing]) / (7 + iRing);
        if (iRing > 0) segment += 6;
        if (iRing > 3) segment += 6;
        m_segmentMap[i] = segment;
      }
 
      // set fractions of active channels
      setActiveHapds();
    }

normalize()

void normalize ( unsigned  timeStamp )

overridevirtual

Normalize accumulated hits (e.g.

transform to rates)

Parameters

timeStamp

time stamp

Implements HitRateBase.

Definition at line 78 of file ARICHHitRateCounter.cc.

    {
      // copy buffer element
      m_rates = m_buffer[timeStamp];
 
      if (not m_rates.valid) return;
 
      // normalize
      m_rates.normalize();
 
      // correct rates for masked-out channels
      if (m_channelMask.hasChanged()) setActiveHapds();
 
      for (int iSegment = 0; iSegment < 18; iSegment++) {
        double nHapds = m_activeHapds[iSegment];
        if (nHapds > 0) m_rates.segmentRates[iSegment] /= nHapds;
        else m_rates.segmentRates[iSegment] = 0;
      }
      m_rates.averageRate /= m_activeTotal;
    }

setActiveHapds()

void setActiveHapds ( )

private

Sets number of active hapds in each segment.

Definition at line 99 of file ARICHHitRateCounter.cc.

    {
      for (auto& nactive : m_activeHapds) nactive = 0;
 
      if (not m_channelMask.isValid()) {
        for (unsigned imod = 1; imod < 421; imod++)  m_activeHapds[m_segmentMap[imod - 1]] += 1.;
        m_activeTotal = 420;
        B2WARNING("ARICHHitRateCounter: no valid channel mask - all HAPDs set to active");
        return;
      }
 
      int nactiveTotal = 0;
      for (unsigned imod = 1; imod < 421; imod++) {
        int nactive = 0;
        for (unsigned ichn = 0; ichn < 144; ichn++) {
          if (m_channelMask->isActive(imod, ichn)) nactive++;
        }
        nactiveTotal += nactive;
        m_activeHapds[m_segmentMap[imod - 1]] += (float)nactive / 144.;
      }
      m_activeTotal = (float)nactiveTotal / 144.;
    }

Member Data Documentation

m_activeHapds

double m_activeHapds[18] = {0}

private

number of active HAPDS in each segment

Definition at line 107 of file ARICHHitRateCounter.h.

m_activeTotal

double m_activeTotal = 0

private

total number of active HAPDS

Definition at line 108 of file ARICHHitRateCounter.h.

m_buffer

std::map<unsigned, TreeStruct> m_buffer

private

average rates in time stamps

Definition at line 97 of file ARICHHitRateCounter.h.

m_channelMask

DBObjPtr<ARICHChannelMask> m_channelMask

private

channel mask

Definition at line 103 of file ARICHHitRateCounter.h.

m_hits

StoreArray<ARICHHit> m_hits

private

collection of digits

Definition at line 100 of file ARICHHitRateCounter.h.

m_modulesInfo

DBObjPtr<ARICHModulesInfo> m_modulesInfo

private

HAPD modules info.

Definition at line 104 of file ARICHHitRateCounter.h.

m_rates

TreeStruct m_rates

private

tree variables

Definition at line 94 of file ARICHHitRateCounter.h.

m_segmentMap

int m_segmentMap[420] = {0}

private

mapping from module ID to segments

Definition at line 109 of file ARICHHitRateCounter.h.

---

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

• background/modules/BeamBkgHitRateMonitor/include/ARICHHitRateCounter.h
• background/modules/BeamBkgHitRateMonitor/src/ARICHHitRateCounter.cc