remove calculation of geometric depths when not needed for particle e…

…ngines
e0404 · Dec 26, 2023 · 8ee2b25 · 8ee2b25
1 parent 3d218c1
commit 8ee2b25
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/doseCalc/+DoseEngines/matRad_ParticleFineSamplingPencilBeamEngine.m b/doseCalc/+DoseEngines/matRad_ParticleFineSamplingPencilBeamEngine.m
@@ -78,8 +78,8 @@ function setDefaults(this)
                 currBeam.ixRadDepths, ...
                 lateralRayCutOff);
 
-            %ray.geoDepths = currBeam.geoDepths{1}(ix); %we don't need that
-            ray.radDepths = currBeam.radDepths{1}(ix);
+            %ray.geoDepths = cellfun(@(rD) rD(ix),currBeam.geoDepths,'UniformOutput',false); %Not needed for this engine
+            ray.radDepths = cellfun(@(rD) rD(ix),currBeam.radDepths,'UniformOutput',false);
             ray.ix = currBeam.ixRadDepths(ix);
             ray.subIxVdoseGrid = currBeam.subIxVdoseGrid(ix);
         end

diff --git a/doseCalc/+DoseEngines/matRad_ParticlePencilBeamEngineAbstract.m b/doseCalc/+DoseEngines/matRad_ParticlePencilBeamEngineAbstract.m
@@ -237,7 +237,7 @@ function chooseLateralModel(this)
                 currBeam.ixRadDepths, ...
                 lateralRayCutOff);
 
-            ray.geoDepths = cellfun(@(rD) rD(ix),currBeam.geoDepths,'UniformOutput',false);
+            %ray.geoDepths = cellfun(@(rD) rD(ix),currBeam.geoDepths,'UniformOutput',false); %usually not needed for particle beams
             ray.radDepths = cellfun(@(rD) rD(ix),currBeam.radDepths,'UniformOutput',false);
             ray.ix = currBeam.ixRadDepths(ix);
             ray.subIxVdoseGrid = currBeam.subIxVdoseGrid(ix);