public void validate_static_method_getpositionatcircleperimeter()
{
Random rng = new MathNet.Numerics.Random.MersenneTwister(0);
var pos = SourceToolbox.GetPositionAtCirclePerimeter(_position, _outerRadius, rng);
Assert.Less(Math.Abs(pos.X - _tp[140]), ACCEPTABLE_PRECISION);
Assert.Less(Math.Abs(pos.Y - _tp[141]), ACCEPTABLE_PRECISION);
Assert.Less(Math.Abs(pos.Z - _tp[142]), ACCEPTABLE_PRECISION);
}
/// <summary>
/// Implements Get next photon
/// </summary>
/// <param name="tissue">tissue</param>
/// <returns>photon</returns>
public Photon GetNextPhoton(ITissue tissue)
{
double curved = 2 * Math.PI * _fiberRadius * _fiberHeightZ * _curvedSurfaceEfficiency;
double bottom = Math.PI * _fiberRadius * _fiberRadius * _bottomSurfaceEfficiency;
Direction finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom(
SourceDefaults.DefaultHalfPolarAngleRange.Clone(),
SourceDefaults.DefaultAzimuthalAngleRange.Clone(),
Rng);
Position finalPosition = SourceDefaults.DefaultPosition.Clone();
if (_fiberRadius > 0.0)
{
if (Rng.NextDouble() > bottom / (curved + bottom)) //Consider
{ /* Curved surface */
// To utilize the final direction given above, we can assume a tube
// parallel to the y-axis. We can rotate it about the x-axis by pi/2
// to compute the new direction.
SourceToolbox.UpdateDirectionAfterRotatingAroundXAxis(-0.5 * Math.PI, finalDirection);
//Sample tube perimeter first to compute x and y coordinates
finalPosition = SourceToolbox.GetPositionAtCirclePerimeter(finalPosition,
_fiberRadius,
Rng);
//Sample tube height to compute z coordinate
finalPosition.Z = _fiberHeightZ * (Rng.NextDouble() - 0.5);
}
else
{ /* Bottom Surface */
//Shift finalPosition by _fiberHeightZ / 2
finalPosition = new Position(0.0, 0.0, _fiberHeightZ * 0.5);
//Sample the bottom face to find x, y coordinates of the emission
finalPosition = SourceToolbox.GetPositionInACircleRandomFlat(finalPosition,
0.0,
_fiberRadius,
Rng);
}
}
//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;
}
/// <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(
SourceDefaults.DefaultHalfPolarAngleRange.Clone(),
SourceDefaults.DefaultAzimuthalAngleRange.Clone(),
Rng);
Position finalPosition = SourceDefaults.DefaultPosition.Clone();
if (_tubeRadius > 0.0)
{
// To utilize the final direction given above, we can assume a tube
// parallel to the y-axis. We can rotate it about the x-axis by pi/2
// to compute the new direction.
SourceToolbox.UpdateDirectionAfterRotatingAroundXAxis(0.5 * Math.PI, finalDirection);
//Sample tube perimeter first to compute x and y coordinates
finalPosition = SourceToolbox.GetPositionAtCirclePerimeter(finalPosition,
_tubeRadius,
Rng);
//Sample tube height to compute z coordinate
finalPosition.Z = _tubeHeightZ * (Rng.NextDouble() - 0.5);
}
//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);
}