I made an estemator for the star star correlation function but I did something wrong. I searched for a star before running the correleation and never normalized taking into accound that the first point chosen was always a star. Need to fix this. I also got some order parameters runing: number of stars, number of flipable plaquets, numer of vertical dimers, number of horizontal dimers.