Applying fit-and-count to mitigate background¶
Sometimes the sideband subtraction technique might not be suitable, for example, if peaking backgrounds are present in the sample of interest. In such cases, a fit-and-count approach must be employed instead, making use of the knowledge of the components in the distribution to model the signal and each background.
In this tutorial we'll use the RooFit
implementation of the fitting.
As before, we start with a fairly standard set of imports:
import matplotlib.pyplot as plt
import mplhep as hep
import numpy as np
import ROOT as R
from triggercalib import HltEff
from triggercalib.utils.helpers import tgraph_to_np
hep.style.use("LHCb2")
R.EnableImplicitMT(8)
Welcome to JupyROOT 6.30/04
Considering again the example of the Hlt1TrackMVA
TOS efficiency in $B^+\to J/\psi\left(\mu\mu\right)K^+$ 2024 MC simulation from the previous tutorial, we will use the $J/\psi K^+$ invariant mass to separate signal and background. We define an observable and PDF, in this case using RooWorkspace (but in principle this can be an existing observable/PDF):
model_description = {
"signal" : { # Name of component
"model" : "CrystalBall", # Name of PDF to be used
"params" : { # Parameters to be used
"mu" : [5279.4, 5275, 5285],
"sigma" : [8.75, 6, 12],
"alphal" : [1.4, 1.4, 1.4],
"nl" : [7, 7, 7],
"alphar" : [1.4, 1.4, 1.4],
"nr" : [7, 7, 7],
},
"yield" : [10000, 0, 1000000] # Associated yield
},
"combinatorial" : {
"model": "Exponential",
"params" : {
"exponent" : [-0.001, -0.01, -0.0001]
},
"yield": [1000, 0, 1000000]
}
}
def build_roofit_pdf(ws, model_description):
model_build_string = "SUM::pdf("
for component_name, component_description in model_description.items():
component_build_string = f"{component_description['model']}::{component_name}_pdf(B_DTF_Jpsi_MASS"
for param_name, (param_val, param_min, param_max) in component_description['params'].items():
component_build_string += f", {component_name}_{param_name}[{param_val}, {param_min}, {param_max}]"
component_build_string += ")"
print(f"Constructing '{component_name}_pdf' using description '{component_build_string}'")
ws.factory(component_build_string)
if "yield" in component_description:
yield_val, yield_min, yield_max = component_description["yield"]
model_build_string += f", {component_name}_yield[{yield_val}, {yield_min}, {yield_max}]*{component_name}_pdf"
model_build_string += ")"
print(f"Constructing 'pdf' using description '{model_build_string}'")
pdf = ws.factory(model_build_string)
return ws, pdf
ws = R.RooWorkspace("ws")
observable = ws.factory("B_DTF_Jpsi_MASS[5200, 5500]")
ws, pdf = build_roofit_pdf(ws, model_description)
Constructing 'signal_pdf' using description 'CrystalBall::signal_pdf(B_DTF_Jpsi_MASS, signal_mu[5279.4, 5275, 5285], signal_sigma[8.75, 6, 12], signal_alphal[1.4, 1.4, 1.4], signal_nl[7, 7, 7], signal_alphar[1.4, 1.4, 1.4], signal_nr[7, 7, 7])' Constructing 'combinatorial_pdf' using description 'Exponential::combinatorial_pdf(B_DTF_Jpsi_MASS, combinatorial_exponent[-0.001, -0.01, -0.0001])' Constructing 'pdf' using description 'SUM::pdf(, signal_yield[10000, 0, 1000000]*signal_pdf, combinatorial_yield[1000, 0, 1000000]*combinatorial_pdf)'
Now that we have an observable and a PDF, we can pass these to HltEff
in place of our sideband from the previous tutorial, opting for a more coarse binning scheme to reduce the number of fits to run:
binning = {
"B_PT" : {
"bins" : [
2e3, 7e3, 15e3, 25e3
]
}
}
hlt_eff = HltEff(
"simple_example",
"root://eoslhcb.cern.ch//eos/lhcb/wg/rta/WP4/TriggerCalib/Bu2JpsiK_Jpsi2MuMu_block1_ntuple.root:Tuple/DecayTree",
tos="Hlt1TrackMVA",
tis=["Hlt1TrackMVA", "Hlt1TwoTrackMVA"],
particle="B",
binning=binning,
observable=observable,
pdf=pdf,
fit_kwargs={
"NumCPU" : 8,
},
output_path="fit-and-count/",
plots=True,
trig_effs=False, # Avoids running fit for 'trig' category (only interested in TOS and TISTOS here)
)
================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0019479496497944945 +/- 0.0001398536284415126 combinatorial_yield: 8939.730030214443 +/- 118.95428179855662 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.58383661268 +/- 0.05657418232794953 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.498370952601476 +/- 0.05159964918657067 signal_yield: 34447.09229525895 +/- 199.14301459686612 Covariance quality: 3 Fit status: 0 Minimum value: -224003.3205400036 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0021460920243351633 +/- 0.00027751217600144947 combinatorial_yield: 2505.695241256127 +/- 67.54951098201764 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.616635768485 +/- 0.07350563085174144 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 9.074356636331393 +/- 0.06517521521318237 signal_yield: 22331.18114545947 +/- 156.1666988607285 Covariance quality: 3 Fit status: 0 Minimum value: -120805.78413112763 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0018860132144080273 +/- 0.000844114542738975 combinatorial_yield: 325.73417220874876 +/- 27.457159529284326 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.583751028043 +/- 0.188220890602679 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 10.537991340554916 +/- 0.16849046606365903 signal_yield: 4555.204043022127 +/- 70.59292477509416 Covariance quality: 3 Fit status: 0 Minimum value: -15723.029343819799 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0023363354777844263 +/- 9.676374823723967e-05 combinatorial_yield: 19443.354076134186 +/- 182.0959750891834 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.659574112436 +/- 0.026263711439696635 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.430851030049693 +/- 0.023261100111540145 signal_yield: 151066.3910102272 +/- 405.9316423664859 Covariance quality: 3 Fit status: 0 Minimum value: -1160911.529720019 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.002570180143229027 +/- 0.0001556359633215537 combinatorial_yield: 8309.289442341917 +/- 128.17279229784253 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.563236534862 +/- 0.0313177604939483 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 9.083662249873004 +/- 0.027336050548275814 signal_yield: 121353.85714670461 +/- 359.8160055667904 Covariance quality: 3 Fit status: 0 Minimum value: -860942.5477918121 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.002503613167499662 +/- 0.0005164491111190777 combinatorial_yield: 909.7706073343844 +/- 47.85608590543768 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.7173220264185 +/- 0.08797370271167892 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 10.308258013188594 +/- 0.0775224508349508 signal_yield: 19681.57555123606 +/- 145.12700433012833 Covariance quality: 3 Fit status: 0 Minimum value: -97929.82790149421 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0017220917869007226 +/- 0.0002387810696026885 combinatorial_yield: 3165.831841811986 +/- 72.17777111336818 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.586740330383 +/- 0.07685231653385927 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.54846704100971 +/- 0.06879488729582661 signal_yield: 18343.344102697196 +/- 142.78374251456444 Covariance quality: 3 Fit status: 0 Minimum value: -99409.06509904448 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0025832400222510197 +/- 0.00032815957720996555 combinatorial_yield: 1809.548117928356 +/- 58.505975219811376 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.609737614319 +/- 0.07833497507499487 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 9.097931462042459 +/- 0.06910815907775092 signal_yield: 19631.474984597022 +/- 145.7561685602086 Covariance quality: 3 Fit status: 0 Minimum value: -102277.69740542985 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0019559038530267622 +/- 0.0009253254542206356 combinatorial_yield: 276.27707140758145 +/- 25.710241902352806 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.5694736695 +/- 0.1944470010785153 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 10.591972534167349 +/- 0.17348899565549125 signal_yield: 4297.630609758679 +/- 68.42773488975581 Covariance quality: 3 Fit status: 0 Minimum value: -14514.169464213992 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0022786326009264627 +/- 5.896338930000021e-05 combinatorial_yield: 50924.5958217767 +/- 288.27164453997466 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.604744121437 +/- 0.01954010923145688 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.395964322098694 +/- 0.017543274778422102 signal_yield: 275660.1020342683 +/- 554.8282581963576 Covariance quality: 3 Fit status: 0 Minimum value: -2395066.2159738475 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.002306157660256243 +/- 0.00013511357357394536 combinatorial_yield: 10896.688192016874 +/- 144.59211318826328 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.563213449513 +/- 0.02964911401431891 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 9.067755680338443 +/- 0.026014897073991072 signal_yield: 135538.46355305746 +/- 381.5083314955409 Covariance quality: 3 Fit status: 0 Minimum value: -984891.3221755793 ================================================ ================================================ Fit performed with RooFit from ROOT 6.30/04 Initial parameters: combinatorial_exponent: -0.001 +/- 0.0 combinatorial_yield: 1000.0 +/- 0.0 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.4 +/- 0.0 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 8.75 +/- 0.0 signal_yield: 10000.0 +/- 0.0 Final parameters: combinatorial_exponent: -0.0024835852920021993 +/- 0.0004863488979532124 combinatorial_yield: 1015.1951798644934 +/- 50.12386147989241 signal_alphal: 1.4 +/- 0.0 signal_alphar: 1.4 +/- 0.0 signal_mu: 5279.709818886293 +/- 0.08607936397402227 signal_nl: 7.0 +/- 0.0 signal_nr: 7.0 +/- 0.0 signal_sigma: 10.308850099973359 +/- 0.07600485111794164 signal_yield: 20594.426591022944 +/- 148.62652858290858 Covariance quality: 3 Fit status: 0 Minimum value: -103618.57044304948 ================================================ [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 265048 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -224003.320540003595 Edm = 3.35333771719563317e-05 Nfcn = 154 combinatorial_exponent = -0.00194795 +/- 0.00013981 (limited) combinatorial_yield = 8939.73 +/- 118.91 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.58 +/- 0.0565708 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 8.49837 +/- 0.05158 (limited) signal_yield = 34447.1 +/- 199.111 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_1) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_1) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_1) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_1) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 100540 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -120805.784131127628 Edm = 2.2122307114079877e-06 Nfcn = 148 combinatorial_exponent = -0.00214609 +/- 0.000277398 (limited) combinatorial_yield = 2505.7 +/- 67.5147 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.62 +/- 0.0735076 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 9.07436 +/- 0.0651565 (limited) signal_yield = 22331.2 +/- 156.149 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_2) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_2) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_2) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_2) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 17653 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -15723.0293438197987 Edm = 4.87513202232002984e-06 Nfcn = 145 combinatorial_exponent = -0.00188601 +/- 0.000843077 (limited) combinatorial_yield = 325.734 +/- 27.4197 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.58 +/- 0.188223 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 10.538 +/- 0.168366 (limited) signal_yield = 4555.2 +/- 70.5767 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_3) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_3) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_3) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tis_pdf_bin_3) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 692882 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -1160911.52972001908 Edm = 0.000199340056770105407 Nfcn = 153 combinatorial_exponent = -0.00233634 +/- 9.67431e-05 (limited) combinatorial_yield = 19443.4 +/- 182.036 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.66 +/- 0.0262636 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 8.43085 +/- 0.0232559 (limited) signal_yield = 151066 +/- 405.895 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_1) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_1) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_1) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_1) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 391847 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -860942.547791812103 Edm = 0.000685838985485316753 Nfcn = 207 combinatorial_exponent = -0.00257018 +/- 0.000155644 (limited) combinatorial_yield = 8309.29 +/- 128.153 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.56 +/- 0.0313189 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 9.08366 +/- 0.027329 (limited) signal_yield = 121354 +/- 359.801 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_2) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_2) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_2) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_2) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 51529 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -97929.8279014942091 Edm = 5.21908757081951673e-05 Nfcn = 145 combinatorial_exponent = -0.00250361 +/- 0.000516196 (limited) combinatorial_yield = 909.771 +/- 47.821 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.72 +/- 0.0879671 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 10.3083 +/- 0.0774744 (limited) signal_yield = 19681.6 +/- 145.112 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_3) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_3) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_3) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tos_pdf_bin_3) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 101073 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -99409.0650990444847 Edm = 0.000409099201074414884 Nfcn = 119 combinatorial_exponent = -0.00172209 +/- 0.000238665 (limited) combinatorial_yield = 3165.83 +/- 72.142 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.59 +/- 0.0768543 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 8.54847 +/- 0.0687777 (limited) signal_yield = 18343.3 +/- 142.763 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_1) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_1) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_1) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_1) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 81286 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -102277.697405429848 Edm = 1.21401104681273255e-05 Nfcn = 132 combinatorial_exponent = -0.00258324 +/- 0.000328115 (limited) combinatorial_yield = 1809.55 +/- 58.4895 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.61 +/- 0.0783411 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 9.09793 +/- 0.069101 (limited) signal_yield = 19631.5 +/- 145.747 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_2) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_2) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_2) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_2) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 16050 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -14514.1694642139919 Edm = 5.65412107022127061e-06 Nfcn = 146 combinatorial_exponent = -0.0019559 +/- 0.000924176 (limited) combinatorial_yield = 276.277 +/- 25.6745 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.57 +/- 0.194434 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 10.592 +/- 0.173346 (limited) signal_yield = 4297.63 +/- 68.4127 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_3) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_3) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_3) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(tistos_pdf_bin_3) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 1737058 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -2395066.21597384755 Edm = 2.74087802244680335e-05 Nfcn = 206 combinatorial_exponent = -0.00227863 +/- 5.8949e-05 (limited) combinatorial_yield = 50924.6 +/- 288.161 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.6 +/- 0.0195398 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 8.39596 +/- 0.0175378 (limited) signal_yield = 275660 +/- 554.75 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_1) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_1) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_1) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_1) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 472177 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -984891.322175579262 Edm = 1.64288719896598195e-05 Nfcn = 201 combinatorial_exponent = -0.00230616 +/- 0.000135144 (limited) combinatorial_yield = 10896.7 +/- 144.592 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.56 +/- 0.0296493 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 9.06776 +/- 0.0260078 (limited) signal_yield = 135538 +/- 381.497 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_2) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_2) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_2) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_2) indirectly selected PDF components: () [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Eval -- RooAbsTestStatistic::initMPMode: started 8 remote server process. [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: activating const optimization [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer Minuit2Minimizer: Minimize with max-calls 4500 convergence for edm < 1 strategy 1 [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#0] WARNING:DataHandling -- RooDataSet.from_numpy() Ignored 56174 out-of-range events [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) Minuit2Minimizer : Valid minimum - status = 0 FVAL = -103618.570443049481 Edm = 0.000467281445584336519 Nfcn = 133 combinatorial_exponent = -0.00248359 +/- 0.000486117 (limited) combinatorial_yield = 1015.2 +/- 50.0894 (limited) signal_alphal = 1.4 (fixed) signal_alphar = 1.4 (fixed) signal_mu = 5279.71 +/- 0.0860778 (limited) signal_nl = 7 (fixed) signal_nr = 7 (fixed) signal_sigma = 10.3089 +/- 0.0759689 (limited) signal_yield = 20594.4 +/- 148.611 (limited) [#1] INFO:Minimization -- RooAbsMinimizerFcn::optimizeConstantTerms: set of constant parameters changed, rerunning const optimizer [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- The following expressions will be evaluated in cache-and-track mode: (signal_pdf,combinatorial_pdf) [#1] INFO:Minimization -- RooAbsMinimizerFcn::setOptimizeConst: deactivating const optimization [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_3) directly selected PDF components: (signal_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_3) indirectly selected PDF components: () [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_3) directly selected PDF components: (combinatorial_pdf) [#1] INFO:Plotting -- RooAbsPdf::plotOn(sel_pdf_bin_3) indirectly selected PDF components: ()
Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -195262.8459 Edm = 23300.54015 NCalls = 21 Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -195262.8459 Edm : 23300.54015 Internal parameters: [ 0.9582415885 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -181.2837747 -188709.0707 -232.4906351 -777.3608554 -262848.1493] Internal covariance matrix: [[ 0.0032699471 0 0 0 0] [ 0 8.094866e-07 0 0 0] [ 0 0 0.0002559031 0 0] [ 0 0 0 0.00051101983 0] [ 0 0 0 0 9.2554826e-07]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -195262.8459 Edm = 23300.54015 NCalls = 21 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -223758.3051 Edm = 98.39807611 NCalls = 41 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -223912.1374 Edm = 13.17776887 NCalls = 54 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -224000.0944 Edm = 1.747880525 NCalls = 71 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -224002.811 Edm = 1.070006864 NCalls = 83 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -224003.2692 Edm = 0.04680936367 NCalls = 95 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -224003.319 Edm = 0.001101787322 NCalls = 107 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -224003.3205 Edm = 4.124648135e-05 NCalls = 119 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -224003.3205 Edm = 3.353337717e-05 NCalls = 154 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tis_bin_1.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -114374.5852 Edm = 6066.301035 NCalls = 19 Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -114374.5852 Edm : 6066.301035 Internal parameters: [ 0.9582415885 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 38.76585677 -44542.77664 -199.4251497 -357.7953218 -123656.8477] Internal covariance matrix: [[ 0.034065822 0 0 0 0] [ 0 2.1146555e-06 0 0 0] [ 0 0 0.00041633947 0 0] [ 0 0 0 0.00074650378 0] [ 0 0 0 0 1.3018288e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -114374.5852 Edm = 6066.301035 NCalls = 19 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -120618.0855 Edm = 125.8061389 NCalls = 32 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -120742.4345 Edm = 70.07966221 NCalls = 46 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -120798.8684 Edm = 4.690242208 NCalls = 61 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -120804.2359 Edm = 2.341996276 NCalls = 73 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -120805.5764 Edm = 0.1513368182 NCalls = 85 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -120805.7833 Edm = 0.0007891701116 NCalls = 97 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -120805.7841 Edm = 2.67989797e-06 NCalls = 109 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -120805.7841 Edm = 2.212230711e-06 NCalls = 148 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tis_bin_2.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -13489.41817 Edm = 2530.042313 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -13500.60228 Edm : 2547.177755 Internal parameters: [ 0.2336312704 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.8422405223 17669.37904 -49.92449445 -251.4373404 55320.95507] Internal covariance matrix: [[ 0.033172914 0 0 0 0] [ 0 3.7006701e-06 0 0 0] [ 0 0 0.0023678147 0 0] [ 0 0 0 0.0033388038 0] [ 0 0 0 0 2.8807692e-06]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -13500.60228 Edm : 2547.177755 Internal parameters: [ 0.2336312704 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.8422405223 17669.37904 -49.92449445 -251.4373404 55320.95507] Internal covariance matrix: [[ 0.033172914 0 0 0 0] [ 0 3.7006701e-06 0 0 0] [ 0 0 0.0023678147 0 0] [ 0 0 0 0.0033388038 0] [ 0 0 0 0 2.8807692e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -13500.60228 Edm = 2547.177755 NCalls = 38 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -15706.34952 Edm = 15.78237936 NCalls = 50 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -15719.92666 Edm = 1.518067518 NCalls = 61 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -15722.71402 Edm = 0.1705294222 NCalls = 73 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -15723.0042 Edm = 0.01701132143 NCalls = 86 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -15723.02711 Edm = 0.001934639248 NCalls = 98 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -15723.02934 Edm = 5.301517747e-06 NCalls = 110 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -15723.02934 Edm = 4.875132022e-06 NCalls = 145 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tis_bin_3.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -852705.3968 Edm = 165178.6774 NCalls = 23 Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -852705.3968 Edm : 165178.6774 Internal parameters: [ 0.9582415885 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 131.9124624 -514838.9199 -1675.802508 777.3853418 -1425191.012] Internal covariance matrix: [[ 0.0030928629 0 0 0 0] [ 0 4.8861382e-07 0 0 0] [ 0 0 5.8245433e-05 0 0] [ 0 0 0 0.00011968512 0] [ 0 0 0 0 2.6138314e-07]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -852705.3968 Edm = 165178.6774 NCalls = 23 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -1153624.483 Edm = 2497.265797 NCalls = 44 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -1154144.235 Edm = 4323.220374 NCalls = 58 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -1156631.809 Edm = 874.2364272 NCalls = 73 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -1160386.776 Edm = 222.3281845 NCalls = 85 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -1160789.23 Edm = 80.94119766 NCalls = 97 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -1160911.38 Edm = 0.1515169761 NCalls = 109 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -1160911.53 Edm = 0.000161139365 NCalls = 120 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -1160911.53 Edm = 0.0001993400568 NCalls = 153 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tos_bin_1.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -659103.9583 Edm = 120940.2044 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -659342.2806 Edm : 125628.3483 Internal parameters: [ 0.2433437453 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.2576314985 -287257.3858 -957.3102763 -689.7931051 -1090363.912] Internal covariance matrix: [[ 0.001443161 0 0 0 0] [ 0 1.3393063e-06 0 0 0] [ 0 0 7.7321211e-05 0 0] [ 0 0 0 0.00014050076 0] [ 0 0 0 0 3.2960094e-07]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -659342.2806 Edm : 125628.3483 Internal parameters: [ 0.2433437453 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.2576314985 -287257.3858 -957.3102763 -689.7931051 -1090363.912] Internal covariance matrix: [[ 0.001443161 0 0 0 0] [ 0 1.3393063e-06 0 0 0] [ 0 0 7.7321211e-05 0 0] [ 0 0 0 0.00014050076 0] [ 0 0 0 0 3.2960094e-07]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -659342.2806 Edm = 125628.3483 NCalls = 40 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -837079.9478 Edm = 7864.124449 NCalls = 60 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -852915.9375 Edm = 2959.082479 NCalls = 78 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -860670.5729 Edm = 193.8871554 NCalls = 97 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -860880.2172 Edm = 65.14698084 NCalls = 109 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -860931.4713 Edm = 10.98531758 NCalls = 124 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -860942.3347 Edm = 0.1874791275 NCalls = 136 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -860942.5284 Edm = 0.01439179088 NCalls = 148 Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -860942.5445 Edm = 0.002249713767 NCalls = 160 Info in <Minuit2>: VariableMetricBuilder 9 - FCN = -860942.5478 Edm = 0.0003241545295 NCalls = 172 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 10 - FCN = -860942.5478 Edm = 0.0006858389855 NCalls = 207 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tos_bin_2.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -94202.47372 Edm = 3280.003223 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -94299.87929 Edm : 3405.218803 Internal parameters: [ -0.1532079556 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -3.467473302 -18026.10853 -309.3187552 -847.5584023 -89754.68264] Internal covariance matrix: [[ 0.009391545 0 0 0 0] [ 0 3.7497835e-06 0 0 0] [ 0 0 0.00053508646 0 0] [ 0 0 0 0.00078189115 0] [ 0 0 0 0 1.4634528e-06]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -94299.87929 Edm : 3405.218803 Internal parameters: [ -0.1532079556 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -3.467473302 -18026.10853 -309.3187552 -847.5584023 -89754.68264] Internal covariance matrix: [[ 0.009391545 0 0 0 0] [ 0 3.7497835e-06 0 0 0] [ 0 0 0.00053508646 0 0] [ 0 0 0 0.00078189115 0] [ 0 0 0 0 1.4634528e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -94299.87929 Edm = 3405.218803 NCalls = 39 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -97602.48086 Edm = 215.0726336 NCalls = 50 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -97884.9596 Edm = 22.31748554 NCalls = 62 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -97929.21305 Edm = 0.291928463 NCalls = 74 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -97929.80661 Edm = 0.01537319615 NCalls = 86 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -97929.82608 Edm = 0.001326997477 NCalls = 98 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -97929.8279 Edm = 2.950113863e-05 NCalls = 110 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -97929.8279 Edm = 5.219087571e-05 NCalls = 145 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tos_bin_3.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -95284.97634 Edm = 4053.142113 NCalls = 21 Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -95284.97634 Edm : 4053.142113 Internal parameters: [ 0.9582415885 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -40.76469293 -53058.0733 -144.8945474 -127.4041997 -87892.19118] Internal covariance matrix: [[ 0.010449099 0 0 0 0] [ 0 1.7218953e-06 0 0 0] [ 0 0 0.00048259444 0 0] [ 0 0 0 0.00096386193 0] [ 0 0 0 0 1.4656263e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -95284.97634 Edm = 4053.142113 NCalls = 21 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -99353.19424 Edm = 38.20056556 NCalls = 34 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -99404.81111 Edm = 3.127179361 NCalls = 48 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -99408.20575 Edm = 0.4331690586 NCalls = 62 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -99408.94635 Edm = 0.09469268657 NCalls = 74 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -99409.0651 Edm = 0.0004066515064 NCalls = 86 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -99409.0651 Edm = 0.0004090992011 NCalls = 119 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tistos_bin_1.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -98369.16877 Edm = 3051.883565 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -98390.28579 Edm : 3774.7178 Internal parameters: [ 0.3565745391 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -18.96837595 -29377.16326 -189.0841562 -216.6510617 -94638.78195] Internal covariance matrix: [[ 0.0077708578 0 0 0 0] [ 0 2.5735494e-06 0 0 0] [ 0 0 0.00047747404 0 0] [ 0 0 0 0.0008615879 0] [ 0 0 0 0 1.4310887e-06]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -98390.28579 Edm : 3774.7178 Internal parameters: [ 0.3565745391 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -18.96837595 -29377.16326 -189.0841562 -216.6510617 -94638.78195] Internal covariance matrix: [[ 0.0077708578 0 0 0 0] [ 0 2.5735494e-06 0 0 0] [ 0 0 0.00047747404 0 0] [ 0 0 0 0.0008615879 0] [ 0 0 0 0 1.4310887e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -98390.28579 Edm = 3774.7178 NCalls = 35 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -102147.5178 Edm = 88.92012008 NCalls = 46 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -102260.9081 Edm = 10.60499564 NCalls = 58 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -102276.9537 Edm = 0.5933993749 NCalls = 71 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -102277.6895 Edm = 0.00724483067 NCalls = 83 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -102277.6974 Edm = 1.434422605e-05 NCalls = 95 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -102277.6974 Edm = 1.214011047e-05 NCalls = 132 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tistos_bin_2.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -12018.10923 Edm = 2840.590598 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -12031.45088 Edm : 2855.000088 Internal parameters: [ 0.1420850049 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.8647156488 18980.40492 -47.23122664 -238.1702703 57962.85615] Internal covariance matrix: [[ 0.035617239 0 0 0 0] [ 0 3.7776867e-06 0 0 0] [ 0 0 0.0025243233 0 0] [ 0 0 0 0.0035286462 0] [ 0 0 0 0 2.932782e-06]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -12031.45088 Edm : 2855.000088 Internal parameters: [ 0.1420850049 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.8647156488 18980.40492 -47.23122664 -238.1702703 57962.85615] Internal covariance matrix: [[ 0.035617239 0 0 0 0] [ 0 3.7776867e-06 0 0 0] [ 0 0 0.0025243233 0 0] [ 0 0 0 0.0035286462 0] [ 0 0 0 0 2.932782e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -12031.45088 Edm = 2855.000088 NCalls = 39 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -14495.97753 Edm = 16.79243992 NCalls = 51 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -14510.62433 Edm = 1.649179459 NCalls = 62 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -14513.86336 Edm = 0.1567519937 NCalls = 74 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -14514.14442 Edm = 0.01615591762 NCalls = 87 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -14514.16721 Edm = 0.002040139052 NCalls = 99 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -14514.16946 Edm = 5.912137628e-06 NCalls = 111 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -14514.16946 Edm = 5.65412107e-06 NCalls = 146 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/tistos_bin_3.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -1599328.068 Edm = 349005.3813 NCalls = 25 Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -1599328.068 Edm : 349005.3813 Internal parameters: [ 0.9582415885 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ -348.1695901 -1285316.134 -2231.711112 -1241.726468 -2735865.2] Internal covariance matrix: [[ 0.00078109197 0 0 0 0] [ 0 1.8360926e-07 0 0 0] [ 0 0 3.1750023e-05 0 0] [ 0 0 0 6.5415922e-05 0] [ 0 0 0 0 1.4593781e-07]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -1599328.068 Edm = 349005.3813 NCalls = 25 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -2385210.24 Edm = 2275.716207 NCalls = 47 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -2389238.657 Edm = 1150.467739 NCalls = 60 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -2392100.183 Edm = 464.0210497 NCalls = 77 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -2392928.935 Edm = 260.2725711 NCalls = 89 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -2393892.514 Edm = 151.74881 NCalls = 101 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -2394584.546 Edm = 410.1559323 NCalls = 117 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -2395042.014 Edm = 25.54699849 NCalls = 130 Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -2395065.676 Edm = 0.7400024314 NCalls = 141 Info in <Minuit2>: VariableMetricBuilder 9 - FCN = -2395066.209 Edm = 0.006880492446 NCalls = 153 Info in <Minuit2>: VariableMetricBuilder 10 - FCN = -2395066.216 Edm = 2.133968214e-05 NCalls = 165 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 11 - FCN = -2395066.216 Edm = 2.740878022e-05 NCalls = 206 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/sel_bin_1.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -740909.0639 Edm = 139159.386 NCalls = 19 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -741068.839 Edm : 143696.9764 Internal parameters: [ 0.4073489348 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 8.381697189 -354918.0527 -1028.55819 -1038.071543 -1235834.43] Internal covariance matrix: [[ 0.0013240359 0 0 0 0] [ 0 9.6865019e-07 0 0 0] [ 0 0 6.9006675e-05 0 0] [ 0 0 0 0.0001255469 0] [ 0 0 0 0 2.9631746e-07]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -741068.839 Edm : 143696.9764 Internal parameters: [ 0.4073489348 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 8.381697189 -354918.0527 -1028.55819 -1038.071543 -1235834.43] Internal covariance matrix: [[ 0.0013240359 0 0 0 0] [ 0 9.6865019e-07 0 0 0] [ 0 0 6.9006675e-05 0 0] [ 0 0 0 0.0001255469 0] [ 0 0 0 0 2.9631746e-07]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -741068.839 Edm = 143696.9764 NCalls = 36 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -964134.112 Edm = 6733.463902 NCalls = 55 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -975511.0818 Edm = 3001.4222 NCalls = 72 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -982235.1038 Edm = 915.8312904 NCalls = 87 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -984158.9461 Edm = 837.9367547 NCalls = 100 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -984592.4981 Edm = 265.5531778 NCalls = 112 Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -984887.1457 Edm = 1.803460879 NCalls = 124 Info in <Minuit2>: VariableMetricBuilder 7 - FCN = -984890.5502 Edm = 0.4743475909 NCalls = 136 Info in <Minuit2>: VariableMetricBuilder 8 - FCN = -984891.3212 Edm = 0.0008595989124 NCalls = 148 Info in <Minuit2>: VariableMetricBuilder 9 - FCN = -984891.3222 Edm = 1.165419655e-06 NCalls = 160 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 10 - FCN = -984891.3222 Edm = 1.642887199e-05 NCalls = 201 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/sel_bin_2.pdf has been created Info in <Minuit2>: MnSeedGenerator Computing seed using NumericalGradient calculator Info in <Minuit2>: MnSeedGenerator Initial state: FCN = -99228.9531 Edm = 3892.616021 NCalls = 21 Info in <Minuit2>: NegativeG2LineSearch Doing a NegativeG2LineSearch since one of the G2 component is negative Info in <Minuit2>: MnSeedGenerator Negative G2 found - new state: Minimum value : -99324.40296 Edm : 4046.629084 Internal parameters: [ -0.09115714119 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.4502505244 -21446.97451 -312.347599 -899.474704 -98826.67162] Internal covariance matrix: [[ 0.0085547221 0 0 0 0] [ 0 3.6215633e-06 0 0 0] [ 0 0 0.00051136493 0 0] [ 0 0 0 0.00074754589 0] [ 0 0 0 0 1.4197205e-06]]] Info in <Minuit2>: MnSeedGenerator Initial state Minimum value : -99324.40296 Edm : 4046.629084 Internal parameters: [ -0.09115714119 -1.507540228 -0.1202898824 -0.08343008661 -1.370461484] Internal gradient : [ 0.4502505244 -21446.97451 -312.347599 -899.474704 -98826.67162] Internal covariance matrix: [[ 0.0085547221 0 0 0 0] [ 0 3.6215633e-06 0 0 0] [ 0 0 0.00051136493 0 0] [ 0 0 0 0.00074754589 0] [ 0 0 0 0 1.4197205e-06]]] Info in <Minuit2>: VariableMetricBuilder Start iterating until Edm is < 0.001 with call limit = 4500 Info in <Minuit2>: VariableMetricBuilder 0 - FCN = -99324.40296 Edm = 4046.629084 NCalls = 39 Info in <Minuit2>: VariableMetricBuilder 1 - FCN = -103242.8149 Edm = 241.9385823 NCalls = 50 Info in <Minuit2>: VariableMetricBuilder 2 - FCN = -103567.1193 Edm = 25.31359978 NCalls = 64 Info in <Minuit2>: VariableMetricBuilder 3 - FCN = -103617.8372 Edm = 0.3460395859 NCalls = 76 Info in <Minuit2>: VariableMetricBuilder 4 - FCN = -103618.5616 Edm = 0.005950694694 NCalls = 88 Info in <Minuit2>: VariableMetricBuilder 5 - FCN = -103618.5704 Edm = 0.0003289427628 NCalls = 100 Info in <Minuit2>: VariableMetricBuilder After Hessian Info in <Minuit2>: VariableMetricBuilder 6 - FCN = -103618.5704 Edm = 0.0004672814456 NCalls = 133 Info in <Minuit2>: Minuit2Minimizer::Hesse Using max-calls 4500 Info in <Minuit2>: Minuit2Minimizer::Hesse Hesse is valid - matrix is accurate Info in <TCanvas::Print>: pdf file fit-and-count/plots/sel_bin_3.pdf has been created
Excellent, our fits have run and we need now only to plot the efficiency, which we can do exactly as in the first tutorial. To check the individual fits, we can find the fit results and plots of each mass fit in fit-and-count/
.
hist = hlt_eff.get_eff("tos_efficiency_B_PT")
midpoints, values, xerrors, yerrors = tgraph_to_np(hist, xscale=1e-3)
xmin = midpoints[0] - xerrors[0][0]
xmax = midpoints[-1] + xerrors[1][-1]
# Plot the efficiency #
plt.figure(figsize=(12,10))
plt.gca()
plt.plot((xmin, xmax), (1,1), color='k', ls='dashed', lw=2)
plt.errorbar(
x=midpoints, y=values, xerr=xerrors, yerr=yerrors,
color='r', elinewidth=2, ls="none", marker='.', markersize=8,
)
plt.xlim(xmin, xmax)
plt.xlabel(r"Transverse momentum, $p_T\left(B^+\right)$ / $\mathrm{GeV}c^{-1}$")
plt.ylim(0, 1.1)
plt.ylabel(r"TOS efficiency, $\varepsilon_\mathrm{TOS}$")
hep.lhcb.label(loc=0, rlabel=r"$B^+\to J/\psi\left(\mu\mu\right)K^+$ 2024 MC")
plt.savefig("tos_efficiencies.pdf")
plt.show()