LmToProjDataWithRandomRejection.h

Go to the documentation of this file.

//
//
/*
    Copyright (C) 2003- 2012, Hammersmith Imanet Ltd
    Copyright (C) 2019, National Physical Laboratory
    Copyright (C) 2019, 2021, University College London
    This file is part of STIR.

    SPDX-License-Identifier: Apache-2.0

    See STIR/LICENSE.txt for details
*/

#ifndef __stir_listmode_LmToProjDataWithRandomRejection_H__
#define __stir_listmode_LmToProjDataWithRandomRejection_H__

#include "stir/listmode/LmToProjData.h"
#include <boost/random/uniform_01.hpp>
#include <boost/random/mersenne_twister.hpp>

START_NAMESPACE_STIR

template <typename LmToProjDataT>
class LmToProjDataWithRandomRejection : public LmToProjDataT
{

public:
  LmToProjDataWithRandomRejection(const char* const par_filename);

  LmToProjDataWithRandomRejection(const char* const par_filename, const unsigned int seed);

  // void set_seed(const unsigned int seed);
  float set_reject_if_above(const float);

  Succeeded set_up() override;

protected:

  void start_new_time_frame(const unsigned int new_frame_num) override;

  void get_bin_from_event(Bin& bin, const ListEvent&) const override;

  // \name parsing variables

  unsigned int seed;
  float reject_if_above;

private:
  typedef LmToProjDataT base_type;
  typedef boost::mt19937 random_generator_type;
  random_generator_type random_generator;

  void set_defaults() override;
  void initialise_keymap() override;
  bool post_processing() override;
};

END_NAMESPACE_STIR

#endif