cschwan/compare_old_vs_new_ckm.sh

## compare_old_vs_new_ckm.sh
#!/bin/bash

set -euo pipefail

grids=(
    ATLAS_WM_7TEV.pineappl.lz4
    ATLAS_WP_7TEV.pineappl.lz4
    LHCB_WM_7TEV.pineappl.lz4
    LHCB_WM_8TEV.pineappl.lz4
    LHCB_WP_7TEV.pineappl.lz4
    LHCB_WP_8TEV.pineappl.lz4
)

old_grids=(
    ATLASWZRAP11CC/ATLASWZRAP11-WM_ylept.pineappl.lz4
    ATLASWZRAP11CC/ATLASWZRAP11-WP_ylept.pineappl.lz4
    LHCBWZMU7TEV/LHCBWZMU7TEV_WM_leptrap.pineappl.lz4
    LHCBWZMU8TEV/LHCBWZMU8TEV_WM_leptrap.pineappl.lz4
    LHCBWZMU7TEV/LHCBWZMU7TEV_WP_leptrap.pineappl.lz4
    LHCBWZMU8TEV/LHCBWZMU8TEV_WP_leptrap.pineappl.lz4
)

test_pdf_set=NNPDF31_nlo_as_0118

for (( i = 0; i != ${#grids[@]}; ++i )); do
    old="${old_grids[$i]}"
    new="${grids[$i]}"
    new_ckm="${new%.pineappl.lz4}_CKM.pineappl.lz4"

    echo ">>> ${old} vs. ${new}"
    echo ">>> CKM (OLD) vs. NO CKM NEW"
    pineappl --silence-lhapdf diff ../applgrids/"${old}" "${new}" "${test_pdf_set}" \
        --orders1=a2,as1a2 --orders2=a2,as1a2 --scale2=1000 --ignore-lumis --ignore-bin-limits
    echo ">>> CKM (OLD) vs. CKM NEW"
    pineappl --silence-lhapdf diff ../applgrids/"${old}" "${new_ckm}" "${test_pdf_set}" \
        --orders1=a2,as1a2 --orders2=a2,as1a2 --scale2=1000 --ignore-lumis --ignore-bin-limits

    echo
    echo
done

## comparison_ckm_no_ckm.txt
>>> ATLASWZRAP11CC/ATLASWZRAP11-WM_ylept.pineappl.lz4 vs. ATLAS_WM_7TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b     x1                O(as^0 a^2)                       O(as^1 a^2)
--+----+----+-----------+-----------+---------+-----------+-----------+---------
 0    0 0.21 4.2240403e5 4.2418021e5 -4.187e-3 7.4550308e3 7.5751404e3 -1.586e-2
 1 0.21 0.42 4.2003285e5 4.2202888e5 -4.730e-3 7.6190059e3 7.9119070e3 -3.702e-2
 2 0.42 0.63 4.1563075e5 4.1764429e5 -4.821e-3 8.1040021e3 8.6391043e3 -6.194e-2
 3 0.63 0.84 4.0928859e5 4.1124265e5 -4.752e-3 8.4422888e3 8.9988383e3 -6.185e-2
 4 0.84 1.05 4.0079290e5 4.0284389e5 -5.091e-3 9.0842830e3 9.4793358e3 -4.168e-2
 5 1.05 1.37 3.8741291e5 3.8958316e5 -5.571e-3 9.7103988e3 1.0151402e4 -4.344e-2
 6 1.37 1.52 3.7401340e5 3.7623615e5 -5.908e-3 1.0039762e4 1.0131505e4 -9.055e-3
 7 1.52 1.74 3.6219912e5 3.6462541e5 -6.654e-3 9.9360130e3 1.0403881e4 -4.497e-2
 8 1.74 1.95 3.4804228e5 3.5053483e5 -7.111e-3 9.5605621e3 9.8270974e3 -2.712e-2
 9 1.95 2.18 3.3289241e5 3.3517621e5 -6.814e-3 8.6304707e3 9.1482968e3 -5.660e-2
10 2.18  2.5 3.1245793e5 3.1473844e5 -7.246e-3 6.8339401e3 7.1154335e3 -3.956e-2
>>> CKM (OLD) vs. CKM NEW
b     x1               O(as^0 a^2)                       O(as^1 a^2)
--+----+----+-----------+-----------+--------+-----------+-----------+---------
 0    0 0.21 4.2240403e5 4.2185924e5 1.291e-3 7.4550308e3 7.2050963e3  3.469e-2
 1 0.21 0.42 4.2003285e5 4.1969088e5 8.148e-4 7.6190059e3 7.5383742e3  1.070e-2
 2 0.42 0.63 4.1563075e5 4.1527346e5 8.604e-4 8.1040021e3 8.2594206e3 -1.882e-2
 3 0.63 0.84 4.0928859e5 4.0882421e5 1.136e-3 8.4422888e3 8.6120819e3 -1.972e-2
 4 0.84 1.05 4.0079290e5 4.0036730e5 1.063e-3 9.0842830e3 9.0828089e3  1.623e-4
 5 1.05 1.37 3.8741291e5 3.8702648e5 9.985e-4 9.7103988e3 9.7388763e3 -2.924e-3
 6 1.37 1.52 3.7401340e5 3.7361073e5 1.078e-3 1.0039762e4 9.7056095e3  3.443e-2
 7 1.52 1.74 3.6219912e5 3.6195622e5 6.711e-4 9.9360130e3 9.9679405e3 -3.203e-3
 8 1.74 1.95 3.4804228e5 3.4783128e5 6.066e-4 9.5605621e3 9.3823231e3  1.900e-2
 9 1.95 2.18 3.3289241e5 3.3246832e5 1.276e-3 8.6304707e3 8.6971424e3 -7.666e-3
10 2.18  2.5 3.1245793e5 3.1207931e5 1.213e-3 6.8339401e3 6.6674362e3  2.497e-2


>>> ATLASWZRAP11CC/ATLASWZRAP11-WP_ylept.pineappl.lz4 vs. ATLAS_WP_7TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b     x1                O(as^0 a^2)                        O(as^1 a^2)
--+----+----+-----------+-----------+---------+------------+------------+---------
 0    0 0.21 5.7135378e5 5.7331872e5 -3.427e-3 -4.7525156e3 -4.4496358e3  6.807e-2
 1 0.21 0.42 5.7284708e5 5.7495947e5 -3.674e-3 -4.6408670e3 -4.3703932e3  6.189e-2
 2 0.42 0.63 5.7554050e5 5.7735379e5 -3.141e-3 -4.3420155e3 -4.6433339e3 -6.489e-2
 3 0.63 0.84 5.7889298e5 5.8143980e5 -4.380e-3 -4.0250183e3 -3.7056226e3  8.619e-2
 4 0.84 1.05 5.8341570e5 5.8516948e5 -2.997e-3 -3.1930355e3 -3.2423913e3 -1.522e-2
 5 1.05 1.37 5.8829101e5 5.9012841e5 -3.114e-3 -1.9993090e3 -1.7503055e3  1.423e-1
 6 1.37 1.52 5.9204449e5 5.9433809e5 -3.859e-3 -1.6054598e2  1.4259565e3  -1.113e0
 7 1.52 1.74 5.9249892e5 5.9433031e5 -3.081e-3  2.6661527e3  2.4117860e3  1.055e-1
 8 1.74 1.95 5.8785760e5 5.8958097e5 -2.923e-3  6.7750134e3  7.5921717e3 -1.076e-1
 9 1.95 2.18 5.7427968e5 5.7609457e5 -3.150e-3  1.2635183e4  1.1698721e4  8.005e-2
10 2.18  2.5 5.3635711e5 5.3784689e5 -2.770e-3  2.1430029e4  2.2286612e4 -3.843e-2
>>> CKM (OLD) vs. CKM NEW
b     x1               O(as^0 a^2)                        O(as^1 a^2)
--+----+----+-----------+-----------+--------+------------+------------+---------
 0    0 0.21 5.7135378e5 5.7061922e5 1.287e-3 -4.7525156e3 -4.7508919e3  3.418e-4
 1 0.21 0.42 5.7284708e5 5.7225623e5 1.032e-3 -4.6408670e3 -4.6741803e3 -7.127e-3
 2 0.42 0.63 5.7554050e5 5.7464748e5 1.554e-3 -4.3420155e3 -4.9507204e3 -1.230e-1
 3 0.63 0.84 5.7889298e5 5.7872675e5 2.872e-4 -4.0250183e3 -4.0226965e3  5.772e-4
 4 0.84 1.05 5.8341570e5 5.8245769e5 1.645e-3 -3.1930355e3 -3.5697458e3 -1.055e-1
 5 1.05 1.37 5.8829101e5 5.8743106e5 1.464e-3 -1.9993090e3 -2.0939761e3 -4.521e-2
 6 1.37 1.52 5.9204449e5 5.9166973e5 6.334e-4 -1.6054598e2  1.0621365e3  -1.151e0
 7 1.52 1.74 5.9249892e5 5.9171330e5 1.328e-3  2.6661527e3  2.0395420e3  3.072e-1
 8 1.74 1.95 5.8785760e5 5.8705863e5 1.361e-3  6.7750134e3  7.2035828e3 -5.949e-2
 9 1.95 2.18 5.7427968e5 5.7371818e5 9.787e-4  1.2635183e4  1.1305429e4  1.176e-1
10 2.18  2.5 5.3635711e5 5.3574500e5 1.143e-3  2.1430029e4  2.1894617e4 -2.122e-2


>>> LHCBWZMU7TEV/LHCBWZMU7TEV_WM_leptrap.pineappl.lz4 vs. LHCB_WM_7TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 3.7292336e5 3.7532595e5 -6.401e-3 6.4700672e4 6.5332710e4 -9.674e-3
1 2.25  2.5 3.5347999e5 3.5608848e5 -7.325e-3 5.8824761e4 5.9946156e4 -1.871e-2
2  2.5 2.75 3.3090524e5 3.3406844e5 -9.469e-3 5.3846944e4 5.4627884e4 -1.430e-2
3 2.75    3 3.0423388e5 3.0718892e5 -9.620e-3 4.6850349e4 4.7985867e4 -2.366e-2
4    3 3.25 2.7190663e5 2.7375224e5 -6.742e-3 4.0099659e4 4.0713216e4 -1.507e-2
5 3.25  3.5 2.3113503e5 2.3297956e5 -7.917e-3 3.2796720e4 3.2999245e4 -6.137e-3
6  3.5    4 1.6028275e5 1.6169805e5 -8.753e-3 2.1119242e4 2.1665733e4 -2.522e-2
7    4  4.5 7.1251353e4 7.1892674e4 -8.921e-3 9.2469580e3 9.4199433e3 -1.836e-2
>>> CKM (OLD) vs. CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 3.7292336e5 3.7228082e5  1.726e-3 6.4700672e4 6.4556375e4  2.235e-3
1 2.25  2.5 3.5347999e5 3.5307300e5  1.153e-3 5.8824761e4 5.9169740e4 -5.830e-3
2  2.5 2.75 3.3090524e5 3.3115068e5 -7.412e-4 5.3846944e4 5.3859539e4 -2.339e-4
3 2.75    3 3.0423388e5 3.0445236e5 -7.176e-4 4.6850349e4 4.7246805e4 -8.391e-3
4    3 3.25 2.7190663e5 2.7128834e5  2.279e-3 4.0099659e4 4.0027055e4  1.814e-3
5 3.25  3.5 2.3113503e5 2.3087760e5  1.115e-3 3.2796720e4 3.2391407e4  1.251e-2
6  3.5    4 1.6028275e5 1.6025570e5  1.688e-4 2.1119242e4 2.1218833e4 -4.694e-3
7    4  4.5 7.1251353e4 7.1284320e4 -4.625e-4 9.2469580e3 9.2071618e3  4.322e-3


>>> LHCBWZMU8TEV/LHCBWZMU8TEV_WM_leptrap.pineappl.lz4 vs. LHCB_WM_8TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 4.2643373e5 4.2995543e5 -8.191e-3 7.3094364e4 7.3808255e4 -9.672e-3
1 2.25  2.5 4.0518638e5 4.0840319e5 -7.877e-3 6.7780367e4 6.8793814e4 -1.473e-2
2  2.5 2.75 3.8177585e5 3.8495435e5 -8.257e-3 6.1521771e4 6.2615569e4 -1.747e-2
3 2.75    3 3.5475650e5 3.5742431e5 -7.464e-3 5.5254814e4 5.5702896e4 -8.044e-3
4    3 3.25 3.2045474e5 3.2372236e5 -1.009e-2 4.7085807e4 4.7987677e4 -1.879e-2
5 3.25  3.5 2.7989006e5 2.8234548e5 -8.697e-3 3.8787358e4 3.9715001e4 -2.336e-2
6  3.5    4 2.0482249e5 2.0668879e5 -9.030e-3 2.6447254e4 2.7031239e4 -2.160e-2
7    4  4.5 1.0172352e5 1.0275727e5 -1.006e-2 1.2542814e4 1.2754997e4 -1.664e-2
>>> CKM (OLD) vs. CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 4.2643373e5 4.2679976e5 -8.576e-4 7.3094364e4 7.3014109e4  1.099e-3
1 2.25  2.5 4.0518638e5 4.0524546e5 -1.458e-4 6.7780367e4 6.7986601e4 -3.033e-3
2  2.5 2.75 3.8177585e5 3.8185789e5 -2.149e-4 6.1521771e4 6.1812200e4 -4.699e-3
3 2.75    3 3.5475650e5 3.5446883e5  8.115e-4 5.5254814e4 5.4918844e4  6.118e-3
4    3 3.25 3.2045474e5 3.2099971e5 -1.698e-3 4.7085807e4 4.7244464e4 -3.358e-3
5 3.25  3.5 2.7989006e5 2.7995161e5 -2.199e-4 3.8787358e4 3.9036955e4 -6.394e-3
6  3.5    4 2.0482249e5 2.0494152e5 -5.808e-4 2.6447254e4 2.6504573e4 -2.163e-3
7    4  4.5 1.0172352e5 1.0191846e5 -1.913e-3 1.2542814e4 1.2472070e4  5.672e-3


>>> LHCBWZMU7TEV/LHCBWZMU7TEV_WP_leptrap.pineappl.lz4 vs. LHCB_WP_7TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 6.5227719e5 6.5524066e5 -4.523e-3 1.0987679e5 1.1243687e5 -2.277e-2
1 2.25  2.5 5.9761510e5 6.0023547e5 -4.366e-3 1.0237238e5 1.0140698e5  9.520e-3
2  2.5 2.75 5.1669345e5 5.1937942e5 -5.171e-3 8.9601247e4 9.0379386e4 -8.610e-3
3 2.75    3 4.1745539e5 4.1841432e5 -2.292e-3 7.3703721e4 7.4337163e4 -8.521e-3
4    3 3.25 3.0906194e5 3.1041740e5 -4.367e-3 5.6379044e4 5.6700595e4 -5.671e-3
5 3.25  3.5 2.0870669e5 2.0961121e5 -4.315e-3 3.9340453e4 3.9651322e4 -7.840e-3
6  3.5    4 9.7024044e4 9.7472812e4 -4.604e-3 1.9818646e4 1.9888208e4 -3.498e-3
7    4  4.5 2.2300258e4 2.2380799e4 -3.599e-3 5.3028838e3 5.3460922e3 -8.082e-3
>>> CKM (OLD) vs. CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 6.5227719e5 6.5259204e5 -4.825e-4 1.0987679e5 1.1173478e5 -1.663e-2
1 2.25  2.5 5.9761510e5 5.9789562e5 -4.692e-4 1.0237238e5 1.0077814e5  1.582e-2
2  2.5 2.75 5.1669345e5 5.1742466e5 -1.413e-3 8.9601247e4 8.9823591e4 -2.475e-3
3 2.75    3 4.1745539e5 4.1688964e5  1.357e-3 7.3703721e4 7.3885035e4 -2.454e-3
4    3 3.25 3.0906194e5 3.0931927e5 -8.319e-4 5.6379044e4 5.6359854e4  3.405e-4
5 3.25  3.5 2.0870669e5 2.0888939e5 -8.746e-4 3.9340453e4 3.9416548e4 -1.931e-3
6  3.5    4 9.7024044e4 9.7143797e4 -1.233e-3 1.9818646e4 1.9774752e4  2.220e-3
7    4  4.5 2.2300258e4 2.2295006e4  2.356e-4 5.3028838e3 5.3159836e3 -2.464e-3


>>> LHCBWZMU8TEV/LHCBWZMU8TEV_WP_leptrap.pineappl.lz4 vs. LHCB_WP_8TEV.pineappl.lz4
>>> CKM (OLD) vs. NO CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 7.4136186e5 7.4379370e5 -3.270e-3 1.2274280e5 1.2291552e5 -1.405e-3
1 2.25  2.5 6.9378357e5 6.9718568e5 -4.880e-3 1.1516651e5 1.1609118e5 -7.965e-3
2  2.5 2.75 6.1967049e5 6.2233439e5 -4.280e-3 1.0422790e5 1.0479152e5 -5.378e-3
3 2.75    3 5.1925652e5 5.2114837e5 -3.630e-3 8.8554912e4 8.9146035e4 -6.631e-3
4    3 3.25 4.0296557e5 4.0414296e5 -2.913e-3 7.0854643e4 7.0828726e4  3.659e-4
5 3.25  3.5 2.8581701e5 2.8723414e5 -4.934e-3 5.1623949e4 5.1813883e4 -3.666e-3
6  3.5    4 1.4463541e5 1.4513268e5 -3.426e-3 2.7691674e4 2.7939505e4 -8.870e-3
7    4  4.5 3.7961159e4 3.8224812e4 -6.897e-3 8.5203092e3 8.6102220e3 -1.044e-2
>>> CKM (OLD) vs. CKM NEW
b    x1                O(as^0 a^2)                       O(as^1 a^2)
-+----+----+-----------+-----------+---------+-----------+-----------+---------
0    2 2.25 7.4136186e5 7.4101268e5  4.712e-4 1.2274280e5 1.2220099e5  4.434e-3
1 2.25  2.5 6.9378357e5 6.9466340e5 -1.267e-3 1.1516651e5 1.1541773e5 -2.177e-3
2  2.5 2.75 6.1967049e5 6.2015658e5 -7.838e-4 1.0422790e5 1.0418630e5  3.993e-4
3 2.75    3 5.1925652e5 5.1938033e5 -2.384e-4 8.8554912e4 8.8633230e4 -8.836e-4
4    3 3.25 4.0296557e5 4.0280962e5  3.872e-4 7.0854643e4 7.0423659e4  6.120e-3
5 3.25  3.5 2.8581701e5 2.8631019e5 -1.723e-3 5.1623949e4 5.1520616e4  2.006e-3
6  3.5    4 1.4463541e5 1.4467739e5 -2.901e-4 2.7691674e4 2.7785390e4 -3.373e-3
7    4  4.5 3.7961159e4 3.8098184e4 -3.597e-3 8.5203092e3 8.5649675e3 -5.214e-3


## patch_grids_with_ckm_matrix.sh
#!/bin/bash

set -euo pipefail

grids=(
    ATLAS_WM_7TEV.pineappl.lz4
    ATLAS_WP_7TEV.pineappl.lz4
    LHCB_WM_7TEV.pineappl.lz4
    LHCB_WM_8TEV.pineappl.lz4
    LHCB_WP_7TEV.pineappl.lz4
    LHCB_WP_8TEV.pineappl.lz4
)

test_pdf_set=NNPDF31_nlo_as_0118

wm_channels_1="l    entry        entry
-+------------+------------
0 1 × ( 1, -2) 1 × ( 3, -4)
1 1 × (22, -4) 1 × (22, -2)
2 1 × (21, -4) 1 × (21, -2)
3 1 × ( 1, 22) 1 × ( 3, 22)
4 1 × ( 1, 21) 1 × ( 3, 21)"

wm_channels_2="l    entry        entry
-+------------+------------
0 1 × ( 1, -2) 1 × ( 3, -4)
1 1 × (21, -4) 1 × (21, -2)
2 1 × (22, -4) 1 × (22, -2)
3 1 × ( 1, 21) 1 × ( 3, 21)
4 1 × ( 1, 22) 1 × ( 3, 22)"

wp_channels_1="l    entry        entry
-+------------+------------
0 1 × ( 2, -1) 1 × ( 4, -3)
1 1 × (21, -3) 1 × (21, -1)
2 1 × (22, -3) 1 × (22, -1)
3 1 × ( 2, 21) 1 × ( 4, 21)
4 1 × ( 2, 22) 1 × ( 4, 22)"

wp_channels_2="l    entry        entry
-+------------+------------
0 1 × ( 2, -1) 1 × ( 4, -3)
1 1 × (22, -3) 1 × (22, -1)
2 1 × (21, -3) 1 × (21, -1)
3 1 × ( 2, 22) 1 × ( 4, 22)
4 1 × ( 2, 21) 1 × ( 4, 21)"

wm_channel_du="0.9490461561 * ( 1, -2) + 0.05089536 * ( 1, -4) + 0.050940490000000005 * ( 3, -2) + 0.9473907556 * ( 3, -4) + 0.0000128881 * ( 5, -2) + 0.0017222500000000003 * ( 5, -4)"
wm_channel_gu="1.0000083656 * (21, -4) + 0.9999995342 * (21, -2)"
wm_channel_dg="0.9999415161 * ( 1, 21) + 0.9983312456 * ( 3, 21) + 0.0017351381000000003 * ( 5, 21)"
wp_channel_ud="0.9490461561 * ( 2, -1) + 0.050940490000000005 * (2, -3) + 0.0000128881 * (2, -5) + 0.05089536 * (4, -1) + 0.9473907556 * (4, -3) + 0.0017222500000000003 * (4, -5)"
wp_channel_gd="0.9999415161 * (21, -1) + 0.9983312456 * (21, -3) + 0.0017351381000000003 * (21, -5)"
wp_channel_ug="0.9999995342 * ( 2, 21) + 1.0000083656 * ( 4, 21)"

for grid in "${grids[@]}"; do
    output="${grid%.pineappl.lz4}_CKM.pineappl.lz4"

    if [[ -f ${output} ]]; then
        echo "!!! ${output} already exists"
    else
        echo ">>> patching ${output}"

        # we could also patch the photon-quark channels, but they're not used in PDF fits right now

        if diff <(echo "${wm_channels_1}") <(pineappl read --lumis "${grid}") >/dev/null; then
            pineappl write \
                --rewrite-channel 0 "${wm_channel_du}" \
                --rewrite-channel 2 "${wm_channel_gu}" \
                --rewrite-channel 4 "${wm_channel_dg}" \
                "${grid}" "${output}"
        elif diff <(echo "${wm_channels_2}") <(pineappl read --lumis "${grid}") >/dev/null; then
            pineappl write \
                --rewrite-channel 0 "${wm_channel_du}" \
                --rewrite-channel 1 "${wm_channel_gu}" \
                --rewrite-channel 3 "${wm_channel_dg}" \
                "${grid}" "${output}"
        elif diff <(echo "${wp_channels_1}") <(pineappl read --lumis "${grid}") >/dev/null; then
            pineappl write \
                --rewrite-channel 0 "${wp_channel_ud}" \
                --rewrite-channel 1 "${wp_channel_gd}" \
                --rewrite-channel 3 "${wp_channel_ug}" \
                "${grid}" "${output}"
        elif diff <(echo "${wp_channels_2}") <(pineappl read --lumis "${grid}") >/dev/null; then
            pineappl write \
                --rewrite-channel 0 "${wp_channel_ud}" \
                --rewrite-channel 2 "${wp_channel_gd}" \
                --rewrite-channel 4 "${wp_channel_ug}" \
                "${grid}" "${output}"
        else
            echo ">>> failed to patch ${grid}"
        fi

        pineappl --silence-lhapdf diff "${grid}" "${output}" "${test_pdf_set}" --orders1=a2,as1a2 --orders2=a2,as1a2 --ignore-lumis
    fi
done
	#!/bin/bash

	set -euo pipefail

	grids=(
	ATLAS_WM_7TEV.pineappl.lz4
	ATLAS_WP_7TEV.pineappl.lz4
	LHCB_WM_7TEV.pineappl.lz4
	LHCB_WM_8TEV.pineappl.lz4
	LHCB_WP_7TEV.pineappl.lz4
	LHCB_WP_8TEV.pineappl.lz4
	)

	old_grids=(
	ATLASWZRAP11CC/ATLASWZRAP11-WM_ylept.pineappl.lz4
	ATLASWZRAP11CC/ATLASWZRAP11-WP_ylept.pineappl.lz4
	LHCBWZMU7TEV/LHCBWZMU7TEV_WM_leptrap.pineappl.lz4
	LHCBWZMU8TEV/LHCBWZMU8TEV_WM_leptrap.pineappl.lz4
	LHCBWZMU7TEV/LHCBWZMU7TEV_WP_leptrap.pineappl.lz4
	LHCBWZMU8TEV/LHCBWZMU8TEV_WP_leptrap.pineappl.lz4
	)

	test_pdf_set=NNPDF31_nlo_as_0118

	for (( i = 0; i != ${#grids[@]}; ++i )); do
	old="${old_grids[$i]}"
	new="${grids[$i]}"
	new_ckm="${new%.pineappl.lz4}_CKM.pineappl.lz4"

	echo ">>> ${old} vs. ${new}"
	echo ">>> CKM (OLD) vs. NO CKM NEW"
	pineappl --silence-lhapdf diff ../applgrids/"${old}" "${new}" "${test_pdf_set}" \
	--orders1=a2,as1a2 --orders2=a2,as1a2 --scale2=1000 --ignore-lumis --ignore-bin-limits
	echo ">>> CKM (OLD) vs. CKM NEW"
	pineappl --silence-lhapdf diff ../applgrids/"${old}" "${new_ckm}" "${test_pdf_set}" \
	--orders1=a2,as1a2 --orders2=a2,as1a2 --scale2=1000 --ignore-lumis --ignore-bin-limits

	echo
	echo
	done
	>>> ATLASWZRAP11CC/ATLASWZRAP11-WM_ylept.pineappl.lz4 vs. ATLAS_WM_7TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	--+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 0 0.21 4.2240403e5 4.2418021e5 -4.187e-3 7.4550308e3 7.5751404e3 -1.586e-2
	1 0.21 0.42 4.2003285e5 4.2202888e5 -4.730e-3 7.6190059e3 7.9119070e3 -3.702e-2
	2 0.42 0.63 4.1563075e5 4.1764429e5 -4.821e-3 8.1040021e3 8.6391043e3 -6.194e-2
	3 0.63 0.84 4.0928859e5 4.1124265e5 -4.752e-3 8.4422888e3 8.9988383e3 -6.185e-2
	4 0.84 1.05 4.0079290e5 4.0284389e5 -5.091e-3 9.0842830e3 9.4793358e3 -4.168e-2
	5 1.05 1.37 3.8741291e5 3.8958316e5 -5.571e-3 9.7103988e3 1.0151402e4 -4.344e-2
	6 1.37 1.52 3.7401340e5 3.7623615e5 -5.908e-3 1.0039762e4 1.0131505e4 -9.055e-3
	7 1.52 1.74 3.6219912e5 3.6462541e5 -6.654e-3 9.9360130e3 1.0403881e4 -4.497e-2
	8 1.74 1.95 3.4804228e5 3.5053483e5 -7.111e-3 9.5605621e3 9.8270974e3 -2.712e-2
	9 1.95 2.18 3.3289241e5 3.3517621e5 -6.814e-3 8.6304707e3 9.1482968e3 -5.660e-2
	10 2.18 2.5 3.1245793e5 3.1473844e5 -7.246e-3 6.8339401e3 7.1154335e3 -3.956e-2
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	--+----+----+-----------+-----------+--------+-----------+-----------+---------
	0 0 0.21 4.2240403e5 4.2185924e5 1.291e-3 7.4550308e3 7.2050963e3 3.469e-2
	1 0.21 0.42 4.2003285e5 4.1969088e5 8.148e-4 7.6190059e3 7.5383742e3 1.070e-2
	2 0.42 0.63 4.1563075e5 4.1527346e5 8.604e-4 8.1040021e3 8.2594206e3 -1.882e-2
	3 0.63 0.84 4.0928859e5 4.0882421e5 1.136e-3 8.4422888e3 8.6120819e3 -1.972e-2
	4 0.84 1.05 4.0079290e5 4.0036730e5 1.063e-3 9.0842830e3 9.0828089e3 1.623e-4
	5 1.05 1.37 3.8741291e5 3.8702648e5 9.985e-4 9.7103988e3 9.7388763e3 -2.924e-3
	6 1.37 1.52 3.7401340e5 3.7361073e5 1.078e-3 1.0039762e4 9.7056095e3 3.443e-2
	7 1.52 1.74 3.6219912e5 3.6195622e5 6.711e-4 9.9360130e3 9.9679405e3 -3.203e-3
	8 1.74 1.95 3.4804228e5 3.4783128e5 6.066e-4 9.5605621e3 9.3823231e3 1.900e-2
	9 1.95 2.18 3.3289241e5 3.3246832e5 1.276e-3 8.6304707e3 8.6971424e3 -7.666e-3
	10 2.18 2.5 3.1245793e5 3.1207931e5 1.213e-3 6.8339401e3 6.6674362e3 2.497e-2


	>>> ATLASWZRAP11CC/ATLASWZRAP11-WP_ylept.pineappl.lz4 vs. ATLAS_WP_7TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	--+----+----+-----------+-----------+---------+------------+------------+---------
	0 0 0.21 5.7135378e5 5.7331872e5 -3.427e-3 -4.7525156e3 -4.4496358e3 6.807e-2
	1 0.21 0.42 5.7284708e5 5.7495947e5 -3.674e-3 -4.6408670e3 -4.3703932e3 6.189e-2
	2 0.42 0.63 5.7554050e5 5.7735379e5 -3.141e-3 -4.3420155e3 -4.6433339e3 -6.489e-2
	3 0.63 0.84 5.7889298e5 5.8143980e5 -4.380e-3 -4.0250183e3 -3.7056226e3 8.619e-2
	4 0.84 1.05 5.8341570e5 5.8516948e5 -2.997e-3 -3.1930355e3 -3.2423913e3 -1.522e-2
	5 1.05 1.37 5.8829101e5 5.9012841e5 -3.114e-3 -1.9993090e3 -1.7503055e3 1.423e-1
	6 1.37 1.52 5.9204449e5 5.9433809e5 -3.859e-3 -1.6054598e2 1.4259565e3 -1.113e0
	7 1.52 1.74 5.9249892e5 5.9433031e5 -3.081e-3 2.6661527e3 2.4117860e3 1.055e-1
	8 1.74 1.95 5.8785760e5 5.8958097e5 -2.923e-3 6.7750134e3 7.5921717e3 -1.076e-1
	9 1.95 2.18 5.7427968e5 5.7609457e5 -3.150e-3 1.2635183e4 1.1698721e4 8.005e-2
	10 2.18 2.5 5.3635711e5 5.3784689e5 -2.770e-3 2.1430029e4 2.2286612e4 -3.843e-2
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	--+----+----+-----------+-----------+--------+------------+------------+---------
	0 0 0.21 5.7135378e5 5.7061922e5 1.287e-3 -4.7525156e3 -4.7508919e3 3.418e-4
	1 0.21 0.42 5.7284708e5 5.7225623e5 1.032e-3 -4.6408670e3 -4.6741803e3 -7.127e-3
	2 0.42 0.63 5.7554050e5 5.7464748e5 1.554e-3 -4.3420155e3 -4.9507204e3 -1.230e-1
	3 0.63 0.84 5.7889298e5 5.7872675e5 2.872e-4 -4.0250183e3 -4.0226965e3 5.772e-4
	4 0.84 1.05 5.8341570e5 5.8245769e5 1.645e-3 -3.1930355e3 -3.5697458e3 -1.055e-1
	5 1.05 1.37 5.8829101e5 5.8743106e5 1.464e-3 -1.9993090e3 -2.0939761e3 -4.521e-2
	6 1.37 1.52 5.9204449e5 5.9166973e5 6.334e-4 -1.6054598e2 1.0621365e3 -1.151e0
	7 1.52 1.74 5.9249892e5 5.9171330e5 1.328e-3 2.6661527e3 2.0395420e3 3.072e-1
	8 1.74 1.95 5.8785760e5 5.8705863e5 1.361e-3 6.7750134e3 7.2035828e3 -5.949e-2
	9 1.95 2.18 5.7427968e5 5.7371818e5 9.787e-4 1.2635183e4 1.1305429e4 1.176e-1
	10 2.18 2.5 5.3635711e5 5.3574500e5 1.143e-3 2.1430029e4 2.1894617e4 -2.122e-2


	>>> LHCBWZMU7TEV/LHCBWZMU7TEV_WM_leptrap.pineappl.lz4 vs. LHCB_WM_7TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 3.7292336e5 3.7532595e5 -6.401e-3 6.4700672e4 6.5332710e4 -9.674e-3
	1 2.25 2.5 3.5347999e5 3.5608848e5 -7.325e-3 5.8824761e4 5.9946156e4 -1.871e-2
	2 2.5 2.75 3.3090524e5 3.3406844e5 -9.469e-3 5.3846944e4 5.4627884e4 -1.430e-2
	3 2.75 3 3.0423388e5 3.0718892e5 -9.620e-3 4.6850349e4 4.7985867e4 -2.366e-2
	4 3 3.25 2.7190663e5 2.7375224e5 -6.742e-3 4.0099659e4 4.0713216e4 -1.507e-2
	5 3.25 3.5 2.3113503e5 2.3297956e5 -7.917e-3 3.2796720e4 3.2999245e4 -6.137e-3
	6 3.5 4 1.6028275e5 1.6169805e5 -8.753e-3 2.1119242e4 2.1665733e4 -2.522e-2
	7 4 4.5 7.1251353e4 7.1892674e4 -8.921e-3 9.2469580e3 9.4199433e3 -1.836e-2
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 3.7292336e5 3.7228082e5 1.726e-3 6.4700672e4 6.4556375e4 2.235e-3
	1 2.25 2.5 3.5347999e5 3.5307300e5 1.153e-3 5.8824761e4 5.9169740e4 -5.830e-3
	2 2.5 2.75 3.3090524e5 3.3115068e5 -7.412e-4 5.3846944e4 5.3859539e4 -2.339e-4
	3 2.75 3 3.0423388e5 3.0445236e5 -7.176e-4 4.6850349e4 4.7246805e4 -8.391e-3
	4 3 3.25 2.7190663e5 2.7128834e5 2.279e-3 4.0099659e4 4.0027055e4 1.814e-3
	5 3.25 3.5 2.3113503e5 2.3087760e5 1.115e-3 3.2796720e4 3.2391407e4 1.251e-2
	6 3.5 4 1.6028275e5 1.6025570e5 1.688e-4 2.1119242e4 2.1218833e4 -4.694e-3
	7 4 4.5 7.1251353e4 7.1284320e4 -4.625e-4 9.2469580e3 9.2071618e3 4.322e-3


	>>> LHCBWZMU8TEV/LHCBWZMU8TEV_WM_leptrap.pineappl.lz4 vs. LHCB_WM_8TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 4.2643373e5 4.2995543e5 -8.191e-3 7.3094364e4 7.3808255e4 -9.672e-3
	1 2.25 2.5 4.0518638e5 4.0840319e5 -7.877e-3 6.7780367e4 6.8793814e4 -1.473e-2
	2 2.5 2.75 3.8177585e5 3.8495435e5 -8.257e-3 6.1521771e4 6.2615569e4 -1.747e-2
	3 2.75 3 3.5475650e5 3.5742431e5 -7.464e-3 5.5254814e4 5.5702896e4 -8.044e-3
	4 3 3.25 3.2045474e5 3.2372236e5 -1.009e-2 4.7085807e4 4.7987677e4 -1.879e-2
	5 3.25 3.5 2.7989006e5 2.8234548e5 -8.697e-3 3.8787358e4 3.9715001e4 -2.336e-2
	6 3.5 4 2.0482249e5 2.0668879e5 -9.030e-3 2.6447254e4 2.7031239e4 -2.160e-2
	7 4 4.5 1.0172352e5 1.0275727e5 -1.006e-2 1.2542814e4 1.2754997e4 -1.664e-2
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 4.2643373e5 4.2679976e5 -8.576e-4 7.3094364e4 7.3014109e4 1.099e-3
	1 2.25 2.5 4.0518638e5 4.0524546e5 -1.458e-4 6.7780367e4 6.7986601e4 -3.033e-3
	2 2.5 2.75 3.8177585e5 3.8185789e5 -2.149e-4 6.1521771e4 6.1812200e4 -4.699e-3
	3 2.75 3 3.5475650e5 3.5446883e5 8.115e-4 5.5254814e4 5.4918844e4 6.118e-3
	4 3 3.25 3.2045474e5 3.2099971e5 -1.698e-3 4.7085807e4 4.7244464e4 -3.358e-3
	5 3.25 3.5 2.7989006e5 2.7995161e5 -2.199e-4 3.8787358e4 3.9036955e4 -6.394e-3
	6 3.5 4 2.0482249e5 2.0494152e5 -5.808e-4 2.6447254e4 2.6504573e4 -2.163e-3
	7 4 4.5 1.0172352e5 1.0191846e5 -1.913e-3 1.2542814e4 1.2472070e4 5.672e-3


	>>> LHCBWZMU7TEV/LHCBWZMU7TEV_WP_leptrap.pineappl.lz4 vs. LHCB_WP_7TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 6.5227719e5 6.5524066e5 -4.523e-3 1.0987679e5 1.1243687e5 -2.277e-2
	1 2.25 2.5 5.9761510e5 6.0023547e5 -4.366e-3 1.0237238e5 1.0140698e5 9.520e-3
	2 2.5 2.75 5.1669345e5 5.1937942e5 -5.171e-3 8.9601247e4 9.0379386e4 -8.610e-3
	3 2.75 3 4.1745539e5 4.1841432e5 -2.292e-3 7.3703721e4 7.4337163e4 -8.521e-3
	4 3 3.25 3.0906194e5 3.1041740e5 -4.367e-3 5.6379044e4 5.6700595e4 -5.671e-3
	5 3.25 3.5 2.0870669e5 2.0961121e5 -4.315e-3 3.9340453e4 3.9651322e4 -7.840e-3
	6 3.5 4 9.7024044e4 9.7472812e4 -4.604e-3 1.9818646e4 1.9888208e4 -3.498e-3
	7 4 4.5 2.2300258e4 2.2380799e4 -3.599e-3 5.3028838e3 5.3460922e3 -8.082e-3
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 6.5227719e5 6.5259204e5 -4.825e-4 1.0987679e5 1.1173478e5 -1.663e-2
	1 2.25 2.5 5.9761510e5 5.9789562e5 -4.692e-4 1.0237238e5 1.0077814e5 1.582e-2
	2 2.5 2.75 5.1669345e5 5.1742466e5 -1.413e-3 8.9601247e4 8.9823591e4 -2.475e-3
	3 2.75 3 4.1745539e5 4.1688964e5 1.357e-3 7.3703721e4 7.3885035e4 -2.454e-3
	4 3 3.25 3.0906194e5 3.0931927e5 -8.319e-4 5.6379044e4 5.6359854e4 3.405e-4
	5 3.25 3.5 2.0870669e5 2.0888939e5 -8.746e-4 3.9340453e4 3.9416548e4 -1.931e-3
	6 3.5 4 9.7024044e4 9.7143797e4 -1.233e-3 1.9818646e4 1.9774752e4 2.220e-3
	7 4 4.5 2.2300258e4 2.2295006e4 2.356e-4 5.3028838e3 5.3159836e3 -2.464e-3


	>>> LHCBWZMU8TEV/LHCBWZMU8TEV_WP_leptrap.pineappl.lz4 vs. LHCB_WP_8TEV.pineappl.lz4
	>>> CKM (OLD) vs. NO CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 7.4136186e5 7.4379370e5 -3.270e-3 1.2274280e5 1.2291552e5 -1.405e-3
	1 2.25 2.5 6.9378357e5 6.9718568e5 -4.880e-3 1.1516651e5 1.1609118e5 -7.965e-3
	2 2.5 2.75 6.1967049e5 6.2233439e5 -4.280e-3 1.0422790e5 1.0479152e5 -5.378e-3
	3 2.75 3 5.1925652e5 5.2114837e5 -3.630e-3 8.8554912e4 8.9146035e4 -6.631e-3
	4 3 3.25 4.0296557e5 4.0414296e5 -2.913e-3 7.0854643e4 7.0828726e4 3.659e-4
	5 3.25 3.5 2.8581701e5 2.8723414e5 -4.934e-3 5.1623949e4 5.1813883e4 -3.666e-3
	6 3.5 4 1.4463541e5 1.4513268e5 -3.426e-3 2.7691674e4 2.7939505e4 -8.870e-3
	7 4 4.5 3.7961159e4 3.8224812e4 -6.897e-3 8.5203092e3 8.6102220e3 -1.044e-2
	>>> CKM (OLD) vs. CKM NEW
	b x1 O(as^0 a^2) O(as^1 a^2)
	-+----+----+-----------+-----------+---------+-----------+-----------+---------
	0 2 2.25 7.4136186e5 7.4101268e5 4.712e-4 1.2274280e5 1.2220099e5 4.434e-3
	1 2.25 2.5 6.9378357e5 6.9466340e5 -1.267e-3 1.1516651e5 1.1541773e5 -2.177e-3
	2 2.5 2.75 6.1967049e5 6.2015658e5 -7.838e-4 1.0422790e5 1.0418630e5 3.993e-4
	3 2.75 3 5.1925652e5 5.1938033e5 -2.384e-4 8.8554912e4 8.8633230e4 -8.836e-4
	4 3 3.25 4.0296557e5 4.0280962e5 3.872e-4 7.0854643e4 7.0423659e4 6.120e-3
	5 3.25 3.5 2.8581701e5 2.8631019e5 -1.723e-3 5.1623949e4 5.1520616e4 2.006e-3
	6 3.5 4 1.4463541e5 1.4467739e5 -2.901e-4 2.7691674e4 2.7785390e4 -3.373e-3
	7 4 4.5 3.7961159e4 3.8098184e4 -3.597e-3 8.5203092e3 8.5649675e3 -5.214e-3