@@ -34,7 +34,7 @@ Now that you've got a grid, you can perform a convolution with a PDF set:
3434We chose to use the default CT18 PDF set for this tutorial, because it's the
3535shortest to type. If you get an error that reads
3636
37- error: Invalid value for '<PDFSETS >...': The PDF set `CT18NNLO` was not found
37+ error: Invalid value for '<CONV_FUNS >...': The PDF set `CT18NNLO` was not found
3838
3939install the PDF set with LHAPDF, or use a different PDF set—the numbers won't
4040matter for the sake of the tutorial. If the command was successful, you should
@@ -122,22 +122,25 @@ following:
122122
123123 Convolutes a PineAPPL grid with a PDF set
124124
125- Usage: pineappl convolve [OPTIONS] <INPUT> <PDFSETS >...
125+ Usage: pineappl convolve [OPTIONS] <INPUT> <CONV_FUNS >...
126126
127127 Arguments:
128- <INPUT> Path of the input grid
129- <PDFSETS >... LHAPDF id(s) or name of the PDF set(s)
128+ <INPUT> Path of the input grid
129+ <CONV_FUNS >... LHAPDF id(s) or name of the PDF set(s)
130130
131131 Options:
132132 -b, --bins <BINS> Selects a subset of bins
133133 -i, --integrated Show integrated numbers (without bin widths) instead of differential ones
134134 -o, --orders <ORDERS> Select orders manually
135+ --xir <XIR> Set the variation of the renormalization scale [default: 1.0]
136+ --xif <XIF> Set the variation of the factorization scale [default: 1.0]
137+ --xia <XIA> Set the variation of the fragmentation scale [default: 1.0]
135138 --digits-abs <ABS> Set the number of fractional digits shown for absolute numbers [default: 7]
136139 --digits-rel <REL> Set the number of fractional digits shown for relative numbers [default: 2]
137140 -h, --help Print help
138141
139142This explains that ` pineappl convolve ` needs at least two arguments, the first
140- being the grid file, denoted as ` <INPUT> ` and a second argument ` <PDFSETS > ` ,
143+ being the grid file, denoted as ` <INPUT> ` and a second argument ` <CONV_FUNS > ` ,
141144which determines the PDF set. Note that this argument has three dots, ` ... ` ,
142145meaning that you're allowed to pass multiple PDF sets, in which case ` pineappl `
143146will perform the convolution with each PDF set, such that you can compare them
@@ -379,32 +382,33 @@ install the set first) as the PDF set:
379382
380383Let's calculate the scale and PDF uncertainties for our grid:
381384
382- pineappl uncert --pdf --scale-env=7 LHCB_WP_7TEV.pineappl.lz4 NNPDF31_nnlo_as_0118_luxqed
385+ pineappl uncert --conv-fun --scale-env=7 LHCB_WP_7TEV.pineappl.lz4 NNPDF31_nnlo_as_0118_luxqed
383386
384387This will show a table very similar to ` pineappl convolve ` :
385388
386- b etal dsig/detal PDF central PDF 7pt-svar (env)
387- [] [pb] [%] [%]
388- -+----+----+-----------+-----------+-----+ ----+-------+-------
389- 0 2 2.25 7.7651327e2 7.7650499e2 -1.01 1.01 -3.69 2.64
390- 1 2.25 2.5 7.1011428e2 7.1008027e2 -1.01 1.01 -3.67 2.69
391- 2 2.5 2.75 6.1683947e2 6.1679433e2 -1.01 1.01 -3.67 2.76
392- 3 2.75 3 4.9791036e2 4.9786461e2 -1.03 1.03 -3.64 2.78
393- 4 3 3.25 3.7016249e2 3.7012355e2 -1.07 1.07 -3.60 2.80
394- 5 3.25 3.5 2.5055318e2 2.5052366e2 -1.13 1.13 -3.55 2.81
395- 6 3.5 4 1.1746882e2 1.1745148e2 -1.33 1.33 -3.48 2.80
396- 7 4 4.5 2.8023753e1 2.8018010e1 -4.05 4.05 -3.40 2.74
389+ b etal dsig/detal NNPDF31_nnlo_as_0118_luxqed 7pt-svar (env)
390+ [] [pb] [pb] [%] [%] [%]
391+ -+----+----+-----------+-----------+--------+---- ----+-------+-------
392+ 0 2 2.25 7.7651327e2 7.7650499e2 -1.01 1.01 -3.69 2.64
393+ 1 2.25 2.5 7.1011428e2 7.1008027e2 -1.01 1.01 -3.67 2.69
394+ 2 2.5 2.75 6.1683947e2 6.1679433e2 -1.01 1.01 -3.67 2.76
395+ 3 2.75 3 4.9791036e2 4.9786461e2 -1.03 1.03 -3.64 2.78
396+ 4 3 3.25 3.7016249e2 3.7012355e2 -1.07 1.07 -3.60 2.80
397+ 5 3.25 3.5 2.5055318e2 2.5052366e2 -1.13 1.13 -3.55 2.81
398+ 6 3.5 4 1.1746882e2 1.1745148e2 -1.33 1.33 -3.48 2.80
399+ 7 4 4.5 2.8023753e1 2.8018010e1 -4.05 4.05 -3.40 2.74
397400
398401The first three columns are exactly the one that ` pineappl convolve ` shows. The
399- next columns are the PDF central predictions, and negative and positive PDF
400- uncertainties. These uncertainties are calculated using LHAPDF, so ` pineappl `
401- always uses the correct algorithm no matter what type of PDF sets you use:
402- Hessian, Monte Carlo, etc. Note that we've chosen a PDF set with Monte Carlo
403- replicas which means that the results shown in ` dsig/detal ` and ` PDF central `
404- aren't exactly the same. The first number is calculated from the zeroth replica
405- of the set and the second number is the average of all replicas (except the
406- zeroth). You'll notice that the PDF uncertainties are symmetric, which will not
407- necessarily be the case with Hessian PDF sets, for instance.
402+ next columns are the convolution function's central predictions, and negative
403+ and positive uncertainties. These uncertainties are calculated using LHAPDF, so
404+ ` pineappl ` always uses the correct algorithm no matter what type of convolution
405+ functions you use: Hessian, Monte Carlo, etc. Note that we've chosen a PDF set
406+ with Monte Carlo replicas which means that the results shown in ` dsig/detal ` and
407+ the central prediction aren't exactly the same. The first number is calculated
408+ from the zeroth replica of the set and the second number is the average of all
409+ replicas (except the zeroth). You'll notice that the PDF uncertainties are
410+ symmetric, which will not necessarily be the case with Hessian PDF sets, for
411+ instance.
408412
409413The two remaining columns show the scale-variation uncertainty which is
410414estimated from taking the envelope of a 7-point scale variation. By changing
0 commit comments