Paul, I can't check right now but I remember the method has been discussed several times in the forum, using the search function should work.
I think it is also presented by Rogelio on his website (or a similar one).
basically:
-delinearize your luminance
-apply a strong HDRWT, may be two passes, to flat aggressively the image
-invoke the starmask tool, even with the default parameters, it gets all the small stars.
this works for the small medium stars that are usually an issue with the starmask tool.
you can then produce a second starmask for the largest stars using the SM tool or working with wavelets.
last you mix the two mask in Pixelmath with the max(a,b) function.
Andrea