  Spacecraft missions provide the unique opportunity to study the properties of
collisionless shocks utilising in situ measurements. In the past years, several
diagnostics have been developed to address key shock parameters using time
series of magnetic field (and plasma) data collected by a single spacecraft
crossing a shock front. A critical aspect of such diagnostics is the averaging
process involved in the evaluation of upstream-downstream quantities. In this
work, we discuss several of these techniques, with a particular focus on the
shock obliquity (defined as the angle between the upstream magnetic field and
the shock normal vector) estimation. We introduce a systematic variation of the
upstream/downstream averaging windows, yielding to an ensemble of shock
parameters, a useful tool to address the robustness of their estimation. This
approach is first tested with a synthetic shock, compliant with the
Rankine-Hugoniot jump conditions for a shock, including the presence of noise
and disturbances. We then employ self-consistent, hybrid kinetic shock
simulations to apply the diagnostics to virtual spacecraft crossing the shock
front at various stages of its evolution, highlighting the role of
shock-induced fluctuations in the parameters estimation. This approach has the
strong advantage of retaining some important properties of collisionless shock
while being able to set a known, nominal set of shock parameters. Finally, two
recent observations of interplanetary shocks from the Solar Orbiter spacecraft
are presented, and the approach is also tested on an interplanetary shock
measured by the four spacecraft of the Magnetospheric Multiscale (MMS) mission.
All the Python software developed and used for the diagnostics (SerPyShock) is
made available for the public, including an example of parameter estimation for
a shock wave recently observed in-situ by the Solar Orbiter spacecraft.
