gambit is hosted by Hepforge, IPPP Durham
GAMBIT  v1.5.0-2191-ga4742ac
a Global And Modular Bsm Inference Tool
ColliderBit_LEP.cpp
Go to the documentation of this file.
1 // GAMBIT: Global and Modular BSM Inference Tool
2 // *********************************************
25 
26 #include <cmath>
27 #include <string>
28 #include <iostream>
29 #include <fstream>
30 #include <memory>
31 #include <numeric>
32 #include <sstream>
33 #include <vector>
34 
40 
41 //#define COLLIDERBIT_DEBUG
42 
43 namespace Gambit
44 {
45 
46  namespace ColliderBit
47  {
48 
49 
50  // *** Limits from e+e- colliders ***
51 
52 
53 
55  double limit_LLike(double x, double x95, double sigma) {
63  static double p95 = 1.;
64  using std::erf;
65  using std::sqrt;
66 
67  if (p95 < 1.01)
68  {
69  for (int i=0; i<20000; i++)
70  {
71  static double step = 0.1;
72  if (0.5 * (1 - erf(p95 + step)) > 0.05) p95 += step;
73  else step /= 10.;
74  }
75  }
76 
77  double result = 0.5 * (1.0 - erf(p95 + (x - x95) / sigma / sqrt(2.)));
78 
79  if (result < 0.0 or Utils::isnan(result))
80  {
81  cout << "result: " << result << endl;
82  cout << "x: " << x << endl;
83  cout << "x95: " << x95 << endl;
84  cout << "sigma: " << sigma << endl;
85  cout << "p95: " << p95 << endl;
86  cout << "(x - x95) / sigma / sqrt(2.): " << (x - x95) / sigma / sqrt(2.) << endl;
87  cout << "erf(p95 + (x - x95) / sigma / sqrt(2.)): " << erf(p95 + (x - x95) / sigma / sqrt(2.)) << endl;
88  ColliderBit_error().raise(LOCAL_INFO, "Suspicious results in limit_LLike!");
89  }
90 
91  return (result == 0.0 ? -1e10 : log(result));
92  }
93 
94 
96  bool is_xsec_sane(const triplet<double>& xsecWithError)
97  {
98  double xsec = xsecWithError.central;
99  double dxsec_upper = xsecWithError.upper - xsecWithError.central;
100  double dxsec_lower = xsecWithError.central - xsecWithError.lower;
101  if (xsec < 0.0 or dxsec_upper < 0.0 or dxsec_lower < 0.0
102  or Utils::isnan(xsec) or Utils::isnan(dxsec_upper) or Utils::isnan(dxsec_lower))
103  {
104  cout << "xsec: " << xsec << endl;
105  cout << "dxsec_upper: " << dxsec_upper << endl;
106  cout << "dxsec_lower: " << dxsec_lower << endl;
107  return false;
108  }
109  return true;
110  }
111 
112 
113 
117  {
119  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
120  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
121  get_sigma_ee_ll(result, 208.0, 1, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
122  if (!is_xsec_sane(result))
123  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
124  }
126  {
128  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
129  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
130  get_sigma_ee_ll(result, 208.0, 1, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
131  if (!is_xsec_sane(result))
132  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
133  }
135  {
137  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
138  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
139  get_sigma_ee_ll(result, 208.0, 1, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
140  if (!is_xsec_sane(result))
141  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
142  }
144  {
145  result = *Pipes::LEP208_SLHA1_convention_xsec_serselbar::Dep::LEP208_xsec_selserbar;
146  if (!is_xsec_sane(result))
147  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
148  }
150  {
152  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
153  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
154  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
155  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
156  get_sigma_ee_ll(result, 208.0, 1, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
157  if (!is_xsec_sane(result))
158  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
159  }
161  {
163  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
164  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
165  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
166  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
167  get_sigma_ee_ll(result, 208.0, 1, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
168  if (!is_xsec_sane(result))
169  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
170  }
172  {
174  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
175  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
176  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
177  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
178  get_sigma_ee_ll(result, 208.0, 1, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
179  if (!is_xsec_sane(result))
180  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
181  }
183  {
184  result = *Pipes::LEP208_SLHA1_convention_xsec_se2se1bar::Dep::LEP208_xsec_se1se2bar;
185  if (!is_xsec_sane(result))
186  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
187  }
189 
193  {
195  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
196  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
197  get_sigma_ee_ll(result, 208.0, 2, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
198  if (!is_xsec_sane(result))
199  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
200  }
202  {
204  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
205  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
206  get_sigma_ee_ll(result, 208.0, 2, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
207  if (!is_xsec_sane(result))
208  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
209  }
211  {
213  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
214  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
215  get_sigma_ee_ll(result, 208.0, 2, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
216  if (!is_xsec_sane(result))
217  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
218  }
220  {
221  result = *Pipes::LEP208_SLHA1_convention_xsec_smursmulbar::Dep::LEP208_xsec_smulsmurbar;
222  if (!is_xsec_sane(result))
223  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
224  }
226  {
228  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
229  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
230  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
231  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
232  get_sigma_ee_ll(result, 208.0, 2, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
233  if (!is_xsec_sane(result))
234  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
235  }
237  {
239  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
240  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
241  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
242  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
243  get_sigma_ee_ll(result, 208.0, 2, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
244  if (!is_xsec_sane(result))
245  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
246  }
248  {
250  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
251  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
252  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
253  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
254  get_sigma_ee_ll(result, 208.0, 2, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
255  if (!is_xsec_sane(result))
256  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
257  }
259  {
260  result = *Pipes::LEP208_SLHA1_convention_xsec_smu2smu1bar::Dep::LEP208_xsec_smu1smu2bar;
261  if (!is_xsec_sane(result))
262  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
263  }
265 
269  {
271  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
272  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
273  get_sigma_ee_ll(result, 208.0, 3, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
274  if (!is_xsec_sane(result))
275  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
276  }
278  {
280  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
281  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
282  get_sigma_ee_ll(result, 208.0, 3, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
283  if (!is_xsec_sane(result))
284  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
285  }
287  {
289  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
290  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
291  get_sigma_ee_ll(result, 208.0, 3, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
292  if (!is_xsec_sane(result))
293  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
294  }
296  {
297  result = *Pipes::LEP208_SLHA1_convention_xsec_staurstaulbar::Dep::LEP208_xsec_staulstaurbar;
298  if (!is_xsec_sane(result))
299  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
300  }
302  {
304  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
305  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
306  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
307  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
308  get_sigma_ee_ll(result, 208.0, 3, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
309  if (!is_xsec_sane(result))
310  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
311  }
313  {
315  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
316  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
317  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
318  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
319  get_sigma_ee_ll(result, 208.0, 3, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
320  if (!is_xsec_sane(result))
321  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
322  }
324  {
326  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
327  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
328  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
329  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
330  get_sigma_ee_ll(result, 208.0, 3, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
331  if (!is_xsec_sane(result))
332  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
333  }
335  {
336  result = *Pipes::LEP208_SLHA1_convention_xsec_stau2stau1bar::Dep::LEP208_xsec_stau1stau2bar;
337  if (!is_xsec_sane(result))
338  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
339  }
341 
342  // @brief \f[ee \to \chi_1\chi_1\f] pair production cross-section at 207 GeV
345  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
346  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
347  get_sigma_ee_chi00(result, 207.0, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
348  if (!is_xsec_sane(result)) {
349  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
350  }
351  }
352 
356  {
358  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
359  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
360  get_sigma_ee_chi00(result, 208.0, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
361  if (!is_xsec_sane(result))
362  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
363  }
365  {
367  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
368  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
369  get_sigma_ee_chi00(result, 208.0, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
370  if (!is_xsec_sane(result))
371  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
372  }
374  {
376  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
377  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
378  get_sigma_ee_chi00(result, 208.0, 1, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
379  if (!is_xsec_sane(result))
380  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
381  }
383  {
385  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
386  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
387  get_sigma_ee_chi00(result, 208.0, 1, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
388  if (!is_xsec_sane(result))
389  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
390  }
392  {
394  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
395  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
396  get_sigma_ee_chi00(result, 208.0, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
397  if (!is_xsec_sane(result))
398  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
399  }
401  {
403  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
404  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
405  get_sigma_ee_chi00(result, 208.0, 2, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
406  if (!is_xsec_sane(result))
407  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
408  }
410  {
412  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
413  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
414  get_sigma_ee_chi00(result, 208.0, 2, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
415  if (!is_xsec_sane(result))
416  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
417  }
419  {
421  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
422  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
423  get_sigma_ee_chi00(result, 208.0, 3, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
424  if (!is_xsec_sane(result))
425  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
426  }
428  {
430  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
431  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
432  get_sigma_ee_chi00(result, 208.0, 3, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
433  if (!is_xsec_sane(result))
434  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
435  }
437  {
439  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
440  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
441  get_sigma_ee_chi00(result, 208.0, 4, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
442  if (!is_xsec_sane(result))
443  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
444  }
446 
447 
451  {
453  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
454  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
455  get_sigma_ee_chipm(result, 208.0, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
456  if (!is_xsec_sane(result))
457  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
458  }
460  {
462  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
463  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
464  get_sigma_ee_chipm(result, 208.0, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
465  if (!is_xsec_sane(result))
466  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
467  }
469  {
471  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
472  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
473  get_sigma_ee_chipm(result, 208.0, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
474  if (!is_xsec_sane(result))
475  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
476  }
478  {
479  result = *Pipes::LEP208_SLHA1_convention_xsec_chipm_21::Dep::LEP208_xsec_chipm_12;
480  if (!is_xsec_sane(result))
481  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
482  }
484 
485 
489  {
491  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
492  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
493  get_sigma_ee_ll(result, 205.0, 1, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
494  if (!is_xsec_sane(result))
495  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
496  }
498  {
500  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
501  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
502  get_sigma_ee_ll(result, 205.0, 1, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
503  if (!is_xsec_sane(result))
504  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
505  }
507  {
509  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
510  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
511  get_sigma_ee_ll(result, 205.0, 1, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
512  if (!is_xsec_sane(result))
513  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
514  }
516  {
517  result = *Pipes::LEP205_SLHA1_convention_xsec_serselbar::Dep::LEP205_xsec_selserbar;
518  if (!is_xsec_sane(result))
519  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
520  }
522  {
524  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
525  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
526  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
527  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
528  get_sigma_ee_ll(result, 205.0, 1, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
529  if (!is_xsec_sane(result))
530  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
531  }
533  {
535  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
536  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
537  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
538  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
539  get_sigma_ee_ll(result, 205.0, 1, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
540  if (!is_xsec_sane(result))
541  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
542  }
544  {
546  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
547  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
548  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
549  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
550  get_sigma_ee_ll(result, 205.0, 1, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
551  if (!is_xsec_sane(result))
552  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
553  }
555  {
556  result = *Pipes::LEP205_SLHA1_convention_xsec_se2se1bar::Dep::LEP205_xsec_se1se2bar;
557  if (!is_xsec_sane(result))
558  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
559  }
561 
565  {
567  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
568  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
569  get_sigma_ee_ll(result, 205.0, 2, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
570  if (!is_xsec_sane(result))
571  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
572  }
574  {
576  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
577  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
578  get_sigma_ee_ll(result, 205.0, 2, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
579  if (!is_xsec_sane(result))
580  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
581  }
583  {
585  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
586  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
587  get_sigma_ee_ll(result, 205.0, 2, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
588  if (!is_xsec_sane(result))
589  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
590  }
592  {
593  result = *Pipes::LEP205_SLHA1_convention_xsec_smursmulbar::Dep::LEP205_xsec_smulsmurbar;
594  if (!is_xsec_sane(result))
595  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
596  }
598  {
600  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
601  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
602  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
603  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
604  get_sigma_ee_ll(result, 205.0, 2, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
605  if (!is_xsec_sane(result))
606  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
607  }
609  {
611  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
612  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
613  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
614  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
615  get_sigma_ee_ll(result, 205.0, 2, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
616  if (!is_xsec_sane(result))
617  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
618  }
620  {
622  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
623  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
624  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
625  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
626  get_sigma_ee_ll(result, 205.0, 2, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
627  if (!is_xsec_sane(result))
628  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
629  }
631  {
632  result = *Pipes::LEP205_SLHA1_convention_xsec_smu2smu1bar::Dep::LEP205_xsec_smu1smu2bar;
633  if (!is_xsec_sane(result))
634  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
635  }
637 
641  {
643  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
644  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
645  get_sigma_ee_ll(result, 205.0, 3, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
646  if (!is_xsec_sane(result))
647  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
648  }
650  {
652  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
653  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
654  get_sigma_ee_ll(result, 205.0, 3, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
655  if (!is_xsec_sane(result))
656  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
657  }
659  {
661  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
662  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
663  get_sigma_ee_ll(result, 205.0, 3, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
664  if (!is_xsec_sane(result))
665  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
666  }
668  {
669  result = *Pipes::LEP205_SLHA1_convention_xsec_staurstaulbar::Dep::LEP205_xsec_staulstaurbar;
670  if (!is_xsec_sane(result))
671  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
672  }
674  {
676  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
677  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
678  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
679  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
680  get_sigma_ee_ll(result, 205.0, 3, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
681  if (!is_xsec_sane(result))
682  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
683  }
685  {
687  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
688  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
689  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
690  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
691  get_sigma_ee_ll(result, 205.0, 3, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
692  if (!is_xsec_sane(result))
693  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
694  }
696  {
698  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
699  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
700  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
701  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
702  get_sigma_ee_ll(result, 205.0, 3, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
703  if (!is_xsec_sane(result))
704  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
705  }
707  {
708  result = *Pipes::LEP205_SLHA1_convention_xsec_stau2stau1bar::Dep::LEP205_xsec_stau1stau2bar;
709  if (!is_xsec_sane(result))
710  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
711  }
713 
714 
718  {
720  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
721  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
722  get_sigma_ee_chi00(result, 205.0, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
723  if (!is_xsec_sane(result))
724  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
725  }
727  {
729  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
730  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
731  get_sigma_ee_chi00(result, 205.0, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
732  if (!is_xsec_sane(result))
733  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
734  }
736  {
738  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
739  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
740  get_sigma_ee_chi00(result, 205.0, 1, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
741  if (!is_xsec_sane(result))
742  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
743  }
745  {
747  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
748  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
749  get_sigma_ee_chi00(result, 205.0, 1, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
750  if (!is_xsec_sane(result))
751  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
752  }
754  {
756  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
757  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
758  get_sigma_ee_chi00(result, 205.0, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
759  if (!is_xsec_sane(result))
760  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
761  }
763  {
765  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
766  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
767  get_sigma_ee_chi00(result, 205.0, 2, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
768  if (!is_xsec_sane(result))
769  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
770  }
772  {
774  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
775  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
776  get_sigma_ee_chi00(result, 205.0, 2, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
777  if (!is_xsec_sane(result))
778  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
779  }
781  {
783  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
784  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
785  get_sigma_ee_chi00(result, 205.0, 3, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
786  if (!is_xsec_sane(result))
787  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
788  }
790  {
792  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
793  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
794  get_sigma_ee_chi00(result, 205.0, 3, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
795  if (!is_xsec_sane(result))
796  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
797  }
799  {
801  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
802  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
803  get_sigma_ee_chi00(result, 205.0, 4, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
804  if (!is_xsec_sane(result))
805  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
806  }
808 
809 
813  {
815  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
816  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
817  get_sigma_ee_chipm(result, 205.0, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
818  if (!is_xsec_sane(result))
819  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
820  }
822  {
824  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
825  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
826  get_sigma_ee_chipm(result, 205.0, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
827  if (!is_xsec_sane(result))
828  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
829  }
831  {
833  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
834  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
835  get_sigma_ee_chipm(result, 205.0, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
836  if (!is_xsec_sane(result))
837  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
838  }
840  {
841  result = *Pipes::LEP205_SLHA1_convention_xsec_chipm_21::Dep::LEP205_xsec_chipm_12;
842  if (!is_xsec_sane(result))
843  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
844  }
845 
849  {
851  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
852  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
853  get_sigma_ee_ll(result, 188.6, 1, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
854  if (!is_xsec_sane(result))
855  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
856  }
858  {
860  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
861  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
862  get_sigma_ee_ll(result, 188.6, 1, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
863  if (!is_xsec_sane(result))
864  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
865  }
867  {
869  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
870  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
871  get_sigma_ee_ll(result, 188.6, 1, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
872  if (!is_xsec_sane(result))
873  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
874  }
876  {
877  result = *Pipes::LEP188_SLHA1_convention_xsec_serselbar::Dep::LEP188_xsec_selserbar;
878  if (!is_xsec_sane(result))
879  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
880  }
882  {
884  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
885  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
886  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
887  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
888  get_sigma_ee_ll(result, 188.6, 1, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
889  if (!is_xsec_sane(result))
890  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
891  }
893  {
895  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
896  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
897  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
898  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
899  get_sigma_ee_ll(result, 188.6, 1, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
900  if (!is_xsec_sane(result))
901  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
902  }
904  {
906  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
907  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
908  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
909  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
910  get_sigma_ee_ll(result, 188.6, 1, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
911  if (!is_xsec_sane(result))
912  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
913  }
915  {
916  result = *Pipes::LEP188_SLHA1_convention_xsec_se2se1bar::Dep::LEP188_xsec_se1se2bar;
917  if (!is_xsec_sane(result))
918  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
919  }
921 
925  {
927  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
928  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
929  get_sigma_ee_ll(result, 188.6, 2, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
930  if (!is_xsec_sane(result))
931  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
932  }
934  {
936  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
937  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
938  get_sigma_ee_ll(result, 188.6, 2, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
939  if (!is_xsec_sane(result))
940  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
941  }
943  {
945  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
946  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
947  get_sigma_ee_ll(result, 188.6, 2, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
948  if (!is_xsec_sane(result))
949  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
950  }
952  {
953  result = *Pipes::LEP188_SLHA1_convention_xsec_smursmulbar::Dep::LEP188_xsec_smulsmurbar;
954  if (!is_xsec_sane(result))
955  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
956  }
958  {
960  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
961  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
962  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
963  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
964  get_sigma_ee_ll(result, 188.6, 2, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
965  if (!is_xsec_sane(result))
966  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
967  }
969  {
971  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
972  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
973  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
974  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
975  get_sigma_ee_ll(result, 188.6, 2, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
976  if (!is_xsec_sane(result))
977  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
978  }
980  {
982  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
983  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
984  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
985  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
986  get_sigma_ee_ll(result, 188.6, 2, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
987  if (!is_xsec_sane(result))
988  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
989  }
991  {
992  result = *Pipes::LEP188_SLHA1_convention_xsec_smu2smu1bar::Dep::LEP188_xsec_smu1smu2bar;
993  if (!is_xsec_sane(result))
994  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
995  }
997 
1001  {
1003  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1004  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1005  get_sigma_ee_ll(result, 188.6, 3, 1, 1, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
1006  if (!is_xsec_sane(result))
1007  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1008  }
1010  {
1012  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1013  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1014  get_sigma_ee_ll(result, 188.6, 3, 1, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
1015  if (!is_xsec_sane(result))
1016  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1017  }
1019  {
1021  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1022  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1023  get_sigma_ee_ll(result, 188.6, 3, 2, 2, tol, tol, pt_error, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, true);
1024  if (!is_xsec_sane(result))
1025  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1026  }
1028  {
1029  result = *Pipes::LEP188_SLHA1_convention_xsec_staurstaulbar::Dep::LEP188_xsec_staulstaurbar;
1030  if (!is_xsec_sane(result))
1031  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1032  }
1034  {
1036  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1037  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1038  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
1039  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
1040  get_sigma_ee_ll(result, 188.6, 3, 1, 1, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
1041  if (!is_xsec_sane(result))
1042  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1043  }
1045  {
1047  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1048  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1049  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
1050  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
1051  get_sigma_ee_ll(result, 188.6, 3, 1, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
1052  if (!is_xsec_sane(result))
1053  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1054  }
1056  {
1058  const static double gtol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1059  const static bool gpt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1060  const static double ftol = runOptions->getValueOrDef<double>(1e-2, "family_mixing_tolerance");
1061  const static bool fpt_error = runOptions->getValueOrDef<bool>(true, "family_mixing_tolerance_invalidates_point_only");
1062  get_sigma_ee_ll(result, 188.6, 3, 2, 2, gtol, ftol, gpt_error, fpt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV, false);
1063  if (!is_xsec_sane(result))
1064  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1065  }
1067  {
1068  result = *Pipes::LEP188_SLHA1_convention_xsec_stau2stau1bar::Dep::LEP188_xsec_stau1stau2bar;
1069  if (!is_xsec_sane(result))
1070  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1071  }
1073 
1074 
1078  {
1080  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1081  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1082  get_sigma_ee_chi00(result, 188.6, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1083  if (!is_xsec_sane(result))
1084  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1085  }
1087  {
1089  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1090  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1091  get_sigma_ee_chi00(result, 188.6, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1092  if (!is_xsec_sane(result))
1093  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1094  }
1096  {
1098  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1099  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1100  get_sigma_ee_chi00(result, 188.6, 1, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1101  if (!is_xsec_sane(result))
1102  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1103  }
1105  {
1107  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1108  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1109  get_sigma_ee_chi00(result, 188.6, 1, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1110  if (!is_xsec_sane(result))
1111  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1112  }
1114  {
1116  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1117  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1118  get_sigma_ee_chi00(result, 188.6, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1119  if (!is_xsec_sane(result))
1120  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1121  }
1123  {
1125  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1126  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1127  get_sigma_ee_chi00(result, 188.6, 2, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1128  if (!is_xsec_sane(result))
1129  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1130  }
1132  {
1134  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1135  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1136  get_sigma_ee_chi00(result, 188.6, 2, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1137  if (!is_xsec_sane(result))
1138  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1139  }
1141  {
1143  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1144  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1145  get_sigma_ee_chi00(result, 188.6, 3, 3, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1146  if (!is_xsec_sane(result))
1147  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1148  }
1150  {
1152  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1153  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1154  get_sigma_ee_chi00(result, 188.6, 3, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1155  if (!is_xsec_sane(result))
1156  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1157  }
1159  {
1161  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1162  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1163  get_sigma_ee_chi00(result, 188.6, 4, 4, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1164  if (!is_xsec_sane(result))
1165  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1166  }
1168 
1169 
1173  {
1175  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1176  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1177  get_sigma_ee_chipm(result, 188.6, 1, 1, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1178  if (!is_xsec_sane(result))
1179  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1180  }
1182  {
1184  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1185  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1186  get_sigma_ee_chipm(result, 188.6, 1, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1187  if (!is_xsec_sane(result))
1188  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1189  }
1191  {
1193  const static double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
1194  const static bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
1195  get_sigma_ee_chipm(result, 188.6, 2, 2, tol, pt_error, *Dep::MSSM_spectrum, Dep::Z_decay_rates->width_in_GeV);
1196  if (!is_xsec_sane(result))
1197  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1198  }
1200  {
1201  result = *Pipes::LEP188_SLHA1_convention_xsec_chipm_21::Dep::LEP188_xsec_chipm_12;
1202  if (!is_xsec_sane(result))
1203  ColliderBit_error().raise(LOCAL_INFO, "Non-physical LEP cross section!");
1204  }
1206 
1207 
1211  {
1213  using std::pow;
1214  using std::log;
1215 
1216  const Spectrum& spec = *Dep::MSSM_spectrum;
1217 
1218  double max_mixing;
1219  const SubSpectrum& mssm = spec.get_HE();
1220  str sel_string = slhahelp::mass_es_from_gauge_es("~e_L", max_mixing, mssm);
1221  str ser_string = slhahelp::mass_es_from_gauge_es("~e_R", max_mixing, mssm);
1222  const double mass_seL=spec.get(Par::Pole_Mass,sel_string);
1223  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1224  const double mass_seR = spec.get(Par::Pole_Mass,ser_string);
1225  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1226  triplet<double> xsecWithError;
1227  double xsecLimit;
1228 
1229  static const ALEPHSelectronLimitAt208GeV limitContainer;
1230  // #ifdef COLLIDERBIT_DEBUG
1231  //static bool dumped=false;
1232  //if(!dumped)
1233  //{
1234  // limitContainer.dumpPlotData(45., 115., 0., 100., mZ, "lepLimitPlanev2/ALEPHSelectronLimitAt208GeV.dump");
1235  // dumped=true;
1236  //}
1237  // #endif
1238 
1239  result = 0;
1240  // Due to the nature of the analysis details of the model independent limit in
1241  // the paper, the best we can do is to try these two processes individually:
1242 
1243  // se_L, se_L
1244  xsecLimit = limitContainer.limitAverage(mass_seL, mass_neut1, mZ);
1245  xsecWithError = *Dep::LEP208_xsec_selselbar;
1246  xsecWithError.upper *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1247  xsecWithError.central *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1248  xsecWithError.lower *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1249 
1250  if (xsecWithError.central < xsecLimit)
1251  {
1252  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1253  }
1254  else
1255  {
1256  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1257  }
1258 
1259  // se_R, se_R
1260  xsecLimit = limitContainer.limitAverage(mass_seR, mass_neut1, mZ);
1261 
1262  xsecWithError = *Dep::LEP208_xsec_serserbar;
1263  xsecWithError.upper *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1264  xsecWithError.central *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1265  xsecWithError.lower *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1266 
1267  if (xsecWithError.central < xsecLimit)
1268  {
1269  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1270  }
1271  else
1272  {
1273  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1274  }
1275 
1276  }
1277 
1278  void ALEPH_Smuon_Conservative_LLike(double& result)
1279  {
1280  using namespace Pipes::ALEPH_Smuon_Conservative_LLike;
1281  using std::pow;
1282  using std::log;
1283 
1284  const Spectrum& spec = *Dep::MSSM_spectrum;
1285 
1286  double max_mixing;
1287  const SubSpectrum& mssm = spec.get_HE();
1288  str smul_string = slhahelp::mass_es_from_gauge_es("~mu_L", max_mixing, mssm);
1289  str smur_string = slhahelp::mass_es_from_gauge_es("~mu_R", max_mixing, mssm);
1290  const double mass_smuL=spec.get(Par::Pole_Mass,smul_string);
1291  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1292  const double mass_smuR = spec.get(Par::Pole_Mass,smur_string);
1293  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1294  triplet<double> xsecWithError;
1295  double xsecLimit;
1296 
1297  static const ALEPHSmuonLimitAt208GeV limitContainer;
1298  // #ifdef COLLIDERBIT_DEBUG
1299  // static bool dumped=false;
1300  // if(!dumped)
1301  // {
1302  // limitContainer.dumpPlotData(45., 115., 0., 100., mZ, "lepLimitPlanev2/ALEPHSmuonLimitAt208GeV.dump");
1303  // dumped=true;
1304  // }
1305  // #endif
1306 
1307  result = 0;
1308  // Due to the nature of the analysis details of the model independent limit in
1309  // the paper, the best we can do is to try these two processes individually:
1310 
1311  // smu_L, smu_L
1312  xsecLimit = limitContainer.limitAverage(mass_smuL, mass_neut1, mZ);
1313  xsecWithError = *Dep::LEP208_xsec_smulsmulbar;
1314  xsecWithError.upper *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1315  xsecWithError.central *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1316  xsecWithError.lower *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1317 
1318  if (xsecWithError.central < xsecLimit)
1319  {
1320  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1321  }
1322  else
1323  {
1324  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1325  }
1326 
1327  // smu_R, smu_R
1328  xsecLimit = limitContainer.limitAverage(mass_smuR, mass_neut1, mZ);
1329 
1330  xsecWithError = *Dep::LEP208_xsec_smursmurbar;
1331  xsecWithError.upper *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1332  xsecWithError.central *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1333  xsecWithError.lower *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1334 
1335  if (xsecWithError.central < xsecLimit)
1336  {
1337  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1338  }
1339  else
1340  {
1341  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1342  }
1343 
1344  }
1345 
1346  void ALEPH_Stau_Conservative_LLike(double& result)
1347  {
1348  using namespace Pipes::ALEPH_Stau_Conservative_LLike;
1349  using std::pow;
1350  using std::log;
1351 
1352  const Spectrum& spec = *Dep::MSSM_spectrum;
1353 
1354  const SubSpectrum& mssm = spec.get_HE();
1355  static const double tol = runOptions->getValueOrDef<double>(1e-5, "family_mixing_tolerance");
1356  static const bool pterror = runOptions->getValueOrDef<bool>(false, "family_mixing_tolerance_invalidates_point_only");
1357  str stau1_string = slhahelp::mass_es_closest_to_family("~tau_1", mssm,tol,LOCAL_INFO,pterror);
1358  str stau2_string = slhahelp::mass_es_closest_to_family("~tau_2", mssm,tol,LOCAL_INFO,pterror);
1359  const double mass_stau1=spec.get(Par::Pole_Mass,stau1_string);
1360  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1361  const double mass_stau2 = spec.get(Par::Pole_Mass,stau2_string);
1362  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1363  triplet<double> xsecWithError;
1364  double xsecLimit;
1365 
1366  static const ALEPHStauLimitAt208GeV limitContainer;
1367  // #ifdef COLLIDERBIT_DEBUG
1368  // static bool dumped=false;
1369  // if(!dumped)
1370  // {
1371  // limitContainer.dumpPlotData(45., 115., 0., 100., mZ, "lepLimitPlanev2/ALEPHStauLimitAt208GeV.dump");
1372  // dumped=true;
1373  // }
1374  // #endif
1375 
1376  result = 0;
1377  // Due to the nature of the analysis details of the model independent limit in
1378  // the paper, the best we can do is to try these two processes individually:
1379 
1380  // stau_1, stau_1
1381  xsecLimit = limitContainer.limitAverage(mass_stau1, mass_neut1, mZ);
1382 
1383  xsecWithError = *Dep::LEP208_xsec_stau1stau1bar;
1384  xsecWithError.upper *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1385  xsecWithError.central *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1386  xsecWithError.lower *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1387 
1388  if (xsecWithError.central < xsecLimit)
1389  {
1390  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1391  }
1392  else
1393  {
1394  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1395  }
1396 
1397  // stau_2, stau_2
1398  xsecLimit = limitContainer.limitAverage(mass_stau2, mass_neut1, mZ);
1399 
1400  xsecWithError = *Dep::LEP208_xsec_stau2stau2bar;
1401  xsecWithError.upper *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1402  xsecWithError.central *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1403  xsecWithError.lower *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1404 
1405  if (xsecWithError.central < xsecLimit)
1406  {
1407  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1408  }
1409  else
1410  {
1411  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1412  }
1413 
1414  }
1415 
1416  void L3_Selectron_Conservative_LLike(double& result)
1417  {
1419  using std::pow;
1420  using std::log;
1421 
1422  const Spectrum& spec = *Dep::MSSM_spectrum;
1423 
1424  double max_mixing;
1425  const SubSpectrum& mssm = spec.get_HE();
1426  str sel_string = slhahelp::mass_es_from_gauge_es("~e_L", max_mixing, mssm);
1427  str ser_string = slhahelp::mass_es_from_gauge_es("~e_R", max_mixing, mssm);
1428  const double mass_seL=spec.get(Par::Pole_Mass,sel_string);
1429  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1430  const double mass_seR = spec.get(Par::Pole_Mass,ser_string);
1431  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1432  triplet<double> xsecWithError;
1433  double xsecLimit;
1434 
1435  static const L3SelectronLimitAt205GeV limitContainer;
1436  // #ifdef COLLIDERBIT_DEBUG
1437  static bool dumped=false;
1438  if(!dumped)
1439  {
1440  limitContainer.dumpPlotData(45., 104., 0., 100., mZ, "lepLimitPlanev2/L3SelectronLimitAt205GeV.dump",200);
1441  dumped=true;
1442  }
1443  // #endif
1444 
1445  result = 0;
1446  // Due to the nature of the analysis details of the model independent limit in
1447  // the paper, the best we can do is to try these two processes individually:
1448 
1449  // se_L, se_L
1450  xsecLimit = limitContainer.limitAverage(mass_seL, mass_neut1, mZ);
1451 
1452  xsecWithError = *Dep::LEP205_xsec_selselbar;
1453  xsecWithError.upper *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1454  xsecWithError.central *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1455  xsecWithError.lower *= pow(Dep::selectron_l_decay_rates->BF("~chi0_1", "e-"), 2);
1456 
1457  if (xsecWithError.central < xsecLimit)
1458  {
1459  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1460  }
1461  else
1462  {
1463  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1464  }
1465 
1466  // se_R, se_R
1467  xsecLimit = limitContainer.limitAverage(mass_seR, mass_neut1, mZ);
1468 
1469  xsecWithError = *Dep::LEP205_xsec_serserbar;
1470  xsecWithError.upper *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1471  xsecWithError.central *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1472  xsecWithError.lower *= pow(Dep::selectron_r_decay_rates->BF("~chi0_1", "e-"), 2);
1473 
1474  if (xsecWithError.central < xsecLimit)
1475  {
1476  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1477  }
1478  else
1479  {
1480  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1481  }
1482 
1483  }
1484 
1485  void L3_Smuon_Conservative_LLike(double& result)
1486  {
1487  using namespace Pipes::L3_Smuon_Conservative_LLike;
1488  using std::pow;
1489  using std::log;
1490 
1491  const Spectrum& spec = *Dep::MSSM_spectrum;
1492 
1493  double max_mixing;
1494  const SubSpectrum& mssm = spec.get_HE();
1495  str smul_string = slhahelp::mass_es_from_gauge_es("~mu_L", max_mixing, mssm);
1496  str smur_string = slhahelp::mass_es_from_gauge_es("~mu_R", max_mixing, mssm);
1497  const double mass_smuL=spec.get(Par::Pole_Mass,smul_string);
1498  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1499  const double mass_smuR = spec.get(Par::Pole_Mass,smur_string);
1500  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1501  triplet<double> xsecWithError;
1502  double xsecLimit;
1503 
1504  static const L3SmuonLimitAt205GeV limitContainer;
1505  // #ifdef COLLIDERBIT_DEBUG
1506  // static bool dumped=false;
1507  // if(!dumped)
1508  // {
1509  // limitContainer.dumpPlotData(45., 115., 0., 100., mZ, "lepLimitPlanev2/L3SmuonLimitAt205GeV.dump");
1510  // dumped=true;
1511  // }
1512  // #endif
1513 
1514  result = 0;
1515  // Due to the nature of the analysis details of the model independent limit in
1516  // the paper, the best we can do is to try these two processes individually:
1517 
1518  // smu_L, smu_L
1519  xsecLimit = limitContainer.limitAverage(mass_smuL, mass_neut1, mZ);
1520 
1521  xsecWithError = *Dep::LEP205_xsec_smulsmulbar;
1522  xsecWithError.upper *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1523  xsecWithError.central *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1524  xsecWithError.lower *= pow(Dep::smuon_l_decay_rates->BF("~chi0_1", "mu-"), 2);
1525 
1526  if (xsecWithError.central < xsecLimit)
1527  {
1528  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1529  }
1530  else
1531  {
1532  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1533  }
1534 
1535  // smu_R, smu_R
1536  xsecLimit = limitContainer.limitAverage(mass_smuR, mass_neut1, mZ);
1537 
1538  xsecWithError = *Dep::LEP205_xsec_smursmurbar;
1539  xsecWithError.upper *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1540  xsecWithError.central *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1541  xsecWithError.lower *= pow(Dep::smuon_r_decay_rates->BF("~chi0_1", "mu-"), 2);
1542 
1543  if (xsecWithError.central < xsecLimit)
1544  {
1545  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1546  }
1547  else
1548  {
1549  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1550  }
1551 
1552  }
1553 
1554  void L3_Stau_Conservative_LLike(double& result)
1555  {
1556  using namespace Pipes::L3_Stau_Conservative_LLike;
1557  using std::pow;
1558  using std::log;
1559 
1560  const Spectrum& spec = *Dep::MSSM_spectrum;
1561 
1562  const SubSpectrum& mssm = spec.get_HE();
1563  static const double tol = runOptions->getValueOrDef<double>(1e-5, "family_mixing_tolerance");
1564  static const bool pterror = runOptions->getValueOrDef<bool>(false, "family_mixing_tolerance_invalidates_point_only");
1565  str stau1_string = slhahelp::mass_es_closest_to_family("~tau_1", mssm,tol,LOCAL_INFO,pterror);
1566  str stau2_string = slhahelp::mass_es_closest_to_family("~tau_2", mssm,tol,LOCAL_INFO,pterror);
1567  const double mass_stau1=spec.get(Par::Pole_Mass,stau1_string);
1568  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1569  const double mass_stau2 = spec.get(Par::Pole_Mass,stau2_string);
1570  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1571  triplet<double> xsecWithError;
1572  double xsecLimit;
1573 
1574  static const L3StauLimitAt205GeV limitContainer;
1575  // #ifdef COLLIDERBIT_DEBUG
1576  // static bool dumped=false;
1577  // if(!dumped)
1578  // {
1579  // limitContainer.dumpPlotData(45., 115., 0., 100., mZ, "lepLimitPlanev2/L3StauLimitAt205GeV.dump");
1580  // dumped=true;
1581  // }
1582  // #endif
1583 
1584  result = 0;
1585  // Due to the nature of the analysis details of the model independent limit in
1586  // the paper, the best we can do is to try these two processes individually:
1587 
1588  // stau_1, stau_1
1589  xsecLimit = limitContainer.limitAverage(mass_stau1, mass_neut1, mZ);
1590 
1591  xsecWithError = *Dep::LEP205_xsec_stau1stau1bar;
1592  xsecWithError.upper *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1593  xsecWithError.central *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1594  xsecWithError.lower *= pow(Dep::stau_1_decay_rates->BF("~chi0_1", "tau-"), 2);
1595 
1596  if (xsecWithError.central < xsecLimit)
1597  {
1598  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1599  }
1600  else
1601  {
1602  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1603  }
1604 
1605  // stau_2, stau_2
1606  xsecLimit = limitContainer.limitAverage(mass_stau2, mass_neut1, mZ);
1607 
1608  xsecWithError = *Dep::LEP205_xsec_stau2stau2bar;
1609  xsecWithError.upper *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1610  xsecWithError.central *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1611  xsecWithError.lower *= pow(Dep::stau_2_decay_rates->BF("~chi0_1", "tau-"), 2);
1612 
1613  if (xsecWithError.central < xsecLimit)
1614  {
1615  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1616  }
1617  else
1618  {
1619  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1620  }
1621 
1622  }
1624 
1628  {
1630  using std::pow;
1631  using std::log;
1632 
1633  const Spectrum& spec = *Dep::MSSM_spectrum;
1634 
1635  const DecayTable& decays = *Dep::decay_rates;
1636  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1637  const double mass_neut2 = spec.get(Par::Pole_Mass,1000023, 0);
1638  const double mass_neut3 = spec.get(Par::Pole_Mass,1000025, 0);
1639  const double mass_neut4 = spec.get(Par::Pole_Mass,1000035, 0);
1640  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1641  triplet<double> xsecWithError;
1642  double xsecLimit, totalBR;
1643 
1644  static const L3NeutralinoAllChannelsLimitAt188pt6GeV limitContainer;
1645  // #ifdef COLLIDERBIT_DEBUG
1646  // static bool dumped=false;
1647  // if(!dumped)
1648  // {
1649  // limitContainer.dumpPlotData(0., 200., 0., 100., mZ, "lepLimitPlanev2/L3NeutralinoAllChannelsLimitAt188pt6GeV.dump");
1650  // dumped=true;
1651  // }
1652  // #endif
1653 
1654  result = 0;
1655  // Due to the nature of the analysis details of the model independent limit in
1656  // the paper, the best we can do is to try these processes individually:
1657 
1658  // neut2, neut1
1659  xsecLimit = limitContainer.limitAverage(mass_neut2, mass_neut1, mZ);
1660 
1661  xsecWithError = *Dep::LEP188_xsec_chi00_12;
1662  // Total up all channels which look like Z* decays
1663  totalBR = 0;
1664  totalBR += decays.at("~chi0_2").BF("~chi0_1", "Z0");
1665  totalBR += decays.at("~chi0_2").BF("~chi0_1", "ubar", "u");
1666  totalBR += decays.at("~chi0_2").BF("~chi0_1", "dbar", "d");
1667  totalBR += decays.at("~chi0_2").BF("~chi0_1", "cbar", "c");
1668  totalBR += decays.at("~chi0_2").BF("~chi0_1", "sbar", "s");
1669  totalBR += decays.at("~chi0_2").BF("~chi0_1", "bbar", "b");
1670  totalBR += decays.at("~chi0_2").BF("~chi0_1", "e+", "e-");
1671  totalBR += decays.at("~chi0_2").BF("~chi0_1", "mu+", "mu-");
1672  totalBR += decays.at("~chi0_2").BF("~chi0_1", "tau+", "tau-");
1673  totalBR += decays.at("~chi0_2").BF("~chi0_1", "nubar_e", "nu_e");
1674  totalBR += decays.at("~chi0_2").BF("~chi0_1", "nubar_mu", "nu_mu");
1675  totalBR += decays.at("~chi0_2").BF("~chi0_1", "nubar_tau", "nu_tau");
1676  xsecWithError.upper *= totalBR;
1677  xsecWithError.central *= totalBR;
1678  xsecWithError.lower *= totalBR;
1679 
1680  if (xsecWithError.central < xsecLimit)
1681  {
1682  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1683  }
1684  else
1685  {
1686  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1687  }
1688 
1689  // neut3, neut1
1690  xsecLimit = limitContainer.limitAverage(mass_neut3, mass_neut1, mZ);
1691 
1692  xsecWithError = *Dep::LEP188_xsec_chi00_13;
1693  // Total up all channels which look like Z* decays
1694  totalBR = 0;
1695  totalBR += decays.at("~chi0_3").BF("~chi0_1", "Z0");
1696  totalBR += decays.at("~chi0_3").BF("~chi0_1", "ubar", "u");
1697  totalBR += decays.at("~chi0_3").BF("~chi0_1", "dbar", "d");
1698  totalBR += decays.at("~chi0_3").BF("~chi0_1", "cbar", "c");
1699  totalBR += decays.at("~chi0_3").BF("~chi0_1", "sbar", "s");
1700  totalBR += decays.at("~chi0_3").BF("~chi0_1", "bbar", "b");
1701  totalBR += decays.at("~chi0_3").BF("~chi0_1", "e+", "e-");
1702  totalBR += decays.at("~chi0_3").BF("~chi0_1", "mu+", "mu-");
1703  totalBR += decays.at("~chi0_3").BF("~chi0_1", "tau+", "tau-");
1704  totalBR += decays.at("~chi0_3").BF("~chi0_1", "nubar_e", "nu_e");
1705  totalBR += decays.at("~chi0_3").BF("~chi0_1", "nubar_mu", "nu_mu");
1706  totalBR += decays.at("~chi0_3").BF("~chi0_1", "nubar_tau", "nu_tau");
1707  xsecWithError.upper *= totalBR;
1708  xsecWithError.central *= totalBR;
1709  xsecWithError.lower *= totalBR;
1710 
1711  if (xsecWithError.central < xsecLimit)
1712  {
1713  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1714  }
1715  else
1716  {
1717  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1718  }
1719 
1720  // neut4, neut1
1721  xsecLimit = limitContainer.limitAverage(mass_neut4, mass_neut1, mZ);
1722 
1723  xsecWithError = *Dep::LEP188_xsec_chi00_14;
1724  // Total up all channels which look like Z* decays
1725  totalBR = 0;
1726  totalBR += decays.at("~chi0_4").BF("~chi0_1", "Z0");
1727  totalBR += decays.at("~chi0_4").BF("~chi0_1", "ubar", "u");
1728  totalBR += decays.at("~chi0_4").BF("~chi0_1", "dbar", "d");
1729  totalBR += decays.at("~chi0_4").BF("~chi0_1", "cbar", "c");
1730  totalBR += decays.at("~chi0_4").BF("~chi0_1", "sbar", "s");
1731  totalBR += decays.at("~chi0_4").BF("~chi0_1", "bbar", "b");
1732  totalBR += decays.at("~chi0_4").BF("~chi0_1", "e+", "e-");
1733  totalBR += decays.at("~chi0_4").BF("~chi0_1", "mu+", "mu-");
1734  totalBR += decays.at("~chi0_4").BF("~chi0_1", "tau+", "tau-");
1735  totalBR += decays.at("~chi0_4").BF("~chi0_1", "nubar_e", "nu_e");
1736  totalBR += decays.at("~chi0_4").BF("~chi0_1", "nubar_mu", "nu_mu");
1737  totalBR += decays.at("~chi0_4").BF("~chi0_1", "nubar_tau", "nu_tau");
1738  xsecWithError.upper *= totalBR;
1739  xsecWithError.central *= totalBR;
1740  xsecWithError.lower *= totalBR;
1741 
1742  if (xsecWithError.central < xsecLimit)
1743  {
1744  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1745  }
1746  else
1747  {
1748  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1749  }
1750 
1751  }
1752 
1754  {
1756  using std::pow;
1757  using std::log;
1758 
1759  const Spectrum& spec = *Dep::MSSM_spectrum;
1760 
1761  const DecayTable& decays = *Dep::decay_rates;
1762  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1763  const double mass_neut2 = spec.get(Par::Pole_Mass,1000023, 0);
1764  const double mass_neut3 = spec.get(Par::Pole_Mass,1000025, 0);
1765  const double mass_neut4 = spec.get(Par::Pole_Mass,1000035, 0);
1766  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1767  triplet<double> xsecWithError;
1768  double xsecLimit, totalBR;
1769 
1770  static const L3NeutralinoLeptonicLimitAt188pt6GeV limitContainer;
1771  // #ifdef COLLIDERBIT_DEBUG
1772  // static bool dumped=false;
1773  // if(!dumped)
1774  // {
1775  // limitContainer.dumpPlotData(0., 200., 0., 100., mZ, "lepLimitPlanev2/L3NeutralinoLeptonicLimitAt188pt6GeV.dump");
1776  // dumped=true;
1777  // }
1778  // #endif
1779 
1780  result = 0;
1781  // Due to the nature of the analysis details of the model independent limit in
1782  // the paper, the best we can do is to try these processes individually:
1783 
1784  // neut2, neut1
1785  xsecLimit = limitContainer.limitAverage(mass_neut2, mass_neut1, mZ);
1786 
1787  xsecWithError = *Dep::LEP188_xsec_chi00_12;
1788  // Total up all channels which look like leptonic Z* decays
1789  // Total up the leptonic Z decays first...
1790  totalBR = 0;
1791  totalBR += decays.at("Z0").BF("e+", "e-");
1792  totalBR += decays.at("Z0").BF("mu+", "mu-");
1793  totalBR += decays.at("Z0").BF("tau+", "tau-");
1794  totalBR = decays.at("~chi0_2").BF("~chi0_1", "Z0") * totalBR;
1795 
1796  totalBR += decays.at("~chi0_2").BF("~chi0_1", "e+", "e-");
1797  totalBR += decays.at("~chi0_2").BF("~chi0_1", "mu+", "mu-");
1798  totalBR += decays.at("~chi0_2").BF("~chi0_1", "tau+", "tau-");
1799  xsecWithError.upper *= totalBR;
1800  xsecWithError.central *= totalBR;
1801  xsecWithError.lower *= totalBR;
1802 
1803  if (xsecWithError.central < xsecLimit)
1804  {
1805  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1806  }
1807  else
1808  {
1809  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1810  }
1811 
1812  // neut3, neut1
1813  xsecLimit = limitContainer.limitAverage(mass_neut3, mass_neut1, mZ);
1814 
1815  xsecWithError = *Dep::LEP188_xsec_chi00_13;
1816  // Total up all channels which look like leptonic Z* decays
1817  // Total up the leptonic Z decays first...
1818  totalBR = 0;
1819  totalBR += decays.at("Z0").BF("e+", "e-");
1820  totalBR += decays.at("Z0").BF("mu+", "mu-");
1821  totalBR += decays.at("Z0").BF("tau+", "tau-");
1822  totalBR = decays.at("~chi0_3").BF("~chi0_1", "Z0") * totalBR;
1823 
1824  totalBR += decays.at("~chi0_3").BF("~chi0_1", "e+", "e-");
1825  totalBR += decays.at("~chi0_3").BF("~chi0_1", "mu+", "mu-");
1826  totalBR += decays.at("~chi0_3").BF("~chi0_1", "tau+", "tau-");
1827  xsecWithError.upper *= totalBR;
1828  xsecWithError.central *= totalBR;
1829  xsecWithError.lower *= totalBR;
1830 
1831  if (xsecWithError.central < xsecLimit)
1832  {
1833  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1834  }
1835  else
1836  {
1837  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1838  }
1839 
1840  // neut4, neut1
1841  xsecLimit = limitContainer.limitAverage(mass_neut4, mass_neut1, mZ);
1842 
1843  xsecWithError = *Dep::LEP188_xsec_chi00_14;
1844  // Total up all channels which look like leptonic Z* decays
1845  // Total up the leptonic Z decays first...
1846  totalBR = 0;
1847  totalBR += decays.at("Z0").BF("e+", "e-");
1848  totalBR += decays.at("Z0").BF("mu+", "mu-");
1849  totalBR += decays.at("Z0").BF("tau+", "tau-");
1850  totalBR = decays.at("~chi0_4").BF("~chi0_1", "Z0") * totalBR;
1851 
1852  totalBR += decays.at("~chi0_4").BF("~chi0_1", "e+", "e-");
1853  totalBR += decays.at("~chi0_4").BF("~chi0_1", "mu+", "mu-");
1854  totalBR += decays.at("~chi0_4").BF("~chi0_1", "tau+", "tau-");
1855  xsecWithError.upper *= totalBR;
1856  xsecWithError.central *= totalBR;
1857  xsecWithError.lower *= totalBR;
1858 
1859  if (xsecWithError.central < xsecLimit)
1860  {
1861  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1862  }
1863  else
1864  {
1865  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1866  }
1867 
1868  }
1869 
1871  {
1873  using std::pow;
1874  using std::log;
1875 
1876  const Spectrum& spec = *Dep::MSSM_spectrum;
1877 
1878  const DecayTable& decays = *Dep::decay_rates;
1879  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1880  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
1881  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
1882  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1883  triplet<double> xsecWithError;
1884  double xsecLimit, totalBR;
1885 
1886  static const L3CharginoAllChannelsLimitAt188pt6GeV limitContainer;
1887  // #ifdef COLLIDERBIT_DEBUG
1888  // static bool dumped=false;
1889  // if(!dumped)
1890  // {
1891  // limitContainer.dumpPlotData(45., 100., 0., 100., mZ, "lepLimitPlanev2/L3CharginoAllChannelsLimitAt188pt6GeV.dump");
1892  // dumped=true;
1893  // }
1894  // #endif
1895 
1896  result = 0;
1897  // Due to the nature of the analysis details of the model independent limit in
1898  // the paper, the best we can do is to try these processes individually:
1899 
1900  // char1, neut1
1901  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
1902 
1903  xsecWithError = *Dep::LEP188_xsec_chipm_11;
1904  // Total up all channels which look like W* decays
1905  totalBR = 0;
1906  totalBR += decays.at("~chi+_1").BF("~chi0_1", "W+");
1907  totalBR += decays.at("~chi+_1").BF("~chi0_1", "u", "dbar");
1908  totalBR += decays.at("~chi+_1").BF("~chi0_1", "c", "sbar");
1909  totalBR += decays.at("~chi+_1").BF("~chi0_1", "e+", "nu_e");
1910  totalBR += decays.at("~chi+_1").BF("~chi0_1", "mu+", "nu_mu");
1911  totalBR += decays.at("~chi+_1").BF("~chi0_1", "tau+", "nu_tau");
1912  xsecWithError.upper *= pow(totalBR, 2);
1913  xsecWithError.central *= pow(totalBR, 2);
1914  xsecWithError.lower *= pow(totalBR, 2);
1915 
1916  if (xsecWithError.central < xsecLimit)
1917  {
1918  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1919  }
1920  else
1921  {
1922  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1923  }
1924 
1925  // char2, neut1
1926  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
1927 
1928  xsecWithError = *Dep::LEP188_xsec_chipm_22;
1929  // Total up all channels which look like W* decays
1930  totalBR = 0;
1931  totalBR += decays.at("~chi+_2").BF("~chi0_1", "W+");
1932  totalBR += decays.at("~chi+_2").BF("~chi0_1", "u", "dbar");
1933  totalBR += decays.at("~chi+_2").BF("~chi0_1", "c", "sbar");
1934  totalBR += decays.at("~chi+_2").BF("~chi0_1", "e+", "nu_e");
1935  totalBR += decays.at("~chi+_2").BF("~chi0_1", "mu+", "nu_mu");
1936  totalBR += decays.at("~chi+_2").BF("~chi0_1", "tau+", "nu_tau");
1937  xsecWithError.upper *= pow(totalBR, 2);
1938  xsecWithError.central *= pow(totalBR, 2);
1939  xsecWithError.lower *= pow(totalBR, 2);
1940 
1941  if (xsecWithError.central < xsecLimit)
1942  {
1943  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
1944  }
1945  else
1946  {
1947  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
1948  }
1949 
1950  }
1951 
1953  {
1955  using std::pow;
1956  using std::log;
1957 
1958  const Spectrum& spec = *Dep::MSSM_spectrum;
1959 
1960  const DecayTable& decays = *Dep::decay_rates;
1961  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
1962  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
1963  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
1964  const double mZ = spec.get(Par::Pole_Mass,23, 0);
1965  triplet<double> xsecWithError;
1966  double xsecLimit, totalBR;
1967 
1968  static const L3CharginoLeptonicLimitAt188pt6GeV limitContainer;
1969  // #ifdef COLLIDERBIT_DEBUG
1970  // static bool dumped=false;
1971  // if(!dumped)
1972  // {
1973  // limitContainer.dumpPlotData(45., 100., 0., 100., mZ, "lepLimitPlanev2/L3CharginoLeptonicLimitAt188pt6GeV.dump");
1974  // dumped=true;
1975  // }
1976  // #endif
1977 
1978  result = 0;
1979  // Due to the nature of the analysis details of the model independent limit in
1980  // the paper, the best we can do is to try these processes individually:
1981 
1982  // char1, neut1
1983  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
1984 
1985  xsecWithError = *Dep::LEP188_xsec_chipm_11;
1986  // Total up all channels which look like leptonic W* decays
1987  // Total up the leptonic W decays first...
1988  totalBR = 0;
1989  totalBR += decays.at("W+").BF("e+", "nu_e");
1990  totalBR += decays.at("W+").BF("mu+", "nu_mu");
1991  totalBR += decays.at("W+").BF("tau+", "nu_tau");
1992  totalBR = decays.at("~chi+_1").BF("~chi0_1", "W+") * totalBR;
1993 
1994  totalBR += decays.at("~chi+_1").BF("~chi0_1", "e+", "nu_e");
1995  totalBR += decays.at("~chi+_1").BF("~chi0_1", "mu+", "nu_mu");
1996  totalBR += decays.at("~chi+_1").BF("~chi0_1", "tau+", "nu_tau");
1997  xsecWithError.upper *= pow(totalBR, 2);
1998  xsecWithError.central *= pow(totalBR, 2);
1999  xsecWithError.lower *= pow(totalBR, 2);
2000 
2001  if (xsecWithError.central < xsecLimit)
2002  {
2003  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2004  }
2005  else
2006  {
2007  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2008  }
2009 
2010  // char2, neut1
2011  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
2012 
2013  xsecWithError = *Dep::LEP188_xsec_chipm_22;
2014  // Total up all channels which look like leptonic W* decays
2015  // Total up the leptonic W decays first...
2016  totalBR = 0;
2017  totalBR += decays.at("W+").BF("e+", "nu_e");
2018  totalBR += decays.at("W+").BF("mu+", "nu_mu");
2019  totalBR += decays.at("W+").BF("tau+", "nu_tau");
2020  totalBR = decays.at("~chi+_2").BF("~chi0_1", "W+") * totalBR;
2021 
2022  totalBR += decays.at("~chi+_2").BF("~chi0_1", "e+", "nu_e");
2023  totalBR += decays.at("~chi+_2").BF("~chi0_1", "mu+", "nu_mu");
2024  totalBR += decays.at("~chi+_2").BF("~chi0_1", "tau+", "nu_tau");
2025  xsecWithError.upper *= pow(totalBR, 2);
2026  xsecWithError.central *= pow(totalBR, 2);
2027  xsecWithError.lower *= pow(totalBR, 2);
2028 
2029  if (xsecWithError.central < xsecLimit)
2030  {
2031  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2032  }
2033  else
2034  {
2035  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2036  }
2037 
2038  }
2039 
2041  {
2043  using std::pow;
2044  using std::log;
2045 
2046  const Spectrum& spec = *Dep::MSSM_spectrum;
2047 
2048  const DecayTable& decays = *Dep::decay_rates;
2049  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2050  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
2051  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
2052  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2053  triplet<double> xsecWithError;
2054  double xsecLimit, totalBR;
2055 
2056  static const OPALCharginoHadronicLimitAt208GeV limitContainer;
2057  // #ifdef COLLIDERBIT_DEBUG
2058  // static bool dumped=false;
2059  // if(!dumped)
2060  // {
2061  // limitContainer.dumpPlotData(75., 105., 0., 105., mZ, "lepLimitPlanev2/OPALCharginoHadronicLimitAt208GeV.dump");
2062  // dumped=true;
2063  // }
2064  // #endif
2065 
2066  result = 0;
2067  // Due to the nature of the analysis details of the model independent limit in
2068  // the paper, the best we can do is to try these processes individually:
2069 
2070  // char1, neut1
2071  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
2072 
2073  xsecWithError = *Dep::LEP208_xsec_chipm_11;
2074  // Total up all channels which look like hadronic W* decays
2075  // Total up the hadronic W decays first...
2076  totalBR = decays.at("W+").BF("hadron", "hadron");
2077  totalBR = decays.at("~chi+_1").BF("~chi0_1", "W+") * totalBR;
2078 
2079  totalBR += decays.at("~chi+_1").BF("~chi0_1", "u", "dbar");
2080  totalBR += decays.at("~chi+_1").BF("~chi0_1", "c", "sbar");
2081  xsecWithError.upper *= pow(totalBR, 2);
2082  xsecWithError.central *= pow(totalBR, 2);
2083  xsecWithError.lower *= pow(totalBR, 2);
2084 
2085  if (xsecWithError.central < xsecLimit)
2086  {
2087  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2088  }
2089  else
2090  {
2091  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2092  }
2093 
2094  // char2, neut1
2095  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
2096 
2097  xsecWithError = *Dep::LEP208_xsec_chipm_22;
2098  // Total up all channels which look like hadronic W* decays
2099  // Total up the hadronic W decays first...
2100  totalBR = decays.at("W+").BF("hadron", "hadron");
2101  totalBR = decays.at("~chi+_2").BF("~chi0_1", "W+") * totalBR;
2102 
2103  totalBR += decays.at("~chi+_2").BF("~chi0_1", "u", "dbar");
2104  totalBR += decays.at("~chi+_2").BF("~chi0_1", "c", "sbar");
2105  xsecWithError.upper *= pow(totalBR, 2);
2106  xsecWithError.central *= pow(totalBR, 2);
2107  xsecWithError.lower *= pow(totalBR, 2);
2108 
2109  if (xsecWithError.central < xsecLimit)
2110  {
2111  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2112  }
2113  else
2114  {
2115  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2116  }
2117 
2118  }
2119 
2121  {
2123  using std::pow;
2124  using std::log;
2125  static const double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
2126  static const bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
2127 
2128  const Spectrum& spec = *Dep::MSSM_spectrum;
2129 
2130  const SubSpectrum& mssm = spec.get_HE();
2131  const DecayTable& decays = *Dep::decay_rates;
2132  const str snue = slhahelp::mass_es_from_gauge_es("~nu_e_L", mssm, tol, LOCAL_INFO, pt_error);
2133  const str snumu = slhahelp::mass_es_from_gauge_es("~nu_mu_L", mssm, tol, LOCAL_INFO, pt_error);
2134  const str snutau = slhahelp::mass_es_from_gauge_es("~nu_tau_L", mssm, tol, LOCAL_INFO, pt_error);
2135  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2136  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
2137  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
2138  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2139  triplet<double> xsecWithError;
2140  double xsecLimit, totalBR;
2141 
2142  static const OPALCharginoSemiLeptonicLimitAt208GeV limitContainer;
2143  // #ifdef COLLIDERBIT_DEBUG
2144  // static bool dumped=false;
2145  // if(!dumped)
2146  // {
2147  // limitContainer.dumpPlotData(75., 105., 0., 105., mZ, "lepLimitPlanev2/OPALCharginoSemiLeptonicLimitAt208GeV.dump");
2148  // dumped=true;
2149  // }
2150  // #endif
2151 
2152  result = 0;
2153  // Due to the nature of the analysis details of the model independent limit in
2154  // the paper, the best we can do is to try these processes individually:
2155 
2156  // char1, neut1
2157  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
2158 
2159  xsecWithError = *Dep::LEP208_xsec_chipm_11;
2160  // Total up all channels which look like leptonic W* decays
2161  // Total up the leptonic W decays first...
2162  totalBR = 0;
2163  totalBR += decays.at("W+").BF("e+", "nu_e");
2164  totalBR += decays.at("W+").BF("mu+", "nu_mu");
2165  totalBR += decays.at("W+").BF("tau+", "nu_tau");
2166  totalBR = decays.at("~chi+_1").BF("~chi0_1", "W+") * totalBR;
2167 
2168  totalBR += decays.at("~chi+_1").BF("~chi0_1", "e+", "nu_e");
2169  totalBR += decays.at("~chi+_1").BF("~chi0_1", "mu+", "nu_mu");
2170  totalBR += decays.at("~chi+_1").BF("~chi0_1", "tau+", "nu_tau");
2171  totalBR += decays.at("~chi+_1").BF(snue, "e+")
2172  * decays.at(snue).BF("~chi0_1", "nu_e");
2173  totalBR += decays.at("~chi+_1").BF(snumu, "mu+")
2174  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2175  totalBR += decays.at("~chi+_1").BF(snutau, "tau+")
2176  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2177  xsecWithError.upper *= totalBR;
2178  xsecWithError.central *= totalBR;
2179  xsecWithError.lower *= totalBR;
2180 
2181  // ALSO, total up all channels which look like hadronic W* decays
2182  // Total up the hadronic W decays first...
2183  totalBR = decays.at("W+").BF("hadron", "hadron");
2184  totalBR = decays.at("~chi+_1").BF("~chi0_1", "W+") * totalBR;
2185 
2186  totalBR += decays.at("~chi+_1").BF("~chi0_1", "u", "dbar");
2187  totalBR += decays.at("~chi+_1").BF("~chi0_1", "c", "sbar");
2188  xsecWithError.upper *= totalBR;
2189  xsecWithError.central *= totalBR;
2190  xsecWithError.lower *= totalBR;
2191 
2192  if (xsecWithError.central < xsecLimit)
2193  {
2194  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2195  }
2196  else
2197  {
2198  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2199  }
2200 
2201  // char2, neut1
2202  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
2203 
2204  xsecWithError = *Dep::LEP208_xsec_chipm_22;
2205  // Total up all channels which look like leptonic W* decays
2206  // Total up the leptonic W decays first...
2207  totalBR = 0;
2208  totalBR += decays.at("W+").BF("e+", "nu_e");
2209  totalBR += decays.at("W+").BF("mu+", "nu_mu");
2210  totalBR += decays.at("W+").BF("tau+", "nu_tau");
2211  totalBR = decays.at("~chi+_2").BF("~chi0_1", "W+") * totalBR;
2212 
2213  totalBR += decays.at("~chi+_2").BF("~chi0_1", "e+", "nu_e");
2214  totalBR += decays.at("~chi+_2").BF("~chi0_1", "mu+", "nu_mu");
2215  totalBR += decays.at("~chi+_2").BF("~chi0_1", "tau+", "nu_tau");
2216  totalBR += decays.at("~chi+_2").BF(snue, "e+")
2217  * decays.at(snue).BF("~chi0_1", "nu_e");
2218  totalBR += decays.at("~chi+_2").BF(snumu, "mu+")
2219  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2220  totalBR += decays.at("~chi+_2").BF(snutau, "tau+")
2221  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2222  xsecWithError.upper *= totalBR;
2223  xsecWithError.central *= totalBR;
2224  xsecWithError.lower *= totalBR;
2225 
2226  // ALSO, total up all channels which look like hadronic W* decays
2227  // Total up the hadronic W decays first...
2228  totalBR = decays.at("W+").BF("hadron", "hadron");
2229  totalBR = decays.at("~chi+_2").BF("~chi0_1", "W+") * totalBR;
2230 
2231  totalBR += decays.at("~chi+_2").BF("~chi0_1", "u", "dbar");
2232  totalBR += decays.at("~chi+_2").BF("~chi0_1", "c", "sbar");
2233  xsecWithError.upper *= totalBR;
2234  xsecWithError.central *= totalBR;
2235  xsecWithError.lower *= totalBR;
2236 
2237  if (xsecWithError.central < xsecLimit)
2238  {
2239  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2240  }
2241  else
2242  {
2243  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2244  }
2245 
2246  }
2247 
2249  {
2251  static const double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
2252  static const bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
2253 
2254  using std::pow;
2255  using std::log;
2256 
2257  const Spectrum& spec = *Dep::MSSM_spectrum;
2258 
2259  const SubSpectrum& mssm = spec.get_HE();
2260  const DecayTable& decays = *Dep::decay_rates;
2261  const str snue = slhahelp::mass_es_from_gauge_es("~nu_e_L", mssm, tol, LOCAL_INFO, pt_error);
2262  const str snumu = slhahelp::mass_es_from_gauge_es("~nu_mu_L", mssm, tol, LOCAL_INFO, pt_error);
2263  const str snutau = slhahelp::mass_es_from_gauge_es("~nu_tau_L", mssm, tol, LOCAL_INFO, pt_error);
2264  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2265  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
2266  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
2267  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2268  triplet<double> xsecWithError;
2269  double xsecLimit, totalBR;
2270 
2271  static const OPALCharginoLeptonicLimitAt208GeV limitContainer;
2272  // #ifdef COLLIDERBIT_DEBUG
2273  // static bool dumped=false;
2274  // if(!dumped)
2275  // {
2276  // limitContainer.dumpPlotData(75., 105., 0., 105., mZ, "lepLimitPlanev2/OPALCharginoLeptonicLimitAt208GeV.dump");
2277  // dumped=true;
2278  // }
2279  // #endif
2280 
2281  result = 0;
2282  // Due to the nature of the analysis details of the model independent limit in
2283  // the paper, the best we can do is to try these processes individually:
2284 
2285  // char1, neut1
2286  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
2287 
2288  xsecWithError = *Dep::LEP208_xsec_chipm_11;
2289  // Total up all channels which look like leptonic W* decays
2290  // Total up the leptonic W decays first...
2291  totalBR = 0;
2292  totalBR += decays.at("W+").BF("e+", "nu_e");
2293  totalBR += decays.at("W+").BF("mu+", "nu_mu");
2294  totalBR += decays.at("W+").BF("tau+", "nu_tau");
2295  totalBR = decays.at("~chi+_1").BF("~chi0_1", "W+") * totalBR;
2296 
2297  totalBR += decays.at("~chi+_1").BF("~chi0_1", "e+", "nu_e");
2298  totalBR += decays.at("~chi+_1").BF("~chi0_1", "mu+", "nu_mu");
2299  totalBR += decays.at("~chi+_1").BF("~chi0_1", "tau+", "nu_tau");
2300  totalBR += decays.at("~chi+_1").BF(snue, "e+")
2301  * decays.at(snue).BF("~chi0_1", "nu_e");
2302  totalBR += decays.at("~chi+_1").BF(snumu, "mu+")
2303  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2304  totalBR += decays.at("~chi+_1").BF(snutau, "tau+")
2305  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2306  xsecWithError.upper *= pow(totalBR, 2);
2307  xsecWithError.central *= pow(totalBR, 2);
2308  xsecWithError.lower *= pow(totalBR, 2);
2309 
2310  if (xsecWithError.central < xsecLimit)
2311  {
2312  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2313  }
2314  else
2315  {
2316  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2317  }
2318 
2319  // char2, neut1
2320  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
2321 
2322  xsecWithError = *Dep::LEP208_xsec_chipm_22;
2323  // Total up all channels which look like leptonic W* decays
2324  // Total up the leptonic W decays first...
2325  totalBR = 0;
2326  totalBR += decays.at("W+").BF("e+", "nu_e");
2327  totalBR += decays.at("W+").BF("mu+", "nu_mu");
2328  totalBR += decays.at("W+").BF("tau+", "nu_tau");
2329  totalBR = decays.at("~chi+_2").BF("~chi0_1", "W+") * totalBR;
2330 
2331  totalBR += decays.at("~chi+_2").BF("~chi0_1", "e+", "nu_e");
2332  totalBR += decays.at("~chi+_2").BF("~chi0_1", "mu+", "nu_mu");
2333  totalBR += decays.at("~chi+_2").BF("~chi0_1", "tau+", "nu_tau");
2334  totalBR += decays.at("~chi+_2").BF(snue, "e+")
2335  * decays.at(snue).BF("~chi0_1", "nu_e");
2336  totalBR += decays.at("~chi+_2").BF(snumu, "mu+")
2337  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2338  totalBR += decays.at("~chi+_2").BF(snutau, "tau+")
2339  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2340  xsecWithError.upper *= pow(totalBR, 2);
2341  xsecWithError.central *= pow(totalBR, 2);
2342  xsecWithError.lower *= pow(totalBR, 2);
2343 
2344  if (xsecWithError.central < xsecLimit)
2345  {
2346  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2347  }
2348  else
2349  {
2350  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2351  }
2352 
2353  }
2354 
2355  // OPAL limit on degenerate chargino--neutralino scenario at 208 GeV
2356  // Sensitive to mass differences between 320 MeV and 5 GeV
2357  // Based on hep-ex/0210043
2359  {
2361 
2362  const Spectrum& spec = *Dep::MSSM_spectrum;
2363  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2364  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
2365  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2366  triplet<double> xsecWithError;
2367  double xsecLimit;
2368 
2369  static const OPALDegenerateCharginoLimitAt208GeV limitContainer;
2370  // #ifdef COLLIDERBIT_DEBUG
2371  // static bool dumped=false;
2372  // if(!dumped)
2373  // {
2374  // limitContainer.dumpPlotData(45.0, 95., 0.320, 5., mZ, "lepLimitPlanev2/OPALDegenerateCharginoLimitAt208GeV.dump");
2375  //
2376  // dumped=true;
2377  // }
2378  // #endif
2379 
2380  result = 0;
2381 
2382  // char1, neut1
2383  xsecLimit = limitContainer.limitAverage(mass_char1, mass_char1-abs(mass_neut1), mZ);
2384  xsecWithError = *Dep::LEP208_xsec_chipm_11;
2385 
2386  if (xsecWithError.central < xsecLimit)
2387  {
2388  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2389  }
2390  else
2391  {
2392  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2393  }
2394 
2395  }
2396 
2397 
2399  {
2401  static const double tol = runOptions->getValueOrDef<double>(1e-2, "gauge_mixing_tolerance");
2402  static const bool pt_error = runOptions->getValueOrDef<bool>(true, "gauge_mixing_tolerance_invalidates_point_only");
2403 
2404  const Spectrum& spec = *Dep::MSSM_spectrum;
2405 
2406  const SubSpectrum& mssm = spec.get_HE();
2407  const DecayTable& decays = *Dep::decay_rates;
2408  const str snue = slhahelp::mass_es_from_gauge_es("~nu_e_L", mssm, tol, LOCAL_INFO, pt_error);
2409  const str snumu = slhahelp::mass_es_from_gauge_es("~nu_mu_L", mssm, tol, LOCAL_INFO, pt_error);
2410  const str snutau = slhahelp::mass_es_from_gauge_es("~nu_tau_L", mssm, tol, LOCAL_INFO, pt_error);
2411  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2412  const double mass_char1 = spec.get(Par::Pole_Mass,1000024, 0);
2413  const double mass_char2 = spec.get(Par::Pole_Mass,1000037, 0);
2414  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2415  triplet<double> xsecWithError;
2416  double xsecLimit, totalBR;
2417 
2418  static const OPALCharginoAllChannelsLimitAt208GeV limitContainer;
2419  // #ifdef COLLIDERBIT_DEBUG
2420  // static bool dumped=false;
2421  // if(!dumped)
2422  // {
2423  // limitContainer.dumpPlotData(75., 105., 0., 105., mZ, "lepLimitPlanev2/OPALCharginoAllChannelsLimitAt208GeV.dump");
2424  // dumped=true;
2425  // }
2426  // #endif
2427 
2428  result = 0;
2429  // Due to the nature of the analysis details of the model independent limit in
2430  // the paper, the best we can do is to try these processes individually:
2431 
2432  // char1, neut1
2433  xsecLimit = limitContainer.limitAverage(mass_char1, mass_neut1, mZ);
2434 
2435  xsecWithError = *Dep::LEP208_xsec_chipm_11;
2436  // Total up all channels which look like W* decays
2437  totalBR = 0;
2438  totalBR += decays.at("~chi+_1").BF("~chi0_1", "W+");
2439  totalBR += decays.at("~chi+_1").BF("~chi0_1", "u", "dbar");
2440  totalBR += decays.at("~chi+_1").BF("~chi0_1", "c", "sbar");
2441  totalBR += decays.at("~chi+_1").BF("~chi0_1", "e+", "nu_e");
2442  totalBR += decays.at("~chi+_1").BF("~chi0_1", "mu+", "nu_mu");
2443  totalBR += decays.at("~chi+_1").BF("~chi0_1", "tau+", "nu_tau");
2444  totalBR += decays.at("~chi+_1").BF(snue, "e+")
2445  * decays.at(snue).BF("~chi0_1", "nu_e");
2446  totalBR += decays.at("~chi+_1").BF(snumu, "mu+")
2447  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2448  totalBR += decays.at("~chi+_1").BF(snutau, "tau+")
2449  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2450  xsecWithError.upper *= pow(totalBR, 2);
2451  xsecWithError.central *= pow(totalBR, 2);
2452  xsecWithError.lower *= pow(totalBR, 2);
2453 
2454  if (xsecWithError.central < xsecLimit)
2455  {
2456  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2457  }
2458  else
2459  {
2460  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2461  }
2462 
2463  // char2, neut1
2464  xsecLimit = limitContainer.limitAverage(mass_char2, mass_neut1, mZ);
2465 
2466  xsecWithError = *Dep::LEP208_xsec_chipm_22;
2467  // Total up all channels which look like W* decays
2468  totalBR = 0;
2469  totalBR += decays.at("~chi+_2").BF("~chi0_1", "W+");
2470  totalBR += decays.at("~chi+_2").BF("~chi0_1", "u", "dbar");
2471  totalBR += decays.at("~chi+_2").BF("~chi0_1", "c", "sbar");
2472  totalBR += decays.at("~chi+_2").BF("~chi0_1", "e+", "nu_e");
2473  totalBR += decays.at("~chi+_2").BF("~chi0_1", "mu+", "nu_mu");
2474  totalBR += decays.at("~chi+_2").BF("~chi0_1", "tau+", "nu_tau");
2475  totalBR += decays.at("~chi+_2").BF(snue, "e+")
2476  * decays.at(snue).BF("~chi0_1", "nu_e");
2477  totalBR += decays.at("~chi+_2").BF(snumu, "mu+")
2478  * decays.at(snumu).BF("~chi0_1", "nu_mu");
2479  totalBR += decays.at("~chi+_2").BF(snutau, "tau+")
2480  * decays.at(snutau).BF("~chi0_1", "nu_tau");
2481  xsecWithError.upper *= pow(totalBR, 2);
2482  xsecWithError.central *= pow(totalBR, 2);
2483  xsecWithError.lower *= pow(totalBR, 2);
2484 
2485  if (xsecWithError.central < xsecLimit)
2486  {
2487  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2488  }
2489  else
2490  {
2491  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2492  }
2493 
2494  }
2495 
2497  {
2499  using std::pow;
2500  using std::log;
2501 
2502  const Spectrum& spec = *Dep::MSSM_spectrum;
2503 
2504  const DecayTable& decays = *Dep::decay_rates;
2505  const double mass_neut1 = spec.get(Par::Pole_Mass,1000022, 0);
2506  const double mass_neut2 = spec.get(Par::Pole_Mass,1000023, 0);
2507  const double mass_neut3 = spec.get(Par::Pole_Mass,1000025, 0);
2508  const double mass_neut4 = spec.get(Par::Pole_Mass,1000035, 0);
2509  const double mZ = spec.get(Par::Pole_Mass,23, 0);
2510  triplet<double> xsecWithError;
2511  double xsecLimit, totalBR;
2512 
2513  static const OPALNeutralinoHadronicLimitAt208GeV limitContainer;
2514  // #ifdef COLLIDERBIT_DEBUG
2515  // static bool dumped=false;
2516  // if(!dumped)
2517  // {
2518  // limitContainer.dumpPlotData(0., 200., 0., 100., mZ, "lepLimitPlanev2/OPALNeutralinoHadronicLimitAt208GeV.dump");
2519  // dumped=true;
2520  // }
2521  // #endif
2522 
2523  result = 0;
2524  // Due to the nature of the analysis details of the model independent limit in
2525  // the paper, the best we can do is to try these processes individually:
2526 
2527  // neut2, neut1
2528  xsecLimit = limitContainer.limitAverage(mass_neut2, mass_neut1, mZ);
2529 
2530  xsecWithError = *Dep::LEP208_xsec_chi00_12;
2531  // Total up all channels which look like Z* decays
2532  totalBR = decays.at("Z0").BF("hadron", "hadron");
2533  totalBR = decays.at("~chi0_2").BF("~chi0_1", "Z0") * totalBR;
2534  totalBR += decays.at("~chi0_2").BF("~chi0_1", "ubar", "u");
2535  totalBR += decays.at("~chi0_2").BF("~chi0_1", "dbar", "d");
2536  totalBR += decays.at("~chi0_2").BF("~chi0_1", "cbar", "c");
2537  totalBR += decays.at("~chi0_2").BF("~chi0_1", "sbar", "s");
2538  totalBR += decays.at("~chi0_2").BF("~chi0_1", "bbar", "b");
2539  xsecWithError.upper *= totalBR;
2540  xsecWithError.central *= totalBR;
2541  xsecWithError.lower *= totalBR;
2542 
2543  if (xsecWithError.central < xsecLimit)
2544  {
2545  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2546  }
2547  else
2548  {
2549  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2550  }
2551 
2552  // neut3, neut1
2553  xsecLimit = limitContainer.limitAverage(mass_neut3, mass_neut1, mZ);
2554 
2555  xsecWithError = *Dep::LEP208_xsec_chi00_13;
2556  // Total up all channels which look like Z* decays
2557  totalBR = decays.at("Z0").BF("hadron", "hadron");
2558  totalBR = decays.at("~chi0_3").BF("~chi0_1", "Z0") * totalBR;
2559  totalBR += decays.at("~chi0_3").BF("~chi0_1", "ubar", "u");
2560  totalBR += decays.at("~chi0_3").BF("~chi0_1", "dbar", "d");
2561  totalBR += decays.at("~chi0_3").BF("~chi0_1", "cbar", "c");
2562  totalBR += decays.at("~chi0_3").BF("~chi0_1", "sbar", "s");
2563  totalBR += decays.at("~chi0_3").BF("~chi0_1", "bbar", "b");
2564  xsecWithError.upper *= totalBR;
2565  xsecWithError.central *= totalBR;
2566  xsecWithError.lower *= totalBR;
2567 
2568  if (xsecWithError.central < xsecLimit)
2569  {
2570  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2571  }
2572  else
2573  {
2574  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2575  }
2576 
2577  // neut4, neut1
2578  xsecLimit = limitContainer.limitAverage(mass_neut4, mass_neut1, mZ);
2579 
2580  xsecWithError = *Dep::LEP208_xsec_chi00_14;
2581  // Total up all channels which look like Z* decays
2582  totalBR = decays.at("Z0").BF("hadron", "hadron");
2583  totalBR = decays.at("~chi0_4").BF("~chi0_1", "Z0") * totalBR;
2584  totalBR += decays.at("~chi0_4").BF("~chi0_1", "ubar", "u");
2585  totalBR += decays.at("~chi0_4").BF("~chi0_1", "dbar", "d");
2586  totalBR += decays.at("~chi0_4").BF("~chi0_1", "cbar", "c");
2587  totalBR += decays.at("~chi0_4").BF("~chi0_1", "sbar", "s");
2588  totalBR += decays.at("~chi0_4").BF("~chi0_1", "bbar", "b");
2589  xsecWithError.upper *= totalBR;
2590  xsecWithError.central *= totalBR;
2591  xsecWithError.lower *= totalBR;
2592 
2593  if (xsecWithError.central < xsecLimit)
2594  {
2595  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.upper - xsecWithError.central);
2596  }
2597  else
2598  {
2599  result += limit_LLike(xsecWithError.central, xsecLimit, xsecWithError.central - xsecWithError.lower);
2600  }
2601 
2602  }
2603 
2605 
2606  void L3_Gravitino_LLike(double& result) {
2619  // Unpack neutralino & gravitino mass
2620  using namespace Pipes::L3_Gravitino_LLike;
2621  const Spectrum& spectrum = *Dep::MSSM_spectrum;
2622  const double m_chi = spectrum.get(Par::Pole_Mass, 1000022, 0);
2623  const double m_gravitino = spectrum.get(Par::Pole_Mass, 1000039, 0);
2624 
2625  // Calculate relevant branching ratio
2627  const auto BF = decay_rates.at("~chi0_1").BF("gamma", "~G");
2628 
2629  // Production cross section of two lightest neutralinos at 207 GeV
2630  const auto production_xsec = *Dep::LEP207_xsec_chi00_11;
2631 
2632  // Make product of cross section and branching ratio squared
2633  triplet<double> xsec;
2634  xsec.upper = production_xsec.upper * pow(BF, 2);
2635  xsec.central = production_xsec.central * pow(BF, 2);
2636  xsec.lower = production_xsec.lower * pow(BF, 2);
2637 
2638  // Construct object for fetching limit (do this once only, hence static)
2639  const std::string fig6c = GAMBIT_DIR "/ColliderBit/data/scraped_fig6c.dat";
2640  static auto L3Gravitino = ImageLimit(fig6c, 0., 103., 0., 103.);
2641  const double limit = L3Gravitino.get_limit(m_chi, m_gravitino);
2642 
2643  // Resulting log-likelihood, taking into account theoretical uncertainty
2644  result = limit_LLike(xsec.central, limit, xsec.upper - xsec.central);
2645  }
2646 
2647  } // namespace ColliderBit
2648 } // namespace Gambit
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chi00_22
void LEP188_SLHA1_convention_xsec_smu1smu1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_staulstaulbar
double get(const Par::Tags partype, const std::string &mass) const
Definition: spectrum.cpp:249
void LEP208_SLHA1_convention_xsec_chi00_44(triplet< double > &result)
Rollcall header for ColliderBit module.
void ALEPH_Selectron_Conservative_LLike(double &result)
LEP Slepton Log-Likelihoods.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_stau2stau2bar
void LEP205_SLHA1_convention_xsec_chipm_21(triplet< double > &result)
Read a limit from an image, represented by a grid of values at every pixel.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry ALEPH_Selectron_Conservative_LLike
A class to contain the limit data from OPAL_EPJC35_2004_1, figure 9b.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chi00_11
void LEP208_SLHA1_convention_xsec_stau1stau2bar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_stau1stau1bar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_smulsmulbar(triplet< double > &result)
ee –> smuon pair production cross-sections at 205 GeV
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chi00_24
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_selselbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_smu2smu2bar
double BF(const std::vector< std::pair< int, int > > &) const
Retrieve branching fraction for decay to a given final state.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_stau1stau1bar
void LEP205_SLHA1_convention_xsec_se1se1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chipm_11
void LEP208_SLHA1_convention_xsec_staurstaulbar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_chi00_23(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_staulstaurbar(triplet< double > &result)
A class to contain the limit data from OPAL_EPJC35_2004_1, figure 5b.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry L3_Neutralino_All_Channels_Conservative_LLike
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_stau1stau1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable L3_Chargino_All_Channels_Conservative_LLike
void LEP188_SLHA1_convention_xsec_smulsmulbar(triplet< double > &result)
ee –> smuon pair production cross-sections at 188.6 GeV
void ALEPH_Stau_Conservative_LLike(double &result)
void LEP205_SLHA1_convention_xsec_chipm_12(triplet< double > &result)
void L3_Smuon_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_xsec_chipm_22
void LEP188_SLHA1_convention_xsec_smu2smu1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_stau2stau2bar
void LEP208_SLHA1_convention_xsec_chipm_11(triplet< double > &result)
ee –> chargino pair production cross-sections at 208 GeV
Entry & at(std::pair< int, int >)
Get entry in decay table for a give particle, throwing an error if particle is absent.
void LEP205_SLHA1_convention_xsec_chi00_12(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_stau1stau1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_smu1smu1bar
A class to contain the limit data from L3PLB_472_2000_420, figure 3b.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_staurstaurbar
A class to contain the limit data from OPAL, hep-ex/0210043, figure 5a (in colour) ...
void LEP188_SLHA1_convention_xsec_smu1smu2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chi00_24
void LEP208_SLHA1_convention_xsec_chipm_12(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_smursmurbar
void LEP188_SLHA1_convention_xsec_smu2smu2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_smu1smu1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chipm_11
void LEP208_SLHA1_convention_xsec_stau1stau1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_smulsmulbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_smu2smu2bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry stau_2_decay_rates
void LEP188_SLHA1_convention_xsec_selselbar(triplet< double > &result)
ee –> selectron pair production cross-sections at 188.6 GeV
void LEP205_SLHA1_convention_xsec_smulsmurbar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_chipm_11(triplet< double > &result)
ee –> chargino pair production cross-sections at 188.6 GeV
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_chipm_22
void OPAL_Degenerate_Chargino_Conservative_LLike(double &result)
void L3_Chargino_Leptonic_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_smulsmulbar
#define LOCAL_INFO
Definition: local_info.hpp:34
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chipm_22
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_smulsmulbar
void LEP208_SLHA1_convention_xsec_stau2stau2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_xsec_chi00_14
A class to contain the limit data from L3_PLB580_2004_37, figure 2a.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_smursmurbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chipm_22
DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_se1se1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_se2se2bar
A class to contain the limit data from OPAL_EPJC35_2004_1, figure 7b.
void LEP208_SLHA1_convention_xsec_chi00_33(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_staulstaulbar(triplet< double > &result)
ee –> stau pair production cross-sections at 205 GeV
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_staurstaurbar
void LEP188_SLHA1_convention_xsec_chi00_11(triplet< double > &result)
ee –> neutralino pair production cross-sections at 188.6 GeV
void LEP208_SLHA1_convention_xsec_smursmulbar(triplet< double > &result)
double limitAverage(double x, double y, double mZ) const
Two-pi averaging interpolator to find limits between limit curves.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_smursmurbar
double limit_LLike(double x, double x95, double sigma)
LEP limit likelihood function.
void LEP205_SLHA1_convention_xsec_stau2stau1bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry selectron_r_decay_rates
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_smu1smu1bar
void LEP208_SLHA1_convention_xsec_selselbar(triplet< double > &result)
ee –> selectron pair production cross-sections at 208 GeV
void LEP205_SLHA1_convention_xsec_se2se2bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_selserbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry selectron_l_decay_rates
DecayTable::Entry LEP208_SLHA1_convention_xsec_serserbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_smursmurbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable DecayTable DecayTable DecayTable triplet< double > L3_Gravitino_LLike
A class to contain the limit data from L3_PLB580_2004_37, figure 2b.
A class to contain the limit data from OPAL_EPJC35_2004_1, figure 8b.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable DecayTable OPAL_Chargino_Hadronic_Conservative_LLike
void LEP205_SLHA1_convention_xsec_chi00_24(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_serserbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_staulstaulbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chi00_34
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_xsec_chi00_13
void L3_Stau_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chi00_24
void LEP208_SLHA1_convention_xsec_se2se2bar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_chipm_22(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_serserbar
Sparticle production cross-section calculators for LEP.
void LEP205_SLHA1_convention_xsec_serserbar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_stau2stau1bar(triplet< double > &result)
A class to contain the limit data from L3PLB_472_2000_420, figure 3a.
void OPAL_Chargino_Leptonic_Conservative_LLike(double &result)
void LEP205_SLHA1_convention_xsec_selserbar(triplet< double > &result)
void get_sigma_ee_chipm(triplet< double > &result, const double sqrts, const int chi_plus, const int chi_minus, const double tol, const bool pt_error, const Spectrum &spec, const double gammaZ)
Retrieve the production cross-section at an e+e- collider for chargino pairs.
Routines to help translate between SLHA2 sfermions and SLHA1 (or similar) sfermions.
void LEP208_SLHA1_convention_xsec_chi00_14(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_smulsmulbar(triplet< double > &result)
ee –> smuon pair production cross-sections at 208 GeV
void LEP188_SLHA1_convention_xsec_staulstaulbar(triplet< double > &result)
ee –> stau pair production cross-sections at 188.6 GeV
void LEP205_SLHA1_convention_xsec_chipm_11(triplet< double > &result)
ee –> chargino pair production cross-sections at 205 GeV
void L3_Neutralino_All_Channels_Conservative_LLike(double &result)
LEP Gaugino Log-Likelihoods.
void LEP188_SLHA1_convention_xsec_stau1stau2bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_se1se1bar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_se2se1bar(triplet< double > &result)
void get_sigma_ee_ll(triplet< double > &result, const double sqrts, const int generation, const int l_chirality, const int lbar_chirality, const double gtol, const double ftol, const bool gpt_error, const bool fpt_error, const Spectrum &spec, const double gammaZ, const bool l_are_gauge_es)
High-level cross section routines.
void LEP188_SLHA1_convention_xsec_serserbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_stau1stau1bar
void LEP188_SLHA1_convention_xsec_staurstaurbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chi00_13
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_xsec_chipm_11
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chi00_34
void LEP188_SLHA1_convention_xsec_chi00_12(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_se1se2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable DecayTable DecayTable OPAL_Chargino_Leptonic_Conservative_LLike
void LEP188_SLHA1_convention_xsec_chi00_24(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_smulsmulbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chi00_11
LEP208_SLHA1_convention_xsec_selselbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chi00_34
void L3_Neutralino_Leptonic_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_se1se1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_chipm_11
void OPAL_Chargino_Hadronic_Conservative_LLike(double &result)
void L3_Gravitino_LLike(double &result)
const double sigma
Definition: SM_Z.hpp:43
void LEP208_SLHA1_convention_xsec_chipm_22(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_chi00_22(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_smursmurbar(triplet< double > &result)
A class to contain the limit data from OPAL_EPJC35_2004_1, figure 6b.
void LEP208_SLHA1_convention_xsec_chi00_24(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_stau2stau2bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_chipm_11
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_smulsmulbar
void LEP188_SLHA1_convention_xsec_chi00_23(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry stau_1_decay_rates
void LEP188_SLHA1_convention_xsec_chi00_14(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_smu2smu2bar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_stau1stau2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_selselbar
void LEP188_SLHA1_convention_xsec_chi00_13(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chi00_22
void LEP188_SLHA1_convention_xsec_chipm_21(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_stau2stau1bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_chipm_22(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_stau2stau2bar
A class to contain the limit data from ALEPH_PLB526_2002_206, figure 3a.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_chi00_13
void dumpPlotData(double xlow, double xhigh, double ylow, double yhigh, double mZ, std::string filename, int ngrid=100) const
Dump limit average data into a file for average debugging.
void LEP208_SLHA1_convention_xsec_chipm_21(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_staulstaurbar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_se2se1bar(triplet< double > &result)
A class to contain the limit data from L3_PLB580_2004_37, figure 2c.
void LEP208_SLHA1_convention_xsec_chi00_11(triplet< double > &result)
ee –> neutralino pair production cross-sections at 208 GeV
Header file that includes all GAMBIT headers required for a module source file.
void LEP188_SLHA1_convention_xsec_chi00_22(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_serserbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_se2se2bar
void LEP208_SLHA1_convention_xsec_chi00_34(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_chi00_34(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chi00_13
void LEP205_SLHA1_convention_xsec_chi00_13(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_staulstaulbar(triplet< double > &result)
ee –> stau pair production cross-sections at 208 GeV
void OPAL_Chargino_All_Channels_Conservative_LLike(double &result)
std::string str
Shorthand for a standard string.
Definition: Analysis.hpp:35
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_smursmurbar
void LEP208_SLHA1_convention_xsec_selserbar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_selselbar(triplet< double > &result)
ee –> selectron pair production cross-sections at 205 GeV
A class to contain the limit data from ALEPH_PLB526_2002_206, figure 3c.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chi00_22
void LEP205_SLHA1_convention_xsec_smursmulbar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_se1se2bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_smulsmurbar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_smu2smu2bar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_staurstaurbar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_se2se2bar(triplet< double > &result)
Virtual base class for interacting with spectrum generator output.
Definition: subspectrum.hpp:87
void LEP205_SLHA1_convention_xsec_smursmurbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_xsec_serserbar
void LEP208_SLHA1_convention_xsec_smu2smu1bar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_smursmurbar(triplet< double > &result)
str mass_es_closest_to_family(str familystate, const SubSpectrum &mssm)
identify the mass eigenstate corresponding to family state takes string and returns only requested st...
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_stau1stau1bar
void L3_Selectron_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_xsec_stau2stau2bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_staulstaulbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_se2se2bar
GAMBIT native decay table class.
Definition: decay_table.hpp:35
void LEP205_SLHA1_convention_xsec_chi00_11(triplet< double > &result)
ee –> neutralino pair production cross-sections at 205 GeV
void LEP188_SLHA1_convention_xsec_chi00_34(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_smu1smu2bar(triplet< double > &result)
A class to contain the limit data from L3PLB_472_2000_420, figure 2b.
void LEP188_SLHA1_convention_xsec_staurstaulbar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_se2se1bar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_staurstaulbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_xsec_chi00_12
void LEP208_SLHA1_convention_xsec_chi00_12(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry smuon_l_decay_rates
double pow(const double &a)
Outputs a^i.
void LEP205_SLHA1_convention_xsec_chi00_23(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable DecayTable DecayTable DecayTable OPAL_Degenerate_Chargino_Conservative_LLike
void LEP205_SLHA1_convention_xsec_stau2stau2bar(triplet< double > &result)
void LEP207_SLHA1_convention_xsec_chi00_11(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP205_SLHA1_convention_xsec_staurstaurbar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry L3_Smuon_Conservative_LLike
bool is_xsec_sane(const triplet< double > &xsecWithError)
LEP limit debugging function.
void LEP205_SLHA1_convention_xsec_staurstaurbar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_serselbar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_smu2smu1bar(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_se1se1bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_smursmulbar(triplet< double > &result)
void OPAL_Chargino_SemiLeptonic_Conservative_LLike(double &result)
void LEP205_SLHA1_convention_xsec_staulstaurbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chipm_22
void get_sigma_ee_chi00(triplet< double > &result, const double sqrts, const int chi_first, const int chi_second, const double tol, const bool pt_error, const Spectrum &spec, const double gammaZ)
Retrieve the production cross-section at an e+e- collider for neutralino pairs.
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry ALEPH_Stau_Conservative_LLike
void LEP208_SLHA1_convention_xsec_chi00_13(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable DecayTable DecayTable DecayTable triplet< double > LEP207_xsec_chi00_11
DecayTable::Entry LEP208_xsec_serserbar
void LEP208_SLHA1_convention_xsec_smu1smu1bar(triplet< double > &result)
str mass_es_from_gauge_es(str gauge_es, double &max_mixing, std::vector< double > &gauge_composition, const SubSpectrum &mssm)
indentifies the state with largest gauge_es content also fills largest max_mixing and full gauge_comp...
A class to contain the limit data from L3PLB_472_2000_420, figure 2a.
void LEP205_SLHA1_convention_xsec_chi00_14(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry decay_rates
void LEP208_SLHA1_convention_xsec_serselbar(triplet< double > &result)
void LEP205_SLHA1_convention_xsec_chi00_44(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_serselbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_chi00_11
void LEP188_SLHA1_convention_xsec_stau2stau2bar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_se1se1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_selselbar
void LEP188_SLHA1_convention_xsec_chipm_12(triplet< double > &result)
A class to contain the limit data from ALEPH_PLB526_2002_206, figure 3b.
TODO: see if we can use this one:
Definition: Analysis.hpp:33
void LEP205_SLHA1_convention_xsec_chi00_33(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_chi00_33(triplet< double > &result)
void ALEPH_Smuon_Conservative_LLike(double &result)
SubSpectrum & get_HE()
Definition: spectrum.cpp:225
void LEP205_SLHA1_convention_xsec_smu1smu1bar(triplet< double > &result)
"Standard Model" (low-energy) plus high-energy model container class
Definition: spectrum.hpp:110
void LEP208_SLHA1_convention_xsec_chi00_22(triplet< double > &result)
void LEP208_SLHA1_convention_xsec_smu1smu2bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_se1se2bar(triplet< double > &result)
void LEP188_SLHA1_convention_xsec_chi00_44(triplet< double > &result)
void OPAL_Neutralino_Hadronic_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_chi00_13
void L3_Chargino_All_Channels_Conservative_LLike(double &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP188_SLHA1_convention_xsec_stau1stau1bar
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry LEP208_SLHA1_convention_xsec_smu2smu2bar
void LEP208_SLHA1_convention_xsec_smulsmurbar(triplet< double > &result)
DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry DecayTable::Entry smuon_r_decay_rates