ProjDataInfoGenericNoArcCorr.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) 2017, ETH Zurich, Institute of Particle Physics and Astrophysics
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/

#ifndef __stir_ProjDataInfoGenericNoArcCorr_H__
#define __stir_ProjDataInfoGenericNoArcCorr_H__

#include "stir/ProjDataInfoGeneric.h"
#include "stir/GeometryBlocksOnCylindrical.h"
#include "stir/DetectionPositionPair.h"
#include "stir/VectorWithOffset.h"
#include "stir/CartesianCoordinate3D.h"

START_NAMESPACE_STIR

class Succeeded;
class ProjDataInfoGenericNoArcCorr : public ProjDataInfoGeneric
{
private:
  typedef ProjDataInfoGeneric base_type;
#ifdef SWIG
  // SWIG needs this typedef to be public
 public:
#endif
  typedef ProjDataInfoGenericNoArcCorr self_type;

public:
  ProjDataInfoGenericNoArcCorr();

  ProjDataInfoGenericNoArcCorr(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);

  ProjDataInfo* clone() const override;

  bool operator==(const self_type&) const;


  inline float get_s(const Bin&) const override;

  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&) const;


  void get_all_det_pos_pairs_for_bin(std::vector<DetectionPositionPair<>>&, const Bin&) 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;

public:

  inline void get_det_pair_for_bin(int& det1_num, int& ring1_num, int& det2_num, int& ring2_num, const Bin&) const;


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



  void find_cartesian_coordinates_of_detection(CartesianCoordinate3D<float>& coord_1,
                                               CartesianCoordinate3D<float>& coord_2,
                                               const Bin& bin) const override;

  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;

private:
  // 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/ProjDataInfoGenericNoArcCorr.inl"

#endif