TOPHitRateCounter.h

/**************************************************************************
 * basf2 (Belle II Analysis Software Framework)                           *
 * Author: The Belle II Collaboration                                     *
 *                                                                        *
 * See git log for contributors and copyright holders.                    *
 * This file is licensed under LGPL-3.0, see LICENSE.md.                  *
 **************************************************************************/
 
#pragma once
 
#include <background/modules/BeamBkgHitRateMonitor/HitRateBase.h>
#include <framework/datastore/StoreArray.h>
#include <framework/database/DBObjPtr.h>
#include <top/dataobjects/TOPDigit.h>
#include <top/dbobjects/TOPCalChannelMask.h>
#include <TTree.h>
#include <TH1F.h>
#include <map>
 
namespace Belle2 {
  namespace Background {
 
    class TOPHitRateCounter: public HitRateBase {
 
    public:
 
      struct TreeStruct {
        float slotRates[16] = {0}; 
        float averageRate = 0; 
        int numEvents = 0; 
        bool valid = false;  
        void normalize()
        {
          if (numEvents == 0) return;
          for (auto& slotRate : slotRates) slotRate /= numEvents;
          averageRate /= numEvents;
        }
 
      };
 
      TOPHitRateCounter(double timeOffset, double timeWindow):
        m_timeOffset(timeOffset), m_timeWindow(timeWindow)
      {}
 
      virtual void initialize(TTree* tree) override;
 
      virtual void clear() override;
 
      virtual void accumulate(unsigned timeStamp) override;
 
      virtual void normalize(unsigned timeStamp) override;
 
    private:
 
      void setActiveFractions();
 
      // class parameters: to be set via constructor or setters
      double m_timeOffset = 0; 
      double m_timeWindow = 0; 
      // tree structure
      TreeStruct m_rates; 
      // buffer
      std::map<unsigned, TreeStruct> m_buffer; 
      // collections
      StoreArray<TOPDigit> m_digits;  
      // DB payloads
      DBObjPtr<TOPCalChannelMask> m_channelMask; 
      // control histograms: to check if time offset and window size are correctly set
      TH1F* m_hits = 0; 
      TH1F* m_hitsInWindow = 0; 
      // other
      double m_activeFractions[16] = {0}; 
      double m_activeTotal = 0; 
    };
 
  } // Background namespace
} // Belle2 namespace