ECLHitRateCounter.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2019 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric, Andrea Fodor                               *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <background/modules/BeamBkgHitRateMonitor/HitRateBase.h>
#include <framework/datastore/StoreArray.h>
#include <ecl/dataobjects/ECLDigit.h>
#include <ecl/dataobjects/ECLDsp.h>
#include <ecl/geometry/ECLGeometryPar.h>
#include <TTree.h>
 
namespace Belle2 {
  namespace Background {
 
    class ECLHitRateCounter: public HitRateBase {
 
    public:
 
      struct TreeStruct {
 
        float averageDspBkgRate[16] = {0}; 
        int numEvents = 0; 
        bool validDspRate = false; 
        void normalize()
        {
          if (numEvents == 0) return;
 
          for (int i = 0; i < 16; i++) {
            averageDspBkgRate[i] /= numEvents;
          }
        }
 
 
      };
 
      ECLHitRateCounter()
      {}
 
      virtual void initialize(TTree* tree) override;
 
      virtual void clear() override;
 
      virtual void accumulate(unsigned timeStamp) override;
 
      virtual void normalize(unsigned timeStamp) override;
 
    private:
 
      // class parameters: to be set via constructor or setters
 
      //functions
      void segmentECL();
 
      int findECLSegment(int cellid)
      {
        return m_segmentMap.find(cellid)->second;
      }
 
      void findElectronicsNoise();
 
      // tree structure
      TreeStruct m_rates; 
      // buffer
      std::map<unsigned, TreeStruct> m_buffer; 
      // collections
      StoreArray<ECLDigit> m_digits;  
      StoreArray<ECLDsp> m_dsps; 
      std::vector<float> m_ADCtoEnergy; 
      std::vector<float> m_waveformNoise; 
      // other
      Belle2::ECL::ECLGeometryPar* m_geometry{nullptr}; 
      std::map<int, int> m_segmentMap; 
      int m_crystalsInSegment[16] = {0}; 
    };
  }
}