[FEATURE] .viz.pulls for normfactors

[andrzejnovak]

.visualize.pulls Should only show the value (maybe just numeric) for unconstrained params

[alexander-held]

Thanks for raising this! Unconstrained parameters do not fit very well into this plot. I like the idea of just showing the best-fit value and uncertainty as text, alternatively we could filter out such parameters automatically. Currently this can be done manually via the exclude kwarg, but that is tedious.

Once conceptual difficulty is that the function currently does not have access to the information about which parameter is unconstrained. We could either add that to the FitResults container (might be convenient to have for other purposes as well), or hand over the model (or just the relevant piece of information) in the pulls() function call.

Another thing to consider would be staterror parameters, they do not fit particularly well here either since they have a different central value and varying constraint term widths.

[andrzejnovak]

Hah indeed, I was using exclude for staterror, but still wanted an easy way to "see" the interesting params. Indeed adding rgx/wildcards parsing to exclude would be one improvement and maybe exclude_by_type to do so by type (perhaps defaulting to ['staterror', 'normfactor']

[alexander-held]

There currently is a hardcoded exclusion for parameters called staterror_*:

cabinetry/src/cabinetry/visualize/__init__.py

Lines 482 to 483 in fce7208

	# exclude staterror parameters from pull plot (they are centered at 1)
	exclude_set.update([label for label in labels_np if label[0:10] == "staterror_"])

which is not very elegant and will fail if users call them differently (which they are free to do).

I think it would be great to do both, as you suggest:

• support exclusion by type with sensible defaults (for which we should probably save the types as metadata in FitResults),
• more flexibility in exclusion via wildcards or regex.

We're currently using fnmatch to support wildcards in other places such as

cabinetry/src/cabinetry/route.py

Line 180 in fce7208

region_matches = fnmatch.fnmatch(region_name, processor["region"])

which probably has a lower barrier to entry but is also less powerful than regex. Either could work here I think.

[vaustrup]

Hi, just stumbled across this issue. I have an implementation for the exclusion by parameter type in my fork of cabinetry. If this is still of interest, I could open a PR?

Edit: Ah, nevermind, I see it being discussed in #342 already.