[1]} provided the analytical formulae to calculate the ESD of
the one-halo term from a NFW profile [2]}
that is specified by two free parameters, the halo mass \(M_{\rm h}\) 
and the concentration. We adopt their formula for the ESD with an additional
parameter, \(R_{\rm off}\) , that specifies the projected off-center distance.
Following the model proposed by [3]},
we describe \(R_{\rm off}\)  by a two-dimensional Gaussian
distribution with mean equal to zero and dispersion given by \(\sigma _{\rm off}\) .
To model the two-halo term, we first use
CAMBhttps://camb.info/ (Code for Anisotropies from
Microwave Background) [4]}
and the mcfithttps://github.com/eelregit/mcfit/
package [5]} to obtain the matter correlation
function, \(\xi _{\rm mm}(r)\) . We then use the halo bias model of
[6]} to obtain the bias factor,
\(b_{\rm h}(M_{\rm h})\) , and to calculate the halo-matter cross-correlation
function, \(\xi _{\rm hm}(r)=b_{\rm h}\xi _{\rm mm}(r)\) .
The projected two-halo term is obtained directly from
\(\xi _{\rm hm}\)  [7]}. Finally, the stellar component
is modelled as a point source and the stellar mass parameter is fixed
as the mean value of \(M_*\)  of the galaxy sample.
We refer the reader to [8]} for a detailed description
about the modelling of the three components.
