Dear roostate experts,
I’m trying to make a plot of CL scan wrt poi. I start working with the StandardHypoTestInvDemo.C which is in the roostate tutorial. However, my fit could not converge at any point of the poi while it does give a converged fit result of upper limit when I set the CL to 95%. I’m thinking if there’s any fitting config that I can tune with to make the fit successful? One thing to notice is that there’s more than 500 NPs in my workspace, mostly are gamma NPs. Would the fit fails because of too many NPs?
Attached at the end of this email you can find my input parameters for the StandardHypoTestInvDemo function and the error message. My workspace and the full log file can be found at
/afs/cern.ch/user/w/wangwe/public/CLscan/
Thank you for your help!
Regards,
Wei
Input parameters:
const char * wsName = "combined",
const char * modelSBName = "ModelConfig",
const char * modelBName = "",
const char * dataName = "obsData",
int calculatorType = 2,
int testStatType = 3,
bool useCLs = true ,
int npoints = 5,
double poimin = 0,
double poimax = 1.2,
int ntoys=1000,
bool useNumberCounting = false,
const char * nuisPriorName = 0
Error message (point 0.3 for example):
[#0] PROGRESS:Eval – Running for SigXsecOverSM = 0.3
----> Doing a re-scan first
----> trying with strategy = 1
----> trying with improve
[#0] ERROR:Minization – AsymptoticCalculator: failure in fitting for qmu or qmuA 0.3 return a dummy result
[#0] PROGRESS:Eval – P values for SigXsecOverSM = 0.3
CLs = nan +/- 0
CLb = nan +/- 0
CLsplusb = nan +/- 0
[#0] WARNING:Eval – HypoTestInverter - Skip invalid result for point SigXsecOverSM = 0.3
[#1] INFO:ObjectHandling – RooWorkspace::saveSnaphot(combined) replacing previous snapshot with name ModelConfig__snapshot