Implement Marginalize() for templates

[mjlarson]

Need to add this in at some point. For a full map convolution, it's trivial: you just sum over each healpix ring. It's less obvious to me how to handle it if we're just evaluating single points.