A solution for performing maximum likelihood estimation on models built from histogram templates.
Project description
Histimator
A solution for performing maximum likelihood estimation on models built from histogram templates.
Free software: GNU General Public License v3
Documentation: https://histimator.readthedocs.io.
Features
TODO
Usage
the histimator core directory has a file called Models containing the core HistiModel class.
the model is initialised:
from histimator.models import HistiModel model = HistiModel("model name")
Each channel is defined as:
from histimodel.Channel import HistiChannel SR = HistiChannel("SignalRegion")
data can be added to the channels as:
SR.SetData([list of data points])
any number of samples are defined as:
from histimator.models import HistiSample sig = HistiSample("Signal") bkg = HistiSample("Background")
each of which needs a histogram:
sig.SetHisto(numpy.histogram) bkg.SetHisto(numpy.histogram)
currently the only parameters available are an overal normalisation on these templates. this is given with a name an initial value (default 1) and a range (default [0.1,10]). Currently no implementation is actually in place to tell Minuit about this range…:
sig.AddNorm("some_norm",1,0,3)
Finally, the samples must be added to the channel and this added to the model.:
SR.AddSample(sig) SR.AddSample(bkg) model.AddChannel(SR)
This model can now be evaluated using probfit Binned Likelihood function:
from iminuit import Minuit from probfit import BinnedLH blh = BinnedLH(model.pdf, data, bins=10, bound=bound, extended=True) m = Minuit(blh, some_norm=0.5, error_some_norm=1.5) m.migrad()
this has various built in plotting functionality.
Credits
This package was created with Cookiecutter and the audreyr/cookiecutter-pypackage project template.
History
0.1.0 (2018-02-16)
First release on PyPI.
Project details
Release history Release notifications | RSS feed
Download files
Download the file for your platform. If you're not sure which to choose, learn more about installing packages.
Source Distribution
Built Distribution
Hashes for Histimator-0.2.1-py2.py3-none-any.whl
Algorithm | Hash digest | |
---|---|---|
SHA256 | 6c941c5a17c15dffc6047e4519758c64c547cbeaf95ee697b8759c01a1d81010 |
|
MD5 | 1225e6ed222ede857f411cac71a67178 |
|
BLAKE2b-256 | 206d96ded101c7aa2a6567f052999043d98fb4511b3e80ed3c505e231046a9b0 |