search

NESTCollaboration / nestpy

Raw Python bindings to the NEST library using pybind11.
GNU General Public License v3.0
13 stars 19 forks source link

PhotonYield giving fluctuating number for Kr83m #91

Open yuema137 opened 2 years ago

yuema137 commented 2 years ago

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 commented 2 years ago

Wouldn't you just simulate ER instead?