def METHOD_NAME(self, xmin=-1., xmax=1., ymin=-1., ymax=1.,
                      zmin=-1., zmax=1, xpmin=-1., xpmax=1.,
                      ypmin=-1., ypmax=1., dpmin=-1., dpmax=1.):
    x = np.random.uniform(xmin, xmax, self.macroparticlenumber)
    y = np.random.uniform(ymin, ymax, self.macroparticlenumber)
    z = np.random.uniform(zmin, zmax, self.macroparticlenumber)
    xp = np.random.uniform(xpmin, xpmax, self.macroparticlenumber)
    yp = np.random.uniform(ypmin, ypmax, self.macroparticlenumber)
    dp = np.random.uniform(dpmin, dpmax, self.macroparticlenumber)
    coords_n_momenta_dict = {
        'x': x, 'y': y, 'z': z,
        'xp': xp, 'yp': yp, 'dp': dp
    }
    return Particles(
        self.macroparticlenumber, self.particlenumber_per_mp, e, m_p,
        self.circumference, self.gamma, coords_n_momenta_dict
    )