TOPNominalQE.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 <top/dbobjects/TOPGeoBase.h>
#include <top/dbobjects/TOPWavelengthFilter.h>
#include <vector>
 
namespace Belle2 {

  class TOPNominalQE : public TOPGeoBase {
  public:

    TOPNominalQE()
    {}

    TOPNominalQE(float lambdaFirst, float lambdaStep, float CE,
                 const std::vector<float>& qe,
                 const std::string& name = "TOPNominalQE"): TOPGeoBase(name),
      m_lambdaFirst(lambdaFirst), m_lambdaStep(lambdaStep), m_CE(CE), m_QE(qe)
    {}

    void set(float lambdaFirst, float lambdaStep, float CE,
             const std::vector<float>& qe,
             const std::string& name)
    {
      m_lambdaFirst = lambdaFirst;
      m_lambdaStep = lambdaStep;
      m_CE = CE;
      m_QE = qe;
      setName(name);
    }

    void applyFilterTransmission(const TOPWavelengthFilter& filter);

    void clear()
    {
      m_lambdaFirst = 0;
      m_lambdaStep = 0;
      m_CE = 0;
      m_QE.clear();
      setName("");
    }

    bool isEmpty() const {return m_QE.empty();}

    const std::vector<float>& getQE() const {return m_QE;}

    double getQE(double lambda) const;

    double getEfficiency(double lambda) const {return getQE(lambda) * m_CE;}

    double getLambdaFirst() const {return m_lambdaFirst;}

    double getLambdaLast() const
    {
      return m_lambdaFirst + (m_QE.size() - 1) * m_lambdaStep;
    }

    double getLambdaStep() const {return m_lambdaStep;}

    double getMinLambda() const;

    double getMaxLambda() const;

    double getCE() const {return m_CE;}

    bool isConsistent() const override;

    void print(const std::string& title = "Nominal quantum efficiency") const override;
 
 
  private:
 
    float m_lambdaFirst = 0; 
    float m_lambdaStep = 0;  
    float m_CE = 0;          
    std::vector<float> m_QE; 
 
    ClassDefOverride(TOPNominalQE, 1); 
 
  };

} // end namespace Belle2