Weighted Average
================
The code further provides an interface to calculate the weighted average of any customized value over the blast surface, where the weight is :math:`dL/d\Omega`. This utility is provided because most observables of GRB afterglows are fundimentally weighted average over the blast surface.

Similar to the way you define your radiation model, you can define it in the Python way or in the C++ way.