##############################################################################
#
# PARAMETER FILE FOR THE PHG BINNING MODULE
#
# RUN NAME:
# CREATED:
# OWNER:
#
# This sample file shows how to set up binning for a 3D PET
# image reconstruction like 3drp.
# See the bin_params example file to see all the binning
# parameters with explanations. Here we show only the
# parameters we are using, and explanations of how we
# set them here.
#
##############################################################################
#
# Energy binning in keV
############################
#
# For PET simulations there are two energy indices, one for
# each photon. Thus the 2 energy bins specified here creates
# 4 options: both photons in the lower bin; the first photon
# in the lower bin, the second in the higher bin; the first
# photon in the higher bin, the second in the lower; both
# photons in the higher bin. (Note as there is no real
# distinction between the first and second photon, the middle
# two options are equivalent. SimSET, however, leaves it to
# the user to combine them.)
#
# As the energy indices are the first given in this file,
# they will be the slowest varying. As a result, the output
# will consist of four 3drp data sets, one for each of the
# energy options given above.
#
INT num_e_bins = 2
REAL min_e = 250.0
REAL max_e = 650.0
#
# Trues/Scatter/Randoms binning options
############################
#
# scatter_random_parameter of 0 histograms scatter events
# together with unscattered events.
#
BOOL accept_randoms = false
INT scatter_random_param = 0
INT min_s = 0
INT max_s = 100
#
# 3D-RP binning parameters
############################
#
# The order these parameters are given will determine which
# index is slowest varying, etc., with the slowest varying being
# the index listed first, the fastest the one listed last. Thus
# phi will be the slowest varying index in the output, theta
# next, then Xr, and Yr will be the fastest.
#
# 'phi'
# Range = [0, pi]
INT num_phi_bins = 100
# 'theta'
# min_theta always = -max_theta
INT num_theta_bins = 20
REAL max_theta = 10.0
# 'Xr'
INT num_xr_bins = 100
REAL min_xr = -26.25
REAL max_xr = 26.25
# 'Yr'
INT num_yr_bins = 20
REAL min_yr = -10.0
REAL max_yr = 10.0
#
# Output files specification
############################
# Should the output from this simulation be added to output from
# an earlier run? (Only works if the count-, weight-, and
# weight-squared- files have the same file names.)
BOOL add_to_existing_img = false
# Specification of weight and weight-squared image data type
# Weight weight_image_type = 2 for four byte reals
# Weight weight_image_type = 3 for 8 byte reals
INT weight_image_type = 2
# Specification of count image data type
# Count count_image_type = 0 for one byte integers
# Count count_image_type = 1 for two byte integers
# Count count_image_type = 2 for four byte integers
INT count_image_type = 2
# Path to the binned values of photon weights
STR weight_image_path = "3drp.weights"
# Path to the binned values of photon weights squared.
STR weight_squared_image_path = "3drp.weights.squared"
# Path to the binned values of absolute photons count
STR count_image_path = "3drp.counts"
#
# Run-time precision
############################
# When the following value is set to false (or omitted), the weight
# and weight-squared values are summed in double precision (eight-byte)
# even if single precision (four-byte) output is requested. The output
# will be more accurate if this variable is set to false, but SimSET
# will require more memory.
#
# We suggest this parameter is left false, however, for some simulations
# SimSET's memory requirements may be so great as to require that it be
# be set true.
BOOL sum_according_to_type = false