EclConfigurationPure.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
 
/* ECL headers. */
#include <ecl/digitization/EclConfiguration.h>
 
class TH1F;
 
namespace Belle2 {
  namespace ECL {
 
    class EclConfigurationPure {
    public:
      static constexpr int    m_nch      = 72 * 16; 
      static constexpr int    m_nsmp     = EclConfiguration::m_nsmp; 
      static constexpr double m_tmin     = -15; 
      static constexpr int    m_ntrg     = EclConfiguration::m_ntrg; 
      static constexpr int    m_nlPure   = EclConfiguration::m_nl * 15; 
      static constexpr int    m_ns       = EclConfiguration::m_ns; 
      static constexpr int    m_ndtPure  = m_ns; 
    private:
      static double m_tickPure;
 
    public:
      static double getTickPure()
      {
        if (m_tickPure < 0) {
          m_tickPure = EclConfiguration::getTick() / EclConfiguration::m_ntrg * 8;
        }
        return m_tickPure;
      }
      static void setTickPure(double newval)
      {
        m_tickPure = newval;
      }
 
      struct signalsamplepure_t {
        double m_sumscale; 
        double m_ft[m_nlPure * m_ns]; 
        double m_ft1[m_nlPure * m_ns]; 
        void InitSample(const TH1F*, const TH1F*);
        double Accumulate(const double a, const double t0, double* s) const;
      };
 
      struct adccountspure_t {
        double total; 
        double c[m_nsmp]; 
        void AddHit(const double a, const double t0, const signalsamplepure_t& q);
      };
 
      struct fitparamspure_t {
        typedef double double_matrix[16][2 * m_ndtPure];
        typedef double fine_array[2 * m_ndtPure]; 
        double invC[16][16]; 
        double_matrix f;     
        double_matrix f1;    
        fine_array c110;
        fine_array c200; 
        fine_array c020; 
        fine_array c101; 
        fine_array c011; 
        double c002;     
        double_matrix c100;
        double_matrix c010; 
        double c001[16];    
      };
 
    };
 
  }
}