nestpy 2.1.1

Python bindings for the NEST noble element simulation

nestpy

Visit the tutorials directory for tutorials on the nestpy calls, maintenance, and benchmark plots.

These are the Python bindings for the NEST library, which provides a direct wrapping of functionality. The library is now pythonic, so be weary of the separate naming conventions for functions/variables from the C++ library.

You do not have to have NEST already installed to use this package.

Installing from PyPI

For 64-bit Linux or Mac systems, instally 'nestpy' should just require running:

pip install nestpy

You can then test that it works by running the example above.

Installing from source

Requirements: You must have CMake>=3.6 and a C++17 compatible compiler (GCC>=4.8) to build.

First, you must check out this repository then simply run the installer:

git clone https://github.com/NESTCollaboration/nestpy
cd nestpy
git submodule update --init --recursive
pip install .

Installing with custom NEST code

Almost all NEST users will want to incorporate some custom code into their workflow, such as custom Detector files or TestSpectra. In order to incorporate that custom code into the nestpy installation, you'll have to copy the files you've edited into the lib/nest/ directory and rerun:

pip install .

In order to create a more efficient workflow, we suggest the user takes the following steps:

1. Fork the official NEST repository into your own public/private one.
2. Make whatever changes to your fork and maintain them with commits.
3. Download or fork nestpy and change the NEST entry in the .gitmodules files to point to your custom fork of NEST.

[submodule "lib/nest"]
	path = lib/nest
	url = https://github.com/NESTCollaboration/nest # change this to any custom fork of NEST if you want to use your own code!
[submodule "lib/gcem"]
	path = lib/gcem
	url = https://github.com/kthohr/gcem.git

Usage

Python bindings to the NEST library:

Simple example covering the basic logic behind NEST

import nestpy

# This is same as C++ NEST with naming
nc = nestpy.NESTcalc(nestpy.detectors.VDetector())

interaction = nestpy.interactions.NR  # NR

E = 10  # keV
print('For an %s keV %s' % (E, interaction))

# Get particle yields
y = nc.GetYields(interaction,
		 E)

print('The photon yield is:', y.PhotonYield)
print('With statistical fluctuations',
      nc.GetQuanta(y).photons)

Examples using a vectorised approach

Get the yields as a function of energy

# Want to use custom yields parameters (e.g. from LZ WS2024)
lz_2024 = nestpy.detectors.LZ_WS2024()
# Get yields as a dataframe 
energy = np.linspace(1, 10, 10000)
nr_yeild = nestpy.helpers.get_yields_df(ne.interactions.NR, energy)
er_devel_yield = ne.helpers.get_yields_df(ne.interactions.beta, energy, nr_parameters = lz_2024.nr_parameters, er_parameters = lz_2024.er_parameters)

Run a vectorised simulation of S1 and S2 parameters

n = int(1e5)

t = nestpy.helpers.run_nest_df(
    # "beta",
    nestpy.interactions.beta,
    lz_2024,
    nestpy.spectra.CH3T(number=n), # Vectorised sampler from energy spectra
    # By default uses random positions in the detector
)

d = nestpy.helpers.run_nest_df(
    nestpy.interactions.NR,
    lz_2024,
    ne.spectra.DD(number=n),
    # nr_yield_params = det.nr_yield_params <- can set custom yield parameters
    # er_yield_params = det.er_yield_params <- can set custom yield parameters
    # width_params = det.width_yield_params <- can set this too
    # s1_mode" = NEST::S1CalculationMode::Hybrid <- Can simulate more detailed calculation or waveforms
    # s2_mode" = NEST::S2CalculationMode::Full <- Can simulate Waveforms with E-trains
)

For more examples on possible calls, please see the tests and tutorials folders.

Support

• Bugs: Please report bugs to the issue tracker on Github such that we can keep track of them and eventually fix them. Please explain how to reproduce the issue (including code) and which system you are running on.
• Help: Help can be provided also via the issue tracker by tagging your issue with 'question'
• Contributing: Please fork this repository then make a pull request. In this pull request, explain the details of your change and include tests.

Technical implementation

This package is a pybind11 wrapper of NEST that uses GitHub Workflows to build binaries using the manylinux Docker image from this page.

• Help from Henry Schreiner, which included a great binding tutorial
• Implementation also based on this
• Implementation of GitHub test and build actions was made possible by Joran Angevaare.

See AUTHORS.md for information on the developers.

Citation

When you use nestpy, please say so in your slides or publications (for publications, see Zenodo link above). You can mention this in addition to how you cite NEST. This is important for us being able to get funding to support this project.