Is it possible to Scan a parameter space using PYHF?

[MohammadMohammadidoust]

Dear All
I am new to pyhf and also new to experimental HEP. I read the pyhf tutorial and have a question. As you mentioned in the section "Using Hep data", to investigate a parameter point, we should create a model using a patch from the given patchset and calculate CLs for this point, and then based on the CLs value we decide to exclude this point or not.
But my question is how to find out about the whole parameter space. the number of parameter points in patchsets is limited and it seems that we should use interpolation methods to find whether each arbitrary point of a parameter space should be excluded or not. could I do this using pyhf?
In other words, I need something like the one you have plotted in the tutorial(ATLAS results plot in the section "using HEPdata" of the tutorial) to flag each point in the parameter space regarding its CLs
can I achieve this using pyhf?

[kratsg]

Hi, this is somewhat out of scope of pyhf. You are right in that you need to perform interpolation (what kind is really up to you; e.g. linear, quartic, ...). You can see what is done in HistFitter using scipy interpolation which relies on making a harvest.json. Once you have this TGraph however, you can make a plot like what's shown in your picture using this gist I made depending on matplotlib and mplhep.