ERROR:Eval -- RooAbsReal::logEvalError RooFit

Roofit and RooStats
1

Hello experts,

I am running some toy generation for my analysis and I am getting the following error message:

[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.841994, !sumCat=()

I have tried adjusting the definition range of the variable tauRaw based on suggestions I found here in the ROOT Forum, but the issue persists. Has anyone encountered a similar problem or has any advice on how to debug this?
Thanks in advance!

2

Welcome to the ROOT Forum!
Maybe @jonas can help, but I think we would need more info, if possible a minimal reproducer or at least a subset of your code

3

Thanks @bellenot!

The problem arises when I include in the total pdf the contribution from a particular background component (“Bs2Dsmunu”).
I can provide the “full” message that I get if this can be of any help:

--- RooAbsArg ---
  Value State: DIRTY
  Shape State: DIRTY
  Attributes:  [SnapShot_ExtRefClone] 
  Address: 0xd7de9d0
  Clients: 
    (0xd7dd2e0,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bd_GenContext "t_resmodel_Bd"
    (0xdc1fb40,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bd_GenContext_conv_exp(-@0/@1)_tauRaw_tau_[lifetime_Bd_GenContext]_GenContext "t_resmodel_Bd convoluted with basis function exp(-@0/@1)_tauRaw_tau"
    (0xdc23a20,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bd_GenContext_conv_exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_dm_[lifetime_Bd_GenContext]_GenContext "t_resmodel_Bd convoluted with basis function exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_dm"
  Servers: 
    (0xd7d65c0,V-) RooRealVar::tauRaw "#font[12]{t}"
    (0xd7df1d0,V-) RooConstVar::SplineAccCoeff0_Bd "v_{1}_Bd"
    (0xd7df7c0,V-) RooConstVar::SplineAccCoeff1_Bd "v_{2}_Bd"
    (0xd7dfe50,V-) RooConstVar::SplineAccCoeff2_Bd "v_{3}_Bd"
    (0xd7e0580,V-) RooConstVar::SplineAccCoeff3_Bd "v_{4}_Bd"
    (0xd7e0cb0,V-) RooConstVar::SplineAccCoeff4_Bd "v_{5}_Bd"
    (0xd7e13e0,V-) RooConstVar::SplineAccCoeff5_Bd "v_{6}_Bd"
    (0xd7e1b10,V-) RooConstVar::SplineAccCoeff6_Bd "v_{7}_Bd"
    (0xd7e2240,V-) RooConstVar::one "1"
    (0xdc18130,V-) RooPolyVar::SplineAccCoeff8_Bd "v_{9}_Bd"
  Proxies: 
    x -> tauRaw
    coefList -> 
      1)  SplineAccCoeff0_Bd
      2)  SplineAccCoeff1_Bd
      3)  SplineAccCoeff2_Bd
      4)  SplineAccCoeff3_Bd
      5)  SplineAccCoeff4_Bd
      6)  SplineAccCoeff5_Bd
      7)  SplineAccCoeff6_Bd
      8)                 one
      9)  SplineAccCoeff8_Bd
--- RooAbsReal ---

  Plot label is "SplineAcceptance_Bd"
Batch data access not initialised.
  1) 0xd7e7810 RooRealVar::              tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xd7ee190 RooCategory::             Mu_Sign = mixed(idx = -1)
  "mixing decision"
  3) 0xd90c2a0 RooRealVar:: k_var_Bd_GenContext = 1 C  L(0 - 1.2)  "k_var_Bd"
  1) 0xe14af70 RooRealVar::              tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xd90c2a0 RooRealVar:: k_var_Bd_GenContext = 1 C  L(0 - 1.2)  "k_var_Bd"
maxEff5 maxCode0
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.03967, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.16299, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.775014, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.7349, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.04168, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.86106, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.858796, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.924034, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.00783, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.806044, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.07258, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.842298, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.04607, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.964468, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.865494, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.834721, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.715263, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.893672, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.881596, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.03491, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.917911, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=0.846778, !sumCat=()
[#0] ERROR:Eval -- RooAbsReal::logEvalError(lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw]) evaluation error, 
 origin       : RooRealIntegral::lifetime_Bs2Dsmunu_GenContext_GenContext_Int[tauRaw]_Norm[Mu_Sign,tauRaw][ Int lifetime_Bs2Dsmunu_GenContext_GenContext_Norm(tauRaw,Mu_Sign) d[Ana](tauRaw) ]
 message      : function value is NAN
 server values: !sumList=(), !intList=(), !anaList=(tauRaw = 7.75), !jacList=(), !facList=(), !func=lifetime_Bs2Dsmunu_GenContext_GenContext=1.11731, !sumCat=()
--- RooAbsArg ---
  Value State: DIRTY
  Shape State: DIRTY
  Attributes:  [SnapShot_ExtRefClone] 
  Address: 0xdb5a1c0
  Clients: 
    (0xdb58900,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bs2Dsmunu_GenContext "t_resmodel_Bs2Dsmunu"
    (0xdb658e0,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bs2Dsmunu_GenContext_conv_exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_Bs2Dsmunu_dm_Bs2Dsmunu_[lifetime_Bs2Dsmunu_GenContext]_GenContext "t_resmodel_Bs2Dsmunu convoluted with basis function exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_Bs2Dsmunu_dm_Bs2Dsmunu"
    (0xe00c0f0,V-) RooDoubleGaussEfficiencyModel::t_resmodel_Bs2Dsmunu_GenContext_conv_exp(-@0/@1)*cosh(@0*@2/2)_tauRaw_tau_Bs2Dsmunu_dgamma_Bs2Dsmunu_[lifetime_Bs2Dsmunu_GenContext]_GenContext "t_resmodel_Bs2Dsmunu convoluted with basis function exp(-@0/@1)*cosh(@0*@2/2)_tauRaw_tau_Bs2Dsmunu_dgamma_Bs2Dsmunu"
  Servers: 
    (0xdbea340,V-) RooRealVar::tauRaw "#font[12]{t}"
    (0xdb5aa20,V-) RooConstVar::SplineAccCoeff0_Bs2Dsmunu "v_{1}_Bs2Dsmunu"
    (0xdb5b010,V-) RooConstVar::SplineAccCoeff1_Bs2Dsmunu "v_{2}_Bs2Dsmunu"
    (0xdb5b690,V-) RooConstVar::SplineAccCoeff2_Bs2Dsmunu "v_{3}_Bs2Dsmunu"
    (0xdb5bde0,V-) RooConstVar::SplineAccCoeff3_Bs2Dsmunu "v_{4}_Bs2Dsmunu"
    (0xdb5c530,V-) RooConstVar::SplineAccCoeff4_Bs2Dsmunu "v_{5}_Bs2Dsmunu"
    (0xdb5cc80,V-) RooConstVar::SplineAccCoeff5_Bs2Dsmunu "v_{6}_Bs2Dsmunu"
    (0xdb5d3d0,V-) RooConstVar::one "1"
    (0xdb5da40,V-) RooPolyVar::SplineAccCoeff6_Bs2Dsmunu "v_{7}_Bs2Dsmunu"
  Proxies: 
    x -> tauRaw
    coefList -> 
      1)  SplineAccCoeff0_Bs2Dsmunu
      2)  SplineAccCoeff1_Bs2Dsmunu
      3)  SplineAccCoeff2_Bs2Dsmunu
      4)  SplineAccCoeff3_Bs2Dsmunu
      5)  SplineAccCoeff4_Bs2Dsmunu
      6)  SplineAccCoeff5_Bs2Dsmunu
      7)                        one
      8)  SplineAccCoeff6_Bs2Dsmunu
--- RooAbsReal ---

  Plot label is "SplineAcceptance_Bs2Dsmunu"
Batch data access not initialised.
  1) 0xe01e900 RooRealVar::                     tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xd868d70 RooCategory::                    Mu_Sign = mixed(idx = -1)
  "mixing decision"
  3) 0xdb431f0 RooRealVar:: k_var_Bs2Dsmunu_GenContext = 1 C  L(0 - 1.2)  "k_var_Bs2Dsmunu"
  1) 0xda5b140 RooRealVar::                     tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xdb431f0 RooRealVar:: k_var_Bs2Dsmunu_GenContext = 1 C  L(0 - 1.2)  "k_var_Bs2Dsmunu"
maxEff5 maxCode0
--- RooAbsArg ---
  Value State: DIRTY
  Shape State: DIRTY
  Attributes:  [SnapShot_ExtRefClone] 
  Address: 0xdfef180
  Clients: 
    (0xdfeda90,V-) RooDoubleGaussEfficiencyModel::t_resmodel_SS_GenContext "t_resmodel_SS"
    (0xdffa700,V-) RooDoubleGaussEfficiencyModel::t_resmodel_SS_GenContext_conv_exp(-@0/@1)_tauRaw_tau_SS_[lifetime_SS_GenContext]_GenContext "t_resmodel_SS convoluted with basis function exp(-@0/@1)_tauRaw_tau_SS"
    (0xdffebf0,V-) RooDoubleGaussEfficiencyModel::t_resmodel_SS_GenContext_conv_exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_SS_dm_SS_[lifetime_SS_GenContext]_GenContext "t_resmodel_SS convoluted with basis function exp(-@0/@1)*cos(@0*@2)_tauRaw_tau_SS_dm_SS"
  Servers: 
    (0xdfe6ef0,V-) RooRealVar::tauRaw "#font[12]{t}"
    (0xdfef980,V-) RooConstVar::SplineAccCoeff0_SS "v_{1}_SS"
    (0xdfeff70,V-) RooConstVar::SplineAccCoeff1_SS "v_{2}_SS"
    (0xdff0560,V-) RooConstVar::SplineAccCoeff2_SS "v_{3}_SS"
    (0xdff0c90,V-) RooConstVar::SplineAccCoeff3_SS "v_{4}_SS"
    (0xdff13c0,V-) RooConstVar::SplineAccCoeff4_SS "v_{5}_SS"
    (0xdff1af0,V-) RooConstVar::SplineAccCoeff5_SS "v_{6}_SS"
    (0xdff2220,V-) RooConstVar::one "1"
    (0xdff2920,V-) RooPolyVar::SplineAccCoeff6_SS "v_{7}_SS"
  Proxies: 
    x -> tauRaw
    coefList -> 
      1)  SplineAccCoeff0_SS
      2)  SplineAccCoeff1_SS
      3)  SplineAccCoeff2_SS
      4)  SplineAccCoeff3_SS
      5)  SplineAccCoeff4_SS
      6)  SplineAccCoeff5_SS
      7)                 one
      8)  SplineAccCoeff6_SS
--- RooAbsReal ---

  Plot label is "SplineAcceptance_SS"
Batch data access not initialised.
  1) 0xd7aa490 RooRealVar::              tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xd7b0c30 RooCategory::             Mu_Sign = mixed(idx = -1)
  "mixing decision"
  3) 0xdf7b9a0 RooRealVar:: k_var_SS_GenContext = 1 C  L(0 - 1.2)  "k_var_SS"
  1) 0xda16da0 RooRealVar::              tauRaw = 7.75  L(-INF - +INF) B(292) // [ps] "#font[12]{t}"
  2) 0xdf7b9a0 RooRealVar:: k_var_SS_GenContext = 1 C  L(0 - 1.2)  "k_var_SS"
maxEff5 maxCode0
[#1] INFO:NumericIntegration -- RooRealIntegral::init(lifetime_Bu_clone_Int[Mu_Sign,tauRaw]) using numeric integrator RooIntegrator1D to calculate Int(tauRaw)
[#1] INFO:NumericIntegration -- RooRealIntegral::init(lifetime_Bu2DD_clone_Int[Mu_Sign,tauRaw]) using numeric integrator RooIntegrator1D to calculate Int(tauRaw)

This is the message I get when I start generating the model with:

for bkg in bkgs:
        bkg_comps[bkg]["signal_%s"%bkg]  = r.RooProdPdf("signal_%s"%bkg, "signal_%s"%bkg, r.RooArgList(bkg_comps[bkg]["lifetime_%s"%bkg], bkg_comps[bkg]["%sPdf"%bkg]) )


for bkg in bkgs[:len(bkgs)]:
        if(bkg=="Bs2Dsmunu"):
            bkg_comps[bkg]["frac_%s"%bkg] = r.RooRealVar("frac_%s"%bkg, "frac_%s"%bkg, 0.009*nEntries)# , 0, 0.072*nEntries )
        else:
            bkg_comps[bkg]["frac_%s"%bkg] = r.RooRealVar("frac_%s"%bkg, "frac_%s"%bkg, eval('frac_'+bkg+'_val'))
            print(f'Fraction value :: {bkg_comps[bkg]["frac_%s"%bkg].getVal()}')
       
        bkgpdf_list.add( bkg_comps[bkg]["signal_%s"%bkg] )
        bkgfrac_list.add( bkg_comps[bkg]["frac_%s"%bkg] )

model = r.RooAddPdf( "model","model", bkgpdf_list, bkgfrac_list )
dataset = model.generate( r.RooArgSet(t,Mcorr,Mu_Sign), tot_events )

Please, let me know if you need any further information. Providing a reproducible is a bit complicated with my code.

4

Your SplineAcceptance_* functions have some RooPolyVars as inputs.

The acceptance functions need to be between 0 and 1 for your full parameter phase space. Is this the case? I’m suspecting that it’s not, because a polynomial is not generally bound between [0, 1].

I’m suspecting that you get pdf values that are <=0 somewhere, because of negative acceptance. RooFit will replace negative pdf values with “NAN”, and then you get errors in the integral like you’re seeing.

What actual RooFit type are your SplineAcceptance_* functions by the way? It comes from some framework built on top of RooFit?

5

Hi, Thanks for your answer.

For that specific background component the RooFit type is the RooCubicSplineFun which is part of the B2DXFitters package.