OSMAPOSLParameters :=

    objective function type:= PoissonLogLikelihoodWithLinearModelForMeanAndProjData
    PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

        input file := input.hs

        projector pair type := Matrix
            Projector Pair Using Matrix Parameters :=
            Matrix type := SPECT UB
            Projection Matrix By Bin SPECT UB Parameters:=
            
                maximum number of sigmas:= 2.0

                ;PSF type of correction { 2D // 3D // Geometrical }
                psf type:= 3D
                ; next 2 parameters define the PSF. They are ignored if psf_type is "Geometrical"
                ; the PSF is modelled as a Gaussian with sigma dependent on the distance from the collimator
                ; sigma_at_depth = collimator_slope * depth_in_cm + collimator sigma 0(cm)
                collimator slope := 0.0163
                collimator sigma 0(cm) := 0.1466

                ;Attenuation correction { Simple // Full // No }
                attenuation type := Simple  
                ;Values in attenuation map in cm-1
                attenuation map := attMapRec.hv

                ;Mask properties { Cylinder // Attenuation Map // Explicit Mask // No}
                mask type := Explicit Mask
                mask file := mask.hv

                keep all views in cache := ${CACHEALLVIEWS3D}

            End Projection Matrix By Bin SPECT UB Parameters:=

    End Projector Pair Using Matrix Parameters :=

    end PoissonLogLikelihoodWithLinearModelForMeanAndProjData Parameters:=

    initial estimate:= init.hv
    output filename prefix := out/OSEM_3DPSF

    number of subsets:= 12
    number of subiterations:= 4
    Save estimates at subiteration intervals:= 4


END :=
