One technique seems to be to calculate Q^2 for all observed reflections, then assume the lowest reflections are (100), (110), (111), (200), ... and see what matches.
If your lattice angles are all 90 deg, then Q^2 = H^2 a*^2 + K^2 b*^2 + L^2 c*^2. For hexagonal this also seems to work quite well. For monoclinic and triclinic this is more of a mess.
I guess in a computer programme one could one by one go through the 14 Bravais lattices, add a few assumptions to each one (like a>c or a<c for tetragonal, as this could invert the order of peaks), and calculate the lattice parameters from the first few reflections.
Note that the space group may lead to systematic absence of certain peaks, so a better critirium for fit/no fit is if there are peaks in the data that have no corresponding HKL.