d2/dc6/PHMicromegasTpcTrackMatching_8h_source.html

// Tell emacs that this is a C++ source

//  -*- C++ -*-.

#ifndef PHMICROMEGASTPCTRACKMATCHING_H

#define PHMICROMEGASTPCTRACKMATCHING_H


#include <fun4all/SubsysReco.h>

#include <tpc/TpcDistortionCorrection.h>

#include <tpc/TpcClusterZCrossingCorrection.h>


#include <trackbase/TrkrDefs.h>


#include <array>

#include <string>

#include <vector>


class ActsGeometry;

class TrkrClusterContainer;

class TrkrClusterIterationMapv1;

class TrackSeedContainer;

class PHCompositeNode;

class PHG4CylinderGeomContainer;

class TrackSeed;

class TrkrCluster;

class TF1;

class TH1;


class PHMicromegasTpcTrackMatching : public SubsysReco

{


  public:

  PHMicromegasTpcTrackMatching(const std::string &name = "PHMicromegasTpcTrackMatching");

  ~PHMicromegasTpcTrackMatching() override = default;


  void set_rphi_search_window_lyr1(const double win){_rphi_search_win[0] = win;}

  void set_z_search_window_lyr1(const double win){_z_search_win[0] = win;}

  void set_rphi_search_window_lyr2(const double win){_rphi_search_win[1] = win;}

  void set_z_search_window_lyr2(const double win){_z_search_win[1] = win;}

  void set_min_tpc_layer(const unsigned int layer){_min_tpc_layer = layer;}

  void set_test_windows_printout(const bool test){_test_windows = test;}

  void SetIteration(int iter){_n_iteration = iter;}


  int InitRun(PHCompositeNode* topNode) override;

  int process_event(PHCompositeNode*) override;

  int End(PHCompositeNode*) override;


  // deprecated calls

  inline void set_sc_calib_mode(const bool) {}

  inline void set_collision_rate(const double) {}


  private:


  int GetNodes(PHCompositeNode* topNode);


  void copyMicromegasClustersToCorrectedMap( );

  Acts::Vector3 getGlobalPosition( TrkrDefs::cluskey key, TrkrCluster* cluster, short int crossing, unsigned int side);


  static constexpr unsigned int _n_mm_layers = 2;


  bool _use_truth_clusters = false;

  TrkrClusterContainer *_cluster_map{nullptr};

  TrkrClusterContainer *_corrected_cluster_map{nullptr};


  TrackSeedContainer *_svtx_seed_map{nullptr};

  TrackSeedContainer *_tpc_track_map{nullptr};

  TrackSeedContainer *_si_track_map{nullptr};


  std::array<double,_n_mm_layers> _rphi_search_win = {0.25, 13.0};


  std::array<double,_n_mm_layers> _z_search_win = {26.0, 0.25};


  // range of TPC layers to use in projection to micromegas

  unsigned int _min_tpc_layer = 38;


  unsigned int _min_mm_layer = 55;


  int _event = -1;


  PHG4CylinderGeomContainer* _geomContainerMicromegas = nullptr;

  TrkrClusterIterationMapv1* _iteration_map = nullptr;

  int _n_iteration = 0;

  //  std::string _track_map_name = "TpcTrackSeedContainer";


  ActsGeometry *_tGeometry = nullptr;


  // crossing z correction

  TpcClusterZCrossingCorrection m_clusterCrossingCorrection;


  // distortion corrections

  TpcDistortionCorrectionContainer* m_dcc_static = nullptr;

  TpcDistortionCorrectionContainer* m_dcc_average = nullptr;

  TpcDistortionCorrectionContainer* m_dcc_fluctuation = nullptr;


  TpcDistortionCorrection m_distortionCorrection;


  bool _test_windows = false;


};


#endif // PHMICROMEGASTPCTRACKMATCHING_H