Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

GetQuanta at zero drift field nuclear recoil #97

Open
Alex-X-Li opened this issue Nov 19, 2022 · 0 comments
Open

GetQuanta at zero drift field nuclear recoil #97

Alex-X-Li opened this issue Nov 19, 2022 · 0 comments

Comments

@Alex-X-Li
Copy link

Alex-X-Li commented Nov 19, 2022

I am just getting start to use NEST, so this might be a dumb question to ask.

I was trying to calculate the quanta (photon and electrons) produce in nuclear recoil with zero drift electric field in your nestpy_tutorial.ipynb. With the code
`

N = 50000
E_mono = 2 # keVnr
E_nr = np.ones(N)*E_mono
field = 0.

photons, electrons = GetQuanta(nestpy.nr, E_nr, field)
print(photons, electrons)
kwargs = {'energy': E_mono, 'drift_field': field}
nest_ly = nestpy.PhotonYield(interaction=nestpy.nr, **kwargs)
nest_qy = nestpy.ElectronYield(interaction=nestpy.nr, **kwargs)

`
And the photon and electrons produce is always 0.

However, for the electron recoil in zero field given by your given example
`

N = 50000
E_mono = 2 # keVee
E_er = np.ones(N)*E_mono
field = 0.

photons, electrons = GetQuanta(nestpy.beta, E_er, field)
kwargs = {'energy': E_mono, 'drift_field': field}
nest_ly = nestpy.PhotonYield(interaction=nestpy.beta, **kwargs)

nest_qy = nestpy.ElectronYield(interaction=nestpy.beta, **kwargs)

`
Shows non-zero photons and electrons production.
4f984d5f-3700-46c7-adc3-e53ae2c429b2

Physically speaking the nuclear recoil should produce photons and electrons as well. What's the reason for getting this result?

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

1 participant