OSMAPOSLParameters :=

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
  End Ray tracing matrix parameters :=
  End Projector Pair Using Matrix Parameters :=

prior type := FilterRootPrior
FilterRootPrior Parameters:=
  penalisation factor := .2
  Filter type:= median
  Median Filter Parameters :=
  mask radius x := 1   
  mask radius y := 1
  mask radius z := 0
  End Median Filter Parameters:=
END FilterRootPrior Parameters:=

end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

output filename prefix := my_test_image_PM_MRP

number of subsets:= 12
start at subset:= 0
number of subiterations:= 6
save estimates at subiteration intervals:= 6


; you normally wouldn't use a filter and a prior, but this is just for testing
inter-iteration filter subiteration interval:= 2
inter-iteration filter type := Chained Data Processor
Chained Data Processor Parameters:=
Data Processor to apply first:= Chained Data Processor
Chained Data Processor Parameters:=
Data Processor to apply first:= Truncate To Cylindrical FOV
Truncate To Cylindrical FOV Parameters:=
strictly_less_than_radius:=1
End Truncate To Cylindrical FOV Parameters:=


Data Processor to apply second:=Separable Cartesian Metz
Separable Cartesian Metz Filter Parameters :=
x-dir filter FWHM (in mm):= 6.0
y-dir filter FWHM (in mm):= 5.0
z-dir filter FWHM (in mm):= 4.0
x-dir filter Metz power:= 2.0
y-dir filter Metz power:= 1.0
z-dir filter Metz power:= 0.0
x-dir maximum kernel size := 13
y-dir maximum kernel size := 13
z-dir maximum kernel size := 31
END Separable Cartesian Metz Filter Parameters :=

END Chained Data Processor Parameters:=

Data Processor to apply second:=Truncate To Cylindrical FOV
Truncate To Cylindrical FOV Parameters:=
strictly_less_than_radius:=0
End Truncate To Cylindrical FOV Parameters:=

END Chained Data Processor Parameters:=

start at subiteration number:=2
enforce initial positivity condition:=0
initial estimate := my_uniform_image.hv

map model := multiplicative
END :=
