OSMAPOSLParameters :=
; test file for OSMAPOSL with a quadratic prior (and ray tracing projection matrix)

objective function type:= PoissonLogLikelihoodWithLinearModelForMeanAndProjData
PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

input file := Utahscat600k_ca_seg4.hs
zero end planes of segment 0:= 1
; if disabled, defaults to maximum segment number in the file
maximum absolute segment number to process := 3

; change to STIR 2.x default for compatibility 
use subset sensitivities:=0
sensitivity filename:= RPTsens_seg3_PM.hv

projector pair type := Matrix
  Projector Pair Using Matrix Parameters :=
  Matrix type := Ray Tracing
  Ray tracing matrix parameters :=
   number of rays in tangential direction to trace for each bin := 2
;   restrict to cylindrical fov := 0
  End Ray tracing matrix parameters :=
  End Projector Pair Using Matrix Parameters :=

prior type := quadratic
  Quadratic Prior Parameters:=
  penalisation factor := 0.5
  ; next defaults to 0, set to 1 for 2D inverse Euclidean weights, 0 for 3D 
  only 2D:= 0
  END Quadratic Prior Parameters:=

end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

output filename prefix := my_test_image_PM_QP
number of subsets:= 12
start at subset:= 0
number of subiterations:= 6
save estimates at subiteration intervals:= 6
enforce initial positivity condition:=1
initial estimate := my_uniform_image_circular.hv

map model := multiplicative

; for compatibility with STIR 1.x
inter-iteration filter subiteration interval:= 1
inter-iteration filter type := Truncate To Cylindrical FOV
Truncate To Cylindrical FOV Parameters:=
End Truncate To Cylindrical FOV Parameters:=
END :=
