##############################################################################
#
# PARAMETER FILE FOR THE PHG SIMULATION
#
# RUN NAME: randPET
# CREATED: Dec 2008
# OWNER: SimSET
#
# This file gives the binning parameters for
# the fastTest simulation randPET.
#
##############################################################################
BOOL accept_randoms = true
# Scatter/Randoms
# For scatter there is no value to ignore the range, if you don't want to
# restrict the ranges just set the max to something huge like 1000
# and follow instructions below
#
# Set scatterRandomParam = 0 to just control the scatter range, i.e., no
# binning by scatter/random state but photons with number_of_scatters
# outside the min_s and max_s range will be rejected.
# Set scatterRandomParam = 1 to bin all non-scattered events/coincidences in
# one index; events with one or more scatters (either photon for PET)
# in a second. Photons are rejected if their number of scatters do not
# fall in the range from min_s to max_s.
# Set scatterRandomParam = 2, (PET-only) to bin by indices computed for each
# photon within min_s to max_s range, indices ranging from
# (num_scatters - min) to (max - min)
# Set scatterRandomParam = 3, (PET-only) is like 2, but photons with
# number_of_scatters greater than max_s are not rejected. Instead
# their index is computed as if number_of_scatters == max_s.
# Set scatterRandomParam = 4, (PET-only) to bin coincidences using
# (blueScatters + pinkScatters) to compute the index. The
# coincidence is rejected if the sum is < min_s or > max_s.
# Set scatterRandomParam = 5, (PET-only) is the same as 4, except that when
# the sum is > max_s the scatterIndex is set to the max rather than
# the coincidence being rejected.
#
# scatterRandomParam from 6-10 (all PET-only) are the same as 1-5 with the
# addition of an extra bin for randoms. The randoms bin is always the last
# bin.
#
# Set scatterRandomParam = 6, (PET-only) to have three indices: one for
# true unscattered coincidences, the second for non-random scattered
# coincidences, the third for random coincidences.
# Set scatterRandomParam = 7, (PET-only) to have the same binning as 2
# except with an extra bin for randoms, i.e., non-random coincidences
# are binned as in 2, but all random coincidences are binned in an
# extra index added at the end of the scatter state array.
# Set scatterRandomParam = 8, (PET-only) to have the same binning as 3
# except with an extra bin for randoms, i.e., non-random coincidences
# are binned as in 3, but all random coincidences are binned in an
# extra index added at the end of the scatter state array.
# Set scatterRandomParam = 9,(PET-only) to have the same binning as 4
# except with an extra bin for randoms, i.e., non-random coincidences
# are binned as in 4, but all random coincidences are binned in an
# extra index added at the end of the scatter state array.
# Set scatterRandomParam = 10, (PET-only) to have the same binning as 5
# except with an extra bin for randoms, i.e., non-random coincidences
# are binned as in 5, but all random coincidences are binned in an
# extra index added at the end of the scatter state array.
INT scatter_param = 6
INT min_s = 0
INT max_s = 100
# Z axis
# There are more 'Z' options below for SSRB and MSRB PVI binning
INT num_z_bins = 1
REAL min_z = -10.0
REAL max_z = 10.0
# Energy bins
# As with scatters and z-position, energy gets binned into two bins, one
# for each coincidence photon, for PET.
#
# For this simulation, we use only one energy bin. Using one bin will
# cause the energy range given by min_e and max_e to be applied. If we
# set the number of energy bins to 0, the energy range would be ignored.
INT num_e1_bins = 1
INT num_e2_bins = 1
REAL min_e = 400.0
REAL max_e = 650.0
# binning parameters for 3d reconstruction (SimSET does not do the
# reconstruction, just does the binning)
############################
# 3DRP 'phi'
# Range = [0, pi]
INT num_phi_bins = 15
# 3DRP 'theta'
# min_theta always = -max_theta
# Specify in degrees
INT num_theta_bins = 3
REAL max_theta = 10.0
# 3DRP 'Xr'
INT num_xr_bins = 5
REAL min_xr = -26.25
REAL max_xr = 26.25
# 3DRP 'Yr'
INT num_yr_bins = 5
REAL min_yr = -10.0
REAL max_yr = 10.0
# MSRB or SSRB?
BOOL do_msrb = false
BOOL do_ssrb = false
# Detector radius
# Typically the inner radius of the detector
REAL detector_radius = 0.0
#############################
# Specification of weight 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
# Determine internal precision, set to false for double otherwise
# it follows the oputput precision specified above
BOOL sum_according_to_type = false
# Specification of count image data type
# Weight weight_image_type = 0 for one byte integers
# Weight weight_image_type = 1 for two byte integers
# Weight weight_image_type = 2 for four byte integers
INT count_image_type = 2
# Flag to indicate whether or not to add to existing images
BOOL add_to_existing_img = false
# Path to the binned values of photon weights
STR weight_image_path = "randPET.wt"
# Path to the binned values of photon weights squared.
STR weight_squared_image_path = "randPET.wtsq"
# Path to the binned values of absolute photons count
STR count_image_path = "randPET.ct"