ECLnOptimal.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 <TObject.h>
#include <framework/logging/Logger.h>
#include <TH2F.h>
 
namespace Belle2 {
//---------------------------------------------------
  class ECLnOptimal: public TObject {
 
//---------------------------------------------------
  public:
 
    ECLnOptimal() :
      m_eUpperBoundariesFwd(),
      m_eUpperBoundariesBrl(),
      m_eUpperBoundariesBwd(),
      m_groupNumber(),
      m_nOptimal(),
      m_peakFracEnergy(),
      m_bias(),
      m_logPeakEnergy()
    { };
 
    //---------------------------------------------------
    //..Getters for the various calibration components
 
    const std::vector<float>& getUpperBoundariesFwd() const {return m_eUpperBoundariesFwd;}
 
    const std::vector<float>& getUpperBoundariesBrl() const {return m_eUpperBoundariesBrl;}
 
    const std::vector<float>& getUpperBoundariesBwd() const {return m_eUpperBoundariesBwd;}
 
    const std::vector<int>& getGroupNumber() const {return m_groupNumber;}
 
    const TH2F getNOptimal() const {return m_nOptimal;}
 
    //..The following have 3 xBins per group for energy interpolation.
    //  Value for nominal energy bin is histogram x bin ix = 3*iGroup + 1;
    //  Value for lower energy bin is ix = 3*iGroup + 2;
    //  Value for higher energy bin is ix = 3*iGroup + 3;
    //  In all three cases, the histogram y bin is iy = iEnergy + 1
 
    const TH2F getPeakFracEnergy() const {return m_peakFracEnergy;}
 
    const TH2F getBias() const {return m_bias;}
 
    const TH2F getLogPeakEnergy() const {return m_logPeakEnergy;}
 
    //---------------------------------------------------
    //..Setters for the various components
 
    void setUpperBoundariesFwd(const std::vector<float>& eUpperBoundariesFwd) {m_eUpperBoundariesFwd = eUpperBoundariesFwd;}
 
    void setUpperBoundariesBrl(const std::vector<float>& eUpperBoundariesBrl) {m_eUpperBoundariesBrl = eUpperBoundariesBrl;}
 
    void setUpperBoundariesBwd(const std::vector<float>& eUpperBoundariesBwd) {m_eUpperBoundariesBwd = eUpperBoundariesBwd;}
 
    void setGroupNumber(const std::vector<int>& groupNumber) {m_groupNumber = groupNumber;}
 
    void setNOptimal(const TH2F& nOptimal) {m_nOptimal = nOptimal;}
 
    void setPeakFracEnergy(const TH2F& peakFracEnergy) {m_peakFracEnergy = peakFracEnergy;}
 
    void setBias(const TH2F& bias) {m_bias = bias;}
 
    void setLogPeakEnergy(const TH2F& logPeakEnergy) {m_logPeakEnergy = logPeakEnergy;}
 
//---------------------------------------------------
  private:
    std::vector<float> m_eUpperBoundariesFwd; 
    std::vector<float> m_eUpperBoundariesBrl; 
    std::vector<float> m_eUpperBoundariesBwd; 
    std::vector<int> m_groupNumber; 
    TH2F m_nOptimal; 
    //..For the following, there are 3 bins in x for each group.
    //  First contains the value for the specified group number and energy point,
    //  corresponding to the optimal number of crystals.
    //  Second is the value, for that number of crystals, for the specified group
    //  but the next lower energy point.
    //  Third is the value for the next higher energy point.
    TH2F m_peakFracEnergy; 
    TH2F m_bias; 
    TH2F m_logPeakEnergy; 
    ClassDef(ECLnOptimal, 1); 
  };
}