PulseHeightGenerator.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2017 - Belle II Collaboration                             *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Marko Staric                                             *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#pragma once
 
#include <math.h>
 
namespace Belle2 {
  namespace TOP {
 
    class PulseHeightGenerator {
 
    public:
 
      PulseHeightGenerator()
      {}
 
      PulseHeightGenerator(double x0, double p1, double p2, double xmax);
 
      double getValue(double x) const
      {
        x /= m_x0;
        return pow(x, m_p1) * exp(-pow(x, m_p2));
      }
 
      double generate() const;
 
    private:
 
      double m_x0 = 0; 
      double m_p1 = 0; 
      double m_p2 = 0; 
      double m_xmax = 0; 
      double m_vPeak = 0; 
    };
 
 
  } // TOP
} // Belle2