search

CMS-HSCP / SUSYBSMAnalysis-HSCP

HSCP analysis code
3 stars 16 forks source link

Recalculate the interpolated stau cross-sections #110

Closed tvami closed 1 year ago

tvami commented 1 year ago

The stau numbers from https://twiki.cern.ch/twiki/bin/view/LHCPhysics/SUSYCrossSections i.e. https://github.com/fuenfundachtzig/xsec/blob/master/json/pp13_stau_LR_NLO%2BNLL_PDF4LHC.json come in 20 GeV blocks while our samples are simulated with the lambda variation which leads to non-20GeV steps. We need to interpolate the curve from the JSON above to the exact numbers

lawrenceleejr commented 1 year ago

Implementing spline-based interpolation of the table numbers. The relevant arrays are

[308.0, 432.0, 557.0, 651.0, 745.0, 871.0, 1029.0]
[0.005617969677491303, 0.0013205993021946, 0.0003982919391389629, 0.00018170273864904894, 8.760649473515753e-05, 3.548748280460506e-05, 1.1684554271190975e-05]

Working on PR to provide interpolation script, and to update the numbers in the limit plot scripts. Plot below demonstrates the interpolation result. The x axis is mass, y axis is cross section. The black points are from the twiki and the red points are the interpolated results above.

test.pdf

lawrenceleejr commented 1 year ago

See here https://github.com/tvami/HSCPbackgroundPred/pull/3 for the PR. I suppose once these new numbers are propagated to the analyzer and that PR merged, this issue can be closed.