adding theory 40_019_013 to scalevariations

Author: jacoterh

validphys2/src/validphys/scalevariations/scalevariationtheoryids.yaml

@@ -478,7 +478,6 @@ n3lo_fhmruvv: &n3lo_fhmruvv
   (175.0): 40_015_003
 
 
-
 scale_variations_for:
   # 4.1 series:
   # NNLO EXA alphas=0.118
@@ -1980,7 +1979,8 @@ scale_variations_for:
       (2, 1): 40_009_009
       (2, 2): 40_009_010
 
-  - theoryid: 40_019_005 # N3LO QCD x NLO QED, ModEv=EXA, mtop=172.5, alphas=0.120
+  # N3LO QCD x NLO QED, ModEv=EXA, mtop=172.5, alphas=0.120
+  - theoryid: 40_019_005
     variations:
 
       # IHOU
@@ -2020,7 +2020,9 @@ scale_variations_for:
       (0.115): 40_019_000
       (0.116): 40_019_001
       (0.117): 40_019_002
+      (0.1176): 40_019_012
       (0.118): 40_019_003
+      (0.1183): 40_019_013
       (0.119): 40_019_004
       (0.120): 40_019_005
       (0.121): 40_019_006
@@ -2029,8 +2031,69 @@ scale_variations_for:
       (0.124): 40_019_009
       (0.125): 40_019_010
       (0.126): 40_019_011
-
+      
       # mtop variations
       (172.5): 40_019_005 # central mtop
       (170.0): 40_020_005
       (175.0): 40_021_005
+
+  # N3LO QCD x NLO QED, ModEv=EXA, mtop=172.5, alphas=0.1183, for FLAG variation
+  - theoryid: 40_019_013
+    variations:
+
+      # IHOU
+      (0, 0, 0, 0, 0, 0, 0): 40_005_000 # central th
+      (1, 0, 0, 0, 0, 0, 0): 40_005_001 # anomalous dimensions variations
+      (2, 0, 0, 0, 0, 0, 0): 40_005_002
+      (0, 1, 0, 0, 0, 0, 0): 40_005_003
+      (0, 2, 0, 0, 0, 0, 0): 40_005_004
+      (0, 0, 1, 0, 0, 0, 0): 40_005_005
+      (0, 0, 2, 0, 0, 0, 0): 40_005_006
+      (0, 0, 0, 1, 0, 0, 0): 40_005_007
+      (0, 0, 0, 2, 0, 0, 0): 40_005_008
+      (0, 0, 0, 0, 1, 0, 0): 40_005_009
+      (0, 0, 0, 0, 2, 0, 0): 40_005_010
+      (0, 0, 0, 0, 0, 1, 0): 40_005_011
+      (0, 0, 0, 0, 0, 2, 0): 40_005_012
+      (0, 0, 0, 0, 0, 0, 1): 40_005_013
+      (0, 0, 0, 0, 0, 0, 2): 40_005_014
+
+      (1, 1): 40_005_000
+      (0.5, 1): 40_005_019 # scale variations
+      (2, 1): 40_005_020
+      (0.5, 0.5): 40_005_021
+      (1, 0.5): 40_005_022
+      (2, 0.5): 40_005_023
+      (0.5, 2): 40_005_024
+      (1, 2): 40_005_025
+      (2, 2): 40_005_026
+      lower bound DIS: 40_005_015 # DIS coeff functions fact variations
+      central DIS: 40_005_000
+      upper bound DIS: 40_005_016
+      # IHOU + MHOU missing prescription
+      (1, 0.5 hadronic): 1025 # As 1018 but DIS from 1000
+      (1, 2 hadronic): 1026 # As 1021 but DIS from 1000
+
+      #alphas variations
+      (0.115): 40_019_000
+      (0.116): 40_019_001
+      (0.117): 40_019_002
+      (0.1176): 40_019_012
+      (0.118): 40_019_003
+      (0.1183): 40_019_013
+      (0.119): 40_019_004
+      (0.120): 40_019_005
+      (0.121): 40_019_006
+      (0.122): 40_019_007
+      (0.123): 40_019_008
+      (0.124): 40_019_009
+      (0.125): 40_019_010
+      (0.126): 40_019_011
+
+      # mtop variations, alphas=0.118,
+      # we keep the same mtop variations for 0.1183 as for 0.118
+      (172.5): 40_019_003 # central mtop
+      (170.0): 40_020_003
+      (175.0): 40_021_003
+
+