TOPPmtQE.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2016 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric, Alessandro Gaz                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <TObject.h>
#include <vector>
 
namespace Belle2 {
  class TOPPmtQE : public TObject {
  public:
    enum {c_NumPmtPixels = 16};
 
    TOPPmtQE()
    {}
 
    TOPPmtQE(const std::string& serialNumber,
             float lambdaFirst, float lambdaStep, float CE0, float CE):
      m_serialNumber(serialNumber), m_lambdaFirst(lambdaFirst), m_lambdaStep(lambdaStep),
      m_CE_noB(CE0), m_CE_withB(CE)
    {}
 
    void setQE(unsigned pmtPixel, const std::vector<float>& qe)
    {
      pmtPixel--;
      if (pmtPixel < c_NumPmtPixels) m_QE[pmtPixel] = qe;
    }
 
    const std::string& getSerialNumber() const {return m_serialNumber;}
 
    const std::vector<float>& getQE(unsigned pmtPixel) const;
 
    double getQE(unsigned pmtPixel, double lambda) const;
 
    const std::vector<float>& getEnvelopeQE() const;
 
    double getEnvelopeQE(double lambda) const;
 
    double getEfficiency(unsigned pmtPixel, double lambda, bool BfieldOn) const;
 
    double getLambdaFirst() const {return m_lambdaFirst;}
 
    double getLambdaLast() const;
 
    double getLambdaLast(unsigned pmtPixel) const;
 
    double getLambdaStep() const {return m_lambdaStep;}
 
    double getCE(bool BfieldOn) const;
 
 
  private:
 
    void setEnvelopeQE() const;
 
    double interpolate(double lambda, const std::vector<float>& QE) const;
 
    std::string m_serialNumber; 
    std::vector<float> m_QE[c_NumPmtPixels]; 
    float m_lambdaFirst = 0; 
    float m_lambdaStep = 0;  
    float m_CE_noB = 0;      
    float m_CE_withB = 0;    
    mutable std::vector<float> m_envelopeQE; 
 
    ClassDef(TOPPmtQE, 3); 
  };
 
} // end namespace Belle2