12 int num = states.size();
15 double mfz = fieldValue.
z();
18 for (std::vector<KinematicState>::const_iterator
i = states.begin();
i != states.end();
i++) {
19 mom[
j] =
i->globalMomentum();
21 a_i[
j] = -
i->particleCharge() * mfz;
25 double pvt2 = pvx * pvx + pvy * pvy;
38 for (
int j = 0;
j != 2; ++
j) {
55 ROOT::Math::SMatrix<double, 2, 7> el_part_d;
56 for (
int j = 0;
j != 2; ++
j) {
65 el_part_d(1, 2) = -1.;
66 el_part_d(0, 3) =
dpos[
j].
y();
67 el_part_d(0, 4) = -
dpos[
j].
x();
75 el_part_d(0, 0) =
mom[
j].
y();
76 el_part_d(0, 1) = -
mom[
j].
x();
77 el_part_d(1, 0) =
mom[
j].
x() * (
mom[
j].
z() * pt2Inverse);
78 el_part_d(1, 1) =
mom[
j].
y() * (
mom[
j].
z() * pt2Inverse);
79 el_part_d(1, 2) = -1.;
80 el_part_d(0, 3) =
dpos[
j].
y();
81 el_part_d(0, 4) = -
dpos[
j].
x();
83 (
mom[
j].
z() * pt2Inverse) -
86 (
mom[
j].
z() * pt2Inverse) -
95 ROOT::Math::SMatrix<double, 2, 3> el_part_e;
96 for (
int j = 0;
j != 2; ++
j) {
110 el_part_e(0, 0) = -
mom[
j].
y();
111 el_part_e(0, 1) =
mom[
j].
x();
112 el_part_e(1, 0) = -
mom[
j].
x() *
mom[
j].
z() * pt2Inverse;
113 el_part_e(1, 1) = -
mom[
j].
y() *
mom[
j].
z() * pt2Inverse;