/// <summary> /// Implements Get next photon /// </summary> /// <param name="tissue">tissue</param> /// <returns>photon</returns> public Photon GetNextPhoton(ITissue tissue) { //Source starts from anywhere in the ellipsoid Position finalPosition = GetFinalPositionFromProfileType(_sourceProfile, _aParameter, _bParameter, _cParameter, Rng); // sample angular distribution Direction finalDirection = GetFinalDirection(); //Find the relevent polar and azimuthal pair for the direction PolarAzimuthalAngles _rotationalAnglesOfPrincipalSourceAxis = SourceToolbox.GetPolarAzimuthalPairFromDirection(_newDirectionOfPrincipalSourceAxis); //Rotation and translation SourceToolbox.UpdateDirectionPositionAfterGivenFlags( ref finalPosition, ref finalDirection, _rotationalAnglesOfPrincipalSourceAxis, _translationFromOrigin, _rotationAndTranslationFlags); var photon = new Photon(finalPosition, finalDirection, tissue, _initialTissueRegionIndex, Rng); return(photon); }
/// <summary> /// Implements Get next photon /// </summary> /// <param name="tissue">tissue</param> /// <returns>photon</returns> public Photon GetNextPhoton(ITissue tissue) { //Source starts from anywhere in the line Position finalPosition = GetFinalPosition(_beamDiameterFWHM, _lineLength, Rng); // sample angular distribution Direction finalDirection = GetFinalDirection(finalPosition); //Find the relevent polar and azimuthal pair for the direction PolarAzimuthalAngles _rotationalAnglesOfPrincipalSourceAxis = SourceToolbox.GetPolarAzimuthalPairFromDirection(_newDirectionOfPrincipalSourceAxis); //Rotation and translation SourceToolbox.UpdateDirectionPositionAfterGivenFlags( ref finalPosition, ref finalDirection, _rotationalAnglesOfPrincipalSourceAxis, _translationFromOrigin, _beamRotationFromInwardNormal, _rotationAndTranslationFlags); var photon = new Photon(finalPosition, finalDirection, 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 emitting side String cSide = SelectEmittingSurface( _cubeLengthX, _cubeWidthY, _cubeHeightZ, Rng); //sample angular distribution Direction finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom( _polarAngleEmissionRange, SourceDefaults.DefaultAzimuthalAngleRange.Clone(), Rng); Position tempPosition = null; Position finalPosition = null; switch (cSide) { case "xpos": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeHeightZ, _cubeWidthY, Rng); finalPosition.X = 0.5 * _cubeLengthX; finalPosition.Y = tempPosition.Y; finalPosition.Z = tempPosition.X; finalDirection = SourceToolbox.UpdateDirectionAfterRotatingAroundXAxis(0.5 * Math.PI, finalDirection); break; case "xneg": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeHeightZ, _cubeWidthY, Rng); finalPosition.X = -0.5 * _cubeLengthX; finalPosition.Y = tempPosition.Y; finalPosition.Z = tempPosition.X; finalDirection = SourceToolbox.UpdateDirectionAfterRotatingAroundXAxis(-0.5 * Math.PI, finalDirection); break; case "ypos": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeLengthX, _cubeHeightZ, Rng); finalPosition.X = tempPosition.X; finalPosition.Y = 0.5 * _cubeWidthY; finalPosition.Z = tempPosition.Y; finalDirection = SourceToolbox.UpdateDirectionAfterRotatingAroundYAxis(0.5 * Math.PI, finalDirection); break; case "yneg": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeLengthX, _cubeHeightZ, Rng); finalPosition.X = tempPosition.X; finalPosition.Y = -0.5 * _cubeWidthY; finalPosition.Z = tempPosition.Y; finalDirection = SourceToolbox.UpdateDirectionAfterRotatingAroundYAxis(-0.5 * Math.PI, finalDirection); break; case "zpos": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeLengthX, _cubeWidthY, Rng); finalPosition.X = tempPosition.X; finalPosition.Y = tempPosition.Y; finalPosition.Z = 0.5 * _cubeHeightZ; break; case "zneg": tempPosition = GetFinalPosition(_beamDiameterFWHM, _cubeLengthX, _cubeWidthY, Rng); finalPosition.X = tempPosition.X; finalPosition.Y = tempPosition.Y; finalPosition.Z = -0.5 * _cubeHeightZ; //reverse signs finalDirection.Uz = -finalDirection.Uz; break; } //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, tissue, _initialTissueRegionIndex, Rng); return(photon); }
/// <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; Position finalPosition; if (_fiberRadius > 0.0) { if (Rng.NextDouble() > bottom / (curved + bottom)) { //sample angular distribution finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom( SourceDefaults.DefaultHalfPolarAngleRange.Clone(), SourceDefaults.DefaultAzimuthalAngleRange.Clone(), Rng); //Translate the photon to _tubeRadius length below the origin. Ring lies on yz plane. finalPosition = new Position(0.0, 0.0, _fiberRadius); //Sample a ring that emits photons outside. SourceToolbox.UpdateDirectionPositionAfterRotatingAroundXAxis( 2.0 * Math.PI * Rng.NextDouble(), ref finalDirection, ref finalPosition); //Ring lies on xy plane. z= 0; SourceToolbox.UpdateDirectionPositionAfterRotatingAroundYAxis( 0.5 * Math.PI, ref finalDirection, ref finalPosition); //Sample tube height finalPosition.Z = _fiberHeightZ * (2.0 * Rng.NextDouble() - 1.0); } else { finalPosition = SourceToolbox.GetPositionInACircleRandomFlat( SourceDefaults.DefaultPosition.Clone(), 0.0, _fiberRadius, Rng); finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom( SourceDefaults.DefaultHalfPolarAngleRange.Clone(), SourceDefaults.DefaultAzimuthalAngleRange.Clone(), Rng); } } else { finalPosition = SourceToolbox.GetPositionInALineRandomFlat( SourceDefaults.DefaultPosition.Clone(), _fiberHeightZ, Rng); finalDirection = SourceToolbox.GetDirectionForGivenPolarAzimuthalAngleRangeRandom( SourceDefaults.DefaultFullPolarAngleRange.Clone(), SourceDefaults.DefaultAzimuthalAngleRange.Clone(), Rng); //Rotate 90degrees around y axis SourceToolbox.UpdateDirectionPositionAfterRotatingAroundYAxis( 0.5 * Math.PI, ref finalDirection, ref finalPosition); } //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, tissue, _initialTissueRegionIndex, Rng); return(photon); }