BinNormalisationFromECAT8.h

Go to the documentation of this file.

/*
  Copyright (C) 2000-2007, Hammersmith Imanet Ltd
  Copyright (C) 2013-2014, 2020, 2023 University College London
 
  Largely a copy of the ECAT7 version.
 
  SPDX-License-Identifier: Apache-2.0
 
  See STIR/LICENSE.txt for details
*/
 
#ifndef __stir_recon_buildblock_BinNormalisationFromECAT8_H__
#define __stir_recon_buildblock_BinNormalisationFromECAT8_H__
 
#include "stir/recon_buildblock/BinNormalisation.h"
#include "stir/recon_buildblock/BinNormalisationWithCalibration.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/ProjData.h"
#include "stir/shared_ptr.h"
#include "stir/ProjDataInfoCylindricalNoArcCorr.h"
#include "stir/data/SinglesRates.h"
#include "stir/Scanner.h"
#include "stir/Array.h"
#include "stir/IO/stir_ecat_common.h"
#include <string>
 
using std::string;
 
START_NAMESPACE_STIR
START_NAMESPACE_ECAT

class BinNormalisationFromECAT8
    : public RegisteredParsingObject<BinNormalisationFromECAT8, BinNormalisation, BinNormalisationWithCalibration>
{
private:
  using base_type = BinNormalisationWithCalibration;
 
public:
  static const char* const registered_name;


  BinNormalisationFromECAT8();

  BinNormalisationFromECAT8(const string& filename);
 
  Succeeded set_up(const shared_ptr<const ExamInfo>& exam_info_sptr, const shared_ptr<const ProjDataInfo>&) override;
  float get_uncalibrated_bin_efficiency(const Bin& bin) const override;
 
  bool use_detector_efficiencies() const;
  bool use_dead_time() const;
  float get_dead_time_efficiency(int ring, float singles_rate) const;
  float get_dead_time_efficiency(const DetectionPosition<>& det_pos, const double start_time, const double end_time) const;
  bool use_geometric_factors() const;
  bool use_crystal_interference_factors() const;
  bool use_axial_effects_factors() const;
 
  void set_singles_rates(const shared_ptr<SinglesRates>& sptr) { singles_rates_ptr = sptr; }
 
private:
  Array<1, float> rng_dtp_array;
  Array<1, float> rng_dtnp_array;
  Array<1, float> axial_t1_array;
  Array<1, float> axial_t2_array;
  Array<1, float> trans_t1_array;
  shared_ptr<SinglesRates> singles_rates_ptr;
  Array<2, float> geometric_factors;
  Array<2, float> efficiency_factors;
  Array<2, float> crystal_interference_factors;
  Array<1, float> axial_effects;
  Array<2, int> sino_index;
  int num_Siemens_sinograms;
 
  shared_ptr<Scanner> scanner_ptr;
  int num_transaxial_crystals_per_block;
  // TODO move to Scanner
  int num_axial_blocks_per_singles_unit;
  shared_ptr<const ProjDataInfo> proj_data_info_ptr;
  shared_ptr<const ProjDataInfo> norm_proj_data_info_sptr;
  ProjDataInfoCylindricalNoArcCorr const* proj_data_info_cyl_ptr;
  shared_ptr<const ProjDataInfoCylindricalNoArcCorr> proj_data_info_cyl_uncompressed_ptr;
  int mash;
  int num_blocks_per_singles_unit;
  float calib_factor, cross_calib_factor;
 
  bool _use_gaps;
  bool _use_detector_efficiencies;
  bool _use_dead_time;
  bool _use_geometric_factors;
  bool _use_crystal_interference_factors;
  bool _use_axial_effects_factors;
  bool _write_components_to_file;
 
  void read_norm_data(const string& filename);

  void construct_sino_lookup_table();
  float find_axial_effects(int ring1, int ring2) const;
  // parsing stuff
  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;
 
  string normalisation_ECAT8_filename;
  std::string singles_filename;
};
 
END_NAMESPACE_ECAT
END_NAMESPACE_STIR
 
#endif