PhotonYield giving fluctuating number for Kr83m

[yuema137]

This is probably related to #84, which (I suppose) was already fixed in v1.5.4. But I still got this fluctuating number problem in v1.5.4 and this causes a trouble when we want to do the field correction with respect to the nominal field.
The minimal example to get this issue is:

for i in range(10):
    print(nestpy.PhotonYield(
                interaction=nestpy.Kr83m,
                energy=41.5,
                drift_field=23))

The ElectronYield function for Kr83m works well and gives a certain number instead of a fluctuating one.
Besides, I'm wondering if it's possible to make deltaT_ns as an input for the PhotonYield and ElectronYield functions for Kr83m (maybe it already does but I didn't know)? Or even better, use deltaT_range_ns=(deltaT_min, deltaT_max) and give the averaged yield accounting for the delay time dependence. This would be useful when we want to compare measurements with different delay time range.

I'm not familiar with this package and am not sure what is causing this problem. Thanks in advance for helping with this and sorry if I understood anything wrong.

[tunnell]

Wouldn't you just simulate ER instead?