The knowledge of the structure of a financial network gives valuable information for the estimation of systemic risk and other important properties. However, since financial data are typically subject to confidentiality, network reconstruction techniques become necessary to infer both the presence of connections (i.e. the topology) and their intensity (i.e. the link weights) from partial information. Recently, various horse races have been conducted to compare the performance of these methods. These comparisons were based on arbitrarily chosen metrics of network similarity. Here we establish a generalised likelihood approach to the rigorous definition, estimation and comparison of methods of reconstruction of weighted networks. The crucial ingredient is the possibility of separately constraining the topology and the link weights to follow certain "tight" empirical relationships that connect them to the observed marginals. We find that a statistically rigorous estimation of the probability of each link weight should incorporate the whole distribution of reconstructed network topologies, and we provide a general procedure to do this. Our results indicate that the best method is obtained by "dressing" the best-performing method available to date with a certain exponential distribution of link weights. The method is fast, unbiased and accurate and reproduces the empirical networks with highest generalised likelihood.Network reconstruction is an active field of research within the broader field of complex networks [20]. Addressing the network reconstruction problem means facing the double challenge represented by the estimation of topology and link weights. The task at hand consists in determining both binary and weighted ensemble distributions, and to understand the interplay between them. Among the methods proposed so far, some assume that the binary and weighted constraints jointly determine the final configuration in terms of both topology and weights while others attribute weights to the binary configuration using a completely separate methodology [17,21]. Amidst the former ones, a special mention is deserved by the Enhanced Configuration Model [15] . This is defined by simultaneously constraining the degrees and the strengths of nodes which jointly affect the estimation of the two sets of quantities, the linkage probabilities and the weight estimates. Since these are jointly determined on the basis of the same information (i.e. constraints), this implies the impossibility to include purely topological additional information. Examples of algorithms belonging to the second group are those iteratively adjusting the link weights (e.g. via the RAS recipe [19]) on top of some previously-determined topological structure, in such a way to satisfy the constraints concerning strengths a posteriori. This approach has encountered critiques in [16]. It is important to notice that this kind of procedure assigns weights deterministically, and therefore the likelihood of observing any real matrix is exactly zero, assuming conti...