1
//==========================================================================
2
// This file has been automatically generated for C++ by
3
// MadGraph5_aMC@NLO v. 2.3.3, 2015-10-25
4
// By the MadGraph5_aMC@NLO Development Team
5
// Visit launchpad.net/madgraph5 and amcatnlo.web.cern.ch
6
//==========================================================================
10
#include "Parameters_sm.h"
15
Parameters_sm::Parameters_sm(const string param_card, const bool verbose):
16
slha(SLHAReader(param_card))
18
setIndependentParameters();
19
setIndependentCouplings();
23
printIndependentParameters();
24
printIndependentCouplings();
29
void Parameters_sm::setIndependentParameters()
34
// Prepare a vector for indices
35
vector<int> indices(2, 0);
36
mdl_WH = slha.get_block_entry("decay", 25, 6.382339e-03);
37
mdl_WW = slha.get_block_entry("decay", 24, 2.047600e+00);
38
mdl_WZ = slha.get_block_entry("decay", 23, 2.441404e+00);
39
mdl_WTau = slha.get_block_entry("decay", 15, 2.270000e-12);
40
mdl_WT = slha.get_block_entry("decay", 6, 1.508336e+00);
41
mdl_ymtau = slha.get_block_entry("yukawa", 15, 1.777000e+00);
42
mdl_ymm = slha.get_block_entry("yukawa", 13, 1.056600e-01);
43
mdl_yme = slha.get_block_entry("yukawa", 11, 5.110000e-04);
44
mdl_ymt = slha.get_block_entry("yukawa", 6, 1.645000e+02);
45
mdl_ymb = slha.get_block_entry("yukawa", 5, 4.200000e+00);
46
mdl_ymc = slha.get_block_entry("yukawa", 4, 1.270000e+00);
47
mdl_etaWS = slha.get_block_entry("wolfenstein", 4, 3.410000e-01);
48
mdl_rhoWS = slha.get_block_entry("wolfenstein", 3, 1.320000e-01);
49
mdl_AWS = slha.get_block_entry("wolfenstein", 2, 8.080000e-01);
50
mdl_lamWS = slha.get_block_entry("wolfenstein", 1, 2.253000e-01);
51
aS = slha.get_block_entry("sminputs", 3, 1.180000e-01);
52
mdl_Gf = slha.get_block_entry("sminputs", 2, 1.166390e-05);
53
aEWM1 = slha.get_block_entry("sminputs", 1, 1.325070e+02);
54
mdl_MH = slha.get_block_entry("mass", 25, 1.250000e+02);
55
mdl_MZ = slha.get_block_entry("mass", 23, 9.118800e+01);
56
mdl_MTA = slha.get_block_entry("mass", 15, 1.777000e+00);
57
mdl_MM = slha.get_block_entry("mass", 13, 1.056600e-01);
58
mdl_Me = slha.get_block_entry("mass", 11, 5.110000e-04);
59
mdl_MT = slha.get_block_entry("mass", 6, 1.720000e+02);
60
mdl_MB = slha.get_block_entry("mass", 5, 4.700000e+00);
61
mdl_MC = slha.get_block_entry("mass", 4, 1.270000e+00);
62
mdl_lamWS__exp__2 = ((mdl_lamWS) * (mdl_lamWS));
63
mdl_CKM1x1 = 1. - mdl_lamWS__exp__2/2.;
64
mdl_CKM1x2 = mdl_lamWS;
65
mdl_complexi = std::complex<double> (0., 1.);
66
mdl_lamWS__exp__3 = ((mdl_lamWS) * (mdl_lamWS) * (mdl_lamWS));
67
mdl_CKM1x3 = mdl_AWS * mdl_lamWS__exp__3 * (-(mdl_etaWS * mdl_complexi) +
69
mdl_CKM2x1 = -mdl_lamWS;
70
mdl_CKM2x2 = 1. - mdl_lamWS__exp__2/2.;
71
mdl_CKM2x3 = mdl_AWS * mdl_lamWS__exp__2;
72
mdl_CKM3x1 = mdl_AWS * mdl_lamWS__exp__3 * (1. - mdl_etaWS * mdl_complexi -
74
mdl_CKM3x2 = -(mdl_AWS * mdl_lamWS__exp__2);
76
mdl_MZ__exp__2 = ((mdl_MZ) * (mdl_MZ));
77
mdl_MZ__exp__4 = ((mdl_MZ) * (mdl_MZ) * (mdl_MZ) * (mdl_MZ));
78
mdl_sqrt__2 = sqrt(2.);
79
mdl_MH__exp__2 = ((mdl_MH) * (mdl_MH));
80
mdl_conjg__CKM1x3 = conj(mdl_CKM1x3);
81
mdl_conjg__CKM2x3 = conj(mdl_CKM2x3);
82
mdl_conjg__CKM3x3 = conj(mdl_CKM3x3);
83
mdl_conjg__CKM2x1 = conj(mdl_CKM2x1);
84
mdl_conjg__CKM3x1 = conj(mdl_CKM3x1);
85
mdl_conjg__CKM2x2 = conj(mdl_CKM2x2);
86
mdl_conjg__CKM3x2 = conj(mdl_CKM3x2);
87
mdl_conjg__CKM1x1 = conj(mdl_CKM1x1);
88
mdl_conjg__CKM1x2 = conj(mdl_CKM1x2);
90
mdl_MW = sqrt(mdl_MZ__exp__2/2. + sqrt(mdl_MZ__exp__4/4. - (mdl_aEW * M_PI *
91
mdl_MZ__exp__2)/(mdl_Gf * mdl_sqrt__2)));
92
mdl_sqrt__aEW = sqrt(mdl_aEW);
93
mdl_ee = 2. * mdl_sqrt__aEW * sqrt(M_PI);
94
mdl_MW__exp__2 = ((mdl_MW) * (mdl_MW));
95
mdl_sw2 = 1. - mdl_MW__exp__2/mdl_MZ__exp__2;
96
mdl_cw = sqrt(1. - mdl_sw2);
97
mdl_sqrt__sw2 = sqrt(mdl_sw2);
98
mdl_sw = mdl_sqrt__sw2;
99
mdl_g1 = mdl_ee/mdl_cw;
100
mdl_gw = mdl_ee/mdl_sw;
101
mdl_vev = (2. * mdl_MW * mdl_sw)/mdl_ee;
102
mdl_vev__exp__2 = ((mdl_vev) * (mdl_vev));
103
mdl_lam = mdl_MH__exp__2/(2. * mdl_vev__exp__2);
104
mdl_yb = (mdl_ymb * mdl_sqrt__2)/mdl_vev;
105
mdl_yc = (mdl_ymc * mdl_sqrt__2)/mdl_vev;
106
mdl_ye = (mdl_yme * mdl_sqrt__2)/mdl_vev;
107
mdl_ym = (mdl_ymm * mdl_sqrt__2)/mdl_vev;
108
mdl_yt = (mdl_ymt * mdl_sqrt__2)/mdl_vev;
109
mdl_ytau = (mdl_ymtau * mdl_sqrt__2)/mdl_vev;
110
mdl_muH = sqrt(mdl_lam * mdl_vev__exp__2);
111
mdl_I1x31 = mdl_yb * mdl_conjg__CKM1x3;
112
mdl_I1x32 = mdl_yb * mdl_conjg__CKM2x3;
113
mdl_I1x33 = mdl_yb * mdl_conjg__CKM3x3;
114
mdl_I2x12 = mdl_yc * mdl_conjg__CKM2x1;
115
mdl_I2x13 = mdl_yt * mdl_conjg__CKM3x1;
116
mdl_I2x22 = mdl_yc * mdl_conjg__CKM2x2;
117
mdl_I2x23 = mdl_yt * mdl_conjg__CKM3x2;
118
mdl_I2x32 = mdl_yc * mdl_conjg__CKM2x3;
119
mdl_I2x33 = mdl_yt * mdl_conjg__CKM3x3;
120
mdl_I3x21 = mdl_CKM2x1 * mdl_yc;
121
mdl_I3x22 = mdl_CKM2x2 * mdl_yc;
122
mdl_I3x23 = mdl_CKM2x3 * mdl_yc;
123
mdl_I3x31 = mdl_CKM3x1 * mdl_yt;
124
mdl_I3x32 = mdl_CKM3x2 * mdl_yt;
125
mdl_I3x33 = mdl_CKM3x3 * mdl_yt;
126
mdl_I4x13 = mdl_CKM1x3 * mdl_yb;
127
mdl_I4x23 = mdl_CKM2x3 * mdl_yb;
128
mdl_I4x33 = mdl_CKM3x3 * mdl_yb;
129
mdl_ee__exp__2 = ((mdl_ee) * (mdl_ee));
130
mdl_sw__exp__2 = ((mdl_sw) * (mdl_sw));
131
mdl_cw__exp__2 = ((mdl_cw) * (mdl_cw));
133
void Parameters_sm::setIndependentCouplings()
137
void Parameters_sm::setDependentParameters()
139
mdl_sqrt__aS = sqrt(aS);
140
G = 2. * mdl_sqrt__aS * sqrt(M_PI);
141
mdl_G__exp__2 = ((G) * (G));
143
void Parameters_sm::setDependentCouplings()
145
GC_11 = mdl_complexi * G;
149
// Routines for printing out parameters
150
void Parameters_sm::printIndependentParameters()
152
cout << "sm model parameters independent of event kinematics:" << endl;
153
cout << setw(20) << "mdl_WH " << "= " << setiosflags(ios::scientific) <<
154
setw(10) << mdl_WH << endl;
155
cout << setw(20) << "mdl_WW " << "= " << setiosflags(ios::scientific) <<
156
setw(10) << mdl_WW << endl;
157
cout << setw(20) << "mdl_WZ " << "= " << setiosflags(ios::scientific) <<
158
setw(10) << mdl_WZ << endl;
159
cout << setw(20) << "mdl_WTau " << "= " << setiosflags(ios::scientific) <<
160
setw(10) << mdl_WTau << endl;
161
cout << setw(20) << "mdl_WT " << "= " << setiosflags(ios::scientific) <<
162
setw(10) << mdl_WT << endl;
163
cout << setw(20) << "mdl_ymtau " << "= " << setiosflags(ios::scientific) <<
164
setw(10) << mdl_ymtau << endl;
165
cout << setw(20) << "mdl_ymm " << "= " << setiosflags(ios::scientific) <<
166
setw(10) << mdl_ymm << endl;
167
cout << setw(20) << "mdl_yme " << "= " << setiosflags(ios::scientific) <<
168
setw(10) << mdl_yme << endl;
169
cout << setw(20) << "mdl_ymt " << "= " << setiosflags(ios::scientific) <<
170
setw(10) << mdl_ymt << endl;
171
cout << setw(20) << "mdl_ymb " << "= " << setiosflags(ios::scientific) <<
172
setw(10) << mdl_ymb << endl;
173
cout << setw(20) << "mdl_ymc " << "= " << setiosflags(ios::scientific) <<
174
setw(10) << mdl_ymc << endl;
175
cout << setw(20) << "mdl_etaWS " << "= " << setiosflags(ios::scientific) <<
176
setw(10) << mdl_etaWS << endl;
177
cout << setw(20) << "mdl_rhoWS " << "= " << setiosflags(ios::scientific) <<
178
setw(10) << mdl_rhoWS << endl;
179
cout << setw(20) << "mdl_AWS " << "= " << setiosflags(ios::scientific) <<
180
setw(10) << mdl_AWS << endl;
181
cout << setw(20) << "mdl_lamWS " << "= " << setiosflags(ios::scientific) <<
182
setw(10) << mdl_lamWS << endl;
183
cout << setw(20) << "aS " << "= " << setiosflags(ios::scientific) <<
184
setw(10) << aS << endl;
185
cout << setw(20) << "mdl_Gf " << "= " << setiosflags(ios::scientific) <<
186
setw(10) << mdl_Gf << endl;
187
cout << setw(20) << "aEWM1 " << "= " << setiosflags(ios::scientific) <<
188
setw(10) << aEWM1 << endl;
189
cout << setw(20) << "mdl_MH " << "= " << setiosflags(ios::scientific) <<
190
setw(10) << mdl_MH << endl;
191
cout << setw(20) << "mdl_MZ " << "= " << setiosflags(ios::scientific) <<
192
setw(10) << mdl_MZ << endl;
193
cout << setw(20) << "mdl_MTA " << "= " << setiosflags(ios::scientific) <<
194
setw(10) << mdl_MTA << endl;
195
cout << setw(20) << "mdl_MM " << "= " << setiosflags(ios::scientific) <<
196
setw(10) << mdl_MM << endl;
197
cout << setw(20) << "mdl_Me " << "= " << setiosflags(ios::scientific) <<
198
setw(10) << mdl_Me << endl;
199
cout << setw(20) << "mdl_MT " << "= " << setiosflags(ios::scientific) <<
200
setw(10) << mdl_MT << endl;
201
cout << setw(20) << "mdl_MB " << "= " << setiosflags(ios::scientific) <<
202
setw(10) << mdl_MB << endl;
203
cout << setw(20) << "mdl_MC " << "= " << setiosflags(ios::scientific) <<
204
setw(10) << mdl_MC << endl;
205
cout << setw(20) << "mdl_lamWS__exp__2 " << "= " <<
206
setiosflags(ios::scientific) << setw(10) << mdl_lamWS__exp__2 << endl;
207
cout << setw(20) << "mdl_CKM1x1 " << "= " << setiosflags(ios::scientific)
208
<< setw(10) << mdl_CKM1x1 << endl;
209
cout << setw(20) << "mdl_CKM1x2 " << "= " << setiosflags(ios::scientific)
210
<< setw(10) << mdl_CKM1x2 << endl;
211
cout << setw(20) << "mdl_complexi " << "= " << setiosflags(ios::scientific)
212
<< setw(10) << mdl_complexi << endl;
213
cout << setw(20) << "mdl_lamWS__exp__3 " << "= " <<
214
setiosflags(ios::scientific) << setw(10) << mdl_lamWS__exp__3 << endl;
215
cout << setw(20) << "mdl_CKM1x3 " << "= " << setiosflags(ios::scientific)
216
<< setw(10) << mdl_CKM1x3 << endl;
217
cout << setw(20) << "mdl_CKM2x1 " << "= " << setiosflags(ios::scientific)
218
<< setw(10) << mdl_CKM2x1 << endl;
219
cout << setw(20) << "mdl_CKM2x2 " << "= " << setiosflags(ios::scientific)
220
<< setw(10) << mdl_CKM2x2 << endl;
221
cout << setw(20) << "mdl_CKM2x3 " << "= " << setiosflags(ios::scientific)
222
<< setw(10) << mdl_CKM2x3 << endl;
223
cout << setw(20) << "mdl_CKM3x1 " << "= " << setiosflags(ios::scientific)
224
<< setw(10) << mdl_CKM3x1 << endl;
225
cout << setw(20) << "mdl_CKM3x2 " << "= " << setiosflags(ios::scientific)
226
<< setw(10) << mdl_CKM3x2 << endl;
227
cout << setw(20) << "mdl_CKM3x3 " << "= " << setiosflags(ios::scientific)
228
<< setw(10) << mdl_CKM3x3 << endl;
229
cout << setw(20) << "mdl_MZ__exp__2 " << "= " <<
230
setiosflags(ios::scientific) << setw(10) << mdl_MZ__exp__2 << endl;
231
cout << setw(20) << "mdl_MZ__exp__4 " << "= " <<
232
setiosflags(ios::scientific) << setw(10) << mdl_MZ__exp__4 << endl;
233
cout << setw(20) << "mdl_sqrt__2 " << "= " << setiosflags(ios::scientific)
234
<< setw(10) << mdl_sqrt__2 << endl;
235
cout << setw(20) << "mdl_MH__exp__2 " << "= " <<
236
setiosflags(ios::scientific) << setw(10) << mdl_MH__exp__2 << endl;
237
cout << setw(20) << "mdl_conjg__CKM1x3 " << "= " <<
238
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM1x3 << endl;
239
cout << setw(20) << "mdl_conjg__CKM2x3 " << "= " <<
240
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM2x3 << endl;
241
cout << setw(20) << "mdl_conjg__CKM3x3 " << "= " <<
242
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM3x3 << endl;
243
cout << setw(20) << "mdl_conjg__CKM2x1 " << "= " <<
244
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM2x1 << endl;
245
cout << setw(20) << "mdl_conjg__CKM3x1 " << "= " <<
246
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM3x1 << endl;
247
cout << setw(20) << "mdl_conjg__CKM2x2 " << "= " <<
248
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM2x2 << endl;
249
cout << setw(20) << "mdl_conjg__CKM3x2 " << "= " <<
250
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM3x2 << endl;
251
cout << setw(20) << "mdl_conjg__CKM1x1 " << "= " <<
252
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM1x1 << endl;
253
cout << setw(20) << "mdl_conjg__CKM1x2 " << "= " <<
254
setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM1x2 << endl;
255
cout << setw(20) << "mdl_aEW " << "= " << setiosflags(ios::scientific) <<
256
setw(10) << mdl_aEW << endl;
257
cout << setw(20) << "mdl_MW " << "= " << setiosflags(ios::scientific) <<
258
setw(10) << mdl_MW << endl;
259
cout << setw(20) << "mdl_sqrt__aEW " << "= " <<
260
setiosflags(ios::scientific) << setw(10) << mdl_sqrt__aEW << endl;
261
cout << setw(20) << "mdl_ee " << "= " << setiosflags(ios::scientific) <<
262
setw(10) << mdl_ee << endl;
263
cout << setw(20) << "mdl_MW__exp__2 " << "= " <<
264
setiosflags(ios::scientific) << setw(10) << mdl_MW__exp__2 << endl;
265
cout << setw(20) << "mdl_sw2 " << "= " << setiosflags(ios::scientific) <<
266
setw(10) << mdl_sw2 << endl;
267
cout << setw(20) << "mdl_cw " << "= " << setiosflags(ios::scientific) <<
268
setw(10) << mdl_cw << endl;
269
cout << setw(20) << "mdl_sqrt__sw2 " << "= " <<
270
setiosflags(ios::scientific) << setw(10) << mdl_sqrt__sw2 << endl;
271
cout << setw(20) << "mdl_sw " << "= " << setiosflags(ios::scientific) <<
272
setw(10) << mdl_sw << endl;
273
cout << setw(20) << "mdl_g1 " << "= " << setiosflags(ios::scientific) <<
274
setw(10) << mdl_g1 << endl;
275
cout << setw(20) << "mdl_gw " << "= " << setiosflags(ios::scientific) <<
276
setw(10) << mdl_gw << endl;
277
cout << setw(20) << "mdl_vev " << "= " << setiosflags(ios::scientific) <<
278
setw(10) << mdl_vev << endl;
279
cout << setw(20) << "mdl_vev__exp__2 " << "= " <<
280
setiosflags(ios::scientific) << setw(10) << mdl_vev__exp__2 << endl;
281
cout << setw(20) << "mdl_lam " << "= " << setiosflags(ios::scientific) <<
282
setw(10) << mdl_lam << endl;
283
cout << setw(20) << "mdl_yb " << "= " << setiosflags(ios::scientific) <<
284
setw(10) << mdl_yb << endl;
285
cout << setw(20) << "mdl_yc " << "= " << setiosflags(ios::scientific) <<
286
setw(10) << mdl_yc << endl;
287
cout << setw(20) << "mdl_ye " << "= " << setiosflags(ios::scientific) <<
288
setw(10) << mdl_ye << endl;
289
cout << setw(20) << "mdl_ym " << "= " << setiosflags(ios::scientific) <<
290
setw(10) << mdl_ym << endl;
291
cout << setw(20) << "mdl_yt " << "= " << setiosflags(ios::scientific) <<
292
setw(10) << mdl_yt << endl;
293
cout << setw(20) << "mdl_ytau " << "= " << setiosflags(ios::scientific) <<
294
setw(10) << mdl_ytau << endl;
295
cout << setw(20) << "mdl_muH " << "= " << setiosflags(ios::scientific) <<
296
setw(10) << mdl_muH << endl;
297
cout << setw(20) << "mdl_I1x31 " << "= " << setiosflags(ios::scientific) <<
298
setw(10) << mdl_I1x31 << endl;
299
cout << setw(20) << "mdl_I1x32 " << "= " << setiosflags(ios::scientific) <<
300
setw(10) << mdl_I1x32 << endl;
301
cout << setw(20) << "mdl_I1x33 " << "= " << setiosflags(ios::scientific) <<
302
setw(10) << mdl_I1x33 << endl;
303
cout << setw(20) << "mdl_I2x12 " << "= " << setiosflags(ios::scientific) <<
304
setw(10) << mdl_I2x12 << endl;
305
cout << setw(20) << "mdl_I2x13 " << "= " << setiosflags(ios::scientific) <<
306
setw(10) << mdl_I2x13 << endl;
307
cout << setw(20) << "mdl_I2x22 " << "= " << setiosflags(ios::scientific) <<
308
setw(10) << mdl_I2x22 << endl;
309
cout << setw(20) << "mdl_I2x23 " << "= " << setiosflags(ios::scientific) <<
310
setw(10) << mdl_I2x23 << endl;
311
cout << setw(20) << "mdl_I2x32 " << "= " << setiosflags(ios::scientific) <<
312
setw(10) << mdl_I2x32 << endl;
313
cout << setw(20) << "mdl_I2x33 " << "= " << setiosflags(ios::scientific) <<
314
setw(10) << mdl_I2x33 << endl;
315
cout << setw(20) << "mdl_I3x21 " << "= " << setiosflags(ios::scientific) <<
316
setw(10) << mdl_I3x21 << endl;
317
cout << setw(20) << "mdl_I3x22 " << "= " << setiosflags(ios::scientific) <<
318
setw(10) << mdl_I3x22 << endl;
319
cout << setw(20) << "mdl_I3x23 " << "= " << setiosflags(ios::scientific) <<
320
setw(10) << mdl_I3x23 << endl;
321
cout << setw(20) << "mdl_I3x31 " << "= " << setiosflags(ios::scientific) <<
322
setw(10) << mdl_I3x31 << endl;
323
cout << setw(20) << "mdl_I3x32 " << "= " << setiosflags(ios::scientific) <<
324
setw(10) << mdl_I3x32 << endl;
325
cout << setw(20) << "mdl_I3x33 " << "= " << setiosflags(ios::scientific) <<
326
setw(10) << mdl_I3x33 << endl;
327
cout << setw(20) << "mdl_I4x13 " << "= " << setiosflags(ios::scientific) <<
328
setw(10) << mdl_I4x13 << endl;
329
cout << setw(20) << "mdl_I4x23 " << "= " << setiosflags(ios::scientific) <<
330
setw(10) << mdl_I4x23 << endl;
331
cout << setw(20) << "mdl_I4x33 " << "= " << setiosflags(ios::scientific) <<
332
setw(10) << mdl_I4x33 << endl;
333
cout << setw(20) << "mdl_ee__exp__2 " << "= " <<
334
setiosflags(ios::scientific) << setw(10) << mdl_ee__exp__2 << endl;
335
cout << setw(20) << "mdl_sw__exp__2 " << "= " <<
336
setiosflags(ios::scientific) << setw(10) << mdl_sw__exp__2 << endl;
337
cout << setw(20) << "mdl_cw__exp__2 " << "= " <<
338
setiosflags(ios::scientific) << setw(10) << mdl_cw__exp__2 << endl;
340
void Parameters_sm::printIndependentCouplings()
342
cout << "sm model couplings independent of event kinematics:" << endl;
345
void Parameters_sm::printDependentParameters()
347
cout << "sm model parameters dependent on event kinematics:" << endl;
348
cout << setw(20) << "mdl_sqrt__aS " << "= " << setiosflags(ios::scientific)
349
<< setw(10) << mdl_sqrt__aS << endl;
350
cout << setw(20) << "G " << "= " << setiosflags(ios::scientific) <<
351
setw(10) << G << endl;
352
cout << setw(20) << "mdl_G__exp__2 " << "= " <<
353
setiosflags(ios::scientific) << setw(10) << mdl_G__exp__2 << endl;
355
void Parameters_sm::printDependentCouplings()
357
cout << "sm model couplings dependent on event kinematics:" << endl;
358
cout << setw(20) << "GC_11 " << "= " << setiosflags(ios::scientific) <<
359
setw(10) << GC_11 << endl;
360
cout << setw(20) << "GC_10 " << "= " << setiosflags(ios::scientific) <<
361
setw(10) << GC_10 << endl;