// ignore DDP error

sqrt( (N1+N2+N3)*k*k + (M1+M2+M3)*p*p )

N1,2,3 = N_gj in bins 1,2,3
M1,2,3 = M_jj in bins 1,2,3

k = N_gj/N_gj_in_gj_mass_template
 
p = N_jj/N_jj_in_jj_mass_template