← Back to branch summary

~sebastien-wertz/mg5amcnlo/standalone_cpp_mem

« back to all changes in this revision

Viewing changes to tests/input_files/IOTestsComparison/IOExportStandaloneCPPMEMTest/MultiProcessTest/src/Parameters_sm.cc

  • Committer: Sébastien Wertz
  • Date: 2015-10-30 14:14:25 UTC
  • Revision ID: swertz@ingrid-ui1-20151030141425-0ymwh93fx6xhcrfs
Finished IOTest, small corrections to exporter

Show diffs side-by-side

added added

removed removed

Lines of Context:
tests/input_files/IOTestsComparison/IOExportStandaloneCPPMEMTest/MultiProcessTest/src/Parameters_sm.cc
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
 
//==========================================================================
7
 
 
8
 
#include <iostream> 
9
 
#include <iomanip> 
10
 
#include "Parameters_sm.h"
11
 
 
12
 
// Initialize static instance
13
 
Parameters_sm * Parameters_sm::instance = 0; 
14
 
 
15
 
// Function to get static instance - only one instance per program
16
 
Parameters_sm * Parameters_sm::getInstance()
17
 
{
18
 
  if (instance == 0)
19
 
    instance = new Parameters_sm(); 
20
 
 
21
 
  return instance; 
22
 
}
23
 
 
24
 
void Parameters_sm::setIndependentParameters(SLHAReader& slha)
25
 
{
26
 
  // Define "zero"
27
 
  zero = 0; 
28
 
  ZERO = 0; 
29
 
  // Prepare a vector for indices
30
 
  vector<int> indices(2, 0); 
31
 
  mdl_WH = slha.get_block_entry("decay", 25, 6.382339e-03); 
32
 
  mdl_WW = slha.get_block_entry("decay", 24, 2.047600e+00); 
33
 
  mdl_WZ = slha.get_block_entry("decay", 23, 2.441404e+00); 
34
 
  mdl_WT = slha.get_block_entry("decay", 6, 1.491500e+00); 
35
 
  mdl_ymtau = slha.get_block_entry("yukawa", 15, 1.777000e+00); 
36
 
  mdl_ymt = slha.get_block_entry("yukawa", 6, 1.730000e+02); 
37
 
  mdl_ymb = slha.get_block_entry("yukawa", 5, 4.700000e+00); 
38
 
  aS = slha.get_block_entry("sminputs", 3, 1.180000e-01); 
39
 
  mdl_Gf = slha.get_block_entry("sminputs", 2, 1.166390e-05); 
40
 
  aEWM1 = slha.get_block_entry("sminputs", 1, 1.325070e+02); 
41
 
  mdl_MH = slha.get_block_entry("mass", 25, 1.250000e+02); 
42
 
  mdl_MZ = slha.get_block_entry("mass", 23, 9.118800e+01); 
43
 
  mdl_MTA = slha.get_block_entry("mass", 15, 1.777000e+00); 
44
 
  mdl_MT = slha.get_block_entry("mass", 6, 1.730000e+02); 
45
 
  mdl_MB = slha.get_block_entry("mass", 5, 4.700000e+00); 
46
 
  mdl_conjg__CKM3x3 = 1.; 
47
 
  mdl_CKM3x3 = 1.; 
48
 
  mdl_conjg__CKM1x1 = 1.; 
49
 
  mdl_complexi = std::complex<double> (0., 1.); 
50
 
  mdl_MZ__exp__2 = ((mdl_MZ) * (mdl_MZ)); 
51
 
  mdl_MZ__exp__4 = ((mdl_MZ) * (mdl_MZ) * (mdl_MZ) * (mdl_MZ)); 
52
 
  mdl_sqrt__2 = sqrt(2.); 
53
 
  mdl_MH__exp__2 = ((mdl_MH) * (mdl_MH)); 
54
 
  mdl_aEW = 1./aEWM1; 
55
 
  mdl_MW = sqrt(mdl_MZ__exp__2/2. + sqrt(mdl_MZ__exp__4/4. - (mdl_aEW * M_PI *
56
 
      mdl_MZ__exp__2)/(mdl_Gf * mdl_sqrt__2)));
57
 
  mdl_sqrt__aEW = sqrt(mdl_aEW); 
58
 
  mdl_ee = 2. * mdl_sqrt__aEW * sqrt(M_PI); 
59
 
  mdl_MW__exp__2 = ((mdl_MW) * (mdl_MW)); 
60
 
  mdl_sw2 = 1. - mdl_MW__exp__2/mdl_MZ__exp__2; 
61
 
  mdl_cw = sqrt(1. - mdl_sw2); 
62
 
  mdl_sqrt__sw2 = sqrt(mdl_sw2); 
63
 
  mdl_sw = mdl_sqrt__sw2; 
64
 
  mdl_g1 = mdl_ee/mdl_cw; 
65
 
  mdl_gw = mdl_ee/mdl_sw; 
66
 
  mdl_vev = (2. * mdl_MW * mdl_sw)/mdl_ee; 
67
 
  mdl_vev__exp__2 = ((mdl_vev) * (mdl_vev)); 
68
 
  mdl_lam = mdl_MH__exp__2/(2. * mdl_vev__exp__2); 
69
 
  mdl_yb = (mdl_ymb * mdl_sqrt__2)/mdl_vev; 
70
 
  mdl_yt = (mdl_ymt * mdl_sqrt__2)/mdl_vev; 
71
 
  mdl_ytau = (mdl_ymtau * mdl_sqrt__2)/mdl_vev; 
72
 
  mdl_muH = sqrt(mdl_lam * mdl_vev__exp__2); 
73
 
  mdl_I1x33 = mdl_yb * mdl_conjg__CKM3x3; 
74
 
  mdl_I2x33 = mdl_yt * mdl_conjg__CKM3x3; 
75
 
  mdl_I3x33 = mdl_CKM3x3 * mdl_yt; 
76
 
  mdl_I4x33 = mdl_CKM3x3 * mdl_yb; 
77
 
  mdl_ee__exp__2 = ((mdl_ee) * (mdl_ee)); 
78
 
  mdl_sw__exp__2 = ((mdl_sw) * (mdl_sw)); 
79
 
  mdl_cw__exp__2 = ((mdl_cw) * (mdl_cw)); 
80
 
}
81
 
void Parameters_sm::setIndependentCouplings()
82
 
{
83
 
 
84
 
}
85
 
void Parameters_sm::setDependentParameters()
86
 
{
87
 
  mdl_sqrt__aS = sqrt(aS); 
88
 
  G = 2. * mdl_sqrt__aS * sqrt(M_PI); 
89
 
  mdl_G__exp__2 = ((G) * (G)); 
90
 
}
91
 
void Parameters_sm::setDependentCouplings()
92
 
{
93
 
  GC_11 = mdl_complexi * G; 
94
 
}
95
 
 
96
 
// Routines for printing out parameters
97
 
void Parameters_sm::printIndependentParameters()
98
 
{
99
 
  cout <<  "sm model parameters independent of event kinematics:" << endl; 
100
 
  cout << setw(20) <<  "mdl_WH " <<  "= " << setiosflags(ios::scientific) <<
101
 
      setw(10) << mdl_WH << endl;
102
 
  cout << setw(20) <<  "mdl_WW " <<  "= " << setiosflags(ios::scientific) <<
103
 
      setw(10) << mdl_WW << endl;
104
 
  cout << setw(20) <<  "mdl_WZ " <<  "= " << setiosflags(ios::scientific) <<
105
 
      setw(10) << mdl_WZ << endl;
106
 
  cout << setw(20) <<  "mdl_WT " <<  "= " << setiosflags(ios::scientific) <<
107
 
      setw(10) << mdl_WT << endl;
108
 
  cout << setw(20) <<  "mdl_ymtau " <<  "= " << setiosflags(ios::scientific) <<
109
 
      setw(10) << mdl_ymtau << endl;
110
 
  cout << setw(20) <<  "mdl_ymt " <<  "= " << setiosflags(ios::scientific) <<
111
 
      setw(10) << mdl_ymt << endl;
112
 
  cout << setw(20) <<  "mdl_ymb " <<  "= " << setiosflags(ios::scientific) <<
113
 
      setw(10) << mdl_ymb << endl;
114
 
  cout << setw(20) <<  "aS " <<  "= " << setiosflags(ios::scientific) <<
115
 
      setw(10) << aS << endl;
116
 
  cout << setw(20) <<  "mdl_Gf " <<  "= " << setiosflags(ios::scientific) <<
117
 
      setw(10) << mdl_Gf << endl;
118
 
  cout << setw(20) <<  "aEWM1 " <<  "= " << setiosflags(ios::scientific) <<
119
 
      setw(10) << aEWM1 << endl;
120
 
  cout << setw(20) <<  "mdl_MH " <<  "= " << setiosflags(ios::scientific) <<
121
 
      setw(10) << mdl_MH << endl;
122
 
  cout << setw(20) <<  "mdl_MZ " <<  "= " << setiosflags(ios::scientific) <<
123
 
      setw(10) << mdl_MZ << endl;
124
 
  cout << setw(20) <<  "mdl_MTA " <<  "= " << setiosflags(ios::scientific) <<
125
 
      setw(10) << mdl_MTA << endl;
126
 
  cout << setw(20) <<  "mdl_MT " <<  "= " << setiosflags(ios::scientific) <<
127
 
      setw(10) << mdl_MT << endl;
128
 
  cout << setw(20) <<  "mdl_MB " <<  "= " << setiosflags(ios::scientific) <<
129
 
      setw(10) << mdl_MB << endl;
130
 
  cout << setw(20) <<  "mdl_conjg__CKM3x3 " <<  "= " <<
131
 
      setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM3x3 << endl;
132
 
  cout << setw(20) <<  "mdl_CKM3x3 " <<  "= " << setiosflags(ios::scientific)
133
 
      << setw(10) << mdl_CKM3x3 << endl;
134
 
  cout << setw(20) <<  "mdl_conjg__CKM1x1 " <<  "= " <<
135
 
      setiosflags(ios::scientific) << setw(10) << mdl_conjg__CKM1x1 << endl;
136
 
  cout << setw(20) <<  "mdl_complexi " <<  "= " << setiosflags(ios::scientific)
137
 
      << setw(10) << mdl_complexi << endl;
138
 
  cout << setw(20) <<  "mdl_MZ__exp__2 " <<  "= " <<
139
 
      setiosflags(ios::scientific) << setw(10) << mdl_MZ__exp__2 << endl;
140
 
  cout << setw(20) <<  "mdl_MZ__exp__4 " <<  "= " <<
141
 
      setiosflags(ios::scientific) << setw(10) << mdl_MZ__exp__4 << endl;
142
 
  cout << setw(20) <<  "mdl_sqrt__2 " <<  "= " << setiosflags(ios::scientific)
143
 
      << setw(10) << mdl_sqrt__2 << endl;
144
 
  cout << setw(20) <<  "mdl_MH__exp__2 " <<  "= " <<
145
 
      setiosflags(ios::scientific) << setw(10) << mdl_MH__exp__2 << endl;
146
 
  cout << setw(20) <<  "mdl_aEW " <<  "= " << setiosflags(ios::scientific) <<
147
 
      setw(10) << mdl_aEW << endl;
148
 
  cout << setw(20) <<  "mdl_MW " <<  "= " << setiosflags(ios::scientific) <<
149
 
      setw(10) << mdl_MW << endl;
150
 
  cout << setw(20) <<  "mdl_sqrt__aEW " <<  "= " <<
151
 
      setiosflags(ios::scientific) << setw(10) << mdl_sqrt__aEW << endl;
152
 
  cout << setw(20) <<  "mdl_ee " <<  "= " << setiosflags(ios::scientific) <<
153
 
      setw(10) << mdl_ee << endl;
154
 
  cout << setw(20) <<  "mdl_MW__exp__2 " <<  "= " <<
155
 
      setiosflags(ios::scientific) << setw(10) << mdl_MW__exp__2 << endl;
156
 
  cout << setw(20) <<  "mdl_sw2 " <<  "= " << setiosflags(ios::scientific) <<
157
 
      setw(10) << mdl_sw2 << endl;
158
 
  cout << setw(20) <<  "mdl_cw " <<  "= " << setiosflags(ios::scientific) <<
159
 
      setw(10) << mdl_cw << endl;
160
 
  cout << setw(20) <<  "mdl_sqrt__sw2 " <<  "= " <<
161
 
      setiosflags(ios::scientific) << setw(10) << mdl_sqrt__sw2 << endl;
162
 
  cout << setw(20) <<  "mdl_sw " <<  "= " << setiosflags(ios::scientific) <<
163
 
      setw(10) << mdl_sw << endl;
164
 
  cout << setw(20) <<  "mdl_g1 " <<  "= " << setiosflags(ios::scientific) <<
165
 
      setw(10) << mdl_g1 << endl;
166
 
  cout << setw(20) <<  "mdl_gw " <<  "= " << setiosflags(ios::scientific) <<
167
 
      setw(10) << mdl_gw << endl;
168
 
  cout << setw(20) <<  "mdl_vev " <<  "= " << setiosflags(ios::scientific) <<
169
 
      setw(10) << mdl_vev << endl;
170
 
  cout << setw(20) <<  "mdl_vev__exp__2 " <<  "= " <<
171
 
      setiosflags(ios::scientific) << setw(10) << mdl_vev__exp__2 << endl;
172
 
  cout << setw(20) <<  "mdl_lam " <<  "= " << setiosflags(ios::scientific) <<
173
 
      setw(10) << mdl_lam << endl;
174
 
  cout << setw(20) <<  "mdl_yb " <<  "= " << setiosflags(ios::scientific) <<
175
 
      setw(10) << mdl_yb << endl;
176
 
  cout << setw(20) <<  "mdl_yt " <<  "= " << setiosflags(ios::scientific) <<
177
 
      setw(10) << mdl_yt << endl;
178
 
  cout << setw(20) <<  "mdl_ytau " <<  "= " << setiosflags(ios::scientific) <<
179
 
      setw(10) << mdl_ytau << endl;
180
 
  cout << setw(20) <<  "mdl_muH " <<  "= " << setiosflags(ios::scientific) <<
181
 
      setw(10) << mdl_muH << endl;
182
 
  cout << setw(20) <<  "mdl_I1x33 " <<  "= " << setiosflags(ios::scientific) <<
183
 
      setw(10) << mdl_I1x33 << endl;
184
 
  cout << setw(20) <<  "mdl_I2x33 " <<  "= " << setiosflags(ios::scientific) <<
185
 
      setw(10) << mdl_I2x33 << endl;
186
 
  cout << setw(20) <<  "mdl_I3x33 " <<  "= " << setiosflags(ios::scientific) <<
187
 
      setw(10) << mdl_I3x33 << endl;
188
 
  cout << setw(20) <<  "mdl_I4x33 " <<  "= " << setiosflags(ios::scientific) <<
189
 
      setw(10) << mdl_I4x33 << endl;
190
 
  cout << setw(20) <<  "mdl_ee__exp__2 " <<  "= " <<
191
 
      setiosflags(ios::scientific) << setw(10) << mdl_ee__exp__2 << endl;
192
 
  cout << setw(20) <<  "mdl_sw__exp__2 " <<  "= " <<
193
 
      setiosflags(ios::scientific) << setw(10) << mdl_sw__exp__2 << endl;
194
 
  cout << setw(20) <<  "mdl_cw__exp__2 " <<  "= " <<
195
 
      setiosflags(ios::scientific) << setw(10) << mdl_cw__exp__2 << endl;
196
 
}
197
 
void Parameters_sm::printIndependentCouplings()
198
 
{
199
 
  cout <<  "sm model couplings independent of event kinematics:" << endl; 
200
 
 
201
 
}
202
 
void Parameters_sm::printDependentParameters()
203
 
{
204
 
  cout <<  "sm model parameters dependent on event kinematics:" << endl; 
205
 
  cout << setw(20) <<  "mdl_sqrt__aS " <<  "= " << setiosflags(ios::scientific)
206
 
      << setw(10) << mdl_sqrt__aS << endl;
207
 
  cout << setw(20) <<  "G " <<  "= " << setiosflags(ios::scientific) <<
208
 
      setw(10) << G << endl;
209
 
  cout << setw(20) <<  "mdl_G__exp__2 " <<  "= " <<
210
 
      setiosflags(ios::scientific) << setw(10) << mdl_G__exp__2 << endl;
211
 
}
212
 
void Parameters_sm::printDependentCouplings()
213
 
{
214
 
  cout <<  "sm model couplings dependent on event kinematics:" << endl; 
215
 
  cout << setw(20) <<  "GC_11 " <<  "= " << setiosflags(ios::scientific) <<
216
 
      setw(10) << GC_11 << endl;
217
 
}
218
 
 
219