  The nitrogen-vacancy (NV) center in diamond is a powerful and versatile
quantum sensor for diverse quantities. In particular, relaxometry (or T1),
allows to detect magnetic noise at the nanoscale. While increasing the number
of NV centers in a nanodiamond allows to collect more signal, a standardized
method to extract information from relaxometry experiments of such NV ensembles
is still missing. In this article, we use T1 relaxation curves acquired at
different concentrations of gadolinium ions to calibrate and optimize the
entire data processing flow, from the acquired raw data to the extracted T1. In
particular, we use a bootstrap to derive a signal to noise ratio (SNR) that can
be quantitatively compared from one method to another. At first, T1 curves are
extracted from photoluminescence pulses. We compare integrating their signal
through an optimized window as performed conventionally, to fitting a known
function on it. Fitting the decaying T1 curves allows to obtain the relevant T1
value. We compared here the three most commonly used fit models that are,
single, bi, and stretched-exponential. We finally investigated the effect of
the bootstrap itself on the precision of the result as well as the use of a
rolling window to allows time-resolution.
