TrgEclDigitizer.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
 
#include <TObject.h>
#include "trg/ecl/TrgEclMapping.h"
#include "trg/ecl/TrgEclDataBase.h"
#include "ecl/dataobjects/ECLHit.h"
#include "trg/ecl/dataobjects/TRGECLBGTCHit.h"
 
namespace Belle2 {

  class TrgEclDigitizer : public TObject {
 
  public:

    TrgEclDigitizer();
    virtual ~TrgEclDigitizer();

    void setup(int);
    void getTCHit(int);
    void digitization01(std::vector<std::vector<double>>&,
                        std::vector<std::vector<double>>&);
    void digitization02(std::vector<std::vector<double>>&,
                        std::vector<std::vector<double>>&);
    void save(int);
    void setWaveform(int wave) { m_SaveTCWaveForm = wave;}
    void setFADC(int fadc) { m_FADC = fadc;}
 
 

    double FADC(int, double);
    double SimplifiedFADC(int, double);
    double interFADC(double);
    double ShapeF(double, double, double, double, double, double, double);
    double ShapeF(double, double);
    double u_max(double, double);
 
  private:

    double m_TimeRange;
    double m_TCEnergy[576][80];
    double m_TCTiming[576][80];
    double m_TCEnergy_tot[576];
    double m_TCTiming_tot[576];
 
    // background TC Energy [GeV]
    // double TCBGEnergy[576][80];
    // background TC Timing [GeV]
    // double TCBGTiming[576][80];

    double m_TCRawEnergy[576][60];
    double m_TCRawTiming[576][60];
    double m_TCRawBkgTag[576][60];

    TrgEclMapping* m_TCMap;
    TrgEclDataBase* m_DataBase;

    std::vector<std::vector<double>> m_MatrixParallel;
    std::vector<std::vector<double>> m_MatrixSerial;
    double m_TCBkgContribution[576][80];
    double m_TCSigContribution[576][80];
    int m_TCBeambkgTag[576][80];
    int m_SaveTCWaveForm;
    double m_WaveForm[576][64];
    int m_FADC;
    int m_BeambkgTag;
  };

} // end namespace Belle2