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) {_waveform = wave;}
    void setFADC(int fadc) {_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 TimeRange;
    double TCEnergy[576][80];
    double TCTiming[576][80];
    double TCEnergy_tot[576];
    double TCTiming_tot[576];
 
    // background TC Energy [GeV]
    // double TCBGEnergy[576][80];
    // background TC Timing [GeV]
    // double TCBGTiming[576][80];
 
    double TCRawEnergy[576][60];
    double TCRawTiming[576][60];
    double TCRawBkgTag[576][60];
 
    TrgEclMapping* _TCMap;
    TrgEclDataBase* _DataBase;
 
    std::vector<std::vector<double>> MatrixParallel;
    std::vector<std::vector<double>> MatrixSerial;
    double TCBkgContribution[576][80];
    double TCSigContribution[576][80];
    int TCBeambkgTag[576][80];
    int _waveform;
    double WaveForm[576][64];
    int _FADC;
    int _BeambkgTag;
  };
} // end namespace Belle2