read in partonic xs

This commit is contained in:
hiro98 2020-05-03 12:26:52 +02:00
parent 61ef7052c0
commit 97083d676e
3 changed files with 11 additions and 1 deletions

View file

@ -377,6 +377,14 @@ We have to convert that to picobarn.
| 2.1137030945166285e-05 | 3.175822429495981e-08 |
That is compatible with sherpa!
#+begin_src jupyter-python :exports both :results raw drawer
sherpa_part, sherpa_part_σ = np.loadtxt('../../runcards/pp_partonic/sherpa_xs')
sherpa_part, sherpa_part_σ # GeV
#+end_src
#+RESULTS:
| 2.11499e-05 | 6.8969e-09 |
We can take some samples as well.
#+begin_src jupyter-python :exports both :results raw drawer

View file

@ -1,5 +1,5 @@
sherpa_xs: Sherpa.yaml
Sherpa -e 0 -a "" 'GENERATE_RESULT_DIRECTORY: false' | sed 's/\x1b\[[0-9;]*m//g' | sed -n 's/.*\s:\s\([0-9]*\.[0-9]*\)\spb\s+-\s(\s\([0-9]*\.[0-9]*e-*[0-9]*\).*/\\\sigma = \\\SI{\1\\\pm \2}{\\\pico\\\barn}/p' > sherpa_xs
Sherpa -e 0 -a "" 'GENERATE_RESULT_DIRECTORY: false' | sed 's/\x1b\[[0-9;]*m//g' | perl -ne 'while(/.*\s:\s([0-9]\.[0-9]+e-?[0-9]+)\spb\s\+-\s\(\s([0-9]\.[0-9]+e-?[0-9]+).*/g){print "$$1\n$$2";}' > sherpa_xs
ROOT_DIR:=$(shell dirname $(realpath $(firstword $(MAKEFILE_LIST))))
analysis: Sherpa.yaml

View file

@ -0,0 +1,2 @@
2.11499e-05
6.8969e-09