def METHOD_NAME(q2, wc_obj, par, B, V):
    lep =  ['e', 'mu', 'tau']
    num = bvnunu_fl_num_summed(q2, wc_obj, par, B, V)
    den = sum([ bvnunu_obs(flavio.physics.bdecays.bvll.observables.dGdq2, q2, wc_obj, par, B, V, 'nu'+nu1, 'nu'+nu2)
        for nu1 in lep for nu2 in lep ])
    if den == 0:
        return 0
    else:
        return num/den