TurnOnHistos.h

Go to the documentation of this file.

// Copyright 2019-2020 CERN and copyright holders of ALICE O2.
// See https://alice-o2.web.cern.ch/copyright for details of the copyright holders.
// All rights not expressly granted are reserved.
//
// This software is distributed under the terms of the GNU General Public
// License v3 (GPL Version 3), copied verbatim in the file "COPYING".
//
// In applying this license CERN does not waive the privileges and immunities
// granted to it by virtue of its status as an Intergovernmental Organization
// or submit itself to any jurisdiction.
 
 
#ifndef PHOS_TURNONHISTOS_H
#define PHOS_TURNONHISTOS_H
 
#include <bitset>
#include <array>
#include "TObject.h"
 
namespace o2
{
 
namespace phos
{
 
class TurnOnHistos
{
 public:
  //class to collect statistics to calculate trigger turn-on curves and trigger bad maps
  static constexpr short NCHANNELS = 3136; 
  static constexpr short NDDL = 14;        
  static constexpr short Npt = 200;        
  static constexpr float dpt = 0.1;        
 
  TurnOnHistos() = default;
 
  TurnOnHistos& operator=(const TurnOnHistos& other) = default;
 
  ~TurnOnHistos() = default;
 
  void merge(TurnOnHistos& other);
 
  void fillTotSp(short ddl, float e)
  {
    short bin = e / dpt;
    if (bin < Npt) {
      mTotSp[ddl][bin]++;
    }
  }
 
  void fillFiredSp(short ddl, float e)
  {
    short bin = e / dpt;
    if (bin < Npt) {
      mTrSp[ddl][bin]++;
    }
  }
 
  void fillFiredMap(const std::bitset<NCHANNELS>& bs)
  {
    for (short i = NCHANNELS; --i;) {
      if (bs[i]) {
        mGoodMap[i]++;
      }
    }
  }
 
  void fillNoisyMap(const std::bitset<NCHANNELS>& bs)
  {
    for (short i = NCHANNELS; --i;) {
      if (bs[i]) {
        mNoisyMap[i]++;
      }
    }
  }
 
  //getters now
  const std::array<float, Npt>& getTotSpectrum(short ddl) const { return mTotSp[ddl]; }
  const std::array<float, Npt>& getTrSpectrum(short ddl) const { return mTrSp[ddl]; }
  const std::array<float, NCHANNELS>& getGoodMap() const { return mGoodMap; }
  const std::array<float, NCHANNELS>& getNoisyMap() const { return mNoisyMap; }
 
 private:
  std::array<float, NCHANNELS> mGoodMap;           
  std::array<float, NCHANNELS> mNoisyMap;          
  std::array<std::array<float, Npt>, NDDL> mTotSp; 
  std::array<std::array<float, Npt>, NDDL> mTrSp;  
 
  ClassDefNV(TurnOnHistos, 1);
};
 
} // namespace phos
 
} // namespace o2
#endif