OSMAPOSL_test_sim_PM_QP.par

OSMAPOSLParameters :=
; test file for OSMAPOSL with a quadratic prior (and ray tracing projection matrix)
objective function type:= PoissonLogLikelihoodWithLinearModelForMeanAndProjData
PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

input file := my_prompts${suffix}.hs
zero end planes of segment 0:= 0 ; segment 0 is has direct and indirect planes
; if disabled, defaults to maximum segment number in the file
;maximum absolute segment number to process := 3

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 := 10
;   do symmetry 90degrees min phi := 0
  End Ray tracing matrix parameters :=
  End Projector Pair Using Matrix Parameters :=

  Bin Normalisation type := from projdata
    Bin Normalisation From ProjData :=
    normalisation projdata filename:= my_acfs${suffix}.hs
    End Bin Normalisation From ProjData:=

additive sinogram := my_additive_sinogram${suffix}.hs

xy output image size (in pixels) := 91
zoom := .5

use subset sensitivities:=1
subset sensitivity filenames:= my_test_sim${suffix}_sens_PM_s${num_subsets}_%d.hv
recompute_sensitivity:=1

prior type := quadratic
  Quadratic Prior Parameters:=
  penalisation factor := 1
  ;weights:={{{0,1,0},{1,0,1},{0,1,0}}}
  END Quadratic Prior Parameters:=

end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

output filename prefix := my_test_sim${suffix}_image_OSMAPOSL_PM_QP
number of subsets:= ${num_subsets}
start at subset:= 0
number of subiterations:= 28
save estimates at subiteration intervals:= ${num_subsets}

;report objective function values interval := 1

END :=