We present the results of 61 new simulations of nonprecessing spinning black hole binaries with mass ratios q=m1/m2 in the range 1/3≤q≤1 and individual spins covering the parameter space −0.85≤α1,2≤0.85. We additionally perform 10 new simulations of nonspinning black hole binaries with mass ratios covering the range 1/6≤q<1. We follow the evolution for typically the last ten orbits before merger down to the formation of the final remnant black hole. This allows for assessment of the accuracy of our previous empirical formulae for relating the binary parameters to the remnant final black hole mass, spin and recoil. We use the new simulation to improve the fit to the above remnant formulae and add a formula for the peak luminosity of gravitational waves, produced around the merger of the two horizons into one. We find excellent agreement (typical errors ∼0.1−0.2%) for the mass and spin, and within ∼5% for the recoil and peak luminosity. These formulae have direct application to parameter estimation techniques applied to LIGO observations of gravitational waves from binary black hole mergers.