Upper limit on production cross section


So, I want to reproduce a plot of the upper limit on the production cross-section as a function of mass from a paper(arXiv:2101.11582v2).

With the given information about Data, Total background, uncertainty in the background, expected signal, and uncertainty in expected signal yields separately in 7 signal region categories. I want to find out the upper limit on cross-section. I have this information for three mass values. First, I want to check at least for these mass points whether I am getting the right result by comparing them with the plot from the paper.

So, how do I do this?

I was looking at some of the exercise. Can I do using this or is there any other way.

However, I think this tutorial has no option to use expected signal events and corresponding uncertainty.

Also, they consider the background a nuisance parameter, but in my case, I already have the background events from the paper. Also, they consider it a simple counting experiment, but how can I do the same when considering seven signal regions?

Is there something where I can do this using RooFit/RooStats? Can this also be done using HistFitter or pyhf? I was also trying out those things, but my results never matched.

Or is there something wrong with my interpretation of the results? Because, I think the upper limit I am getting is on the number of signal events. Then, how do I correlate this to the cross-section? Since my branching ratio is one, I am multiplying the cross section with this upper limit on no. of signal events. Is this wrong? Or should I divide it with (LuminosityAcceptanceEfficiency)? However they have given the information about acceptance times efficiency for all signal regions combined. Then, can I use it here?
Any insights would be of great help.

Thanking you
Debabrata Sahoo

Welcome to the ROOT Forum!
I’m not sure that what you describe is sufficient to provide a concrete answer, but I invite @jonas to have a look nevertheless

I want to reproduce this plot (expected limit )

Using the information:

I have information for expected Signal events for three mass points (0.9TeV, 1.1 and 1.3 TeV). I should be able to find the upper limit on cross-section, at least for these three masses. My question is, how can I find the desired upper limit, considering all seven signal regions combined?

I hope my question is clear.

Dear @jonas, Could you please say something about this?