ProjDataInfoCylindricalNoArcCorr.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2000- 2011-06-24, Hammersmith Imanet Ltd
    Copyright (C) 2011-07-01 - 2011, Kris Thielemans
    Copyright (C) 2016, University of Hull
    Copyright (C) 2015-2018, 2022, 2024, 2026, University College London
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
#ifndef __stir_ProjDataInfoCylindricalNoArcCorr_H__
#define __stir_ProjDataInfoCylindricalNoArcCorr_H__
 
#include "stir/ProjDataInfoCylindrical.h"
#include "stir/DetectionPositionPair.h"
#include "stir/VectorWithOffset.h"
#include "stir/CartesianCoordinate3D.h"
 
START_NAMESPACE_STIR
 
class Succeeded;
class ProjDataInfoCylindricalNoArcCorr : public ProjDataInfoCylindrical
{
private:
  typedef ProjDataInfoCylindrical base_type;
#ifdef STIR_COMPILING_SWIG_WRAPPER
  // SWIG needs this typedef to be public
public:
#endif
  typedef ProjDataInfoCylindricalNoArcCorr self_type;
 
public:
  ProjDataInfoCylindricalNoArcCorr();

  ProjDataInfoCylindricalNoArcCorr(const shared_ptr<Scanner> scanner_ptr,
                                   const float ring_radius,
                                   const float angular_increment,
                                   const VectorWithOffset<int>& num_axial_pos_per_segment,
                                   const VectorWithOffset<int>& min_ring_diff_v,
                                   const VectorWithOffset<int>& max_ring_diff_v,
                                   const int num_views,
                                   const int num_tangential_poss,
                                   const int tof_mash_factor = 0);


  ProjDataInfoCylindricalNoArcCorr(const shared_ptr<Scanner> scanner_ptr,
                                   const VectorWithOffset<int>& num_axial_pos_per_segment,
                                   const VectorWithOffset<int>& min_ring_diff_v,
                                   const VectorWithOffset<int>& max_ring_diff_v,
                                   const int num_views,
                                   const int num_tangential_poss,
                                   const int tof_mash_factor = 0);
 
  self_type* clone() const override;
 
  bool operator==(const self_type&) const;


  inline float get_s(const Bin&) const override;

  inline float get_angular_increment() const;
 
  std::string parameter_info() const override;



  inline bool
  get_view_tangential_pos_num_for_det_num_pair(int& view_num, int& tang_pos_num, const int det1_num, const int det2_num) const;

  inline void
  get_det_num_pair_for_view_tangential_pos_num(int& det1_num, int& det2_num, const int view_num, const int tang_pos_num) const;


  inline Succeeded get_bin_for_det_pos_pair(Bin&, const DetectionPositionPair<>&) const;


  inline void get_det_pos_pair_for_bin(DetectionPositionPair<>&, const Bin&) const;


  unsigned int get_num_det_pos_pairs_for_bin(const Bin&, bool ignore_non_spatial_dimensions = true) const;


  void get_all_det_pos_pairs_for_bin(std::vector<DetectionPositionPair<>>&,
                                     const Bin&,
                                     bool ignore_non_spatial_dimensions = true) const;
 
private:
  // old function, now private. Use get_bin_for_det_pos_pair instead.

  inline Succeeded get_bin_for_det_pair(Bin&,
                                        const int det1_num,
                                        const int ring1_num,
                                        const int det2_num,
                                        const int ring2_num,
                                        const int mashed_timing_pos_num = 0) const;
 
public:

 
  inline void get_det_pair_for_bin(int& det_num1, int& ring_num1, int& det_num2, int& ring_num2, const Bin& bin) const;

 
  Bin get_bin(const LOR<float>&, const double delta_time = 0.0) const override;



  virtual Succeeded find_scanner_coordinates_given_cartesian_coordinates(int& det1,
                                                                         int& det2,
                                                                         int& ring1,
                                                                         int& ring2,
                                                                         const CartesianCoordinate3D<float>& c1,
                                                                         const CartesianCoordinate3D<float>& c2) const;
 
  void find_cartesian_coordinates_of_detection(CartesianCoordinate3D<float>& coord_1,
                                               CartesianCoordinate3D<float>& coord_2,
                                               const Bin& bin) const;
 
protected:
  // give access to protected member find_cartesian_coordinates_given_scanner_coordinates
  friend class ProjDataInfoCylindricalNoArcCorrTests;
  template <typename ProjDataInfoT>
  friend class CListEventScannerWithDiscreteDetectors;
  friend class LmToProjDataWithMCTests;


  virtual void find_cartesian_coordinates_given_scanner_coordinates(CartesianCoordinate3D<float>& coord_1,
                                                                    CartesianCoordinate3D<float>& coord_2,
                                                                    const int Ring_A,
                                                                    const int Ring_B,
                                                                    const int det1,
                                                                    const int det2,
                                                                    const int timing_pos_num) const;
 
public:
  void find_bin_given_cartesian_coordinates_of_detection(Bin& bin,
                                                         const CartesianCoordinate3D<float>& coord_1,
                                                         const CartesianCoordinate3D<float>& coord_2) const;
 
private:
  float ring_radius;
  float angular_increment;

  float get_psi_offset() const;
 
  // used in get_view_tangential_pos_num_for_det_num_pair()
  struct Det1Det2
  {
    int det1_num;
    int det2_num;
  };
  mutable VectorWithOffset<VectorWithOffset<Det1Det2>> uncompressed_view_tangpos_to_det1det2;
  mutable bool uncompressed_view_tangpos_to_det1det2_initialised;
  void initialise_uncompressed_view_tangpos_to_det1det2() const;
 
  // used in get_view_tangential_pos_num_for_det_num_pair()
  // we prestore a lookup-table in terms for unmashed view/tangpos
  struct ViewTangPosSwap
  {
    int view_num;
    int tang_pos_num;
    bool swap_detectors;
  };
  mutable VectorWithOffset<VectorWithOffset<ViewTangPosSwap>> det1det2_to_uncompressed_view_tangpos;
  mutable bool det1det2_to_uncompressed_view_tangpos_initialised;
  void initialise_det1det2_to_uncompressed_view_tangpos() const;

  inline void initialise_uncompressed_view_tangpos_to_det1det2_if_not_done_yet() const;
  inline void initialise_det1det2_to_uncompressed_view_tangpos_if_not_done_yet() const;
 
protected:
  bool blindly_equals(const root_type* const) const override;
};
 
END_NAMESPACE_STIR
 
#include "stir/ProjDataInfoCylindricalNoArcCorr.inl"
 
#endif