4449   const double s_sqrt = sqrt(s);
  4450   const double spin_deg_factor = 1.0;
  4453   double tmax = mandelstam_t[0];
  4454   double tmin = mandelstam_t[1];
  4460           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4463                  2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4466                      (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  4467                       4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  4469                      (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4470                       8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4474                  2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  4475                  1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4477                      (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  4479                      (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  4480                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4481                       2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4484                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4485                       pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
  4489                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  4491                      (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  4492                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4493                       4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4496                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4497                       pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  4498                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  4501                       2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4502                       2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  4504                           (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
  4505               (1. * pow(
a1_mass, 2) - 1. * tmin) +
  4507            (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  4510            (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  4511               (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
  4517                0.5 * 
eta2 * pow(m_rho, 2) + 0.5 * 
eta1 * s - 1. * 
eta2 * s) *
  4521                 (12. + 1. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  4522                  delta * (-8. + 8. * 
C4 * pow(m_rho, 2))) +
  4523             (-4. - 3. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  4524              delta * (8. + 8. * 
C4 * pow(m_rho, 2))) *
  4540                        0.5 * pow(m_rho, 2) + 1. * s)) *
  4542           (0.25 * (-2. + 1. * 
delta) *
  4543            (-8. * 
C4 * pow(m_rho, 4) +
  4545             (-2. - 3. * 
delta) * s +
  4546             pow(m_rho, 2) * (2. + 1. * 
delta + 16. * 
C4 * s)) *
  4550            (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
  4551             8 * 
C4 * pow(m_rho, 6) * (-6 + 
delta - 8 * 
C4 * s) +
  4552             2 * 
delta * pow(m_rho, 2) * s * (-6 + 
delta - 8 * 
C4 * s) +
  4553             pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 * 
C4 * (6 + 
delta) * s +
  4554                              32 * pow(
C4, 2) * pow(s, 2))) *
  4559                 (0.75 * pow(m_rho, 4) - 0.125 * 
delta * pow(m_rho, 4) -
  4560                  1. * 
C4 * pow(m_rho, 6) +
  4562                      (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4564                      (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) -
  4565                  0.25 * pow(m_rho, 2) * s - 0.375 * 
delta * pow(m_rho, 2) * s +
  4566                  2. * 
C4 * pow(m_rho, 4) * s + 0.25 * 
delta * pow(s, 2) -
  4567                  1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4568             eta1 * (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  4570                         (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) +
  4572                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4573                     0.25 * 
delta * pow(s, 2) +
  4574                     1. * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4577           0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  4580                                2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4581                                2. * pow(m_rho, 4) +
  4583                                                   pow(m_rho, 2) - 1. * s) -
  4584                                3. * pow(m_rho, 2) * s + pow(s, 2)) +
  4587                                2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4589                                                   pow(m_rho, 2) - 1. * s) +
  4590                                pow(m_rho, 2) * s + pow(s, 2)) +
  4594                     4. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4595                     2. * pow(m_rho, 4) - 2. * pow(s, 2) +
  4597                                        2. * pow(m_rho, 2) + 2. * s))) *
  4599           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4602                     2. * pow(m_rho, 4) +
  4604                     8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4609                     pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4610                     4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4612                                        4. * pow(m_rho, 2) + 4. * s)) +
  4616                     pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4617                     2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4619                                        4. * pow(m_rho, 2) + 4. * s))) *
  4621           (0.125 * (-2. + 1. * 
delta) *
  4622            (2. + 1. * 
delta - 8. * 
C4 * pow(m_rho, 2)) * pow(tmin, 2)) /
  4624           0.5 * pow(1. * 
eta1 - 1. * 
eta2, 2) *
  4625               (-0.5 + 1. * 
C4 * pow(m_rho, 2)) * pow(tmin, 2) -
  4627            (0.5 - 0.125 * pow(
delta, 2) - 2. * 
C4 * pow(m_rho, 2) +
  4628             1. * 
C4 * 
delta * pow(m_rho, 2)) *
  4631           0.0625 * pow(1. * 
eta1 - 1. * 
eta2, 4) *
  4632               (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
  4634           0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  4636                        1. * pow(m_rho, 2) - 1. * s) +
  4638                        1. * pow(m_rho, 2) + s)) *
  4640           0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) * pow(tmin, 3) -
  4641           0.020833333333333332 * pow(1. * 
eta1 - 1. * 
eta2, 4) * pow(tmin, 3) +
  4642           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4646                     2. * pow(m_rho, 4) +
  4648                     4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
  4652                     2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  4653                     4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4655                                        4. * pow(m_rho, 2) + 4. * s)) +
  4659                     2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  4660                     2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4662                                        4. * pow(m_rho, 2) + 4. * s))) *
  4663               (-1. * pow(m_rho, 2) + s + tmin) -
  4664           0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  4665               (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
  4666                eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
  4667               pow(-1. * pow(m_rho, 2) + s + tmin, 2) +
  4668           0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) *
  4669               pow(-1. * pow(m_rho, 2) + s + tmin, 3) +
  4671               (-1. + 2. * 
C4 * pow(m_rho, 2)) *
  4673                       1. * pow(m_rho, 2) + s + tmin,
  4676            (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  4677                     0.5 * 
C4 * pow(m_rho, 6) +
  4679                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4681                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4682                     0.125 * pow(m_rho, 2) * s -
  4683                     0.1875 * 
delta * pow(m_rho, 2) * s +
  4684                     1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  4685                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4686             eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  4688                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4690                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4691                     0.125 * 
delta * pow(s, 2) +
  4692                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4694             1. * s + 1. * tmin)) /
  4697            (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  4698                     0.5 * 
C4 * pow(m_rho, 6) +
  4700                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4702                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4703                     0.125 * pow(m_rho, 2) * s -
  4704                     0.1875 * 
delta * pow(m_rho, 2) * s +
  4705                     1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  4706                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4707             eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  4709                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4711                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4712                     0.125 * 
delta * pow(s, 2) +
  4713                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4722                      (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  4724                      (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  4725                  pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
  4727                      (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
  4734                         (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  4735                          pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
  4736                          1. * pow(m_rho, 2) * s))) *
  4745                  2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4746                  1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  4748                      (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
  4750                  3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4751                  0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
  4752                  1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  4756                      (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  4757                       1. * pow(m_rho, 2) * s +
  4758                       pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
  4761                  4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  4762                  1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  4767                  4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  4768                  0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
  4769                  1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
  4772                       (1. * pow(m_rho, 2) - 1. * s) * s +
  4773                       pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
  4779                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  4780                2. * pow(m_rho, 2) * s + pow(s, 2) +
  4782                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
  4783           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4789                      (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  4791                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4792                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4793                       2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4796                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4797                       pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
  4800                       1. * pow(m_rho, 4) +
  4802                                          6. * pow(m_rho, 2) - 6. * s) +
  4803                       2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  4804                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
  4808                  2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4811                      (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  4812                       4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  4814                      (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4815                       8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4819                       2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
  4820                       4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
  4825                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  4827                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4828                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4829                       4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4832                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4833                       pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  4834                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  4837                       1. * pow(m_rho, 4) +
  4839                                          6. * pow(m_rho, 2) - 6. * s) +
  4840                       4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  4841                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
  4844                       2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4845                       2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  4847                           (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
  4856                  28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  4857                  20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4858                  10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
  4860                  16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  4861                  12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  4862                  4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
  4863                  6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  4865                  4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
  4867                                     28. * pow(m_rho, 2) + 28. * s) +
  4869                      (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
  4870                       pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
  4871                       36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
  4875                                          4. * pow(m_rho, 2) - 4. * s) +
  4876                       (4. * pow(m_rho, 2) - 4. * s) * s +
  4877                       pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
  4879                      (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
  4880                       18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
  4882                       pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
  4884                           (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
  4885                            12. * pow(s, 2)))) +
  4889                  14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  4890                  7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
  4891                  2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
  4893                      (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  4895                  11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
  4896                  6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  4897                  1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
  4898                  1. * pow(m_rho, 4) * pow(s, 2) -
  4900                  1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  4903                       1. * pow(m_rho, 4) +
  4904                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  4905                       1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4907                                          2. * pow(m_rho, 2) + 2. * s)) +
  4909                      (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
  4910                       21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  4911                       pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  4913                      (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  4914                       pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
  4915                       9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
  4918                           (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
  4923                  14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  4924                  1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4925                  6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
  4927                      (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  4929                  1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  4930                  16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  4931                  7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
  4932                  14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  4933                  9. * pow(m_rho, 4) * pow(s, 2) -
  4935                  5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  4938                       3. * pow(m_rho, 4) +
  4939                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  4940                       5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4942                                          2. * pow(m_rho, 2) + 2. * s)) +
  4944                      (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
  4945                       9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  4946                       pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  4948                      (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  4949                       pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
  4950                       9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
  4953                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  4954                            6. * pow(s, 2))))) *
  4960                pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  4961                2. * pow(m_rho, 2) * s + pow(s, 2) +
  4963                    (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
  4964           0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  4969                         (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
  4970                          2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
  4972                         (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4973                          8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4978                     pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
  4979                     pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
  4981                                        3. * pow(m_rho, 2) + 3. * s) +
  4983                         (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4984                          pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4985                          4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
  4988                     1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4989                     pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
  4991                                        3. * pow(m_rho, 2) + 3. * s) +
  4993                         (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4994                          pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4995                          2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
  4996               log(abs(-1. * pow(
a1_mass, 2) + tmin)) -
  5000                  0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5002                      (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5004                      (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
  5007                         (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5008                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  5009                          0.5 * pow(m_rho, 2) * s) +
  5011                         (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
  5012                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5013            log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
  5018                         (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
  5019                     pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
  5020                     pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
  5023                          pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5024                          (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
  5027                     pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
  5028                     (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
  5029                     pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
  5031                         (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
  5032                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5033            log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
  5035                1. * pow(m_rho, 2) + 1. * s) -
  5036           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5038             0.5 * pow(m_rho, 2) + 0.5 * s) *
  5042                      (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
  5044                      (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
  5045                       4. * pow(m_rho, 2) * s +
  5046                       pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
  5049                       4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5050                       2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
  5051                       2. * pow(m_rho, 2) * pow(s, 2) +
  5053                           (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
  5055                      (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  5056                       pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
  5057                       2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
  5060                           (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
  5065                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5067                                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  5068                                       pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5070                      (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5071                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
  5072                       pow(m_rho, 2) * s) +
  5074                      (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
  5075                       5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
  5077                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  5079                           (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
  5083                  2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5085                      (-2. * pow(m_rho, 2) + 2. * s) +
  5087                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5089                      (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
  5091                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5092                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5093                       3. * pow(m_rho, 2) * s) +
  5095                      (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
  5097                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  5099                           (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
  5100                            2. * pow(s, 2))))) *
  5101            log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
  5104                0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
  5105                0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
  5107                    (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
  5111                      (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5113                      (pow(
pion_mass, 2) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5114                       (-0.25 + 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  5116                      (-1. * 
C4 * pow(m_rho, 6) +
  5118                           (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  5119                       0.25 * 
delta * pow(s, 2) +
  5121                           (-0.25 - 0.375 * 
delta - 1. * 
C4 * s) +
  5122                       pow(m_rho, 4) * (0.75 - 0.125 * 
delta + 2. * 
C4 * s))) +
  5125                      (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  5127                      (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  5129                           (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  5130                       0.25 * 
delta * pow(s, 2) +
  5131                       1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5133                      (pow(
pion_mass, 4) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5134                       s * ((0.25 - 0.125 * 
delta) * pow(m_rho, 2) +
  5135                            (-0.25 + 0.125 * 
delta) * s) +
  5137                           (2. * 
C4 * pow(m_rho, 4) + (0.5 + 0.25 * 
delta) * s +
  5138                            pow(m_rho, 2) * (-1. - 2. * 
C4 * s))))) *
  5139            log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
  5142               log(abs(-1. * pow(
pion_mass, 2) + tmin)) +
  5145             8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  5147             5.999999999999999 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 2) *
  5149             1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
  5150            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5151               (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5153            (0. + 
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
  5154             eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
  5155                     pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
  5156            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5158           (2. * (-2. + 1. * 
delta) *
  5159            (0. + (-0.25 + 0.125 * 
delta) * pow(m_rho, 2) * s +
  5161                                  pow(m_rho, 2) * (1. + 2. * 
C4 * s))) *
  5162            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5165            (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
  5167                         (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5168                          (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  5170                         (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
  5172                     s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5174             eta2 * (0.5 * pow(m_rho, 6) +
  5175                     pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
  5176                     1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
  5179                         (1.5 * pow(m_rho, 4) - 3. * pow(m_rho, 2) * s +
  5182                         (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  5184                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5185            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5188                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5189                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5191                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  5193               log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s +
  5198                 (-1. * pow(m_rho, 2) + 1. * s) +
  5200                 (-0.5 * pow(m_rho, 4) +
  5201                  pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5202                  0.5 * pow(m_rho, 2) * s) +
  5204                 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
  5205                  pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  5206            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5210            (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5211                     (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  5212             eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  5214                     pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  5215            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5217                1. * pow(m_rho, 2) + 1. * s) -
  5220             8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  5222             1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 4) * s +
  5224                 (
C4 * (32. - 16. * 
delta) * pow(m_rho, 6) +
  5227                      (-8. + 24. * 
delta - 10. * pow(
delta, 2) + 32. * 
C4 * s -
  5229                  pow(m_rho, 4) * (-16. + 8. * 
delta - 64. * 
C4 * s +
  5231            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5232               (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5233           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5241                         (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
  5242                          2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
  5244                         (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
  5245                          8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
  5249                     3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5251                         (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
  5253                     1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  5254                     1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
  5257                          1. * pow(m_rho, 2) + s) +
  5259                         (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5260                          4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5261                          pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
  5264                     3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5266                         (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
  5271                          pow(m_rho, 2) + s) +
  5273                         (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5274                          2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5275                          pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
  5279                       2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmin) -
  5281                           (-1. * pow(m_rho, 2) + s + tmin) +
  5282                       pow(-1. * pow(m_rho, 2) + s + tmin, 2))) -
  5285                         (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5287                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  5289                         (pow(
pion_mass, 2) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5290                          (0.25 - 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  5292                         (1. * 
C4 * pow(m_rho, 6) +
  5294                              (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  5295                          0.25 * 
delta * pow(s, 2) +
  5296                          pow(m_rho, 4) * (-0.75 + 0.125 * 
delta - 2. * 
C4 * s) +
  5298                              (0.25 + 0.375 * 
delta + 1. * 
C4 * s))) +
  5300                         (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5302                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5304                         (-0.5 * pow(m_rho, 4) + 1. * 
C4 * pow(m_rho, 6) +
  5306                              (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  5307                          0.25 * 
delta * pow(s, 2) -
  5308                          1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5310                         (pow(
pion_mass, 4) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5311                          s * ((-0.25 + 0.125 * 
delta) * pow(m_rho, 2) +
  5312                               (0.25 - 0.125 * 
delta) * s) +
  5314                              (-2. * 
C4 * pow(m_rho, 4) +
  5315                               (-0.5 - 0.25 * 
delta) * s +
  5316                               pow(m_rho, 2) * (1. + 2. * 
C4 * s))))) *
  5319                            1. * pow(m_rho, 2) + s + tmin,
  5327                  0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5330                           (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
  5333                                          0.5 * pow(m_rho, 2) - 1. * s)) +
  5335                      (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5337                      (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
  5341                                     1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
  5345                       pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  5346                       0.5 * pow(m_rho, 2) * s) +
  5348                      (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5349                       pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
  5350                       0.5 * pow(m_rho, 2) * s) +
  5354                       0.5 * pow(m_rho, 2) * s +
  5355                       pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
  5358                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5359                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5360                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5361                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5363                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5369                     1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5372                     4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  5373                     1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  5375                     4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  5376                     0.5 * pow(m_rho, 4) * pow(s, 2) -
  5378                     1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  5382                          0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
  5385                          (1. * pow(m_rho, 2) - 0.5 * s) * s +
  5386                          pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
  5389                          pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
  5390                          0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
  5392                              (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
  5398                      (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
  5400                  pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
  5401                                     (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
  5403                      (-0.5 * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5406                      (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5409                      (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
  5410                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
  5411                       1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
  5413                                          1. * pow(m_rho, 2) + 1.5 * s)) +
  5416                       pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
  5418                           (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5420                       s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5421                            0.5 * pow(s, 2))))) *
  5424                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5425                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5426                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5427                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5429                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5433                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5434                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5436                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  5437           (0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  5441                      (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  5444                       0.5 * pow(m_rho, 2) + 0.5 * s) +
  5446                      (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
  5447                       pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
  5448                       6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
  5452                       (1. * pow(m_rho, 2) - 1. * s) * s +
  5453                       pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  5455                      (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  5456                       pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
  5457                       1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
  5459                           (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
  5462                       0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  5463                       1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5464                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  5466                                          6. * pow(m_rho, 2) + 6. * s) +
  5468                           (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  5471                           (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  5472                            pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  5473                            6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  5475                      (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
  5476                       3.5 * pow(m_rho, 6) * s +
  5477                       4.5 * pow(m_rho, 4) * pow(s, 2) -
  5478                       2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  5479                       pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  5481                           (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
  5484                           (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
  5485                            7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
  5489                      (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  5492                       0.5 * pow(m_rho, 2) + 0.5 * s) +
  5494                      (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
  5495                       5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
  5496                       pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  5498                      (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  5499                       pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
  5500                       3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
  5503                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
  5506                       0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  5507                       1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5508                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  5510                                          6. * pow(m_rho, 2) + 6. * s) +
  5512                           (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  5515                           (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  5516                            pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  5517                            6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  5520                       pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
  5522                           (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
  5524                           (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5527                           (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
  5528                            4.5 * pow(m_rho, 2) * pow(s, 2))) +
  5531                       pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  5533                           (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
  5536                           (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
  5537                            3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
  5538                       s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
  5539                            0.5 * pow(m_rho, 2) * pow(s, 2) +
  5540                            0.5 * pow(s, 3)))) +
  5545                       1. * pow(m_rho, 2) - 1. * s) +
  5547                      (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
  5549                      (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
  5550                       pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
  5551                       8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5553                      (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  5554                       6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
  5555                       pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
  5557                           (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
  5559                      (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
  5560                       3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
  5561                       1. * pow(m_rho, 2) * pow(s, 3) +
  5562                       pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
  5564                           (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
  5566                           (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
  5567                            6. * pow(m_rho, 2) * pow(s, 2))) +
  5569                      (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
  5570                       pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
  5571                       2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
  5574                           (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
  5577                           (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
  5581                       1. * pow(m_rho, 6) +
  5582                       pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
  5584                                          12. * pow(m_rho, 2) - 12. * s) -
  5585                       3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
  5588                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  5591                           (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
  5592                            12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  5594                                (-24. * pow(m_rho, 2) + 24. * s))))) *
  5597                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5598                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5599                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5600                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5602                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5606                pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  5607                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5609                    (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
  5613          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5616                 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5619                     (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  5620                      4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  5622                     (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  5623                      8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5627                 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5628                 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5630                     (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  5632                     (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5633                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5634                      2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5637                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5638                      pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
  5642                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5644                     (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5645                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5646                      4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5649                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5650                      pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  5651                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  5654                      2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5655                      2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  5657                          (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
  5658              (1. * pow(
a1_mass, 2) - 1. * tmax) +
  5660           (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  5663           (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  5664              (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
  5670               0.5 * 
eta2 * pow(m_rho, 2) + 0.5 * 
eta1 * s - 1. * 
eta2 * s) *
  5674                (12. + 1. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  5675                 delta * (-8. + 8. * 
C4 * pow(m_rho, 2))) +
  5676            (-4. - 3. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  5677             delta * (8. + 8. * 
C4 * pow(m_rho, 2))) *
  5692                       0.5 * pow(m_rho, 2) + 1. * s)) *
  5694          (0.25 * (-2. + 1. * 
delta) *
  5695           (-8. * 
C4 * pow(m_rho, 4) +
  5697            (-2. - 3. * 
delta) * s +
  5698            pow(m_rho, 2) * (2. + 1. * 
delta + 16. * 
C4 * s)) *
  5702           (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
  5703            8 * 
C4 * pow(m_rho, 6) * (-6 + 
delta - 8 * 
C4 * s) +
  5704            2 * 
delta * pow(m_rho, 2) * s * (-6 + 
delta - 8 * 
C4 * s) +
  5705            pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 * 
C4 * (6 + 
delta) * s +
  5706                             32 * pow(
C4, 2) * pow(s, 2))) *
  5711                (0.75 * pow(m_rho, 4) - 0.125 * 
delta * pow(m_rho, 4) -
  5712                 1. * 
C4 * pow(m_rho, 6) +
  5714                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5716                     (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) -
  5717                 0.25 * pow(m_rho, 2) * s - 0.375 * 
delta * pow(m_rho, 2) * s +
  5718                 2. * 
C4 * pow(m_rho, 4) * s + 0.25 * 
delta * pow(s, 2) -
  5719                 1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5720            eta1 * (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  5722                        (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) +
  5724                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5725                    0.25 * 
delta * pow(s, 2) +
  5726                    1. * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5729          0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  5733                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) +
  5735                        (2. * pow(
pion_mass, 2) + pow(m_rho, 2) - 1. * s) -
  5736                    3. * pow(m_rho, 2) * s + pow(s, 2)) +
  5739                               2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  5741                                                  pow(m_rho, 2) - 1. * s) +
  5742                               pow(m_rho, 2) * s + pow(s, 2)) +
  5746                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) -
  5749                                       2. * pow(m_rho, 2) + 2. * s))) *
  5751          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5754                    2. * pow(m_rho, 4) +
  5756                    8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5761                    pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5762                    4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5764                                       4. * pow(m_rho, 2) + 4. * s)) +
  5768                    pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5769                    2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5771                                       4. * pow(m_rho, 2) + 4. * s))) *
  5773          (0.125 * (-2. + 1. * 
delta) *
  5774           (2. + 1. * 
delta - 8. * 
C4 * pow(m_rho, 2)) * pow(tmax, 2)) /
  5776          0.5 * pow(1. * 
eta1 - 1. * 
eta2, 2) *
  5777              (-0.5 + 1. * 
C4 * pow(m_rho, 2)) * pow(tmax, 2) -
  5779           (0.5 - 0.125 * pow(
delta, 2) - 2. * 
C4 * pow(m_rho, 2) +
  5780            1. * 
C4 * 
delta * pow(m_rho, 2)) *
  5783          0.0625 * pow(1. * 
eta1 - 1. * 
eta2, 4) *
  5784              (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
  5786          0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  5788                       1. * pow(m_rho, 2) - 1. * s) +
  5790                       1. * pow(m_rho, 2) + s)) *
  5792          0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) * pow(tmax, 3) -
  5793          0.020833333333333332 * pow(1. * 
eta1 - 1. * 
eta2, 4) * pow(tmax, 3) +
  5794          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5798                    2. * pow(m_rho, 4) +
  5800                    4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
  5804                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5805                    4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5807                                       4. * pow(m_rho, 2) + 4. * s)) +
  5811                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5812                    2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5814                                       4. * pow(m_rho, 2) + 4. * s))) *
  5815              (-1. * pow(m_rho, 2) + s + tmax) -
  5816          0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  5817              (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
  5818               eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
  5819              pow(-1. * pow(m_rho, 2) + s + tmax, 2) +
  5820          0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) *
  5821              pow(-1. * pow(m_rho, 2) + s + tmax, 3) +
  5823              (-1. + 2. * 
C4 * pow(m_rho, 2)) *
  5829                (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  5830                 0.5 * 
C4 * pow(m_rho, 6) +
  5832                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5834                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5835                 0.125 * pow(m_rho, 2) * s - 0.1875 * 
delta * pow(m_rho, 2) * s +
  5836                 1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  5837                 0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5838            eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  5840                        (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5842                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5843                    0.125 * 
delta * pow(s, 2) +
  5844                    0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5846            1. * s + 1. * tmax)) /
  5850                (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  5851                 0.5 * 
C4 * pow(m_rho, 6) +
  5853                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5855                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5856                 0.125 * pow(m_rho, 2) * s - 0.1875 * 
delta * pow(m_rho, 2) * s +
  5857                 1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  5858                 0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5859            eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  5861                        (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5863                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5864                    0.125 * 
delta * pow(s, 2) +
  5865                    0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5874                     (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5876                     (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5877                 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
  5879                     (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
  5886                        (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5887                         pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
  5888                         1. * pow(m_rho, 2) * s))) *
  5897                 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  5899                     (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
  5901                 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5902                 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
  5903                 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5907                     (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5908                      1. * pow(m_rho, 2) * s +
  5909                      pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
  5912                 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5913                 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  5918                 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  5919                 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
  5920                 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
  5923                      (1. * pow(m_rho, 2) - 1. * s) * s +
  5924                      pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
  5930               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5931               2. * pow(m_rho, 2) * s + pow(s, 2) +
  5933                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
  5934          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5940                     (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  5942                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5943                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5944                      2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5947                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5948                      pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
  5951                      1. * pow(m_rho, 4) +
  5953                                         6. * pow(m_rho, 2) - 6. * s) +
  5954                      2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5955                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
  5959                 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5962                     (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  5963                      4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  5965                     (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  5966                      8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5970                      2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
  5971                      4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
  5976                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5978                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5979                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5980                      4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5983                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5984                      pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  5985                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  5988                      1. * pow(m_rho, 4) +
  5990                                         6. * pow(m_rho, 2) - 6. * s) +
  5991                      4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5992                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
  5995                      2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5996                      2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  5998                          (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
  6007                 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6008                 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6009                 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
  6011                 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6012                 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  6013                 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
  6014                 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6016                 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
  6018                     (-56. * pow(
pion_mass, 2) - 28. * pow(m_rho, 2) + 28. * s) +
  6020                     (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
  6021                      pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
  6022                      36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
  6026                                         4. * pow(m_rho, 2) - 4. * s) +
  6027                      (4. * pow(m_rho, 2) - 4. * s) * s +
  6028                      pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
  6030                     (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
  6031                      18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
  6033                      pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
  6035                          (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
  6036                           12. * pow(s, 2)))) +
  6040                 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  6041                 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
  6042                 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
  6044                     (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  6046                 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
  6047                 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  6048                 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
  6049                 1. * pow(m_rho, 4) * pow(s, 2) -
  6051                 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  6054                      1. * pow(m_rho, 4) +
  6055                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  6056                      1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  6058                                         2. * pow(m_rho, 2) + 2. * s)) +
  6060                     (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
  6061                      21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6062                      pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  6064                     (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  6065                      pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
  6066                      9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
  6069                          (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
  6074                 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  6075                 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6076                 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
  6078                     (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  6080                 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6081                 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  6082                 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
  6083                 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6084                 9. * pow(m_rho, 4) * pow(s, 2) -
  6086                 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  6089                      3. * pow(m_rho, 4) +
  6090                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  6091                      5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  6093                                         2. * pow(m_rho, 2) + 2. * s)) +
  6095                     (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
  6096                      9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6097                      pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  6099                     (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  6100                      pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
  6101                      9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
  6104                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  6105                           6. * pow(s, 2))))) *
  6111               pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  6112               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6114                   (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
  6115          0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  6120                        (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
  6121                         2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
  6123                        (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  6124                         8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  6129                    pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
  6130                    pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
  6132                        (-6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) + 3. * s) +
  6134                        (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  6135                         pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  6136                         4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
  6139                    1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  6140                    pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
  6142                        (-6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) + 3. * s) +
  6144                        (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  6145                         pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  6146                         2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
  6147              log(abs(-1. * pow(
a1_mass, 2) + tmax)) -
  6151                 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6153                     (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  6155                     (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
  6158                        (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  6159                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  6160                         0.5 * pow(m_rho, 2) * s) +
  6162                        (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
  6163                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6164           log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
  6169                        (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
  6170                    pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
  6171                    pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
  6174                         pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6175                         (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
  6178                    pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
  6179                    (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
  6180                    pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
  6182                        (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
  6183                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6184           log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
  6186               1. * pow(m_rho, 2) + 1. * s) -
  6187          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  6193                     (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
  6195                     (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
  6196                      4. * pow(m_rho, 2) * s +
  6197                      pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
  6200                      4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  6201                      2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
  6202                      2. * pow(m_rho, 2) * pow(s, 2) +
  6204                          (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
  6206                     (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  6207                      pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
  6208                      2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
  6211                          (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
  6216                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  6218                                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  6219                                      pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6221                     (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6222                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
  6223                      pow(m_rho, 2) * s) +
  6225                     (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
  6226                      5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
  6228                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  6230                          (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
  6234                 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6236                     (-2. * pow(m_rho, 2) + 2. * s) +
  6238                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  6240                     (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
  6242                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  6243                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  6244                      3. * pow(m_rho, 2) * s) +
  6246                     (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
  6248                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  6250                          (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
  6251                           2. * pow(s, 2))))) *
  6252           log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
  6255               0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
  6256               0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
  6258                   (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
  6262                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  6264                     (pow(
pion_mass, 2) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6265                      (-0.25 + 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  6267                     (-1. * 
C4 * pow(m_rho, 6) +
  6269                          (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  6270                      0.25 * 
delta * pow(s, 2) +
  6271                      pow(m_rho, 2) * s * (-0.25 - 0.375 * 
delta - 1. * 
C4 * s) +
  6272                      pow(m_rho, 4) * (0.75 - 0.125 * 
delta + 2. * 
C4 * s))) +
  6275                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  6277                     (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  6279                          (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  6280                      0.25 * 
delta * pow(s, 2) +
  6281                      1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  6283                     (pow(
pion_mass, 4) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6284                      s * ((0.25 - 0.125 * 
delta) * pow(m_rho, 2) +
  6285                           (-0.25 + 0.125 * 
delta) * s) +
  6287                          (2. * 
C4 * pow(m_rho, 4) + (0.5 + 0.25 * 
delta) * s +
  6288                           pow(m_rho, 2) * (-1. - 2. * 
C4 * s))))) *
  6289           log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
  6292              log(abs(-1. * pow(
pion_mass, 2) + tmax)) +
  6295            8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  6297            5.999999999999999 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 2) *
  6299            1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
  6300           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6301              (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6303           (0. + 
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
  6304            eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
  6305                    pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
  6306           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6308          (2. * (-2. + 1. * 
delta) *
  6309           (0. + (-0.25 + 0.125 * 
delta) * pow(m_rho, 2) * s +
  6311                                 pow(m_rho, 2) * (1. + 2. * 
C4 * s))) *
  6312           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6315           (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
  6317                        (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6318                         (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  6320                        (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
  6322                    s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6325                (0.5 * pow(m_rho, 6) +
  6326                 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
  6327                 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
  6329                 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 4) -
  6330                                      3. * pow(m_rho, 2) * s + 1.5 * pow(s, 2)) +
  6332                     (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  6334                      pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6335           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6338               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  6339               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6341                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  6343              log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax)) +
  6347                (-1. * pow(m_rho, 2) + 1. * s) +
  6349                (-0.5 * pow(m_rho, 4) +
  6350                 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6351                 0.5 * pow(m_rho, 2) * s) +
  6353                (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
  6354                 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  6355           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6359           (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6360                    (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  6361            eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  6363                    pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  6364           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6366               1. * pow(m_rho, 2) + 1. * s) -
  6369            8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  6371            1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 4) * s +
  6373                (
C4 * (32. - 16. * 
delta) * pow(m_rho, 6) +
  6376                     (-8. + 24. * 
delta - 10. * pow(
delta, 2) + 32. * 
C4 * s -
  6378                 pow(m_rho, 4) * (-16. + 8. * 
delta - 64. * 
C4 * s +
  6380           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6381              (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6382          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  6390                        (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
  6391                         2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
  6393                        (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
  6394                         8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
  6398                    3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6400                        (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
  6402                    1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  6403                    1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
  6406                         1. * pow(m_rho, 2) + s) +
  6408                        (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6409                         4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  6410                         pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
  6413                    3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6415                        (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
  6420                         pow(m_rho, 2) + s) +
  6422                        (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6423                         2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  6424                         pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
  6428                      2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) -
  6429                      4. * pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) +
  6430                      pow(-1. * pow(m_rho, 2) + s + tmax, 2))) -
  6433                        (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6435                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  6437                        (pow(
pion_mass, 2) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6438                         (0.25 - 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  6440                        (1. * 
C4 * pow(m_rho, 6) +
  6442                             (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  6443                         0.25 * 
delta * pow(s, 2) +
  6444                         pow(m_rho, 4) * (-0.75 + 0.125 * 
delta - 2. * 
C4 * s) +
  6446                             (0.25 + 0.375 * 
delta + 1. * 
C4 * s))) +
  6449                     (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6451                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  6453                     (-0.5 * pow(m_rho, 4) + 1. * 
C4 * pow(m_rho, 6) +
  6455                          (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  6456                      0.25 * 
delta * pow(s, 2) -
  6457                      1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  6459                     (pow(
pion_mass, 4) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6460                      s * ((-0.25 + 0.125 * 
delta) * pow(m_rho, 2) +
  6461                           (0.25 - 0.125 * 
delta) * s) +
  6463                          (-2. * 
C4 * pow(m_rho, 4) + (-0.5 - 0.25 * 
delta) * s +
  6464                           pow(m_rho, 2) * (1. + 2. * 
C4 * s))))) *
  6467                           1. * pow(m_rho, 2) + s + tmax,
  6475                 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6478                          (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
  6481                                         0.5 * pow(m_rho, 2) - 1. * s)) +
  6483                     (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  6485                     (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
  6488                                       1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
  6492                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  6493                         0.5 * pow(m_rho, 2) * s) +
  6495                        (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6496                         pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
  6497                         0.5 * pow(m_rho, 2) * s) +
  6501                         0.5 * pow(m_rho, 2) * s +
  6502                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
  6505                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6506                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6507                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6508                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6510                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6516                    1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6519                    4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6520                    1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  6522                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6523                    0.5 * pow(m_rho, 4) * pow(s, 2) -
  6525                    1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  6529                         0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
  6532                         (1. * pow(m_rho, 2) - 0.5 * s) * s +
  6533                         pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
  6536                         pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
  6537                         0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
  6539                             (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
  6545                     (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
  6547                 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
  6548                                    (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
  6549                 pow(
pion_mass, 4) * (-0.5 * pow(m_rho, 4) +
  6550                                      2. * pow(m_rho, 2) * s - 1.5 * pow(s, 2)) +
  6552                     (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6555                     (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
  6556                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
  6557                      1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
  6559                                         1. * pow(m_rho, 2) + 1.5 * s)) +
  6562                      pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
  6564                          (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  6566                      s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6567                           0.5 * pow(s, 2))))) *
  6570                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6571                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6572                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6573                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6575                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6579               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  6580               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6582                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  6583          (0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  6587                     (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  6590                      0.5 * pow(m_rho, 2) + 0.5 * s) +
  6592                     (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
  6593                      pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
  6594                      6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
  6598                      (1. * pow(m_rho, 2) - 1. * s) * s +
  6599                      pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  6601                     (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  6602                      pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
  6603                      1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
  6605                          (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
  6608                      0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  6609                      1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  6610                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  6612                                         6. * pow(m_rho, 2) + 6. * s) +
  6614                          (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  6617                          (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  6618                           pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  6619                           6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  6621                     (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
  6622                      3.5 * pow(m_rho, 6) * s + 4.5 * pow(m_rho, 4) * pow(s, 2) -
  6623                      2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  6624                      pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  6626                          (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
  6629                          (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
  6630                           7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
  6634                     (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  6637                      0.5 * pow(m_rho, 2) + 0.5 * s) +
  6639                     (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
  6640                      5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
  6641                      pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  6643                     (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  6644                      pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
  6645                      3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
  6648                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
  6651                      0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  6652                      1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  6653                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  6655                                         6. * pow(m_rho, 2) + 6. * s) +
  6657                          (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  6660                          (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  6661                           pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  6662                           6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  6665                      pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
  6667                          (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
  6669                          (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  6672                          (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
  6673                           4.5 * pow(m_rho, 2) * pow(s, 2))) +
  6676                      pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  6678                          (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
  6681                          (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
  6682                           3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
  6683                      s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
  6684                           0.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3)))) +
  6689                      1. * pow(m_rho, 2) - 1. * s) +
  6691                     (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
  6693                     (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
  6694                      pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
  6695                      8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  6697                     (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  6698                      6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
  6699                      pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
  6701                          (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
  6703                     (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
  6704                      3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
  6705                      1. * pow(m_rho, 2) * pow(s, 3) +
  6706                      pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
  6708                          (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
  6710                          (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
  6711                           6. * pow(m_rho, 2) * pow(s, 2))) +
  6713                     (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
  6714                      pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
  6715                      2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
  6718                          (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
  6721                          (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
  6725                      1. * pow(m_rho, 6) +
  6726                      pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
  6728                                         12. * pow(m_rho, 2) - 12. * s) -
  6729                      3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
  6732                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  6735                          (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
  6736                           12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6738                               (-24. * pow(m_rho, 2) + 24. * s))))) *
  6741                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6742                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6743                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6744                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6746                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6750               pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  6751               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6753                   (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
 constexpr double a1_mass
a1 mass in GeV. 
 
static constexpr double C4
 
static constexpr double Pi
Value of . 
 
static constexpr double delta
 
static constexpr double eta1
 
constexpr double pion_mass
Pion mass in GeV. 
 
static constexpr double Const
 
std::array< T, 2 > get_t_range(const T sqrts, const T m1, const T m2, const T m3, const T m4)
Get the range of Mandelstam-t values allowed in a particular 2->2 process, see PDG 2014 booklet...
 
constexpr double gev2_mb
GeV^-2 <-> mb conversion factor. 
 
static constexpr double Gammaa1
 
double cut_off(const double sigma_mb)
Cross section after cut off. 
 
static constexpr double eta2
 
static constexpr double ghat