From 43755f9c304717927ac49cb154cc9a32d33ec18d Mon Sep 17 00:00:00 2001 From: Greg Rischbieter Date: Thu, 23 May 2024 09:59:54 -0400 Subject: [PATCH] Removed Constraints on Yield Generators within helpers.py --- src/nestpy/helpers.py | 9 --------- 1 file changed, 9 deletions(-) diff --git a/src/nestpy/helpers.py b/src/nestpy/helpers.py index 7aab4f9..5d5e382 100644 --- a/src/nestpy/helpers.py +++ b/src/nestpy/helpers.py @@ -10,8 +10,6 @@ Main ingredients for the above steps: 1. np.vectorize, dictionary with yields, field array and energy array. - - Note: some of the yields will crash the plots at too high of energies (way above physical meaning) - so np.nan is returned to keep things running. """ import numpy as np @@ -105,13 +103,6 @@ def GetYieldsVectorized(interaction, yield_type, nc=None, **kwargs): else: interaction_object = interaction - if 'energy' in kwargs.keys(): - if interaction_object == GetInteractionObject('nr') and kwargs['energy'] > 2e2: - return np.nan - if interaction_object == GetInteractionObject('gammaray') and kwargs['energy'] > 3e3: - return np.nan - if interaction_object == GetInteractionObject('beta') and kwargs['energy'] > 3e3: - return np.nan yield_object = nc.GetYields(interaction = interaction_object, **kwargs) # returns the yields for the type of yield we are considering return getattr(yield_object, yield_type)