Our topic concerns a long standing puzzle: the energy of gravitating systems. More precisely we want to consider, for gravitating systems, how to best describe energy-momentum and angular momentum/center-of-mass momentum (CoMM). It is known that these quantities cannot be given by a local density. The modern understanding is that (i) they are quasi-local (associated with a closed 2-surface), (ii) they have no unique formula, (iii) they have no reference frame independent description. In the first part of this work we review some early history, much of it not so well known, on the subject of gravitational energy in Einstein's general relativity (GR), noting especially Noether's contribution. In the second part we review (including some new results) much of our covariant Hamiltonian formalism and apply it to Poincaré gauge theories of gravity (PG), with GR as a special case. The key point is that the Hamiltonian boundary term has two roles, it determines the quasi-local quantities, and, furthermore it determines the boundary conditions for the dynamical variables. Energy-momentum and angular momentum/CoMM are associated with the geometric symmetries under Poincaré transformations. They are best described in a local Poincaré gauge theory. The type of spacetime that naturally has this symmetry is Riemann-Cartan spacetime, with a metric compatible connection having, in general, both curvature and torsion. Thus our expression for the energy-momentum of physical systems is obtained via our covariant Hamiltonian formulation applied to the PG. d One reason was his famous "hole" argument. 61,98 e Einstein and Hilbert had quite different agendas; 112,115,136,141 Hilbert in his Foundation of Physics papers, based on the work of Einstein and Mie, was using his axiomatic method with the objective of finding a unified field theory of gravity and electromagnetism. 15,28,58,108,109,114 July 28, 2015 0:26 World Scientific Review Volume -9.75in x 6.5in 100GR˙Ch4˙150726v1 page 5Gravitational energy for GR and Poincaré gauge theories:a covariant Hamiltonian approach 5 l The sign in this expression is dictated by the condition for positive energy determined by the Hamiltonian using our local Minkowski signature convention: Pµ = (−E/c, p).