this cannot be always correct. It should depend on the forward operator and/or inverse option if the spatial density is not in the forward fields but in an extra structure.