ATLAS Offline Software
Pythia8_ShowerWeights.py
Go to the documentation of this file.
1 # Warnings:
2 # UserHooks Warning: the calculation of uncertainty variations will only be consistent in the absence of any external modifications to the shower branching probabilities via the UserHooks framework. It is therefore strongly advised to avoid combining the automated uncertainty calculations with any such UserHooks modifications.
3 #
4 #Merging Warning: in multi-jet merging approaches, trial showers are used to generate missing Sudakov factor corrections to the hard matrix elements. Currently that framework is not consistently combined with the variations introduced here, so the two should not be used simultaneously. This means shower weights should NOT be used with CKKLW or FxFx merging!
5 #
6 #
7 # Pythia8 shower weights are only available in release 2.26 and later. The
8 # test below checks the Pythia8 release and also verifies that the Pythia_i
9 # tag is recent enought to store the shower weights.
10 
11 if "ShowerWeightNames" in genSeq.Pythia8.__slots__.keys():
12  print ("Initalizing Shower Weights from Pythia8_ShowerWeights.py")
13  if("NNPDF" in evgenConfig.tune):
14  genSeq.Pythia8.Commands += [
15  'UncertaintyBands:doVariations = on',
16  'UncertaintyBands:List = {\
17  Var3cUp isr:muRfac=0.549241,\
18  Var3cDown isr:muRfac=1.960832,\
19  isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
20  isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
21  isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
22  isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
23  isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
24  isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
25  isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
26  isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
27  isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
28  isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
29  isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
30  isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
31  isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
32  isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
33  isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
34  isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
35  isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
36  isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
37  isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
38  isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
39  hardHi fsr:cNS=2.0 isr:cNS=2.0,\
40  hardLo fsr:cNS=-2.0 isr:cNS=-2.0,\
41  isr:PDF:plus isr:PDF:plus=1,\
42  isr:PDF:minus isr:PDF:minus=2\
43  }']
44 
45  genSeq.Pythia8.ShowerWeightNames = [
46  "Var3cUp",
47  "Var3cDown",
48  "isr:muRfac=2.0_fsr:muRfac=2.0",
49  "isr:muRfac=2.0_fsr:muRfac=1.0",
50  "isr:muRfac=2.0_fsr:muRfac=0.5",
51  "isr:muRfac=1.0_fsr:muRfac=2.0",
52  "isr:muRfac=1.0_fsr:muRfac=0.5",
53  "isr:muRfac=0.5_fsr:muRfac=2.0",
54  "isr:muRfac=0.5_fsr:muRfac=1.0",
55  "isr:muRfac=0.5_fsr:muRfac=0.5",
56  "isr:muRfac=1.75_fsr:muRfac=1.0",
57  "isr:muRfac=1.5_fsr:muRfac=1.0",
58  "isr:muRfac=1.25_fsr:muRfac=1.0",
59  "isr:muRfac=0.625_fsr:muRfac=1.0",
60  "isr:muRfac=0.75_fsr:muRfac=1.0",
61  "isr:muRfac=0.875_fsr:muRfac=1.0",
62  "isr:muRfac=1.0_fsr:muRfac=1.75",
63  "isr:muRfac=1.0_fsr:muRfac=1.5",
64  "isr:muRfac=1.0_fsr:muRfac=1.25",
65  "isr:muRfac=1.0_fsr:muRfac=0.625",
66  "isr:muRfac=1.0_fsr:muRfac=0.75",
67  "isr:muRfac=1.0_fsr:muRfac=0.875",
68  "hardHi",
69  "hardLo",
70  "isr:PDF:plus",
71  "isr:PDF:minus"
72  ]
73  else:
74  genSeq.Pythia8.Commands += [
75  'UncertaintyBands:doVariations = on',
76  'UncertaintyBands:List = {\
77  Var3cUp isr:muRfac=0.549241,\
78  Var3cDown isr:muRfac=1.960832,\
79  isr:muRfac=2.0_fsr:muRfac=2.0 isr:muRfac=2.0 fsr:muRfac=2.0,\
80  isr:muRfac=2.0_fsr:muRfac=1.0 isr:muRfac=2.0 fsr:muRfac=1.0,\
81  isr:muRfac=2.0_fsr:muRfac=0.5 isr:muRfac=2.0 fsr:muRfac=0.5,\
82  isr:muRfac=1.0_fsr:muRfac=2.0 isr:muRfac=1.0 fsr:muRfac=2.0,\
83  isr:muRfac=1.0_fsr:muRfac=0.5 isr:muRfac=1.0 fsr:muRfac=0.5,\
84  isr:muRfac=0.5_fsr:muRfac=2.0 isr:muRfac=0.5 fsr:muRfac=2.0,\
85  isr:muRfac=0.5_fsr:muRfac=1.0 isr:muRfac=0.5 fsr:muRfac=1.0,\
86  isr:muRfac=0.5_fsr:muRfac=0.5 isr:muRfac=0.5 fsr:muRfac=0.5,\
87  isr:muRfac=1.75_fsr:muRfac=1.0 isr:muRfac=1.75 fsr:muRfac=1.0,\
88  isr:muRfac=1.5_fsr:muRfac=1.0 isr:muRfac=1.5 fsr:muRfac=1.0,\
89  isr:muRfac=1.25_fsr:muRfac=1.0 isr:muRfac=1.25 fsr:muRfac=1.0,\
90  isr:muRfac=0.625_fsr:muRfac=1.0 isr:muRfac=0.625 fsr:muRfac=1.0,\
91  isr:muRfac=0.75_fsr:muRfac=1.0 isr:muRfac=0.75 fsr:muRfac=1.0,\
92  isr:muRfac=0.875_fsr:muRfac=1.0 isr:muRfac=0.875 fsr:muRfac=1.0,\
93  isr:muRfac=1.0_fsr:muRfac=1.75 isr:muRfac=1.0 fsr:muRfac=1.75,\
94  isr:muRfac=1.0_fsr:muRfac=1.5 isr:muRfac=1.0 fsr:muRfac=1.5,\
95  isr:muRfac=1.0_fsr:muRfac=1.25 isr:muRfac=1.0 fsr:muRfac=1.25,\
96  isr:muRfac=1.0_fsr:muRfac=0.625 isr:muRfac=1.0 fsr:muRfac=0.625,\
97  isr:muRfac=1.0_fsr:muRfac=0.75 isr:muRfac=1.0 fsr:muRfac=0.75,\
98  isr:muRfac=1.0_fsr:muRfac=0.875 isr:muRfac=1.0 fsr:muRfac=0.875,\
99  hardHi fsr:cNS=2.0 isr:cNS=2.0,\
100  hardLo fsr:cNS=-2.0 isr:cNS=-2.0\
101  }']
102 
103  genSeq.Pythia8.ShowerWeightNames = [
104  "Var3cUp",
105  "Var3cDown",
106  "isr:muRfac=2.0_fsr:muRfac=2.0",
107  "isr:muRfac=2.0_fsr:muRfac=1.0",
108  "isr:muRfac=2.0_fsr:muRfac=0.5",
109  "isr:muRfac=1.0_fsr:muRfac=2.0",
110  "isr:muRfac=1.0_fsr:muRfac=0.5",
111  "isr:muRfac=0.5_fsr:muRfac=2.0",
112  "isr:muRfac=0.5_fsr:muRfac=1.0",
113  "isr:muRfac=0.5_fsr:muRfac=0.5",
114  "isr:muRfac=1.75_fsr:muRfac=1.0",
115  "isr:muRfac=1.5_fsr:muRfac=1.0",
116  "isr:muRfac=1.25_fsr:muRfac=1.0",
117  "isr:muRfac=0.625_fsr:muRfac=1.0",
118  "isr:muRfac=0.75_fsr:muRfac=1.0",
119  "isr:muRfac=0.875_fsr:muRfac=1.0",
120  "isr:muRfac=1.0_fsr:muRfac=1.75",
121  "isr:muRfac=1.0_fsr:muRfac=1.5",
122  "isr:muRfac=1.0_fsr:muRfac=1.25",
123  "isr:muRfac=1.0_fsr:muRfac=0.625",
124  "isr:muRfac=1.0_fsr:muRfac=0.75",
125  "isr:muRfac=1.0_fsr:muRfac=0.875",
126  "hardHi",
127  "hardLo"
128  ]
129 
130 
if
if(febId1==febId2)
Definition: LArRodBlockPhysicsV0.cxx:567
Generated on Mon Dec 23 2024 21:16:58 for ATLAS Offline Software by   doxygen 1.8.18