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 <trackbase/ActsTrackingGeometry.h>

#include <trackbase/ActsSurfaceMaps.h>


#include <array>

#include <string>

#include <vector>


class TrkrClusterContainer;

class TrkrClusterIterationMapv1;

class SvtxTrackMap;

class AssocInfoContainer;

class PHCompositeNode;

class PHG4CylinderGeomContainer;

class SvtxTrack;

class TrkrCluster;

class TF1;

class TH1;

class TrkrHitSetContainer;


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 set_sc_calib_mode(const bool mode){_sc_calib_mode = mode;}

  void set_collision_rate(const double rate){_collision_rate = rate;}

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

  void set_track_map_name(const std::string &map_name) { _track_map_name = map_name; }


  int InitRun(PHCompositeNode* topNode) override;

  int process_event(PHCompositeNode*) override;

  int End(PHCompositeNode*) override;


  private:


  int GetNodes(PHCompositeNode* topNode);


  void copyMicromegasClustersToCorrectedMap( );


  static constexpr unsigned int _n_mm_layers = 2;


  bool _use_truth_clusters = false;

  TrkrClusterContainer *_cluster_map{nullptr};

  TrkrClusterContainer *_corrected_cluster_map{nullptr};


  SvtxTrackMap *_track_map{nullptr};

  AssocInfoContainer *_assoc_container{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;


  bool _sc_calib_mode = false;


  double _collision_rate = 50e3;


  double _reference_collision_rate = 50e3;


  int _event = -1;


  PHG4CylinderGeomContainer* _geomContainerMicromegas = nullptr;

  TrkrClusterIterationMapv1* _iteration_map = nullptr;

  int _n_iteration = 0;

  std::string _track_map_name = "SvtxTrackMap";


  ActsTrackingGeometry *_tGeometry = nullptr;

  ActsSurfaceMaps *_surfmaps = nullptr;


  TF1 *fdrphi{nullptr};

  double _par0 = -0.36619;

  double _par1 = 0.00375714;


  bool _test_windows = false;


};


#endif // PHMICROMEGASTPCTRACKMATCHING_H