Hi all, I was trying to use the RooMultiPdf and CMS combine tool for bias studies where it uses RooMinimizer but I always get this error… I checked my fits using different polynomials to data but it all turned out okay (i.e, STATUS=OK). I really don’t know why it throws me this error… Any help is appreciated. Thanks!
WARNING: underflow to nan in pdf_binCat2_obsOnly for bin 0, weight 5 FASTEXIT from pdf_binCat2_obsOnly [#0] WARNING:Minization -- RooMinimizerFcn: Minimized function has error status. Returning maximum FCN so far (5.26398e-05) to force MIGRAD to back out of this region. Error log follows Parameter values: Bern1_p0=nan, Bern2_p0=nan, Bern2_p1=nan, CMS_Ele_ID=nan, CMS_PU_13TeV=nan, CMS_hzg_Pho_CSEV_SF=nan, CMS_hzg_Pho_ID_SF=nan, CMS_hzg_res_EE_Cat2=nan, CMS_hzg_scale_EE_Cat2=nan, CMS_hzg_trig_SF=nan, CMS_mH_BtagSF=nan, Laurent1_l1=nan, Prod_XS_PDF_Alpha=nan, exp1_p1=nan, lumi_13TeV=nan, powerlaw1_p1=nan, r=-150, shapeBkg_bkg_mass_Cat2__norm=nan RooProduct::n_exp_final_binCat2_proc_bkg_mass[ n_exp_binCat2_proc_bkg_mass * shapeBkg_bkg_mass_Cat2__norm ] function value is NAN @ !compRSet=(n_exp_binCat2_proc_bkg_mass = 1,shapeBkg_bkg_mass_Cat2__norm = nan +/- 12.0361), !compCSet=() RooProduct::newmean_EE_Cat2[ mean_EE_Cat2 * CMS_hzg_scale_EE_Cat2 ] function value is NAN @ !compRSet=(mean_EE_Cat2 = 125.826 +/- 16.4588,CMS_hzg_scale_EE_Cat2 = nan +/- 0.00618534), !compCSet=() ProcessNormalization::n_exp_binCat2_proc_hzg[ thetaList=(lumi_13TeV,Prod_XS_PDF_Alpha,CMS_PU_13TeV,CMS_mH_BtagSF,CMS_Ele_ID,CMS_hzg_trig_SF,CMS_hzg_Pho_ID_SF,CMS_hzg_Pho_CSEV_SF) asymmThetaList=() otherFactorList=(r) ] function value is NAN @ thetaList=(lumi_13TeV = nan +/- 0.994458,Prod_XS_PDF_Alpha = nan +/- 0.994458,CMS_PU_13TeV = nan +/- 0.994458,CMS_mH_BtagSF = nan +/- 0.994458,CMS_Ele_ID = nan +/- 0.994458,CMS_hzg_trig_SF = nan +/- 0.994458,CMS_hzg_Pho_ID_SF = nan +/- 0.994458,CMS_hzg_Pho_CSEV_SF = nan +/- 0.994458), asymmThetaList=(), otherFactorList=(r = -150 +/- 68.2483) RooProduct::newsigma_EE_Cat2[ sigma_EE_Cat2 * CMS_hzg_res_EE_Cat2 ] function value is NAN @ !compRSet=(sigma_EE_Cat2 = 11.2476 +/- 12.1663,CMS_hzg_res_EE_Cat2 = nan +/- 0.0791254), !compCSet=() RooFormulaVar::Bern1_p0_sq[ actualVars=(Bern1_p0) formula="@0*@0" ] function value is NAN @ actualVars=(Bern1_p0 = nan +/- 41.8899) RooAbsReal::[ params=(lumi_13TeV,lumi_13TeV_In,Prod_XS_PDF_Alpha,Prod_XS_PDF_Alpha_In,CMS_PU_13TeV,CMS_PU_13TeV_In,CMS_mH_BtagSF,CMS_mH_BtagSF_In,CMS_Ele_ID,CMS_Ele_ID_In,CMS_hzg_trig_SF,CMS_hzg_trig_SF_In,CMS_hzg_Pho_ID_SF,CMS_hzg_Pho_ID_SF_In,CMS_hzg_Pho_CSEV_SF,CMS_hzg_Pho_CSEV_SF_In,CMS_hzg_scale_EE_Cat2,CMS_hzg_scale_EE_Cat2_In,CMS_hzg_res_EE_Cat2,CMS_hzg_res_EE_Cat2_In,Bern1_p0,Bern2_p0,Bern2_p1,Laurent1_l1,alpha2_EE_Cat2,alpha_EE_Cat2,exp1_p1,mean_EE_Cat2,n2_EE_Cat2,n_EE_Cat2,n_exp_binCat2_proc_bkg_mass,powerlaw1_p1,r,shapeBkg_bkg_mass_Cat2__norm,sigma_EE_Cat2) ] function value is NAN @ params=(lumi_13TeV = nan +/- 0.994458,lumi_13TeV_In = 0,Prod_XS_PDF_Alpha = nan +/- 0.994458,Prod_XS_PDF_Alpha_In = 0,CMS_PU_13TeV = nan +/- 0.994458,CMS_PU_13TeV_In = 0,CMS_mH_BtagSF = nan +/- 0.994458,CMS_mH_BtagSF_In = 0,CMS_Ele_ID = nan +/- 0.994458,CMS_Ele_ID_In = 0,CMS_hzg_trig_SF = nan +/- 0.994458,CMS_hzg_trig_SF_In = 0,CMS_hzg_Pho_ID_SF = nan +/- 0.994458,CMS_hzg_Pho_ID_SF_In = 0,CMS_hzg_Pho_CSEV_SF = nan +/- 0.994458,CMS_hzg_Pho_CSEV_SF_In = 0,CMS_hzg_scale_EE_Cat2 = nan +/- 0.00618534,CMS_hzg_scale_EE_Cat2_In = 1,CMS_hzg_res_EE_Cat2 = nan +/- 0.0791254,CMS_hzg_res_EE_Cat2_In = 1,Bern1_p0 = nan +/- 41.8899,Bern2_p0 = nan +/- 41.9154,Bern2_p1 = nan +/- 41.9187,Laurent1_l1 = nan +/- 7.65461,alpha2_EE_Cat2 = 0.894593 +/- 2.2633,alpha_EE_Cat2 = 0.463823 +/- 2.10824,exp1_p1 = nan +/- 2.27154,mean_EE_Cat2 = 125.826 +/- 16.4588,n2_EE_Cat2 = 0.543447 +/- 1.18292,n_EE_Cat2 = 0.439995 +/- 1.04775,n_exp_binCat2_proc_bkg_mass = 1,powerlaw1_p1 = nan +/- 19.5489,r = -150 +/- 68.2483,shapeBkg_bkg_mass_Cat2__norm = nan +/- 12.0361,sigma_EE_Cat2 = 11.2476 +/- 12.1663) RooAbsReal::Cat2[ params=(Bern1_p0,Bern2_p0,Bern2_p1,CMS_Ele_ID,CMS_PU_13TeV,CMS_hzg_Pho_CSEV_SF,CMS_hzg_Pho_ID_SF,CMS_hzg_res_EE_Cat2,CMS_hzg_scale_EE_Cat2,CMS_hzg_trig_SF,CMS_mH_BtagSF,Laurent1_l1,Prod_XS_PDF_Alpha,alpha2_EE_Cat2,alpha_EE_Cat2,exp1_p1,lumi_13TeV,mean_EE_Cat2,n2_EE_Cat2,n_EE_Cat2,n_exp_binCat2_proc_bkg_mass,powerlaw1_p1,r,shapeBkg_bkg_mass_Cat2__norm,sigma_EE_Cat2) ] Number of events is negative or error @ params=(Bern1_p0 = nan +/- 41.8899,Bern2_p0 = nan +/- 41.9154,Bern2_p1 = nan +/- 41.9187,CMS_Ele_ID = nan +/- 0.994458,CMS_PU_13TeV = nan +/- 0.994458,CMS_hzg_Pho_CSEV_SF = nan +/- 0.994458,CMS_hzg_Pho_ID_SF = nan +/- 0.994458,CMS_hzg_res_EE_Cat2 = nan +/- 0.0791254,CMS_hzg_scale_EE_Cat2 = nan +/- 0.00618534,CMS_hzg_trig_SF = nan +/- 0.994458,CMS_mH_BtagSF = nan +/- 0.994458,Laurent1_l1 = nan +/- 7.65461,Prod_XS_PDF_Alpha = nan +/- 0.994458,alpha2_EE_Cat2 = 0.894593 +/- 2.2633,alpha_EE_Cat2 = 0.463823 +/- 2.10824,exp1_p1 = nan +/- 2.27154,lumi_13TeV = nan +/- 0.994458,mean_EE_Cat2 = 125.826 +/- 16.4588,n2_EE_Cat2 = 0.543447 +/- 1.18292,n_EE_Cat2 = 0.439995 +/- 1.04775,n_exp_binCat2_proc_bkg_mass = 1,powerlaw1_p1 = nan +/- 19.5489,r = -150 +/- 68.2483,shapeBkg_bkg_mass_Cat2__norm = nan +/- 12.0361,sigma_EE_Cat2 = 11.2476 +/- 12.1663) RooBernstein::Bern1[ x=CMS_hzg_mass coefList=(Bern1_p0_sq) ] p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=111, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=113, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=115, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=117, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=119, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=121, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=123, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=125, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=127, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=129, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=131, coefList=(Bern1_p0_sq = nan) p.d.f value is Not-a-Number (nan), forcing value to zero @ x=CMS_hzg_mass=133, coefList=(Bern1_p0_sq = nan) ... (remaining 25 messages suppressed) RooRealIntegral::shapeSig_hzg_Cat2_Int[CMS_hzg_mass][ Int shapeSig_hzg_Cat2d[Ana](CMS_hzg_mass) ] function value is NAN @ !sumList=(), !intList=(), !anaList=(CMS_hzg_mass = 111), !jacList=(), !facList=(), !func=shapeSig_hzg_Cat2=99/nan, !sumCat=() RooRealIntegral::Bern1_Int[CMS_hzg_mass][ Int Bern1d[Ana](CMS_hzg_mass) ] function value is NAN @ !sumList=(), !intList=(), !anaList=(CMS_hzg_mass = 111), !jacList=(), !facList=(), !func=Bern1=nan/nan, !sumCat=() ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range... ERROR evaluating range...