cosmo montecarlo cluster stat hdr

function ds_z = cosmo_montecarlo_cluster_stat(ds, nbrhood, varargin)