BinNormalisationFromPETSIRD.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2025, 2026, University Medical Center Groningen
    Copyright (C) 2025, 2026, London South Bank University
    This file is part of STIR.
 
    SPDX-License-Identifier: Apache-2.0
 
    See STIR/LICENSE.txt for details
*/
 
#ifndef __stir_recon_buildblock_BinNormalisationFromPETSIRD_H__
#define __stir_recon_buildblock_BinNormalisationFromPETSIRD_H__
 
#include "stir/recon_buildblock/BinNormalisation.h"
#include "stir/recon_buildblock/BinNormalisationWithCalibration.h"
#include "stir/RegisteredParsingObject.h"
#include "stir/shared_ptr.h"
#include "stir/PETSIRDInfo.h"
#include "stir/Array.h"
 
using std::string;
 
START_NAMESPACE_STIR
 
class BinNormalisationFromPETSIRD
    : public RegisteredParsingObject<BinNormalisationFromPETSIRD, BinNormalisation, BinNormalisationWithCalibration>
{
private:
  using base_type = BinNormalisationWithCalibration;
 
public:
  static const char* const registered_name;
 
  BinNormalisationFromPETSIRD();
 
  BinNormalisationFromPETSIRD(const std::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;
 
  inline bool with_detector_efficiencies() const { return m_with_detector_efficiencies; }
  inline bool use_dead_time() const { return m_use_dead_time; }
  inline bool with_geometric_factors() const { return m_with_geometric_factors; }
 
private:
  bool use_hdf5;
 
  void set_defaults() override;
 
  void initialise_keymap() override;
 
  bool post_processing() override;
 
  void read_norm_data(const string& filename);
 
  string normalisation_filename;

  bool m_with_detector_efficiencies;
  bool m_use_dead_time;
  bool m_with_geometric_factors;
  //   shared_ptr<PETSIRDInfo> petsird_info_sptr;
  shared_ptr<petsird::PETSIRDReaderBase> petsird_data_sptr;
 
  shared_ptr<petsird::ScannerInformation> scanner_info_sptr;
 
  shared_ptr<const PETSIRDInfo> petsird_info_sptr;
  // dead time blocks copied as stored in the PETSIRD file
  std::vector<petsird::DeadTimeTimeBlock> dead_time_blocks;
  float get_dead_time_efficiency(const DetectionPosition<>& det_pos, const double start_time, const double end_time) const;
};
 
END_NAMESPACE_STIR
 
#endif