pyjet

The interface between FastJet and NumPy

These details have not been verified by PyPI

Project links

Homepage

Project description

pyjet allows you to perform jet clustering with FastJet on NumPy arrays. By default pyjet only depends on NumPy and internally uses FastJet’s standalone fjcore release. The interface code is written in Cython that then becomes compiled C++, so it’s fast. Remember that if you use pyjet then you are using FastJet and should cite the papers listed here.

Strict dependencies

Python (2.7+, 3.5+)
Numpy

Getting started

pyjet provides the cluster() function that takes a NumPy array as input and returns a ClusterSequence from which you can access the jets:

from pyjet import cluster
from pyjet.testdata import get_event

vectors = get_event()
sequence = cluster(vectors, R=1.0, p=-1)
jets = sequence.inclusive_jets()  # list of PseudoJets
exclusivejets = sequence.exclusive_jets(3)  # Find the cluster history when there are 3 jets

The input is given in the form of a structured array in numpy. The first four fields of the input array vectors must be either:

np.dtype([('pT', 'f8'), ('eta', 'f8'), ('phi', 'f8'), ('mass', 'f8')])

or if cluster(..., ep=True):

np.dtype([('E', 'f8'), ('px', 'f8'), ('py', 'f8'), ('pz', 'f8')])

Note that the field names of the input array need not match ‘pT’, ‘eta’, ‘phi’, ‘mass’ etc. pyjet only assumes that the first four fields are those quantities. This array may also have additional fields of any type. Additional fields will then become attributes of the PseudoJet objects.

See the examples to get started:

https://github.com/scikit-hep/pyjet/raw/master/examples/jet_areas.png

Standalone Installation

To simply use the built-in FastJet source, from your virtual environment, run:

python -m pip install pyjet

And you’re good to go! If you have a old version of pip (<10), you will need to have Cython and Numpy already installed to build from source - however on most systems, you should get a binary wheel.

Get example.py and run it:

    curl -O https://raw.githubusercontent.com/scikit-hep/pyjet/master/examples/example.py
    python example.py
    jet#          pT        eta        phi       mass  #constit.
    1        983.280     -0.868      2.905     36.457         34
    2        901.745      0.221     -0.252     51.850         34
    3         67.994     -1.194     -0.200     11.984         32
    4         12.465      0.433      0.673      5.461         13
    5          6.568     -2.629      1.133      2.099          9
    6          6.498     -1.828     -2.248      3.309          6

    The 6th jet has the following constituents:
    PseudoJet(pt=0.096, eta=-2.166, phi=-2.271, mass=0.000)
    PseudoJet(pt=2.200, eta=-1.747, phi=-1.972, mass=0.140)
    PseudoJet(pt=1.713, eta=-2.037, phi=-2.469, mass=0.940)
    PseudoJet(pt=0.263, eta=-1.682, phi=-2.564, mass=0.140)
    PseudoJet(pt=1.478, eta=-1.738, phi=-2.343, mass=0.940)
    PseudoJet(pt=0.894, eta=-1.527, phi=-2.250, mass=0.140)

    Get the constituents as an array (pT, eta, phi, mass):
    [( 0.09551261, -2.16560157, -2.27109083,   4.89091390e-06)
     ( 2.19975694, -1.74672746, -1.97178728,   1.39570000e-01)
     ( 1.71301882, -2.03656511, -2.46861524,   9.39570000e-01)
     ( 0.26339374, -1.68243005, -2.56397904,   1.39570000e-01)
     ( 1.47781519, -1.7378898 , -2.34304346,   9.39570000e-01)
     ( 0.89353864, -1.52729244, -2.24973202,   1.39570000e-01)]

    or (E, px, py, pz):
    [( 0.42190436, -0.06155242, -0.07303395, -0.41095089)
     ( 6.50193926, -0.85863306, -2.02526044, -6.11692764)
     ( 6.74203628, -1.33952806, -1.06775374, -6.45273802)
     ( 0.74600384, -0.22066287, -0.1438199 , -0.68386087)
     ( 4.43164941, -1.0311407 , -1.05862485, -4.07096881)
     ( 2.15920027, -0.56111108, -0.69538886, -1.96067711)]

Reclustering the constituents of the hardest jet with the kt algorithm
[PseudoJet(pt=983.280, eta=-0.868, phi=2.905, mass=36.457)]

Go back in the clustering sequence to when there were two jets
PseudoJet(pt=946.493, eta=-0.870, phi=2.908, mass=20.117)
PseudoJet(pt=36.921, eta=-0.800, phi=2.821, mass=4.119)

Ask how many jets there are with a given dcut
There are 9 jets with a dcut of 0.5

Get the jets with the given dcut
1 PseudoJet(pt=308.478, eta=-0.865, phi=2.908, mass=2.119)
2 PseudoJet(pt=256.731, eta=-0.868, phi=2.906, mass=0.140)
3 PseudoJet(pt=142.326, eta=-0.886, phi=2.912, mass=0.829)
4 PseudoJet(pt=135.971, eta=-0.870, phi=2.910, mass=0.140)
5 PseudoJet(pt=91.084, eta=-0.864, phi=2.899, mass=1.530)
6 PseudoJet(pt=30.970, eta=-0.831, phi=2.822, mass=2.124)
7 PseudoJet(pt=7.123, eta=-0.954, phi=2.939, mass=1.017)
8 PseudoJet(pt=5.951, eta=-0.626, phi=2.818, mass=0.748)
9 PseudoJet(pt=4.829, eta=-0.812, phi=3.037, mass=0.384)

Using an External FastJet Installation

To take advantage of the full FastJet library and optimized O(NlnN) kt and anti-kt algorithms you can first build and install FastJet and then install pyjet with the --external-fastjet flag. Before building FastJet you will need to install CGAL and GMP.

On a Debian-based system (Ubuntu):

sudo apt-get install libcgal-dev libcgal11v5 libgmp-dev libgmp10

On an RPM-based system (Fedora):

sudo dnf install gmp.x86_64 gmp-devel.x86_64 CGAL.x86_64 CGAL-devel.x86_64

On Mac OS:

brew install cgal gmp wget

Then run pyjet’s install-fastjet.sh script:

curl -O https://raw.githubusercontent.com/scikit-hep/pyjet/master/install-fastjet.sh
chmod +x install-fastjet.sh
sudo ./install-fastjet.sh

Now install pyjet like:

python -m pip install numpy Cython
python setup.py install --external-fastjet

pyjet will now use the external FastJet installation on your system.

Note on units

The package is indifferent to particular units, which are merely “propagated” through the code. We do recommend that the HEP units be used, as defined in the units module of the hepunits package.

It is worth noting that the azimuthal angle phi is expressed in radians and varies from pi to pi.

Developing

If you want to setup for development:

python3 -m venv .env
source .env/bin/activate
pip install -e .[dev]
pytest

Project details

These details have not been verified by PyPI

Project links

Homepage

Release history Release notifications | RSS feed

1.9.0

Jul 1, 2022

1.8.2

Jan 18, 2021

1.8.1

Jan 12, 2021

1.8.0

Nov 24, 2020

1.7.0

Nov 13, 2020

This version

1.6.0

Dec 28, 2019

1.5.0

May 7, 2019

1.4.0

Jan 10, 2019

1.3.0

Aug 11, 2018

1.2.0

Jul 18, 2018

1.1.1

Apr 4, 2018

1.0.1

Aug 10, 2017

1.0.0

Jun 7, 2017

0.5.1

Jun 6, 2017

0.5.0

Jun 6, 2017

0.4.0

Jun 2, 2017

0.3.0

Jun 1, 2017

0.2.0

May 29, 2017

0.1.0

May 27, 2017

0.0.4

May 24, 2017

0.0.3

May 23, 2017

0.0.2

May 23, 2017

0.0.1

May 23, 2017

Download files

Download the file for your platform. If you're not sure which to choose, learn more about installing packages.

Source Distribution

pyjet-1.6.0.tar.gz (99.9 kB view details)

Uploaded Dec 28, 2019 Source

Built Distributions

pyjet-1.6.0-cp38-cp38-manylinux2010_x86_64.whl (1.6 MB view details)

Uploaded Dec 28, 2019 CPython 3.8 manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp38-cp38-manylinux1_x86_64.whl (1.2 MB view details)

Uploaded Dec 28, 2019 CPython 3.8

pyjet-1.6.0-cp38-cp38-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 3.8

pyjet-1.6.0-cp38-cp38-macosx_10_9_x86_64.whl (317.7 kB view details)

Uploaded Dec 28, 2019 CPython 3.8 macOS 10.9+ x86-64

pyjet-1.6.0-cp37-cp37m-manylinux2010_x86_64.whl (1.5 MB view details)

Uploaded Dec 28, 2019 CPython 3.7m manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp37-cp37m-manylinux1_x86_64.whl (1.2 MB view details)

Uploaded Dec 28, 2019 CPython 3.7m

pyjet-1.6.0-cp37-cp37m-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 3.7m

pyjet-1.6.0-cp37-cp37m-macosx_10_9_x86_64.whl (317.6 kB view details)

Uploaded Dec 28, 2019 CPython 3.7m macOS 10.9+ x86-64

pyjet-1.6.0-cp36-cp36m-manylinux2010_x86_64.whl (1.5 MB view details)

Uploaded Dec 28, 2019 CPython 3.6m manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp36-cp36m-manylinux1_x86_64.whl (1.2 MB view details)

Uploaded Dec 28, 2019 CPython 3.6m

pyjet-1.6.0-cp36-cp36m-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 3.6m

pyjet-1.6.0-cp36-cp36m-macosx_10_9_x86_64.whl (318.7 kB view details)

Uploaded Dec 28, 2019 CPython 3.6m macOS 10.9+ x86-64

pyjet-1.6.0-cp35-cp35m-manylinux2010_x86_64.whl (1.5 MB view details)

Uploaded Dec 28, 2019 CPython 3.5m manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp35-cp35m-manylinux1_x86_64.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 3.5m

pyjet-1.6.0-cp35-cp35m-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 3.5m

pyjet-1.6.0-cp27-cp27mu-manylinux2010_x86_64.whl (1.5 MB view details)

Uploaded Dec 28, 2019 CPython 2.7mu manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp27-cp27mu-manylinux1_x86_64.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 2.7mu

pyjet-1.6.0-cp27-cp27mu-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 2.7mu

pyjet-1.6.0-cp27-cp27m-manylinux2010_x86_64.whl (1.5 MB view details)

Uploaded Dec 28, 2019 CPython 2.7m manylinux: glibc 2.12+ x86-64

pyjet-1.6.0-cp27-cp27m-manylinux1_x86_64.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 2.7m

pyjet-1.6.0-cp27-cp27m-manylinux1_i686.whl (1.1 MB view details)

Uploaded Dec 28, 2019 CPython 2.7m

pyjet-1.6.0-cp27-cp27m-macosx_10_9_x86_64.whl (315.2 kB view details)

Uploaded Dec 28, 2019 CPython 2.7m macOS 10.9+ x86-64

File details

Details for the file pyjet-1.6.0.tar.gz.

File metadata

Download URL: pyjet-1.6.0.tar.gz
Upload date: Dec 28, 2019
Size: 99.9 kB
Tags: Source
Uploaded using Trusted Publishing? No
Uploaded via: twine/3.1.1 pkginfo/1.5.0.1 requests/2.22.0 setuptools/40.6.2 requests-toolbelt/0.9.1 tqdm/4.41.0 CPython/3.6.8

File hashes

Hashes for pyjet-1.6.0.tar.gz
Algorithm	Hash digest
SHA256	`ab6e63f8a8fd73bbd76ef2a384eea69bc1c201f2ce876faa4151ade6c0b20615`
MD5	`c31fefda590ee5885db97c8ef533bd20`
BLAKE2b-256	`18694c9f5e872b865f8a9713b71c95e19b6bbaadf9c725adc8bdf63c407c1cdb`

Algorithm	Hash digest
SHA256	`4c48617d1afd225798f7f92e6237a7cfc4905ce9b7bc148ec44d8d042b393325`
MD5	`56a6d74d8d3e406b9fdb4126dde7fc85`
BLAKE2b-256	`235b0bf380764599615a2f21cb0e57faec5f825ae49446fe001d740fdcb33c36`

Algorithm	Hash digest
SHA256	`486a24791c046e5f126cc35ba5fed732a85db3f97f90781e0aa20fca92bfd9f4`
MD5	`1ea6bd6084a852c89515d4a075ac57d1`
BLAKE2b-256	`6fab6f3e9b04544e834bda72e893d690b344f0a78f74172b24964a21f07ca9c5`

Algorithm	Hash digest
SHA256	`f5ad68445a5e49837d6ae2930f40cd157b451bee38751ba0ec1a639b0dc09f6e`
MD5	`d7bd58c8c8d64e7ddfa9965eb3b9aaca`
BLAKE2b-256	`1b94e1bffb8cb54cdbae3c4fb2d13a946b09cace9f5a4ac5042666fc13c92359`

Algorithm	Hash digest
SHA256	`e7c8a9899f97256440c4b6d9b8165dcbacd8b259d89421a5ee9934123e9f1d64`
MD5	`28f2d30be63e8bb3b1d7138ae362287c`
BLAKE2b-256	`6febb01b3a157249ba0f6543ea5b2dfc2a749baed676b16c28fa51a4d3136eb2`

Algorithm	Hash digest
SHA256	`078177c6f5dfc54277cea80f6c7026c71da7e8381952ddd2dc3b139ac46d4594`
MD5	`7a994790bccfd2ccafe2de67bbc2c209`
BLAKE2b-256	`72958a9a1bc6281dc294710c47c74ea32a021841f00126d17d5f195bd7d7f108`

Algorithm	Hash digest
SHA256	`8d1c5930e08968d0385e14ad257b994e8ccb06a58d79a8d9b71f7171854091ed`
MD5	`d3b8fb0226df51c6eb40712c13cb7b68`
BLAKE2b-256	`95c63cbd1f3cfef66bfdcf838a97c278913f655059e892756a0f5505ec58f2ed`

Algorithm	Hash digest
SHA256	`9ac9b83a5b8916a1db4d5f8148ebb3fc6aa935e0dc1cab82856bb68e0235373a`
MD5	`7ae743b9bbb31e8f88622884dde37c1e`
BLAKE2b-256	`ac09a9c3ff8ff0fdf1b28d81e7cc21c852860cbbdcd3c56bcf63c3ce24f683c8`

Algorithm	Hash digest
SHA256	`85b0d6d32ca309b11d85dd93b65e4c9c664189addd1ac15557a37753e20ee317`
MD5	`364bea3b614f27bfce5536e24de7c7db`
BLAKE2b-256	`6a2f5f0eb751ec76f30d569d0166d19ce200973661fb4b4144265258bc79ebf9`

Algorithm	Hash digest
SHA256	`b70244b0e67014f7fa8299c8d5b0de53deb5a2a14cdbcda94dab47680e413174`
MD5	`e5cff4e15126ac88971f352b46761df8`
BLAKE2b-256	`4f40f734ed845beb91b83115102f04cfda7c301e5e12294c90435e325c3f1ee4`

Algorithm	Hash digest
SHA256	`254446793b839510f8fd2c4e6a6a4b631c59498d44b650107fa00d6f4d434abb`
MD5	`416be0c924b6d331e4a6fee7cc56d9c8`
BLAKE2b-256	`0ccd172a2beae13d489cca042815c6eba4cf84f899f58714d947e11e5d882b4a`

Algorithm	Hash digest
SHA256	`988d6e20d794a11ce6418e06eb4c9a3c0dcb73a839fca88ff084990a4110f42f`
MD5	`48b29646d30d12abb50fa9c773ac81ee`
BLAKE2b-256	`eb3c0cf00a3e03f1c6165604cf3a062a5c2b5d38ec2dc012aa0381ecc4b9b607`

Algorithm	Hash digest
SHA256	`9c31ab673f6064c4652187d0bb13f30a47900d24aa10fba1629b98210348ffbc`
MD5	`4e7d1c7fcb473d123c7e956735e22eb3`
BLAKE2b-256	`a51dcf113c4964c16c056980f0b9acf1faac0fefda57984e5cdf53a1c331ed54`

Algorithm	Hash digest
SHA256	`c40ec44eae21c48137fd63467c7636cf1602683846869e2577b93d626ef9a84b`
MD5	`90cea8b7ea7fe1e37eaa261d67472b81`
BLAKE2b-256	`128d90887cbe59ca56219283e105e91429d827df41e2e99320a7a42df4e0fd83`

Algorithm	Hash digest
SHA256	`fbc4c455c3024b168c46187238f40f842299f57316fb7fc1a305caa9070a7dd9`
MD5	`6d378f83904415e4e560fa060a80818a`
BLAKE2b-256	`24dfcb3e53c175af91574979e21b99947fda6b8abafcdd93ea6ba08d30dd9bce`

Algorithm	Hash digest
SHA256	`e2be39adf39dd3727aab9683606c40ef218b7985c615eb1d47334b960715ff3d`
MD5	`16b379365ed1589feeff3f868a3aee06`
BLAKE2b-256	`6807c22ab62f78c53564d558fc308cf67e512889f5630658fec2683ccc2883b6`

Algorithm	Hash digest
SHA256	`a8e32570ad1af3884c37fc5b34c8583bfacab8e0e5362a0ebd8dc671269033ad`
MD5	`2bba23dc9146012f83c170b0ce111b6f`
BLAKE2b-256	`03e8e106bbc34874c83b0bf8e4534b5c3c6f824716672d41d0758c1c47d51afe`

Algorithm	Hash digest
SHA256	`9cd68f26618ee2f62b972910353e0068775c9fd8fbaf4d6555b0a88fa323da20`
MD5	`8e0c3b1717cc694b5f97508c96393cae`
BLAKE2b-256	`b4fd5c6cce55ecc08ee51115434a39d7cdc49d1985cedbc0e21e8e1dbb315fe4`

Algorithm	Hash digest
SHA256	`c311e34ca8603e358bbb1b9049f2308bfcb0515959fe24df8cd0b014002cfccd`
MD5	`2681cd0f75c7cdf5c76dacca9339755e`
BLAKE2b-256	`38d0f2cc75dce968fa20e804311d0c5bb0cc239ab01017c0946081c00dc53c37`

Algorithm	Hash digest
SHA256	`c4ebe13298b752892fbe0e708083b020b0aba755da1be6497d59acbb1bb99e9c`
MD5	`9d700e7c73ef75675c90e11c2e851700`
BLAKE2b-256	`08407f019068cc67a5b5423a6a3501c603c0aa17ab61fc79dc728411e94eae07`

Algorithm	Hash digest
SHA256	`7faf90d8125080854af45aadb23946111a48c06b293549d36901328c0d583688`
MD5	`5d9d2a8e7984f63c835727135d6b84f9`
BLAKE2b-256	`f156a396790dd40afb5d9728a8fe82f88f00ed71e387524d627d5da0b93d18b9`

pyjet 1.6.0

Navigation

Verified details

Maintainers

Unverified details

Project links

Meta

Classifiers

Project description

Strict dependencies

Getting started

Standalone Installation

Using an External FastJet Installation

Note on units

Developing

Project details

Verified details

Maintainers

Unverified details

Project links

Meta

Classifiers

Release history Release notifications | RSS feed

Download files

Source Distribution

Built Distributions

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

File details

File metadata

File hashes

Algorithm	Hash digest
SHA256	`b6f8dd907d75cd57332ffe8abc23061aa29070a32cca656bd3d470cc28f0d7aa`
MD5	`0c80c0ce1d424ae47a34d04634815164`
BLAKE2b-256	`8ba993bf24119b9dcc66c0dfe7c95e5ac62b13fa05dd9087a5cac9f52a6afe33`

Algorithm	Hash digest
SHA256	`3d4dfdf4dcf2a85b2e1a1794f173bdd2851d8d9ddfde4126a4e6060a8288f9d0`
MD5	`653af8024ec3116ddff983f2e070a226`
BLAKE2b-256	`99d481c3bab0f9ac393f33697ae566b49886668d69bacaa68a9ab3ca34996d59`