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
    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 SWIG
  // 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);

  ProjDataInfo* 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) const override;



  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;

  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;

  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;

  bool blindly_equals(const root_type* const) const override;
};

END_NAMESPACE_STIR

#include "stir/ProjDataInfoCylindricalNoArcCorr.inl"

#endif