JGetRiseTime.hh

Go to the documentation of this file.

#ifndef __JTRIGGER__JGETRISETIME__
#define __JTRIGGER__JGETRISETIME__

#include <vector>
#include <cmath>

#include "JLang/JCC.hh"


/**
 * \file
 *
 * Rise time evaluation.
 * \author mdejong
 */
namespace JTRIGGER {}
namespace JPP { using namespace JTRIGGER; }

namespace JTRIGGER {

  /**
   * Auxiliary class for rise time evaluation.
   *
   *     \f[ \Delta t  =  p_0 \exp(p_1 \sqrt{x} + p_2 x) ~+~ p_3 \f]
   *
   * where
   * \f$\Delta t\f$ (ns) is the time it takes the analogue pulse to rise to the discriminator threshold and
   * \f$x\f$ the time over threshold (ns).
   */
  struct JGetRiseTime :
    public std::vector<double>
  {
    /**
     * Default constructor.
     */
    JGetRiseTime()
    {
      for (int i = 0; i <= 256; ++i) {

        const double tot = i;
        const double t1  = p0() * exp(p1()*sqrt(tot) + p2()*tot) + p3();

        this->push_back(t1);
      }
    }


    /**
     * Get time of analogue signal to rise to discriminator threshold (so-called time slewing).
     *
     * \param  tot            time over threshold [ns]
     * \return                rise time [ns]
     */
    double operator()(const double tot) const
    {
      int i = (int) tot;

      if      (i <  0)
        return *(this->begin());
      else if (i >= (int) this->size())
        return *(this->rbegin());
      else
        return (*this)[i];
    }


    static double p0() { return      7.70824; }
    static double p1() { return   0.00879447; }
    static double p2() { return   -0.0621101; }
    static double p3() { return     -1.90226; }
  };


  /**
   * Function object for rise time evaluation.
   */
  static const JGetRiseTime getRiseTime;
}

#endif