Package: boost-histogram
boost-histogram for PythonPython bindings forBoost::Histogram(source), a C++14 library. This is one of thefastest librariesfor
histogramming, while still providing the power of a full histogram object. Seewhat's new.Other members of the boost-histogram family include:Hist: The first-party analyst-friendly histogram library that extends
boost-histogram with named axes, many new shortcuts including UHI+, plotting
shortcuts, and more.UHI: Specification for Histogram library interop, especially for plotting.mplhep: Plotting extension for matplotlib with support for UHI histograms.histoprint: Histogram display library for the command line with support for UHI.dask-histogram: Dask support for boost-histogram.UsageText intro (click to expand)importboost_histogramasbh# Compose axis however you like; this is a 2D histogramhist=bh.Histogram(bh.axis.Regular(2,0,1),bh.axis.Regular(4,0.0,1.0),)# Filling can be done with arrays, one per dimensionhist.fill([0.3,0.5,0.2],[0.1,0.4,0.9],)# NumPy array view into histogram counts, no overflow binsvalues=hist.values()# Make a new histogram with just the second axis, summing over the first, and# rebinning the second into larger bins:h2=hist[::sum,::bh.rebin(2)]We support theuhiPlottableHistogramprotocol, so boost-histogram/Histhistograms can be plotted via any compatible library, such asmplhep.CheatsheetSimplified list of features (click to expand)Many axis types (all supportmetadata=...)bh.axis.Regular(n, start, stop, ...): Make a regular axis. Options listed below.overflow=False: Turn off overflow binunderflow=False: Turn off underflow bingrowth=True: Turn on growing axis, bins added when out-of-range items addedcircular=True: Turn on wrapping, so that out-of-range values wrap around into the axistransform=bh.axis.transform.Log: Log spacingtransform=bh.axis.transform.Sqrt: Square root spacingtransform=bh.axis.transform.Pow(v): Power spacingSee also the flexibleFunction transformbh.axis.Integer(start, stop, *, underflow=True, overflow=True, growth=False, circular=False): Special high-speed version ofregularfor evenly spaced bins of width 1bh.axis.Variable([start, edge1, edge2, ..., stop], *, underflow=True, overflow=True, circular=False): Uneven bin spacingbh.axis.IntCategory([...], *, growth=False): Integer categoriesbh.axis.StrCategory([...], *, growth=False): String categoriesbh.axis.Boolean(): A True/False axisAxis features:.index(value): The index at a point (or points) on the axis.value(index): The value for a fractional bin (or bins) in the axis.bin(i): The bin edges (continuous axis) or a bin value (discrete axis).centers: The N bin centers (if continuous).edges: The N+1 bin edges (if continuous).extent: The number of bins (including under/overflow).metadata: Anything a user wants to store.traits: The options set on the axis.size: The number of bins (not including under/overflow).widths: The N bin widthsMany storage typesbh.storage.Double(): Doubles for weighted values (default)bh.storage.Int64(): 64-bit unsigned integersbh.storage.Unlimited(): Starts small, but can go up to unlimited precision ints or doubles.bh.storage.AtomicInt64(): Threadsafe filling, experimental. Does not support growing axis in threads.bh.storage.Weight(): Stores a weight and sum of weights squared.bh.storage.Mean(): Accepts a sample and computes the mean of the samples (profile).bh.storage.WeightedMean(): Accepts a sample and a weight. It computes the weighted mean of the samples.Accumulatorsbh.accumulator.Sum: High accuracy sum (Neumaier) - used by the sum method when summing a numerical histogrambh.accumulator.WeightedSum: Tracks a weighted sum and variancebh.accumulator.Mean: Running count, mean, and variance (Welfords's incremental algorithm)bh.accumulator.WeightedMean: Tracks a weighted sum, mean, and variance (West's incremental algorithm)Histogram operationsh.ndim: The number of dimensionsh.size or len(h): The number of bins+: Add two histograms (storages must match types currently)*=: Multiply by a scaler (not all storages) (hist * scalarandscalar * histsupported too)/=: Divide by a scaler (not all storages) (hist / scalarsupported too).kind: Eitherbh.Kind.COUNTorbh.Kind.MEAN, depending on storage.storage_type: Fetch the histogram storage type.sum(flow=False): The total count of all bins.project(ax1, ax2, ...): Project down to listed axis (numbers). Can also reorder axes..to_numpy(flow=False, view=False): Convert to a NumPy style tuple (with or without under/overflow bins).view(flow=False): Get a view on the bin contents (with or without under/overflow bins).values(flow=False): Get a view on the values (counts or means, depending on storage).variances(flow=False): Get the variances if available.counts(flow=False): Get the effective counts for all storage types.reset(): Set counters to 0 (growing axis remain the same size).empty(flow=False): Check to see if the histogram is empty (can check flow bins too if asked).copy(deep=False): Make a copy of a histogram.axes: Get the axes as a tuple-like (all properties of axes are available too).axes[0]: Get the 0th axis.axes.edges: The lower values as a broadcasting-ready array.axes.centers: The centers of the bins broadcasting-ready array.axes.widths: The bin widths as a broadcasting-ready array.axes.metadata: A tuple of the axes metadata.axes.size: A tuple of the axes sizes (size without flow).axes.extent: A tuple of the axes extents (size with flow).axes.bin(*args): Returns the bin edges as a tuple of pairs (continuous axis) or values (describe).axes.index(*args): Returns the bin index at a value for each axis.axes.value(*args): Returns the bin value at an index for each axisIndexing - SupportsUHI IndexingBin content access / settingv = h[b]: Access bin content by index numberv = h[{0:b}]: All actions can be represented byaxis:itemdictionary instead of by position (mostly useful for slicing)Slicing to get histogram or set array of valuesh2 = h[a:b]: Access a slice of a histogram, cut portions go to flow bins if presenth2 = h[:, ...]: Using:and...supported just like NumPyh2 = h[::sum]: Third item in slice is the "action"h[...] = array: Set the bin contents, either include or omit flow binsSpecial accessorsbh.loc(v): Supply value in axis coordinates instead of bin numberbh.underflow: The underflow bin (use empty beginning on slice for slicing instead)bh.overflow: The overflow bin (use empty end on slice for slicing instead)Special actions (third item in slice)sum: Remove axes via projection; if limits are given, use thosebh.rebin(n): Rebin an axisNumPy compatibilitybh.numpyprovides fasterdrop in replacementsfor NumPy histogram functionsHistograms follow the buffer interface, and provide.view()Histograms can be converted to NumPy style output tuple with.to_numpy()DetailsAll objects support copy/deepcopy/pickleFully statically typed, tested with MyPy.InstallationYou can install this library fromPyPIwith pip:python3-mpipinstallboost-histogramAll the normal best-practices for Python apply; Pip should not be very old (Pip
9 is very old), you should be in a virtual environment, etc. Python 3.7+ is
required; for older versions of Python (3.5 and 2.7),0.13will be installed
instead, which is API equivalent to 1.0, but will not be gaining new features.
1.3.x was the last series to support Python 3.6.Binaries available:The easiest way to get boost-histogram is to use a binary wheel, which happens
when you run the above command on a supported platform. Wheels are produced usingcibuildwheel; all common
platforms have wheels provided in boost-histogram:SystemArchPython versionsPyPy versionsManyLinux201464-bit3.7, 3.8, 3.9, 3.10, 3.11, 3.123.7, 3.8, 3.9, 3.10ManyLinux2014ARM643.7, 3.8, 3.9, 3.10, 3.11, 3.123.7, 3.8, 3.9, 3.10MuslLinux_1_164-bit3.7, 3.8, 3.9, 3.10, 3.11, 3.12macOS 10.9+64-bit3.7, 3.8, 3.9, 3.10, 3.11, 3.123.7, 3.8, 3.9, 3.10macOS Universal2Arm643.8, 3.9, 3.10, 3.11, 3.12Windows32 & 64-bit3.7, 3.8, 3.9, 3.10, 3.11, 3.12Windows64-bit3.7, 3.8, 3.9, 3.10manylinux2014: Requires pip 19.3.ARM on Linux is supported. PowerPC or IBM-Z available on request.macOS Universal2 wheels for Apple Silicon and Intel provided for Python 3.8+ (requires Pip 21.0.1 or newer).If you are on a Linux system that is not part of the "many" in manylinux or musl in musllinux, such as ClearLinux, building from source is usually fine, since the compilers on those systems are often quite new. It will just take longer to install when it is using the sdist instead of a wheel. All dependencies are header-only and included.Conda-ForgeThe boost-histogram package is available onconda-forge, as well. All supported variants are available.condainstall-cconda-forgeboost-histogramSource buildsFor a source build, for example from an "SDist" package, the only requirements are a C++14 compatible compiler. The compiler requirements are dictated by Boost.Histogram's C++ requirements: gcc >= 5.5, clang >= 3.8, or msvc >= 14.1. You should have a version of pip less than 2-3 years old (10+).Boost is not required or needed (this only depends on included header-only dependencies). You can install directly from GitHub if you would like.python-mpipinstallgit+https://github.com/scikit-hep/boost-histogram.git@developDevelopingSeeCONTRIBUTING.mdfor details on how to set up a development environment.ContributorsWe would like to acknowledge the contributors that made this project possible (emoji key):Henry Schreiner🚧💻📖Hans Dembinski🚧💻N!no⚠️📖Jim Pivarski🤔Nicholas Smith🐛physicscitizen🐛Chanchal Kumar Maji📖Doug Davis🐛Pierre Grimaud📖Beojan Stanislaus🐛Popinaodude🐛Congqiao Li🐛alexander-held🐛Chris Burr📖Konstantin Gizdov📦🐛Kyle Cranmer📖Aman Goel📖💻Jay Gohil📖This project follows theall-contributorsspecification.Talks and other documentation/tutorial sourcesTheofficial documentation is here, and includes aquickstart.2019-4-15 IRIS-HEP Topical meeting2019-10-17 PyHEP Histogram session-repo with talks and workbook2019-11-7 CHEP2020-07-07 SciPy2020-07-17 PyHEPAcknowledgementsThis library was primarily developed by Henry Schreiner and Hans Dembinski.Support for this work was provided by the National Science Foundation cooperative agreement OAC-1836650 (IRIS-HEP) and OAC-1450377 (DIANA/HEP). Any opinions, findings, conclusions or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of the National Science Foundation.