The peak power (or luminosity) is roughly proportional to $\nu^2$, where
$$\nu = \frac{m_1 m_2}{(m_1+m_2)^2}$$
is the symmetric mass-ratio (This follows from the fact that the luminosity in general is proportional to $\nu^2$.). So, yes, the maximal peak power is reach for equal mass-mergers.
In general, the peak power will depend on the spins (and their orientation), eccentricity, (and if you wish to consider it charges) of the merging binaries. However, there is no way to establish these dependencies analytically. Instead we have to rely on numerical relativity simulations. For example, this paper, studies the relation between mass-ratio, spin, and peak luminosity for spin aligned binaries. They find that the maximal peak luminosity seems to be attained in the equal mass case with both spins maximal and aligned, reaching roughly 0.002 in natural units (about a factor two higher than the non-spinning case). The intuitive explanation for this is that aligned spins allow the inspiral phase to continue to much smaller separations, leading to a higher peak luminosity. Following this intuition any misalignment of the spins should result in a lower peak luminosity. A similar intuition applies to the charges case. Adding similar charges to the black holes delays the merger and should increase the peak luminosity. I'm not aware of any systematic NR campaigns that looked specifically at this, though.
People have only fairly recently started exploring the parameter space of eccentric mergers with NR simulations with no (published) reports on the peak luminosity. However, the dependence of the peak luminosity on eccentricity is likely to be fairly weak, because most inspirals will shed most of their eccentricity before merger. Intuitively, one should expect a slight increase in the peak luminosity though.