TpcDigitizerModule.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.                  *
 **************************************************************************/
 
#ifndef TPCDIGITIZERMODULE_H
#define TPCDIGITIZERMODULE_H
 
#include <framework/core/Module.h>
#include <framework/datastore/StoreArray.h>
#include <beast/microtpc/dataobjects/MicrotpcHit.h>
 
#include <TF1.h>
#include <TVector3.h>
#include <TVector2.h>
 
#include <vector>
#include <map>
#include <tuple>
 
constexpr int MAXSIZE = 26880;
constexpr int MAXtSIZE = 10000;
 
namespace Belle2 {
  namespace microtpc {
    class TpcDigitizerModule : public Module {
 
    public:
 
      TpcDigitizerModule();
 
      virtual ~TpcDigitizerModule();
 
      virtual void initialize() override;
 
      virtual void beginRun() override;
 
      virtual void event() override;
 
      virtual void endRun() override;
 
      virtual void terminate() override;
 
 
    private:
 
      StoreArray<MicrotpcHit> m_microtpcHit; 
      void getXMLData();
 
      //bool Pixelization(int);
      void Pixelization();
 
      /*TLorentzVector Drift(
        double x1, double y1, double z1,
        double st, double sl, double vd
      );*/
      virtual void Drift(double, double, double, double&, double&, double&, double&, double, double, double);
 
      TVector2 GEMGeo1(double x1, double y1);
 
      TVector2 GEMGeo2(double x1, double y1);
 
      Int_t m_phase;
      TF1* fctToT_Calib1;
      TF1* fctToT_Calib2;
      double m_GEMGain1;
      double m_GEMGainRMS1;
      double m_ScaleGain1;
      double m_GEMGain2;
      double m_GEMGainRMS2;
      double m_ScaleGain2;
      double m_GEMpitch;
      int m_PixelThreshold;
      int m_PixelThresholdRMS;
      int m_ChipColumnNb;
      int m_ChipRowNb;
      double m_ChipColumnX;
      double m_ChipRowY;
      int m_PixelTimeBinNb;
      double m_PixelTimeBin;
      double m_TOTA1;
      double m_TOTB1;
      double m_TOTC1;
      double m_TOTQ1;
      double m_TOTA2;
      double m_TOTB2;
      double m_TOTC2;
      double m_TOTQ2;
      double m_z_DG;
      double m_z_TG;
      double m_z_CG;
      double m_Dt_DG;
      double m_Dt_TG;
      double m_Dt_CG;
      double m_Dl_DG;
      double m_Dl_TG;
      double m_Dl_CG;
      double m_v_DG;
      double m_v_TG;
      double m_v_CG;
      double m_Workfct;
      double m_Fanofac;
      double m_GasAbs;
      //int m_dchip[10][80][336][MAXtSIZE];
      std::map<std::tuple<int, int, int>, int> m_dchip;
      std::map<std::tuple<int, int>, int> m_dchip_map;
      std::map<std::tuple<int, int>, int> m_dchip_pdg_map;
      std::map<std::tuple<int, int>, int> m_dchip_trkID_map;
      std::map<std::tuple<int, int>, int> m_dchip_detNb_map;
      int m_LookAtRec;
      int m_nTPC = 0;
      std::vector<TVector3> m_TPCCenter;
      std::vector<float> m_TPCAngleX;
      std::vector<float> m_TPCAngleZ;
 
      double m_lowerTimingCut = 0;
      double m_upperTimingCut = 1000000;
 
      int olddetNb = -1;
 
      int oldtrkID = -1;
 
    };
 
  }
}
 
#endif /* TPCDIGITIZERMODULE_H */