r.m 
% index_maker.m> as specified to get the necessary subset indices. The line |[gaussK, 
% options] = PSFKernel(options);| is only applicable if PSF reconstruction is 
% used as it creates the PSF kernel gaussK.
% 
% To use your own measurement data, replace |load('Cylindrical_PET_example_cylpet_example_sinograms_combined_static_200x168x703_span3.mat','raw_SinM'|) 
% with your own custom load or simply run the code up to that and manually load 
% your data. You should either name your measurement data as