We use a Hamiltonian dynamics to discuss the statistical mechanics of long-lasting quasistationary states particularly relevant for long-range interacting systems. Despite the presence of an anomalous single-particle velocity distribution, we find that the Central Limit Theorem implies the Boltzmann expression in Gibbs' Γ-space. We identify the nonequilibrium sub-manifold of Γ-space characterizing the anomalous behavior and show that by restricting the Boltzmann-Gibbs approach to this sub-manifold we obtain the statistical mechanics of the quasi-stationary states.PACS numbers: 05.70.Ln, In comparison with its equilibrium counterpart, nonequilibrium statistical mechanics does not rely on universal notions, like the ensembles ones, through which one can handle large classes of physical systems [1]. Incomplete (or partial) equilibrium states [2,3] are in this respect a remarkable exception, since in these cases concepts of equilibrium statistical mechanics can be used to describe nonequilibrium situations. Incomplete equilibrium states arise when different parts of the system themselves reach a state of equilibrium long before they equilibrate with each other [2]. The classical understanding on how a system approaches equilibrium is based on the short time-scale collisions mechanism which links any initial condition to the statistical equilibrium. For longrange interacting systems, this picture is not valid anymore since the time-scale for microscopic collisions diverges with the range of the interactions. This implies that the Boltzmann equation must be substituted with other approximations such as the Vlasov or the BalescuLenard equations [4], where the interparticle correlations are negligible or almost negligible and a nonequilibrium initial configuration could stay frozen or almost frozen for a very long time. This applies, e.g., to gravitational systems, Bose-Einstein condensates and plasma physics [5]. Due to the physical relevance of long-range interacting systems and to the privileged position of incomplete equilibrium states in nonequilibrium statistical mechanics, it is important to investigate whether the notion of incomplete equilibrium plays an important role in understanding the nonequilibrium properties of these systems.Recently we showed [6] that nonequilibrium states in which the value of macroscopic quantities remains stationary or quasi-stationary for reasonably long time (quasi-stationary states -QSSs) are important, e.g., for experiments, since they appear even when the long-range system exchanges energy with a thermal bath (TB). Using the same paradigmatic long-range interacting system of Ref.[6], the Hamiltonian Mean Field (HMF) model [7], here we discuss the Gibbs' Γ-space statistical mechan- * Electronic address: baldovin@pd.infn.it, orlandin@pd.infn.it ics description of the QSSs in a canonical ensemble perspective. We identify the nonequilibrium sub-manifold of Γ-space within which the quasi-stationary dynamics is confined and we show that the Boltzmann-Gibbs (BG) approach, restricted t...