TDCRPy 1.11.2__py3-none-any.whl → 2.15.8__py3-none-any.whl

tdcrpy-2.15.8.dist-info/RECORD

@@ -0,0 +1,82 @@
+tdcrpy/TDCRPy.py,sha256=3puNWaySzWsCg08OYzI89rWRy5f2CV2aa5HuC0tVCFM,67965
+tdcrpy/TDCRPy1.py,sha256=QTBZh5B5JWnGB0BQfD-cFmwA9W080OD4sG-aj50-ejo,38106
+tdcrpy/TDCR_model_lib.py,sha256=aaQfweMgml6_yzvEAMmVC0zXcYd2WBxqY-os-nzmktA,138157
+tdcrpy/TDCRoptimize.py,sha256=c2XIGveeLdVYYek4Rg6dygMvVA2xIrIkMb3L-_jUucM,6496
+tdcrpy/__init__.py,sha256=9Djir8dPNchcJVQvhl-oRHEOsoDkiZlkOhWT-eHR7wQ,95
+tdcrpy/config.toml,sha256=UvIV6oUFjkk96c0Z053l14vekVc0eZ2-C0xy8MTs2zQ,1725
+tdcrpy/test2.py,sha256=poLLXJyIaCeqh1VSkwgbi-udvY7lQjxz_YStKjJXGhU,501
+tdcrpy/test_randomGen.py,sha256=CWN8jWJ3VAu6WhkRIR9BHv-e4V4hGx5jSfWBqGE0TT4,2195
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_C-14.txt,sha256=Eh3KaNbfYHakk_uStLu8K1aFciO6_i_rS2yKxGGppDE,8416
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Ca-45.txt,sha256=TymodcK4ttoO1duZuW3RGOwHFwPPzw2ESPc_H_QQN8k,8830
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Co-60.txt,sha256=kxD5E7tk_Gc1Ylg8qCG1r3oB21m7wUT4zBWsmbseiMU,40203
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_H-3.txt,sha256=QJNhdLbIPpm8njOc1YBR3m16x0cG-2hhF1bDVIlzsic,8320
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Ni-63.txt,sha256=xvc6RPEmYOrUGGHOgVHmSajPuRG49fmKxXU-tZjpdLo,8869
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Pm-147.txt,sha256=MvxB_Ax2_AM_PQOLdTdUUt8LNrmh7li9hJk_Ugu7Xe0,8651
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Pu-241.txt,sha256=ZeAn2opKFgvvYzFtan18zqDf5_TG5gNstU2_4fyITgU,9228
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_S-35.txt,sha256=HUdghpHIKukza2_uw62w2COv6i4-1gCA3gbv8qeZc2Y,9016
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Sr-89.txt,sha256=dLdl_yULusR-C3xHqNK7cAs3Z_DgEwL3Hf_a0pSCFFU,8356
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Sr-90.txt,sha256=BcdNCg8azMW-M25dXroAeD-DNgUCGmaL9ay5zdc_Cpo,15248
+tdcrpy/MCNP-MATRIX/Spectra_for_analytical_model/dep_spectrum_Tc-99.txt,sha256=_8KEDdembKT79zM1LsaTiFdPp6bHW-S9f30eXliaTCc,8014
+tdcrpy/MCNP-MATRIX/matrice/fichier/E_depose.txt,sha256=NkeaRgZAiz9fB_sgsf0zJRdDxOEsonjqMU3HSbCSQo8,64129
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-beta-_1_200k.txt,sha256=bAjifmJH_9nRc8MWlqyUfMrKNrHSTu85dZYd2k7Sz_o,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-beta-_2000_10000k.txt,sha256=_K3dDYkwzAYbqwqBPv_lMB94LGcaJPPoLaGhvNmxGxA,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-beta-_200_2000k.txt,sha256=cVGktWekTTSFkGmgG-uIVcbRNNQ5jQf43DVBqDj7Sw0,18978766
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-photon_1_200k.txt,sha256=JqlkXfAP0zx27asFwp9ofxpRUHAVkwfbCq9uJ3wqaEs,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-photon_2000_10000k.txt,sha256=lPNY-a8__0HAvvE5kaC5Mzv5I1yLlV-1YL4PSBwtzJQ,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_10ml-photon_200_2000k.txt,sha256=8bzS9V0apkiGYAVWAzinGs0subLZgvkWoHaYVivSJ1o,18978766
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-beta-_1_200k.txt,sha256=b8J-4hEyr54ImgC_C_M4jZWVOHoM0cSu13A17W7nfYg,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-beta-_2000_10000k.txt,sha256=Y4kf4LyxVGoKQgbF3pgkJUfV0GkB9bGeEmSNAedqXiw,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-beta-_200_2000k.txt,sha256=6lGJUeQkP8sRw2tCZDZlQ3nrtPju6LFMhUQFLi3mBKY,18978766
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-photon_1_200k.txt,sha256=AIl93JKnOv84mS7D1T1ryECV4X1pcNnFEiZvRCgQEZo,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-photon_2000_10000k.txt,sha256=gQSUAK5J9cwzeIRCada-pVHCX8XVK8DIUoXZpLDmoAs,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_13ml-photon_200_2000k.txt,sha256=3j2Y1lBVefHH2B0i-UZYoB07qk2GFhyst4a3ZhkWuKE,18978766
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-beta-_1_200k.txt,sha256=3JEuFXWctT4KmPTDoMmO94uouRS8D8sUA_Yndf1mesQ,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-beta-_2000_10000k.txt,sha256=ojcqnJ8IoSvlpdWzmucVtMHUXR857qK5oVudRGvQSS4,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-beta-_200_2000k.txt,sha256=mqo4-ZdlSn38lJ692YSLsxCYEnzNcS8_dhZpX7LfITI,18978766
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-photon_1_200k.txt,sha256=kxDaTMq-3H_telytHRV-bJXynqSyW9SmUMRVDW6y3yY,4213603
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-photon_2000_10000k.txt,sha256=oKl_dTNoQCPAliK1erzvtO1PHBCRItAZvY0wkBa_qN4,16872466
+tdcrpy/MCNP-MATRIX/matrice/fichier/matrice_16ml-photon_200_2000k.txt,sha256=B1f_TWs6wzreI4mbYYJDOkm2aDik49w4yQ96tR7Awe0,18978766
+tdcrpy/Micelle/2nmfaq01.csv,sha256=nLnOq8HyJGJSHeJ91-ZLRUZPgmaogiHcQk1SicHYaqE,1340
+tdcrpy/Micelle/faq01.csv,sha256=pAERSkEojSzlJW_Lufp_UmCLDP2Xj0d0gI6EgMveC_s,7573
+tdcrpy/Quenching/QuenchEnergyAlpha_0.015.txt,sha256=kaSid3pQuBSwIeJ7CTmx6NsFgr2_A6u5GOgMqJiVvQA,195818
+tdcrpy/Quenching/QuenchEnergyAlpha_1.1e-05.txt,sha256=lozxOyzSB4FAukpuTQUGbtQHqiWldqMXXMgstEC76sc,182522
+tdcrpy/Quenching/QuenchEnergyAlpha_1.2e-05.txt,sha256=pjH7BVj9ohsPjL6hD6nJvL76nIL6HhkcVsmRYskI7dU,182825
+tdcrpy/Quenching/QuenchEnergyAlpha_1.3e-05.txt,sha256=HxecTi-9A-V5niGVz_w73SKHvoQDUrimBun3XlRJtzU,182903
+tdcrpy/Quenching/QuenchEnergyAlpha_1.4e-05.txt,sha256=tP4_1vUzOTYTreTYWA1w7922k2CZtC8aTbEPynC9doM,183085
+tdcrpy/Quenching/QuenchEnergyAlpha_1.5e-05.txt,sha256=QzP2sSQe5P_OMbfgrPbIEHwOmHFzy3M7yS001zFeEwQ,183344
+tdcrpy/Quenching/QuenchEnergyAlpha_1e-05.txt,sha256=39uybKh8rSGlewQ34vo_BIIzT_0WeVfJm9MFmnRzMKI,182418
+tdcrpy/Quenching/QuenchEnergyAlpha_6e-06.txt,sha256=J_bLgnNolr4zce5xPRlCwE5ow8SdacvNtiy-wA7fgI8,183575
+tdcrpy/Quenching/QuenchEnergyAlpha_7e-06.txt,sha256=ixhrw1dF8Y22cGSGTE1RsEPrzUWSb8zURZH3kC4h81g,183332
+tdcrpy/Quenching/QuenchEnergyAlpha_8e-06.txt,sha256=xfu7weEsCAPnDxTjy8ifw7rZdOi5U54CY0_5s78bhW0,182986
+tdcrpy/Quenching/QuenchEnergyAlpha_9e-06.txt,sha256=KdfsVNvmfRy6PWQXOAQwQaZFzVc10BQ9x_hXNGSoL2c,182611
+tdcrpy/Quenching/QuenchEnergyElectron_0.006.txt,sha256=HeoyfYbKI_6YUxNDrynlj00af-_TU6Mcynkt7xThksQ,181564
+tdcrpy/Quenching/QuenchEnergyElectron_0.007.txt,sha256=1quYASu3nJ8Vzhi-58sw1AAwaDHomsXJKDiPDszOr20,181528
+tdcrpy/Quenching/QuenchEnergyElectron_0.008.txt,sha256=cdwGUrm9k3blIOFAyhSsQOPt3doqv8_356WZZt5VdGE,181612
+tdcrpy/Quenching/QuenchEnergyElectron_0.009.txt,sha256=PUHKJihb8ziCniKtVi2mdfC3w6xSZz2ekENx7tC4nEU,181560
+tdcrpy/Quenching/QuenchEnergyElectron_0.01.txt,sha256=BX23Y57xg19TBoFg7J1oU0NFD7uk69ETywyfcyM2j0g,181558
+tdcrpy/Quenching/QuenchEnergyElectron_0.011.txt,sha256=LFVN8d9lnRcHkKSLsSW9LKfVSC437PjBmpPW7W8qSlo,181474
+tdcrpy/Quenching/QuenchEnergyElectron_0.012.txt,sha256=ImvhAp_S5roe2kpQXJs00BN7_tS1DsjorskBLOXW-5I,181477
+tdcrpy/Quenching/QuenchEnergyElectron_0.013.txt,sha256=qkSGhRTPkRaj9tHaxcYMrN1HnJCfZtIP-SCgSiXvWLg,181545
+tdcrpy/Quenching/QuenchEnergyElectron_0.014.txt,sha256=ds0cxeDL7aiHjPRnSIqj5mfcaVcwAPLz2CFz4TXniLQ,181661
+tdcrpy/Quenching/QuenchEnergyElectron_0.015.txt,sha256=AMYTsBO9TVvC1v29k-1SsvY7_Xl1cw7SQyHlgtwxH3Y,181566
+tdcrpy/Quenching/TandataUG.txt,sha256=AkUGz2dAZ6g3bQIIacR5f92aa8gmiklgRbN4Dyz-12k,260012
+tdcrpy/Quenching/alpha_toulene.txt,sha256=rJB9xCLM3GLg91EKfhu8sGwro281ySq6R0G5mXMt3f8,1553
+tdcrpy/Quenching/inputVecteurAlpha.txt,sha256=ee4emplX51vhomC7cGNP_eE5eqS9By5PaSSaa4hh1eg,178652
+tdcrpy/Quenching/inputVecteurElectron.txt,sha256=5NBRBaiB7op1SJ6v952w9gFKqh0x83pn8zYcoUpy33Y,181498
+tdcrpy/decayData/All-nuclides_BetaShape.zip,sha256=k-UOym1MPPzGrOwal2mrlcXi3NSR1vAqzO15qTsW25c,8000329
+tdcrpy/decayData/All-nuclides_Ensdf.zip,sha256=J7m6Uq9UqEvBgazuu31NPDO-43xMh3Ktq6w6yMEWBF0,453398
+tdcrpy/decayData/All-nuclides_PenNuc.zip,sha256=o2DftXAJW6-zpBpi_ozHZiDIzl3XcEcAjFJZWAeHbI8,410519
+tdcrpy/decayData/atom-ENDF-VII0.zip,sha256=oP4Z7qWweWF4j12QAIfvLFLlqeY1RomXwXD2iFgQCqo,1624444
+tdcrpy/decayData/photo-ENDF.zip,sha256=eOiVUt6OKGl0a3pErfaibZzF_yoSzswdUCTbnGNlu7g,249674
+tdcrpy/docs/_build/html/genindex.html,sha256=FopjoEyhTbeh7qOem-A7er-6yCIfdhv1BIzs9dJocCc,2522
+tdcrpy/docs/_build/html/index.html,sha256=KbgBWH_0ixaswXk_lCM_pMk10v4rLs5qdPubj25SUMc,3443
+tdcrpy/docs/_build/html/search.html,sha256=IeTP6CsRdqKWOJttbtskg1oBgdN8norKlzai4Tmk8Kc,2799
+tdcrpy/docs/_build/html/source/modules.html,sha256=Jf-qxVBId0UgpwyvYuyjtMNG-ezPOK1oH_IxRXAWTRs,3601
+tdcrpy/docs/_build/html/source/tdcrpy.html,sha256=-38lHMNFB22p1tWJEeN3yDqfDiCYE304vxDamO1-iRc,3779
+tdcrpy/test/__init__.py,sha256=47DEQpj8HBSa-_TImW-5JCeuQeRkm5NMpJWZG3hSuFU,0
+tdcrpy/test/test_tdcrpy.py,sha256=JINqSEMFoNpptE4f3h6ZzTYW1rBx90KkaoQzltSg-No,4692
+tdcrpy-2.15.8.dist-info/licenses/LICENCE.md,sha256=ZTpWyGU3qv_iwEpgvCijoCuCYpOPpyzJCgOk46WpUKU,1066
+tdcrpy-2.15.8.dist-info/METADATA,sha256=G68HRdhNUqMMynusKChlR44rsKGNXoD7qFWHo-qisY8,45363
+tdcrpy-2.15.8.dist-info/WHEEL,sha256=_zCd3N1l69ArxyTb8rzEoP9TpbYXkqRFSNOD5OuxnTs,91
+tdcrpy-2.15.8.dist-info/top_level.txt,sha256=f4vzFFcKSEnonAACs0ZXuRczmroLLqtPTqXFymU_VU0,14
+tdcrpy-2.15.8.dist-info/RECORD,,

{TDCRPy-1.11.2.dist-info → tdcrpy-2.15.8.dist-info}/WHEEL

@@ -1,5 +1,5 @@
 Wheel-Version: 1.0
-Generator: bdist_wheel (0.43.0)
+Generator: setuptools (80.9.0)
 Root-Is-Purelib: true
 Tag: py3-none-any
 

TDCRPy-1.11.2.dist-info/METADATA

@@ -1,334 +0,0 @@
-Metadata-Version: 2.1
-Name: TDCRPy
-Version: 1.11.2
-Summary: TDCR model
-Home-page: https://pypi.org/project/TDCRPy/
-Author: RomainCoulon (Romain Coulon)
-Author-email: <romain.coulon@bipm.org>
-Project-URL: Documentation, https://github.com/RomainCoulon/TDCRPy/
-Keywords: Python,TDCR,Monte-Carlo,radionuclide,scintillation,counting
-Classifier: Development Status :: 4 - Beta
-Classifier: Intended Audience :: Science/Research
-Classifier: License :: OSI Approved :: MIT License
-Classifier: Natural Language :: English
-Classifier: Natural Language :: French
-Classifier: Programming Language :: Python :: 3
-Classifier: Operating System :: Unix
-Classifier: Operating System :: MacOS :: MacOS X
-Classifier: Operating System :: Microsoft :: Windows
-Classifier: Topic :: Scientific/Engineering :: Physics
-Description-Content-Type: text/markdown
-License-File: LICENCE.md
-Requires-Dist: numpy
-Requires-Dist: tqdm
-Requires-Dist: setuptools
-Requires-Dist: scipy
-Requires-Dist: configparser
-Requires-Dist: importlib.resources
-Requires-Dist: matplotlib
-
-# TDCRPy
-
-`TDCRPy` is a Python code to calculate detection efficiency of a liquid scintillation counter using 3-photomultiplier tubes.
-The calculation is based on the photo-physical model called of the Triple-to-Double-Coincidence-Ratio method (TDCR) and a Monte-Carlo sampling allowing to adress complexe decay schemes and radionuclide mixtures.
-
-The code is developped and maintained by the BIPM (MIT license).
-
-## Installation
-
-`TDCRPy` requires that the following packages are installed in your `Python` environement.
-
-```shell
-pip install importlib.resources configparser numpy tqdm setuptools scipy
-```
-or in `conda` environement:
-
-```shell
-conda install importlib.resources configparser numpy tqdm setuptools scipy
-```
-
-Then, `TDCRPy` can be installed.
-
-```shell
-pip install TDCRPy
-```
-
-To obtain the last version.
-
-```shell
-pip install TDCRPy --upgrade
-```
-
-The module can be imported in your Python code such as.
-
-```python
-import tdcrpy
-```
-
-## Test
-
-To run the unit tests of the package:
-
-```shell
-python -m unittest tdcrpy.test.test_tdcrpy
-```
-
-or (using the coverage package)
-
-```shell
-coverage run -m unittest tdcrpy.test.test_tdcrpy
-coverage report -m
-```
-
-## User guide
-
-At first, the TDCRPy module must be imported.
-
-```python
-import tdcrpy
-```
-### TDCRPy()
-
-The main function of the TDCRPy module is `TDCRPy()` in which the Monte-Carlo Triple-to-Double Coincidence Ratio model is implemented. The computation is made for a given solution containing a radionuclide (or a mixture of radionuclides) `rad`, a given volume of scintillator `V` and a given Birks constant `kB`.
-
-It can operates in two modes:
-*  In `mode="eff"`, it calculates the efficiency of the TDCR system as a function of a value (or a triplet) of the free parameter(s) `L`. The measurement data is not used is this mode.
-*  In `mode="res"`, it calculates the residual of the TDCR model parametrized by a value (or triplet) of the free parameter(s) `L` and the measurement data `TD`, `TAB`, `TBC`, `TAC`. This residual is used in the optimisation procedure `tdcrpy.TDCRoptimize.eff()` that aims to estimate the dection efficiency of the TDCR system (see above).
-
-Also, two configurations can be set:
-
-* In `mode2="sym"`, the symmetry is considered between the 3 photomultiplier tubes (PMTs). Here, the free parameter `L` is a scalar and only the global TDCR value `TD` is used as measurement data.
-* In `mode2="asym"`, a possible asymmetry between the 3 photomultiplier tubes is considered in the model. Here, the free parameter `L` is a triplet corresponding to each PMTs and only the specific TDCR values `TAB`, `TBC` and `TAC` are used as measurement data.
-
-The parmeter `N` sets the number of Monte-Carlo trails used for the estimation. Each MC trial corresponds to a simulated radiactive decay.
-
-TDCRPy() is parametrised by the following parameters:
-*  `L` is the free parameter(s) in keV^-1^. It represents the number of photoelectrons produced by the PMTs per unit of ionisation energy released in the scintillator during a decay without considering the scintillation quenching. It is of type `float` when `mode2="sym"` and of type `tuple` or `list` when `mode2="asym"`.
-*  `TD` (type `float`) is the quotient of the triple coincidence count rate divided by the logic sum of the double coincidence count rate.
-*  `TAB` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel A and B.
-*  `TBC` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel B and C.
-*  `TAC` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel A and C.
-*  `rad` (type `string`) is the list of radionuclides contained in the solution, e.g. `"H-3"`, `"H-3, Co-60"`.
-*  `pmf_1` (type `string`) is the probability of each radionuclide contained in `rad`. If `rad = "H-3, Co-60"` and `pmf_1="0.8, 0.2"`, the solution contains 80 % of tritium and 20 % of cobalt 60.
-*  `N` (type `integer`) is the number of Monte-Carlo trials corresponding each to a simulated nuclear decay in the scintillator. Monte-Carlo calculation is not used in the case of pure beta emitting radionuclides where the analytical model is implemented.
-*  `kB` (type `float`) is the Birks constant caracterizing the scintillator. It is espressed in cm/keV.
-*  `V` (type `float`) is the volume in ml of the scintillator.
-*  `mode` sets the type of output: `mode="res"` to return the residual, `mode="eff"` to return efficiencies.
-*  `mode2` sets whether the TDCR system has to be considered as symetrical by the model `mode2="sym"` or not symetrical `mode2="asym"`.
-*  `Display` (type `boolean`) is an optional parameter set by default, `Display=False`. If `Display=True`, detailed on the simulated decays are displayed.
-*  `barp` (type `boolean`) is an optional parameter set by default, `barp=False`. If `barp=True`, the progression bar of the calculation is displayed.
-
-
-In  `mode="eff"`, `TDCRPy()` returns a `tuple` composed of:
-*  the estimation of the counting efficiency of single events *S*,
-*  the standard uncertainty of the estimation of *S*,
-*  the estimation of the counting efficiency of double coincidences *D*,
-*  the standard uncertainty of the estimation of *D*,
-*  the estimation of the counting efficiency of triple coincidences *T*,
-*  the standard uncertainty of the estimation of *T*,
-
-Example: 
-
-```python
-import tdcrpy
-
-L = 1.5
-TD = 0.977667386529166
-TAB = 0.992232838598821
-TBC = 0.992343419459002
-TAC = 0.99275350064608
-Rad="Co-60"
-pmf_1="1"
-N = 100
-kB =1.0e-5
-V = 10
-mode = "eff"
-mode2 = "sym"
-
-result = tdcrpy.TDCRPy.TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2)
-
-print("efficiency S = ", round(result[0],4), "+/-", round(result[1],4))
-print("efficiency D = ", round(result[2],4), "+/-", round(result[3],4))
-print("efficiency T = ", round(result[4],4), "+/-", round(result[5],4))
-```
-
-```Console
-efficiency S =  0.9875 +/- 0.007
-efficiency D =  0.9903 +/- 0.0072
-efficiency T =  0.9749 +/- 0.0121
-```
-
-In  `mode="res"`, TDCRPy() returns the residuals *R* (type `float`) of the TDCR model for the given value of the free parameter `L`.
-
-Example: 
-
-```python
-import tdcrpy
-
-L = 1.5
-TD = 0.977667386529166
-TAB = 0.992232838598821
-TBC = 0.992343419459002
-TAC = 0.99275350064608
-Rad="Co-60"
-pmf_1="1"
-N = 100
-kB =1.0e-5
-V = 10
-mode = "res"
-mode2 = "sym"
-
-result = tdcrpy.TDCRPy.TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2)
-
-print("R = ",result)
-```
-
-```Console
-R =  9.72238771394384e-06
-```
-
-When a triplet of free parameters is considered, `L` is of type `tuple` and `mode2 = "asym"`. It must be noted that the calculation time is three times higher than in the symetrical model. 
-
-Example: 
-
-```python
-import tdcrpy
-
-L = (1.5, 1.2, 1.4)
-TD = 0.977667386529166
-TAB = 0.992232838598821
-TBC = 0.992343419459002
-TAC = 0.99275350064608
-Rad="Co-60"
-pmf_1="1"
-N = 100
-kB =1.0e-5
-V = 10
-mode = "res"
-mode2 = "asym"
-
-result = tdcrpy.TDCRPy.TDCRPy(L, TD, TAB, TBC, TAC, Rad, pmf_1, N, kB, V, mode, mode2)
-
-print("R = ",result)
-```
-
-```Console
-R =  1.4353680366705997e-05
-```
-
-By default, the Monte-Carlo model is not applied for the following list radionuclides when used solely (not in a mixture): ^3^H, ^14^C, ^35^S, ^45^Ca, ^63^Ni, ^89^Sr, ^90^Sr, ^99^Tc, ^147^Pm, ^241^Pu. In the latter cases, the analytical model is applied and the calculation time does not dependent to the parameter `N`.
-
-### Advance settings
-
-An advanced setting can be configured in the `config.toml` file for functions `TDCRPy.TDCRPy()` and `TDCRoptimize.eff()`. In this file:
-*  By default `Y = True` so that the analytical model is applied for solution containing only pure beta emitting radionuclides. If you would like to apply the MC calculation also for these nuclides, set `Y = False`.
-*  The list of radionuclides for which the analytical model is applied is defined by default such as `radListPureBeta = H-3, C-14, S-35, Ca-45, Ni-63, Sr-89, Sr-90, Tc-99, Pm-147, Pu-241`.
-*  The number of bins to discretize the linear energy space for quenching calculation for the radionuclides listed above has been set to induce an error from numerical approximation below 10^-4^. Thus the parameter `nE = 7000, 1000, 1000, 500, 2000, 500, 200, 500, 1000, 7000`.
-*  In the case of Monte-Carlo calculation, the number of bins to discretize the linear energy space for quenching calculation can be adjusted. `nE_electron` and `nE_alpha` parameters for respectively electrons and alpha particles are respectiveley set by default such as `nE_electron = 1000`
-and `nE_alpha = 1000`. These values ensure an error on quenched energy estimation from numerical approximation below 10^-3^.
-* By default the calculation is set for Ultima-Gold cocktail mixed with a small amount of aqueous solution. You can adapt for a specific scintillator by changing the `density` (default `density=0.96`), the mean charge number `Z` (default `Z=5.2`) and the mean mass number `A` (default `A=11.04`) of the scintillator.
-* To process photoelectric interactions, the atomic concentration of the scintillator is spectified. By default `pH = 0.578772`, `pC = 0.338741`, `pN = 0.000302`, `pO = 0.082022`, `pP = 0.000092`, `pCl = 0.000071`
-* A correction due reverse micelles in the scintillator is implemented. The diameter of the micelles and the fraction of aqueous solution is specified. By defaut `diam_micelle = 2`, `fAq = 0.1`
-* To optimize the speed of the Monte-Carlo calculation, a spline interpolation on precalculated quenched energy is applied. The parameter `depthSpline` (default `depthSpline = 5`) sets the number of bins on each side of the energy point on which the interpolation is applied. The parameter `Einterp` (default `Einterp = 1`) set the energy (in keV) above which the interpolation is applied.
-
-### TDCRoptimize.eff()
-
-An optimization procedure `TDCRoptimize.eff()` is proposed to estimate the counting efficiency of a TDCR system. It approximates the free parameter(s) `L` in order to minimize the residual estimated by `TDCRPy()` for a given set of measurement data. In `mode2 = "sym"`, the *golden section search technique* is applied with a maximum number of iterations of 20. In `mode2 = "asym"`, the *Nelder-Mead technique* is applied with a maximum number of iterations of 20.
-
-`TDCRoptimize.eff()` is parametrised by the following parameters:
-*  `TD` (type `float`) is the quotient of the triple coincidence count rate divided by the logic sum of the double coincidence count rate.
-*  `TAB` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel A and B.
-*  `TBC` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel B and C.
-*  `TAC` (type `float`) is the quotient of the triple coincidence count rate divided by the double coincidence count rate between the channel A and C.
-*  `rad` (type `string`) is the list of radionuclides contained in the solution, e.g. `"H-3"`, `"H-3, Co-60"`.
-*  `pmf_1` (type `string`) is the probability of each radionuclide contained in `rad`. If `rad = "H-3, Co-60"` and `pmf_1="0.8, 0.2"`, the solution contains 80 % of tritium and 20 % of cobalt 60.
-*  `kB` (type `float`) is the Birks constant caracterizing the scintillator. It is espressed in cm/keV.
-*  `V` (type `float`) is the volume in ml of the scintillator.
-*  `mode2` sets whether the TDCR system has to be considered symetrical `mode2="sym"` or not symetrical `mode2="asym"`.
-*  `N` (type `integer`) is the number of Monte-Carlo trials corresponding each, to a simulated nuclear decay in the scintillator. Monte-Carlo calculation is not used in the case of pure beta emitting radionuclides where the analytical model is implemented. It is an optional parameter set by default at `N = 10000`.
-*  `L` (type `float`) initial guess of the free parameter(s) in keV^-1^. It represents the number of photoelectrons produced by the PMTs per unit of ionisation energy released in the scintillator during a decay and in the abscence of scintillation quenching. It is set by default at `L = 1`.
-
-`TDCRoptimize.eff()` returns:
-*  the global free parameter *L*,
-*  the free parameters related to each channel (*L<sub>A</sub>*, *L<sub>B</sub>*, *L<sub>C</sub>*),
-*  the estimation of the counting efficiency of single events *S*,
-*  the standard uncertainty of the estimation of *S*,
-*  the estimation of the counting efficiency of double coincidences *D*,
-*  the standard uncertainty of the estimation of *D*,
-*  the estimation of the counting efficiency of triple coincidences *T*,
-*  the standard uncertainty of the estimation of *T*.
-
-Example with `mode2 = "sym"`: 
-
-```python
-import tdcrpy
-
-TD = 0.977667386529166
-TAB = 0.992232838598821
-TBC = 0.992343419459002
-TAC = 0.99275350064608
-Rad="Co-60"
-pmf_1="1"
-N = 250
-kB =1.0e-5
-V = 10
-mode2 = "sym"
-
-result = tdcrpy.TDCRoptimize.eff(TD, TAB, TBC, TAC, Rad, pmf_1, kB, V, mode2, N=N)
-
-print("Global free parameter = \t", round(result[0],4), " keV-1")
-print("Free parameter (PMT A) = \t", round(result[1][0],4) , " keV-1")
-print("Free parameter (PMT B) = \t", round(result[1][1],4) , " keV-1")
-print("Free parameter (PMT C) = \t", round(result[1][2],4) , " keV-1")
-print("efficiency S = \t", round(result[2],4), "+/-", round(result[3],4))
-print("efficiency D = \t", round(result[4],4), "+/-", round(result[5],4))
-print("efficiency T = \t", round(result[6],4), "+/-", round(result[7],4))
-```
-
-```Console
-Global free parameter = 	 1.3061  keV-1
-Free parameter (PMT A) = 	 1.3061  keV-1
-Free parameter (PMT B) = 	 1.3061  keV-1
-Free parameter (PMT C) = 	 1.3061  keV-1
-efficiency S = 	 0.979 +/- 0.0067
-efficiency D = 	 0.9833 +/- 0.0064
-efficiency T = 	 0.963 +/- 0.0096
-```
-
-Example with `mode2 = "asym"`: 
-
-```python
-import tdcrpy
-
-TD = 0.977667386529166
-TAB = 0.992232838598821
-TBC = 0.992343419459002
-TAC = 0.99275350064608
-Rad="Co-60"
-pmf_1="1"
-N = 250
-kB =1.0e-5
-V = 10
-mode2 = "asym"
-
-result = tdcrpy.TDCRoptimize.eff(TD, TAB, TBC, TAC, Rad, pmf_1, kB, V, mode2, N=N)
-
-print("Global free parameter = \t", round(result[0],4), " keV-1")
-print("Free parameter (PMT A) = \t", round(result[1][0],4) , " keV-1")
-print("Free parameter (PMT B) = \t", round(result[1][1],4) , " keV-1")
-print("Free parameter (PMT C) = \t", round(result[1][2],4) , " keV-1")
-print("efficiency S = \t", round(result[2],4), "+/-", round(result[3],4))
-print("efficiency D = \t", round(result[4],4), "+/-", round(result[5],4))
-print("efficiency T = \t", round(result[6],4), "+/-", round(result[7],4))
-```
-
-```Console
-Global free parameter = 	 1.3061  keV-1
-Free parameter (PMT A) = 	 1.3061  keV-1
-Free parameter (PMT B) = 	 1.3061  keV-1
-Free parameter (PMT C) = 	 1.3061  keV-1
-efficiency S = 	 0.979 +/- 0.0067
-efficiency D = 	 0.9833 +/- 0.0064
-efficiency T = 	 0.963 +/- 0.0096
-```
-