packages/tpc/TpcClusterMover.h

0001 #ifndef TPC_TPCCLUSTERMOVER_H
0002 #define TPC_TPCCLUSTERMOVER_H
0003
0004 /*!
0005  * \file TpcClusterMover.h
0006  * \Moves TPC clusters to the readout TPC surface after distortion corrections
0007  * \author Tony Frawley, May 2022
0008  */
0009 #include <trackbase/TrkrDefs.h>
0010 #include <trackbase_historic/ActsTransformations.h>
0011 #include <vector>
0012
0013 class PHG4TpcCylinderGeomContainer;
0014
0015 class TpcClusterMover
0016 {
0017  public:
0018   //! constructor
0019   TpcClusterMover();
0020
0021   void set_verbosity(int verb) { _verbosity = verb; }
0022
0023   std::vector<std::pair<TrkrDefs::cluskey, Acts::Vector3>> processTrack(const std::vector<std::pair<TrkrDefs::cluskey, Acts::Vector3>>& global_in);
0024
0025   //! Updates the assumed default geometry below to that contained in the
0026   //! cell geo
0027   void initialize_geometry(PHG4TpcCylinderGeomContainer *cellgeo);
0028
0029  private:
0030   int get_circle_circle_intersection(double target_radius, double R, double X0, double Y0, double xclus, double yclus, double &x, double &y);
0031
0032   double _z_start = 0.0;
0033   double _y_start = 0.0;
0034   double _x_start = 0.0;
0035
0036   double _z_proj = 0.0;
0037   double _y_proj = 0.0;
0038   double _x_proj = 0.0;
0039
0040   double layer_radius[48] = {0};
0041   double inner_tpc_min_radius = 30.0;
0042   double mid_tpc_min_radius = 40.0;
0043   double outer_tpc_min_radius = 60.0;
0044   double outer_tpc_max_radius = 76.4;
0045
0046   double inner_tpc_spacing = 0.0;
0047   double mid_tpc_spacing = 0.0;
0048   double outer_tpc_spacing = 0.0;
0049
0050   int _verbosity = 0;
0051 };
0052
0053 #endif