I have a 2D binned fit with different components. Now I want to test the fit using toy data generated directly from the PDFs.
I know that we can generate the toy data using PDF>generate(), but here I have ~5 PDfs of both variables (x and y).
By the way, don’t forget to use the Extended() command argument if you also want to randomize the number of events in the toy dataset according to a Poisson distribution around the nominal sum of yields in the model.
Thanks for the reply. yes, I did this, and working fine.
Now I want to check the distribution of (fitted yield - input yield)/fit uncertainty in each toy dataset.
As far as I know in RooMCStudy, we can do this using mcstudy->plotPull().
Is there any way to do this here?