Hammer  1.0.0
Helicity Amplitude Module for Matrix Element Reweighting
 All Classes Namespaces Files Functions Variables Typedefs Enumerations Enumerator Friends Macros Groups Pages
RateBD1starLepNu.cc
Go to the documentation of this file.
1 ///
2 /// @file RateBD1starLepNu.cc
3 /// @brief \f$ B \rightarrow D_1^* \tau\nu \f$ total rate
4 ///
5 
6 //**** This file is a part of the HAMMER library
7 //**** Copyright (C) 2016 - 2020 The HAMMER Collaboration
8 //**** HAMMER is licensed under version 3 of the GPL; see COPYING for details
9 //**** Please note the MCnet academic guidelines; see GUIDELINES for details
10 
11 // -*- C++ -*-
12 #include <cmath>
13 #include "Hammer/Rates/RateBD1starLepNu.hh"
14 #include "Hammer/IndexLabels.hh"
15 #include "Hammer/Math/Constants.hh"
16 #include "Hammer/Math/MultiDim/SparseContainer.hh"
17 #include "Hammer/Particle.hh"
18 #include "Hammer/Tools/Pdg.hh"
19 
20 using namespace std;
21 
22 namespace Hammer {
23 
24  namespace MD = MultiDimensional;
25 
26  RateBD1starLepNu::RateBD1starLepNu() {
27  //Create tensor rank and dimensions
28  vector<IndexType> dims{{11, 8, 11, 8, _nPoints}};
29  string name {"RateBD1starLepNuQ2"};
30  auto& pdg = PID::instance();
31  addProcessSignature(PID::BPLUS,{-PID::DSSD1STAR,PID::NU_TAU,PID::ANTITAU});
32  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTITAU), pdg.getMass(PID::BPLUS)-pdg.getMass(PID::DSSD1STAR))});
33  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU,FF_BDSSD1STAR,WILSON_BCTAUNU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
34 
35  addProcessSignature(PID::BZERO,{PID::DSSD1STARMINUS,PID::NU_TAU,PID::ANTITAU});
36  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTITAU), pdg.getMass(PID::BZERO)-pdg.getMass(PID::DSSD1STARMINUS))});
37  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU,FF_BDSSD1STAR,WILSON_BCTAUNU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
38 
39  addProcessSignature(PID::BPLUS,{-PID::DSSD1STAR,PID::NU_MU,PID::ANTIMUON});
40  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTIMUON), pdg.getMass(PID::BPLUS)-pdg.getMass(PID::DSSD1STAR))});
41  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCMUNU,FF_BDSSD1STAR,WILSON_BCMUNU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
42 
43  addProcessSignature(PID::BZERO,{PID::DSSD1STARMINUS,PID::NU_MU,PID::ANTIMUON});
44  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTIMUON), pdg.getMass(PID::BZERO)-pdg.getMass(PID::DSSD1STARMINUS))});
45  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCMUNU,FF_BDSSD1STAR,WILSON_BCMUNU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
46 
47  addProcessSignature(PID::BPLUS,{-PID::DSSD1STAR,PID::NU_E,PID::POSITRON});
48  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::POSITRON), pdg.getMass(PID::BPLUS)-pdg.getMass(PID::DSSD1STAR))});
49  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCENU,FF_BDSSD1STAR,WILSON_BCENU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
50 
51  addProcessSignature(PID::BZERO,{PID::DSSD1STARMINUS,PID::NU_E,PID::POSITRON});
52  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::POSITRON), pdg.getMass(PID::BZERO)-pdg.getMass(PID::DSSD1STARMINUS))});
53  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCENU,FF_BDSSD1STAR,WILSON_BCENU_HC,FF_BDSSD1STAR_HC,INTEGRATION_INDEX})});
54 
55  //bs -> cs
56  addProcessSignature(PID::BS,{PID::DSSDS1STARMINUS,PID::NU_TAU,PID::ANTITAU});
57  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTITAU), pdg.getMass(PID::BS)-pdg.getMass(PID::DSSDS1STARMINUS))});
58  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCTAUNU,FF_BSDSSDS1STAR,WILSON_BCTAUNU_HC,FF_BSDSSDS1STAR_HC,INTEGRATION_INDEX})});
59 
60  addProcessSignature(PID::BS,{PID::DSSDS1STARMINUS,PID::NU_MU,PID::ANTIMUON});
61  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::ANTIMUON), pdg.getMass(PID::BS)-pdg.getMass(PID::DSSDS1STARMINUS))});
62  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCMUNU,FF_BSDSSDS1STAR,WILSON_BCMUNU_HC,FF_BSDSSDS1STAR_HC,INTEGRATION_INDEX})});
63 
64  addProcessSignature(PID::BS,{PID::DSSDS1STARMINUS,PID::NU_E,PID::POSITRON});
65  addIntegrationBoundaries({PS::makeQ2Function(pdg.getMass(PID::POSITRON), pdg.getMass(PID::BS)-pdg.getMass(PID::DSSDS1STARMINUS))});
66  addTensor(Tensor{name, MD::makeEmptySparse(dims, {WILSON_BCENU,FF_BSDSSDS1STAR,WILSON_BCENU_HC,FF_BSDSSDS1STAR_HC,INTEGRATION_INDEX})});
67 
68  setSignatureIndex();
69  }
70 
71  Tensor RateBD1starLepNu::evalAtPSPoint(const vector<double>& point) {
72  auto labs = getTensor().labels();
73  labs.pop_back();
74  auto dimensions = getTensor().dims();
75  dimensions.pop_back();
76  Tensor result{"RateBD1starLepNu", MD::makeEmptySparse(dimensions, labs)};
77 
78  const double Mb = masses()[0];
79  const double Mc = masses()[1];
80  const double Mt = masses()[3];
81  const double Mb2 = Mb*Mb;
82  const double Mc2 = Mc*Mc;
83  const double Mt2 = Mt*Mt;
84 
85  //kinematic objects
86  const double Sqq = point[0];
87  // const double Sqq2 = Sqq*Sqq;
88  const double mSqq = Sqq/Mb2;
89  const double w = (Mb2 + Mc2 - Sqq)/(2 * Mb * Mc);
90  const double rC = Mc/Mb;
91  const double rt = Mt/Mb;
92 
93  const double mSqqSq = mSqq*mSqq;
94  const double wSq = w*w;
95  const double rtSq = rt*rt;
96  const double rCSq = rC*rC;
97  // const double rCp1Sq = (rC + 1.)*(rC + 1.);
98  // const double rCm1Sq = (rC - 1.)*(rC - 1.);
99  // const double wp1Sq = (w + 1.)*(w + 1.);
100  // const double wm1Sq = (w - 1.)*(w - 1.);
101  const double wSqm1Sq = (wSq - 1.)*(wSq - 1.);
102  // const double wp1Cu = wp1Sq*(w + 1.);
103  // const double wm1Cu = wm1Sq*(w - 1.);
104  // const double wSqm1Cu = wSqm1Sq*(wSq - 1.);
105 
106  double RateNorm = (GFermi*GFermi*Mb2*Mc*rCSq*(-Mt2 + Sqq)*(1 + rCSq - rtSq - 2*rC*w)*sqrt(w*w - 1.))/(32.*pow(pi,3.)*(1 + rCSq - 2*rC*w));
107 
108  //set non-zero tensor elements
109  result.element({0,1,0,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
110  result.element({0,1,0,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
111  result.element({0,1,0,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
112  result.element({0,1,1,0}) = (rt*(-1 + wSq))/mSqq;
113  result.element({0,1,3,0}) = (rt*(-1 + wSq))/mSqq;
114  result.element({0,1,5,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
115  result.element({0,1,5,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
116  result.element({0,1,5,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
117  result.element({0,1,7,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
118  result.element({0,1,7,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
119  result.element({0,1,7,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
120  result.element({0,1,9,5}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
121  result.element({0,1,9,6}) = (12*rt*(rC - w))/mSqq;
122  result.element({0,1,9,7}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
123  result.element({0,2,0,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
124  result.element({0,2,0,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
125  result.element({0,2,0,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
126  result.element({0,2,1,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
127  result.element({0,2,3,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
128  result.element({0,2,5,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
129  result.element({0,2,5,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
130  result.element({0,2,5,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
131  result.element({0,2,7,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
132  result.element({0,2,7,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
133  result.element({0,2,7,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
134  result.element({0,2,9,5}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
135  result.element({0,2,9,6}) = (-4*rC*rt*(-1 + wSq))/mSqq;
136  result.element({0,2,9,7}) = (4*rC*rt*wSqm1Sq)/mSqq;
137  result.element({0,3,0,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
138  result.element({0,3,0,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
139  result.element({0,3,0,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
140  result.element({0,3,1,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
141  result.element({0,3,3,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
142  result.element({0,3,5,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
143  result.element({0,3,5,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
144  result.element({0,3,5,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
145  result.element({0,3,7,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
146  result.element({0,3,7,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
147  result.element({0,3,7,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
148  result.element({0,3,9,5}) = (-4*rt*w*(-1 + wSq))/mSqq;
149  result.element({0,3,9,6}) = (-4*rt*(-1 + wSq))/mSqq;
150  result.element({0,3,9,7}) = (4*rt*wSqm1Sq)/mSqq;
151  result.element({0,4,0,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
152  result.element({0,4,5,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
153  result.element({0,4,7,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
154  result.element({0,4,9,5}) = (8*rC*rt*(-1 + wSq))/mSqq;
155  result.element({0,4,9,6}) = (8*rt*(-1 + wSq))/mSqq;
156  result.element({1,0,0,1}) = (rt*(-1 + wSq))/mSqq;
157  result.element({1,0,0,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
158  result.element({1,0,0,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
159  result.element({1,0,1,0}) = -1 + wSq;
160  result.element({1,0,3,0}) = -1 + wSq;
161  result.element({1,0,5,1}) = (rt*(-1 + wSq))/mSqq;
162  result.element({1,0,5,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
163  result.element({1,0,5,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
164  result.element({1,0,7,1}) = (rt*(-1 + wSq))/mSqq;
165  result.element({1,0,7,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
166  result.element({1,0,7,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
167  result.element({2,0,2,0}) = -1 + wSq;
168  result.element({2,0,4,0}) = -1 + wSq;
169  result.element({2,0,6,1}) = (rt*(-1 + wSq))/mSqq;
170  result.element({2,0,6,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
171  result.element({2,0,6,3}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
172  result.element({2,0,8,1}) = (rt*(-1 + wSq))/mSqq;
173  result.element({2,0,8,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
174  result.element({2,0,8,3}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
175  result.element({3,0,0,1}) = (rt*(-1 + wSq))/mSqq;
176  result.element({3,0,0,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
177  result.element({3,0,0,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
178  result.element({3,0,1,0}) = -1 + wSq;
179  result.element({3,0,3,0}) = -1 + wSq;
180  result.element({3,0,5,1}) = (rt*(-1 + wSq))/mSqq;
181  result.element({3,0,5,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
182  result.element({3,0,5,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
183  result.element({3,0,7,1}) = (rt*(-1 + wSq))/mSqq;
184  result.element({3,0,7,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
185  result.element({3,0,7,3}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
186  result.element({4,0,2,0}) = -1 + wSq;
187  result.element({4,0,4,0}) = -1 + wSq;
188  result.element({4,0,6,1}) = (rt*(-1 + wSq))/mSqq;
189  result.element({4,0,6,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
190  result.element({4,0,6,3}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
191  result.element({4,0,8,1}) = (rt*(-1 + wSq))/mSqq;
192  result.element({4,0,8,2}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
193  result.element({4,0,8,3}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
194  result.element({5,1,0,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
195  result.element({5,1,0,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
196  result.element({5,1,0,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
197  result.element({5,1,1,0}) = (rt*(-1 + wSq))/mSqq;
198  result.element({5,1,3,0}) = (rt*(-1 + wSq))/mSqq;
199  result.element({5,1,5,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
200  result.element({5,1,5,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
201  result.element({5,1,5,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
202  result.element({5,1,7,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
203  result.element({5,1,7,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
204  result.element({5,1,7,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
205  result.element({5,1,9,5}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
206  result.element({5,1,9,6}) = (12*rt*(rC - w))/mSqq;
207  result.element({5,1,9,7}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
208  result.element({5,2,0,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
209  result.element({5,2,0,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
210  result.element({5,2,0,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
211  result.element({5,2,1,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
212  result.element({5,2,3,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
213  result.element({5,2,5,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
214  result.element({5,2,5,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
215  result.element({5,2,5,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
216  result.element({5,2,7,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
217  result.element({5,2,7,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
218  result.element({5,2,7,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
219  result.element({5,2,9,5}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
220  result.element({5,2,9,6}) = (-4*rC*rt*(-1 + wSq))/mSqq;
221  result.element({5,2,9,7}) = (4*rC*rt*wSqm1Sq)/mSqq;
222  result.element({5,3,0,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
223  result.element({5,3,0,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
224  result.element({5,3,0,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
225  result.element({5,3,1,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
226  result.element({5,3,3,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
227  result.element({5,3,5,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
228  result.element({5,3,5,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
229  result.element({5,3,5,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
230  result.element({5,3,7,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
231  result.element({5,3,7,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
232  result.element({5,3,7,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
233  result.element({5,3,9,5}) = (-4*rt*w*(-1 + wSq))/mSqq;
234  result.element({5,3,9,6}) = (-4*rt*(-1 + wSq))/mSqq;
235  result.element({5,3,9,7}) = (4*rt*wSqm1Sq)/mSqq;
236  result.element({5,4,0,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
237  result.element({5,4,5,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
238  result.element({5,4,7,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
239  result.element({5,4,9,5}) = (-8*rC*rt*(-1 + wSq))/mSqq;
240  result.element({5,4,9,6}) = (-8*rt*(-1 + wSq))/mSqq;
241  result.element({6,1,2,0}) = (rt*(-1 + wSq))/mSqq;
242  result.element({6,1,4,0}) = (rt*(-1 + wSq))/mSqq;
243  result.element({6,1,6,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
244  result.element({6,1,6,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
245  result.element({6,1,6,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
246  result.element({6,1,8,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
247  result.element({6,1,8,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
248  result.element({6,1,8,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
249  result.element({6,1,10,5}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
250  result.element({6,1,10,6}) = (12*rt*(rC - w))/mSqq;
251  result.element({6,1,10,7}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
252  result.element({6,2,2,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
253  result.element({6,2,4,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
254  result.element({6,2,6,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
255  result.element({6,2,6,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
256  result.element({6,2,6,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
257  result.element({6,2,8,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
258  result.element({6,2,8,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
259  result.element({6,2,8,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
260  result.element({6,2,10,5}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
261  result.element({6,2,10,6}) = (-4*rC*rt*(-1 + wSq))/mSqq;
262  result.element({6,2,10,7}) = (4*rC*rt*wSqm1Sq)/mSqq;
263  result.element({6,3,2,0}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
264  result.element({6,3,4,0}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
265  result.element({6,3,6,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
266  result.element({6,3,6,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
267  result.element({6,3,6,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
268  result.element({6,3,8,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
269  result.element({6,3,8,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
270  result.element({6,3,8,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
271  result.element({6,3,10,5}) = (-4*rt*w*(-1 + wSq))/mSqq;
272  result.element({6,3,10,6}) = (-4*rt*(-1 + wSq))/mSqq;
273  result.element({6,3,10,7}) = (4*rt*wSqm1Sq)/mSqq;
274  result.element({6,4,6,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
275  result.element({6,4,8,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
276  result.element({6,4,10,5}) = (8*rC*rt*(-1 + wSq))/mSqq;
277  result.element({6,4,10,6}) = (8*rt*(-1 + wSq))/mSqq;
278  result.element({7,1,0,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
279  result.element({7,1,0,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
280  result.element({7,1,0,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
281  result.element({7,1,1,0}) = (rt*(-1 + wSq))/mSqq;
282  result.element({7,1,3,0}) = (rt*(-1 + wSq))/mSqq;
283  result.element({7,1,5,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
284  result.element({7,1,5,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
285  result.element({7,1,5,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
286  result.element({7,1,7,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
287  result.element({7,1,7,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
288  result.element({7,1,7,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
289  result.element({7,1,9,5}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
290  result.element({7,1,9,6}) = (12*rt*(rC - w))/mSqq;
291  result.element({7,1,9,7}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
292  result.element({7,2,0,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
293  result.element({7,2,0,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
294  result.element({7,2,0,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
295  result.element({7,2,1,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
296  result.element({7,2,3,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
297  result.element({7,2,5,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
298  result.element({7,2,5,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
299  result.element({7,2,5,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
300  result.element({7,2,7,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
301  result.element({7,2,7,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
302  result.element({7,2,7,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
303  result.element({7,2,9,5}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
304  result.element({7,2,9,6}) = (-4*rC*rt*(-1 + wSq))/mSqq;
305  result.element({7,2,9,7}) = (4*rC*rt*wSqm1Sq)/mSqq;
306  result.element({7,3,0,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
307  result.element({7,3,0,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
308  result.element({7,3,0,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
309  result.element({7,3,1,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
310  result.element({7,3,3,0}) = (rt*(-rC + w)*(-1 + wSq))/mSqq;
311  result.element({7,3,5,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
312  result.element({7,3,5,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
313  result.element({7,3,5,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
314  result.element({7,3,7,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
315  result.element({7,3,7,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
316  result.element({7,3,7,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
317  result.element({7,3,9,5}) = (-4*rt*w*(-1 + wSq))/mSqq;
318  result.element({7,3,9,6}) = (-4*rt*(-1 + wSq))/mSqq;
319  result.element({7,3,9,7}) = (4*rt*wSqm1Sq)/mSqq;
320  result.element({7,4,0,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
321  result.element({7,4,5,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
322  result.element({7,4,7,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
323  result.element({7,4,9,5}) = (8*rC*rt*(-1 + wSq))/mSqq;
324  result.element({7,4,9,6}) = (8*rt*(-1 + wSq))/mSqq;
325  result.element({8,1,2,0}) = (rt*(-1 + wSq))/mSqq;
326  result.element({8,1,4,0}) = (rt*(-1 + wSq))/mSqq;
327  result.element({8,1,6,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
328  result.element({8,1,6,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
329  result.element({8,1,6,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
330  result.element({8,1,8,1}) = (4*mSqqSq + 2*mSqq*(rCSq + rtSq - 2*rC*w + wSq) + rtSq*(-3 + rCSq - 2*rC*w + 4*wSq))/(3.*mSqqSq);
331  result.element({8,1,8,2}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
332  result.element({8,1,8,3}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
333  result.element({8,1,10,5}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
334  result.element({8,1,10,6}) = (12*rt*(rC - w))/mSqq;
335  result.element({8,1,10,7}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
336  result.element({8,2,2,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
337  result.element({8,2,4,0}) = -((rt*(-1 + rC*w)*(-1 + wSq))/mSqq);
338  result.element({8,2,6,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
339  result.element({8,2,6,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
340  result.element({8,2,6,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
341  result.element({8,2,8,1}) = -((-1 + wSq)*(2*mSqq*rC*(rC - w) + rtSq*(-3 + rCSq + 2*rC*w)))/(3.*mSqqSq);
342  result.element({8,2,8,2}) = ((-1 + wSq)*(2*mSqq*rCSq*(-1 + wSq) + rtSq*(3 - 6*rC*w + rCSq*(-1 + 4*wSq))))/(3.*mSqqSq);
343  result.element({8,2,8,3}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
344  result.element({8,2,10,5}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
345  result.element({8,2,10,6}) = (-4*rC*rt*(-1 + wSq))/mSqq;
346  result.element({8,2,10,7}) = (4*rC*rt*wSqm1Sq)/mSqq;
347  result.element({8,3,2,0}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
348  result.element({8,3,4,0}) = -((rt*(rC - w)*(-1 + wSq))/mSqq);
349  result.element({8,3,6,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
350  result.element({8,3,6,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
351  result.element({8,3,6,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
352  result.element({8,3,8,1}) = (-2*(mSqq + 2*rtSq)*(rC - w)*(-1 + wSq))/(3.*mSqqSq);
353  result.element({8,3,8,2}) = ((-1 + wSq)*(2*mSqq*rC*(-1 + wSq) + rtSq*(3*w + 3*rCSq*w - 2*rC*(2 + wSq))))/(3.*mSqqSq);
354  result.element({8,3,8,3}) = ((-1 + wSq)*(2*mSqq*(-1 + wSq) + rtSq*(-1 + 3*rCSq - 6*rC*w + 4*wSq)))/(3.*mSqqSq);
355  result.element({8,3,10,5}) = (-4*rt*w*(-1 + wSq))/mSqq;
356  result.element({8,3,10,6}) = (-4*rt*(-1 + wSq))/mSqq;
357  result.element({8,3,10,7}) = (4*rt*wSqm1Sq)/mSqq;
358  result.element({8,4,6,4}) = (-2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
359  result.element({8,4,8,4}) = (2*(2*mSqq + rtSq)*(-1 + wSq))/(3.*mSqq);
360  result.element({8,4,10,5}) = (-8*rC*rt*(-1 + wSq))/mSqq;
361  result.element({8,4,10,6}) = (-8*rt*(-1 + wSq))/mSqq;
362  result.element({9,5,0,1}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
363  result.element({9,5,0,2}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
364  result.element({9,5,0,3}) = (-4*rt*w*(-1 + wSq))/mSqq;
365  result.element({9,5,0,4}) = (8*rC*rt*(-1 + wSq))/mSqq;
366  result.element({9,5,5,1}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
367  result.element({9,5,5,2}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
368  result.element({9,5,5,3}) = (-4*rt*w*(-1 + wSq))/mSqq;
369  result.element({9,5,5,4}) = (-8*rC*rt*(-1 + wSq))/mSqq;
370  result.element({9,5,7,1}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
371  result.element({9,5,7,2}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
372  result.element({9,5,7,3}) = (-4*rt*w*(-1 + wSq))/mSqq;
373  result.element({9,5,7,4}) = (8*rC*rt*(-1 + wSq))/mSqq;
374  result.element({9,5,9,5}) = (16*(mSqqSq*wSq + 4*rtSq*(1 - 2*rC*w + rCSq*(-1 + 2*wSq)) + 2*mSqq*(1 - 2*rC*w + rtSq*wSq + rCSq*(-1 + 2*wSq))))/(3.*mSqqSq);
375  result.element({9,5,9,6}) = (16*(mSqq + 2*rtSq)*(-4*rC + (2 + mSqq)*w + 2*rCSq*w))/(3.*mSqqSq);
376  result.element({9,5,9,7}) = (-16*(mSqq + 2*rtSq)*w*(-1 + wSq))/(3.*mSqq);
377  result.element({9,6,0,1}) = (12*rt*(rC - w))/mSqq;
378  result.element({9,6,0,2}) = (-4*rC*rt*(-1 + wSq))/mSqq;
379  result.element({9,6,0,3}) = (-4*rt*(-1 + wSq))/mSqq;
380  result.element({9,6,0,4}) = (8*rt*(-1 + wSq))/mSqq;
381  result.element({9,6,5,1}) = (12*rt*(rC - w))/mSqq;
382  result.element({9,6,5,2}) = (-4*rC*rt*(-1 + wSq))/mSqq;
383  result.element({9,6,5,3}) = (-4*rt*(-1 + wSq))/mSqq;
384  result.element({9,6,5,4}) = (-8*rt*(-1 + wSq))/mSqq;
385  result.element({9,6,7,1}) = (12*rt*(rC - w))/mSqq;
386  result.element({9,6,7,2}) = (-4*rC*rt*(-1 + wSq))/mSqq;
387  result.element({9,6,7,3}) = (-4*rt*(-1 + wSq))/mSqq;
388  result.element({9,6,7,4}) = (8*rt*(-1 + wSq))/mSqq;
389  result.element({9,6,9,5}) = (16*(mSqq + 2*rtSq)*(-4*rC + (2 + mSqq)*w + 2*rCSq*w))/(3.*mSqqSq);
390  result.element({9,6,9,6}) = (16*(mSqqSq + 4*rtSq*(-1 + rCSq - 2*rC*w + 2*wSq) + 2*mSqq*(-1 + rCSq + rtSq - 2*rC*w + 2*wSq)))/(3.*mSqqSq);
391  result.element({9,6,9,7}) = (-16*(mSqq + 2*rtSq)*(-1 + wSq))/(3.*mSqq);
392  result.element({9,7,0,1}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
393  result.element({9,7,0,2}) = (4*rC*rt*wSqm1Sq)/mSqq;
394  result.element({9,7,0,3}) = (4*rt*wSqm1Sq)/mSqq;
395  result.element({9,7,5,1}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
396  result.element({9,7,5,2}) = (4*rC*rt*wSqm1Sq)/mSqq;
397  result.element({9,7,5,3}) = (4*rt*wSqm1Sq)/mSqq;
398  result.element({9,7,7,1}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
399  result.element({9,7,7,2}) = (4*rC*rt*wSqm1Sq)/mSqq;
400  result.element({9,7,7,3}) = (4*rt*wSqm1Sq)/mSqq;
401  result.element({9,7,9,5}) = (-16*(mSqq + 2*rtSq)*w*(-1 + wSq))/(3.*mSqq);
402  result.element({9,7,9,6}) = (-16*(mSqq + 2*rtSq)*(-1 + wSq))/(3.*mSqq);
403  result.element({9,7,9,7}) = (16*(mSqq + 2*rtSq)*wSqm1Sq)/(3.*mSqq);
404  result.element({10,5,6,1}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
405  result.element({10,5,6,2}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
406  result.element({10,5,6,3}) = (-4*rt*w*(-1 + wSq))/mSqq;
407  result.element({10,5,6,4}) = (8*rC*rt*(-1 + wSq))/mSqq;
408  result.element({10,5,8,1}) = (4*rt*(-2 + 3*rC*w - wSq))/mSqq;
409  result.element({10,5,8,2}) = (-4*rC*rt*w*(-1 + wSq))/mSqq;
410  result.element({10,5,8,3}) = (-4*rt*w*(-1 + wSq))/mSqq;
411  result.element({10,5,8,4}) = (-8*rC*rt*(-1 + wSq))/mSqq;
412  result.element({10,5,10,5}) = (16*(mSqqSq*wSq + 4*rtSq*(1 - 2*rC*w + rCSq*(-1 + 2*wSq)) + 2*mSqq*(1 - 2*rC*w + rtSq*wSq + rCSq*(-1 + 2*wSq))))/(3.*mSqqSq);
413  result.element({10,5,10,6}) = (16*(mSqq + 2*rtSq)*(-4*rC + (2 + mSqq)*w + 2*rCSq*w))/(3.*mSqqSq);
414  result.element({10,5,10,7}) = (-16*(mSqq + 2*rtSq)*w*(-1 + wSq))/(3.*mSqq);
415  result.element({10,6,6,1}) = (12*rt*(rC - w))/mSqq;
416  result.element({10,6,6,2}) = (-4*rC*rt*(-1 + wSq))/mSqq;
417  result.element({10,6,6,3}) = (-4*rt*(-1 + wSq))/mSqq;
418  result.element({10,6,6,4}) = (8*rt*(-1 + wSq))/mSqq;
419  result.element({10,6,8,1}) = (12*rt*(rC - w))/mSqq;
420  result.element({10,6,8,2}) = (-4*rC*rt*(-1 + wSq))/mSqq;
421  result.element({10,6,8,3}) = (-4*rt*(-1 + wSq))/mSqq;
422  result.element({10,6,8,4}) = (-8*rt*(-1 + wSq))/mSqq;
423  result.element({10,6,10,5}) = (16*(mSqq + 2*rtSq)*(-4*rC + (2 + mSqq)*w + 2*rCSq*w))/(3.*mSqqSq);
424  result.element({10,6,10,6}) = (16*(mSqqSq + 4*rtSq*(-1 + rCSq - 2*rC*w + 2*wSq) + 2*mSqq*(-1 + rCSq + rtSq - 2*rC*w + 2*wSq)))/(3.*mSqqSq);
425  result.element({10,6,10,7}) = (-16*(mSqq + 2*rtSq)*(-1 + wSq))/(3.*mSqq);
426  result.element({10,7,6,1}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
427  result.element({10,7,6,2}) = (4*rC*rt*wSqm1Sq)/mSqq;
428  result.element({10,7,6,3}) = (4*rt*wSqm1Sq)/mSqq;
429  result.element({10,7,8,1}) = (-4*rt*(rC - w)*(-1 + wSq))/mSqq;
430  result.element({10,7,8,2}) = (4*rC*rt*wSqm1Sq)/mSqq;
431  result.element({10,7,8,3}) = (4*rt*wSqm1Sq)/mSqq;
432  result.element({10,7,10,5}) = (-16*(mSqq + 2*rtSq)*w*(-1 + wSq))/(3.*mSqq);
433  result.element({10,7,10,6}) = (-16*(mSqq + 2*rtSq)*(-1 + wSq))/(3.*mSqq);
434  result.element({10,7,10,7}) = (16*(mSqq + 2*rtSq)*wSqm1Sq)/(3.*mSqq);
435 
436  result*=(RateNorm);
437 
438  return result;
439 
440  }
441 
442 }
Hammer::MultiDimensional::makeEmptySparse
TensorData makeEmptySparse(const IndexList &dimensions, const LabelsList &labels)
Definition: SparseContainer.cc:287
Hammer::Tensor::labels
LabelsList labels() const
get the labels of all the indices at once
Definition: Tensor.cc:83
Hammer::GFermi
static constexpr double GFermi
Definition: Constants.hh:29
IndexLabels.hh
Tensor indices label definitions.
SparseContainer.hh
Sparse tensor data container.
Hammer::Tensor
Multidimensional tensor class with complex numbers as elements.
Definition: Tensor.hh:33
Constants.hh
Various numerical constants.
Pdg.hh
Hammer particle data class.
Hammer::PS::makeQ2Function
BoundaryFunction makeQ2Function(double qmin, double qmax)
Definition: RateBase.hh:122
Particle.hh
Hammer particle class.
Hammer::pi
static constexpr double pi
Definition: Constants.hh:21