shaperdsp.h

/**************************************************************************
 * BASF2 (Belle Analysis Framework 2)                                     *
 * Copyright(C) 2015-2018 - Belle II Collaboration                        *
 *                                                                        *
 * Author: The Belle II Collaboration                                     *
 * Contributors: Alexei Sibidanov                                         *
 *               Alex Bobrov                                              *
 *                                                                        *
 * This software is provided "as is" without any warranty.                *
 **************************************************************************/
 
#ifndef SHAPERDSP_T_H
#define SHAPERDSP_T_H
 
#include <vector>
#include <utility>
#include <limits>
 
static const double realNaN = std::numeric_limits<double>::quiet_NaN();
 
namespace Belle2 {
  namespace ECL {
    typedef std::pair<double, double> dd_t;
 
    class ShaperDSP_t {
    public:
 
      struct sv123shift_t {
        double t;
        double s0;
        double c0;
        double s1;
        double c1;
        double e0;
        double e1;
        double es;
        double ed;
 
        sv123shift_t() {} 
        sv123shift_t(double _t, const ShaperDSP_t& _p) { init(_t, _p);}  
        void init(double, const ShaperDSP_t&) __attribute__((noinline)); 
        sv123shift_t& operator +=(const sv123shift_t&);      
        sv123shift_t operator +(const sv123shift_t&) const;  
        bool validshift(const sv123shift_t& x) const {return t + x.t > 0;} 
      };
 
      struct shaperdspshift_t: public sv123shift_t {
        double et0; 
        double et1; 
        shaperdspshift_t() {} 
        shaperdspshift_t(double _t, const ShaperDSP_t& _p) { init(_t, _p);} 
        void init(double, const ShaperDSP_t&) __attribute__((noinline)); 
        shaperdspshift_t& operator +=(const shaperdspshift_t&); 
        shaperdspshift_t operator +(const shaperdspshift_t&) const; 
      };
 
    private:
      static const double _defs[]; //
      static constexpr double _filterdt = 0.2;
 
      double _cs0 = realNaN;
      double _cc0 = realNaN;
      double _cs1 = realNaN;
      double _cc1 = realNaN;
      double _ces = realNaN;
      double _ced = realNaN;
      double  _dw0 = realNaN;
      double  _dw1 = realNaN;
      double _dks0 = realNaN;
      double _dks1 = realNaN;
      double   _ds = realNaN;
      double   _dd = realNaN;
      double  _dt0 = realNaN;
      double  _dt1 = realNaN;
 
      double _toff = realNaN;
      double  _w0 = realNaN;
      double  _w1 = realNaN;
 
      double _ccc = realNaN;
 
      sv123shift_t _tp;
 
      sv123shift_t _tm;
 
      shaperdspshift_t _tstride;
 
      shaperdspshift_t  _toffset;
      shaperdspshift_t  _tzero;
 
      void Sv123_init(double t01, double tb1, double t02, double tb2, double td1, double ts1);
      void init(const double*, double) __attribute__((noinline));
      void init(const double*) __attribute__((noinline));
 
      double Sv123(const sv123shift_t&) const;
      dd_t ddSv123(const sv123shift_t&) const;
 
      double Sv123_filtered(const sv123shift_t&) const __attribute__((noinline));
      dd_t ddSv123_filtered(const sv123shift_t&) const __attribute__((noinline));
 
      double ShaperDSP(const shaperdspshift_t&) const;
      dd_t ddShaperDSP(const shaperdspshift_t&) const;
 
    public:
      ShaperDSP_t() { init(_defs); }
      explicit ShaperDSP_t(const std::vector<double>& s) { init(s); }
      ShaperDSP_t(const std::vector<double>& s, double u) { init(s, u); }
 
      ~ShaperDSP_t() {}
      void init(const std::vector<double>& s, double u = 27.7221);
      double operator()(double) const;
      double operator()(double*, double*);
      void settimestride(double);
      void setseedoffset(double);
      void settimeseed(double);
      void nextseed();
      void fillvector(std::vector<double>&) const;
      void fillvector(std::vector<dd_t>&) const;
      void fillvector(double, std::vector<double>&) const;
      void fillvector(double, std::vector<dd_t>&) const;
      void fillarray(int, double*) const;
      void fillarray(int, dd_t*) const;
      void fillarray(double, int, double*) const;
      void fillarray(double, int, dd_t*) const;
    };
  }
}
#endif//SHAPERDSP_T_H