SVDPulseShapeCalibrations.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 <vxd/dataobjects/VxdID.h>
#include <svd/dbobjects/SVDCalibrationsBase.h>
#include <svd/dbobjects/SVDCalibrationsVector.h>
#include <svd/dbobjects/SVDStripCalAmp.h>
#include <framework/database/DBObjPtr.h>
#include <framework/logging/Logger.h>
#include <string>
 
 
#include <math.h>
 
namespace Belle2 {
  class SVDPulseShapeCalibrations {
  public:
    static std::string calAmp_name; 
    typedef SVDCalibrationsBase< SVDCalibrationsVector< SVDStripCalAmp > >
    t_calAmp_payload; 
    SVDPulseShapeCalibrations()
      : m_calAmp_aDBObjPtr(calAmp_name)
    {
      m_calAmp_aDBObjPtr.addCallback([ this ](const std::string&) -> void {
        B2DEBUG(20, "SVDPulseShapeCalibrations: from now on we are using " <<
                this->m_calAmp_aDBObjPtr -> get_uniqueID()); });
    }
 
    inline double getChargeFromADC(
      const Belle2::VxdID& sensorID,
      const bool& isU, const unsigned short& strip,
      const double& pulseADC
    ) const
    {
      return  pulseADC / getGain(sensorID, isU, strip);
    }
 
    inline long int getADCFromCharge(
      const Belle2::VxdID& sensorID,
      const bool& isU, const unsigned short& strip,
      const double& charge) const
    {
      return roundl(charge * getGain(sensorID, isU, strip));
    }
 
    inline unsigned char getCappedADCFromCharge(
      const Belle2::VxdID& sensorID,
      const bool& isU, const unsigned short& strip,
      const double& charge) const
    {
      auto chargeLongInt = getADCFromCharge(sensorID, isU, strip, charge);
      if (chargeLongInt < 0)
        return 0;
      return chargeLongInt > 255 ? 255 : chargeLongInt;
    }
 
    inline float getPeakTime(const VxdID& sensorID, const bool& isU,
                             const unsigned short& strip) const
    {
      return m_calAmp_aDBObjPtr->getReference(sensorID.getLayerNumber(), sensorID.getLadderNumber(),
                                              sensorID.getSensorNumber(), m_calAmp_aDBObjPtr->sideIndex(isU),
                                              strip).peakTime;
    }
 
    inline float getWidth(const VxdID& sensorID, const bool& isU,
                          const unsigned short& strip) const
    {
      return m_calAmp_aDBObjPtr->getReference(sensorID.getLayerNumber(), sensorID.getLadderNumber(), sensorID.getSensorNumber(),
                                              m_calAmp_aDBObjPtr->sideIndex(isU), strip).pulseWidth;
 
    }
 
    TString getUniqueID() { return m_calAmp_aDBObjPtr->get_uniqueID(); }
 
    bool isValid()
    {
      return m_calAmp_aDBObjPtr.isValid();
    }
 
  private:
 
    inline float getGain(const VxdID& sensorID, const bool& isU,
                         const unsigned short& strip) const
    {
      return m_calAmp_aDBObjPtr->getReference(sensorID.getLayerNumber(), sensorID.getLadderNumber(), sensorID.getSensorNumber(),
                                              m_calAmp_aDBObjPtr->sideIndex(isU), strip).gain ;
 
    }
 
  private:
    DBObjPtr< t_calAmp_payload > m_calAmp_aDBObjPtr; 
  };
}