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