/// <summary>
/// Implements Get next photon
/// </summary>
/// <param name="tissue">tissue</param>
/// <returns>photon</returns>
public Photon GetNextPhoton(ITissue tissue)
{
// sample angular distribution
Direction finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom(_polarAngleRangeToDefineSphericalSurface, _azimuthalAngleRangeToDefineSphericalSurface, Rng);
//Source starts from anywhere in the sphere
Position finalPosition = GetFinalPositionFromProfileType(finalDirection, _radius, Rng);
//Lambertian distribution (uniform hemispherical distribution)
PolarAzimuthalAngles polarAzimuthalPair = SourceToolbox.GetPolarAzimuthalPairForGivenAngleRangeRandom(
SourceDefaults.DefaultHalfPolarAngleRange.Clone(),
SourceDefaults.DefaultAzimuthalAngleRange.Clone(),
Rng);
//Avoid updating the finalDirection during following rotation
Position dummyPosition = finalPosition;
//Rotate polar azimutahl angle by polarAzimuthalPair vector
SourceToolbox.UpdateDirectionPositionAfterRotatingByGivenAnglePair(polarAzimuthalPair, ref finalDirection, ref dummyPosition);
//Find the relevent polar and azimuthal pair for the direction
PolarAzimuthalAngles _rotationalAnglesOfPrincipalSourceAxis = SourceToolbox.GetPolarAzimuthalPairFromDirection(_newDirectionOfPrincipalSourceAxis);
//Translation and source rotation
SourceToolbox.UpdateDirectionPositionAfterGivenFlags(
ref finalPosition,
ref finalDirection,
_rotationalAnglesOfPrincipalSourceAxis,
_translationFromOrigin,
_rotationAndTranslationFlags);
var photon = new Photon(finalPosition, finalDirection, 1.0, tissue, _initialTissueRegionIndex, Rng);
return(photon);
}
public void validate_static_method_getpolarazimuthalpairforgivenanglerangerandom()
{
Random rng = new MathNet.Numerics.Random.MersenneTwister(0);
var angPair = SourceToolbox.GetPolarAzimuthalPairForGivenAngleRangeRandom(_polRange, _aziRange, rng);
Assert.Less(Math.Abs(angPair.Theta - _tp[43]), ACCEPTABLE_PRECISION);
Assert.Less(Math.Abs(angPair.Phi - _tp[44]), ACCEPTABLE_PRECISION);
}