Same as for partcile solid angle, take into account * plasma domain sampling (parallelize on R, sum on phi) * LOS blocked (intersection detection)
Same as for partcile solid angle, take into account