Close issue #20, overall performance increases 

Changes sampling to method outlined in issue #20.
Overall reduce calculation time for test case by 40%.
20% faster for Monte Carlo npa. 
Monte Carlo speed depends on the fraction of the overlap of the npa diagnostic volume and beam volume. Smaller fraction = faster calculation compared to old sampling method.