release-05-02-19/doxygen/TpcDigitizerModule_8h_source.html

/**************************************************************************

 * BASF2 (Belle Analysis Framework 2)                                     *

 * Copyright(C) 2014 - Belle II Collaboration                             *

 *                                                                        *

 * Author: The Belle II Collaboration                                     *

 * Contributors: Igal Jaegle                                              *

 *                                                                        *

 * This software is provided "as is" without any warranty.                *

 **************************************************************************/


#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 */