I am trying to plot the profileLL of my model with and without external constraints (gaussian functions) to visualize the impact on the precision of my variable. My model with constraints is a RooSimultaneous pdf called
combModel and the one without constraints is also a RooSimultaneous
combModel_bare consists of (for the moment) two sub-models parametrizing the (signal+background)*constraints and (signal+background) of the two categories respectively.
I am first fitting the
combModel to the combined dataset made of the two categories datasets (Asimov datasets created out of MC data), then create the NLL function and minimize it with the same specifications than the fit. Lastly, I create the profile LL of my variable out of the minimized NLL.
Then I fix the global observables appearing in the
combModel to the expectation and fit the
combModel_bare to another combined dataset made for the unconstrained models (I can also use the dataset from above, there is no difference regarding the problem) and also create the NLL function and minimize it with the same specifications as the second fit. Also here, I construct the profile LL of the variable out of the second minimized NLL.
The issue now is, that I observe different minima of the variable in the profile LL calculations than in the original fits. It seems to me, that the call of the profile likelihood changes the “best fit” value of the fitTo result. Is there something, I can do to use the minimized NLL variable as “best fit” of my POI in the profile LL?
Thank you very much!