text: The Binder cumulant for the average order parameter 
  deﬁned in Eq. (3) was obtained close to the cv peak and is 
  reported in Fig. 3a; its analysis leads to an estimate of the 
  critical temperature of the sample (given by the location 
  of the common crossing point of the diﬀerent curves re- 
  ported in the ﬁgure) of TC (16) = 133.2(5) This value can 
  be considered in a rather good agreement with the exper- 
  imental ordering temperature of Holmium T Ho 
  N , the rel- 
  ative diﬀerence being about 1%. Even such a mismatch 
  between T Ho 
  N and TC (16) could be completely eliminated 
  by slightly adjusting the in-plane coupling constant J0, 
  but, as discussed in Sec. II, we shall preserve the value 
  reported in Refs. 13, and 12 in order to allow for a correct 
  comparison with the results reported in those papers. 
  
  We now move to describe and discuss MC simulation 
  data for thinner samples. A graphical synthesis of the 
  results obtained for n = 8 in reported in Fig. 4a-d. The 
  speciﬁc heat cv, shown in Figs. 4a, reveals very small 
  ﬁnite-size eﬀects, which, however, cannot be unambigu- 
  ously detected for the largest lattice size (L = 64), as 
  they fall comfortably within the error range. Surpris- 
  ingly, the speciﬁc heat maximum is located close to the 
  bulk transition temperature as found for n = 16, and 