ostaphep / ostap Goto Github PK

User friendly and more intuitive interface to ROOT and PyROOT

License: BSD 3-Clause "New" or "Revised" License

CMake 0.08% Python 13.42% Shell 0.01% C++ 15.04% C 71.46% Dockerfile 0.01%

python cern scientific-computing physics statistics root-cern pyroot interactive

ostap's Introduction

Ostap Project

Nowadays ROOT and PyROOT are de-facto standard tools for performing physics analysis. The Ostap project is a community-driven initiative aiming to provide more user friendly and more intuitive interface to ROOT and PyROOT and extending the existing functionality.

Project started in 2009 from the private collections of python functions used in Kali - framework for calibration of LHCb electromagnetic calorimeter. A lot of functionality is picked from Bender project - python based physics analysis environment used in LHCb experiemnt. Till Autumn 2016 the project was a part of LHCb software suit and with great success has been used for preparation of approximately 60 physics papers. A standalone, LHCb independent version, has appeared at start of 2017

Key features include:

Very easy manipulations with ROOT and RooFit objects: histograms, trees, datasets, etc
Very easy interface to RooFit machinery
Extended set of models.PDFs for RooFit
Powerful, pickle-based persistency for object
Interactive ostap analysis environment

Dependencies

mandatory: ROOT version <6.22 or >=6.24, RooFit
highly desirable: numpy
- mandatory for Fast Fourier Transform, used in histogram/function parameterization;
- optional for some other issues, in particular for the prime number treatment;
optional: scipy
- numerical integration (quadratures, cubatures), root finding, minimization;
- ostap offers home-made replacements, but the native methods from scipy are more efficient;
optional: pathos, dill, multiprocess and ppt
- needed for parallel processing;
- ostap offers a multiprocessing-based replacement with reduced functionality;
optional: terminaltables
- nice format of tables - in particular for nice printout for TTree, TChain, RooDataSet, ... ;
- ostap offers a home-made replacement with a bit reduced functionalty.
desirable: [more_itertoos]
- chunked stuff is used in parallelization for spliting jobs into chunks
- ostap offers a home-made replacement
optional (only for python3, 3.6<=version ) : berkeleydb
- python interface to Berkeley DB (libdb needs to be installed!)
optional (only for python3, 3.3<=version<3.10 ) : bsddb3
- python interface to Berkeley DB (libdb needs to be installed!)

Setup

There are several possibilities to start working with Ostap, you can build Ostap on the Linux or lxplus/7, run the docker container and now the Ostap is available on conda and SWAN.

The possible setup options are described here

ostap's People

Contributors

Stargazers

Watchers

Forkers

gitter-badger cpawley pro100tema bopopescu dmitryyugolubkov tatianaovsiannikova chrisburr mohamedelashri trellixvulnteam simonwaldherr

ostap's Issues

Bad interferecne for makePlots and calls for ROC-curves for tmva/chopping

Bad interference for makePlots and calls for ROC-curves for tmva/chopping
certains standard TMVA macros produced segfaults for various ROOT versions

cmp_fit does not from Feb 21

From Feb 21 cmp_fit does not for ROOT 6.27/1

prepare for ZENODO

prepare .zenodo.json file
prepare CITATION.CFF

Confusing output of ROOT.TTree.statVar method

Output prints " ... mean = XXX +- YYY ...", which is a bit confusing, as YYY is an individual uncertainty of single entry, but not an uncertainty of the mean, which is also random value. Better switch to ""...mean=XXX rms=YYY...".

In principle, estimates of uncertainties for mean and rms random valuers (under the Gaussian hypothesis ) are also worth to be printed.

Let me alsothsnk authors for the incredible tool!!

`PhaseSpace3` and `PhaseSpace3s` do not agree

Ostap::Math::PhaseSpace3 and Ostap::Math::PhaeSpace3s do not agree

import ROOT, math
from   ostap.core.pyrouts  import Ostap
from   ostap.logger.logger import getLogger

logger = getLogger ("ChicPi")

PS3 = Ostap.Math.PhaseSpace3  ## <-- HERE 
PS3 = Ostap.Math.PhaseSpace3s ## <-- HERE  
GeV = 1.0
MeV = 0.001 * GeV


m_chic1 = 3510.67 * MeV
m_K     = 493.677 * MeV
m_pi    = 139.570 * MeV
m_p     = 938.272 * MeV
m_Lb    = 5619.60 * MeV


ps_pi = PS3 ( m_chic1 , m_pi , m_p )
ps_K  = PS3 ( m_chic1 , m_K  , m_p )

logger.info ( "PS(pi)                       : %s "  % ps_pi ( m_Lb ) )
logger.info ( "PS(K)                        : %s "  % ps_K  ( m_Lb ) )
r = ps_pi ( m_Lb ) / ps_K ( m_Lb )

theta_C = 13.02 * math.pi / 180.0
tc2     = math.tan (  theta_C ) **2
logger.info ( "PS(pi)/PS(K)                 : %s "  % r    )
logger.info ( "tan^2 theta_C                : %s "  % tc2  )
logger.info ( "PS(pi)/PS(K) * tan^2 theta_C : %s "  % ( r  * tc2 ) )

restore state for `PDF.wilks*`

Restore the initial state for PDF.wilks* and similar methods.

Problem during ostap installation

I have follow the direction to install ostap on my Ubuntu machine. However, when I execute

cmake .. -DCMAKE_INSTALL_PREFIX=..

I have the following error. I am wondering if you could tell me how to fix it. Thank you.

........
........
-- Configuring done (3.8s)
CMake Error at source/CMakeROOT_6_23.cmake:225 (add_dependencies):
The dependency target "ROOT::MathMore" of target "ostapDict" does not
exist.
Call Stack (most recent call first):
source/CMakeLists.txt:14 (include)

CMake Error at source/CMakeROOT_6_23.cmake:194 (target_link_libraries):
Target "ostap" links to:

ROOT::MathMore

but the target was not found. Possible reasons include:

* There is a typo in the target name.
* A find_package call is missing for an IMPORTED target.
* An ALIAS target is missing.

Call Stack (most recent call first):
source/CMakeLists.txt:14 (include)

CMake Error at source/CMakeROOT_6_23.cmake:226 (target_link_libraries):
Target "ostapDict" links to:

ROOT::MathMore

but the target was not found. Possible reasons include:

* There is a typo in the target name.
* A find_package call is missing for an IMPORTED target.
* An ALIAS target is missing.

Call Stack (most recent call first):
source/CMakeLists.txt:14 (include)

-- Generating done (0.0s)
CMake Generate step failed. Build files cannot be regenerated correctly.

thisostap for zsh/ksh/csh/tcsh/fish

Current thisostap.sh works only for bash

versions for zsh/ksh and csh/tcsh are needed
version for fish coulb be nice

Error with ROOT.TMinuit

Seems TMinuit works not correctly as it fails if one tries to follow examples. (Ostap: 1.4.5.0, Python 3.8.0 | packaged by conda-forge | (default, Nov 22 2019, 19:11:38) )

In [1]: m = ROOT.TMinuit(9)
In [2]: help(m)
In [3]: m.addPar ( 'ququ' , 10 , 0.1 )

TypeError Traceback (most recent call last)
/home/adzyuba/miniconda3/envs/ostap-env/bin/ostap in
----> 1 m.addPar ( 'ququ' , 10 , 0.1 )

/home/adzyuba/miniconda3/envs/ostap-env/lib/python3.8/site-packages/ostap/fitting/minuit.py in mn_add_par(self, name, start, step, low, high)
225 ipar = len ( self )
226 ierr = ROOT.Long ( 0 )
--> 227 self.mnparm ( ipar , name , start , step , low , high , ierr )
228 #
229 return ierr

TypeError: void TMinuit::mnparm(int k, TString cnamj, double uk, double wk, double a, double b, int& ierflg) =>
could not convert argument 7 (use ctypes.c_int for pass-by-ref of ints)