Please correct me if I’m wrong, but my understanding is that using the EnergyFunctionFilter with dose_coefficients does give units of [pSv-cm^3/source_particle], since the flux tally is given in [particles-cm/source particle] and the dose_coefficients are in [pSv*cm^2].
So to get [pSv/second] you would need to multiply but the source strength in [particles/second] and also divide by the cell tally volume (or mesh voxel volume for mesh tallies) in [cm^3]. Is this correct?