4459   const double s_sqrt = sqrt(s);
  4460   const double spin_deg_factor = 1.0;
  4463   double tmax = mandelstam_t[0];
  4464   double tmin = mandelstam_t[1];
  4470           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4473                  2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4476                      (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  4477                       4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  4479                      (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4480                       8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4484                  2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  4485                  1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4487                      (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  4489                      (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  4490                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4491                       2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4494                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4495                       pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
  4499                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  4501                      (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  4502                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4503                       4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4506                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4507                       pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  4508                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  4511                       2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4512                       2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  4514                           (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
  4515               (1. * pow(
a1_mass, 2) - 1. * tmin) +
  4517            (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  4520            (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  4521               (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
  4527                0.5 * 
eta2 * pow(m_rho, 2) + 0.5 * 
eta1 * s - 1. * 
eta2 * s) *
  4531                 (12. + 1. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  4532                  delta * (-8. + 8. * 
C4 * pow(m_rho, 2))) +
  4533             (-4. - 3. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  4534              delta * (8. + 8. * 
C4 * pow(m_rho, 2))) *
  4550                        0.5 * pow(m_rho, 2) + 1. * s)) *
  4552           (0.25 * (-2. + 1. * 
delta) *
  4553            (-8. * 
C4 * pow(m_rho, 4) +
  4555             (-2. - 3. * 
delta) * s +
  4556             pow(m_rho, 2) * (2. + 1. * 
delta + 16. * 
C4 * s)) *
  4560            (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
  4561             8 * 
C4 * pow(m_rho, 6) * (-6 + 
delta - 8 * 
C4 * s) +
  4562             2 * 
delta * pow(m_rho, 2) * s * (-6 + 
delta - 8 * 
C4 * s) +
  4563             pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 * 
C4 * (6 + 
delta) * s +
  4564                              32 * pow(
C4, 2) * pow(s, 2))) *
  4569                 (0.75 * pow(m_rho, 4) - 0.125 * 
delta * pow(m_rho, 4) -
  4570                  1. * 
C4 * pow(m_rho, 6) +
  4572                      (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4574                      (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) -
  4575                  0.25 * pow(m_rho, 2) * s - 0.375 * 
delta * pow(m_rho, 2) * s +
  4576                  2. * 
C4 * pow(m_rho, 4) * s + 0.25 * 
delta * pow(s, 2) -
  4577                  1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4578             eta1 * (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  4580                         (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) +
  4582                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4583                     0.25 * 
delta * pow(s, 2) +
  4584                     1. * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4587           0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  4590                                2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4591                                2. * pow(m_rho, 4) +
  4593                                                   pow(m_rho, 2) - 1. * s) -
  4594                                3. * pow(m_rho, 2) * s + pow(s, 2)) +
  4597                                2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4599                                                   pow(m_rho, 2) - 1. * s) +
  4600                                pow(m_rho, 2) * s + pow(s, 2)) +
  4604                     4. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  4605                     2. * pow(m_rho, 4) - 2. * pow(s, 2) +
  4607                                        2. * pow(m_rho, 2) + 2. * s))) *
  4609           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4612                     2. * pow(m_rho, 4) +
  4614                     8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4619                     pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4620                     4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4622                                        4. * pow(m_rho, 2) + 4. * s)) +
  4626                     pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4627                     2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4629                                        4. * pow(m_rho, 2) + 4. * s))) *
  4631           (0.125 * (-2. + 1. * 
delta) *
  4632            (2. + 1. * 
delta - 8. * 
C4 * pow(m_rho, 2)) * pow(tmin, 2)) /
  4634           0.5 * pow(1. * 
eta1 - 1. * 
eta2, 2) *
  4635               (-0.5 + 1. * 
C4 * pow(m_rho, 2)) * pow(tmin, 2) -
  4637            (0.5 - 0.125 * pow(
delta, 2) - 2. * 
C4 * pow(m_rho, 2) +
  4638             1. * 
C4 * 
delta * pow(m_rho, 2)) *
  4641           0.0625 * pow(1. * 
eta1 - 1. * 
eta2, 4) *
  4642               (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
  4644           0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  4646                        1. * pow(m_rho, 2) - 1. * s) +
  4648                        1. * pow(m_rho, 2) + s)) *
  4650           0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) * pow(tmin, 3) -
  4651           0.020833333333333332 * pow(1. * 
eta1 - 1. * 
eta2, 4) * pow(tmin, 3) +
  4652           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4656                     2. * pow(m_rho, 4) +
  4658                     4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
  4662                     2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  4663                     4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4665                                        4. * pow(m_rho, 2) + 4. * s)) +
  4669                     2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  4670                     2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4672                                        4. * pow(m_rho, 2) + 4. * s))) *
  4673               (-1. * pow(m_rho, 2) + s + tmin) -
  4674           0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  4675               (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
  4676                eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
  4677               pow(-1. * pow(m_rho, 2) + s + tmin, 2) +
  4678           0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) *
  4679               pow(-1. * pow(m_rho, 2) + s + tmin, 3) +
  4681               (-1. + 2. * 
C4 * pow(m_rho, 2)) *
  4683                       1. * pow(m_rho, 2) + s + tmin,
  4686            (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  4687                     0.5 * 
C4 * pow(m_rho, 6) +
  4689                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4691                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4692                     0.125 * pow(m_rho, 2) * s -
  4693                     0.1875 * 
delta * pow(m_rho, 2) * s +
  4694                     1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  4695                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4696             eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  4698                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4700                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4701                     0.125 * 
delta * pow(s, 2) +
  4702                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4704             1. * s + 1. * tmin)) /
  4707            (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  4708                     0.5 * 
C4 * pow(m_rho, 6) +
  4710                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4712                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4713                     0.125 * pow(m_rho, 2) * s -
  4714                     0.1875 * 
delta * pow(m_rho, 2) * s +
  4715                     1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  4716                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  4717             eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  4719                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  4721                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  4722                     0.125 * 
delta * pow(s, 2) +
  4723                     0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  4732                      (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  4734                      (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  4735                  pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
  4737                      (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
  4744                         (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  4745                          pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
  4746                          1. * pow(m_rho, 2) * s))) *
  4755                  2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4756                  1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  4758                      (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
  4760                  3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4761                  0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
  4762                  1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  4766                      (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  4767                       1. * pow(m_rho, 2) * s +
  4768                       pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
  4771                  4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  4772                  1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  4777                  4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  4778                  0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
  4779                  1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
  4782                       (1. * pow(m_rho, 2) - 1. * s) * s +
  4783                       pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
  4789                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  4790                2. * pow(m_rho, 2) * s + pow(s, 2) +
  4792                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
  4793           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  4799                      (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  4801                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4802                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  4803                       2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4806                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4807                       pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
  4810                       1. * pow(m_rho, 4) +
  4812                                          6. * pow(m_rho, 2) - 6. * s) +
  4813                       2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  4814                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
  4818                  2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4821                      (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  4822                       4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  4824                      (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4825                       8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4829                       2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
  4830                       4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
  4835                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  4837                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  4838                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  4839                       4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  4842                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4843                       pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  4844                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  4847                       1. * pow(m_rho, 4) +
  4849                                          6. * pow(m_rho, 2) - 6. * s) +
  4850                       4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  4851                       pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
  4854                       2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  4855                       2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  4857                           (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
  4864                 (-2. * pow(
Gammaa1, 4) * pow(a1_mass, 4) +
  4865                  14. * pow(a1_mass, 8) + 14. * pow(
pion_mass, 8) +
  4866                  28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  4867                  20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4868                  10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
  4870                  16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  4871                  12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  4872                  4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
  4873                  6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  4875                  4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
  4876                  pow(a1_mass, 6) * (-56. * pow(
pion_mass, 2) -
  4877                                     28. * pow(m_rho, 2) + 28. * s) +
  4879                      (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
  4880                       pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
  4881                       36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
  4882                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  4883                      (-4. * pow(a1_mass, 4) - 4. * pow(
pion_mass, 4) +
  4884                       pow(a1_mass, 2) * (8. * pow(
pion_mass, 2) +
  4885                                          4. * pow(m_rho, 2) - 4. * s) +
  4886                       (4. * pow(m_rho, 2) - 4. * s) * s +
  4887                       pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
  4889                      (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
  4890                       18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
  4892                       pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
  4894                           (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
  4895                            12. * pow(s, 2)))) +
  4897                 (1. * pow(
Gammaa1, 4) * pow(a1_mass, 4) - 7. * pow(a1_mass, 8) -
  4899                  14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  4900                  7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
  4901                  2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
  4903                      (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  4905                  11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
  4906                  6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  4907                  1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
  4908                  1. * pow(m_rho, 4) * pow(s, 2) -
  4910                  1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  4911                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  4912                      (2. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) +
  4913                       1. * pow(m_rho, 4) +
  4914                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  4915                       1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4916                       pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  4917                                          2. * pow(m_rho, 2) + 2. * s)) +
  4919                      (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
  4920                       21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  4921                       pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  4923                      (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  4924                       pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
  4925                       9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
  4928                           (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
  4931                 (1. * pow(
Gammaa1, 4) * pow(a1_mass, 4) - 7. * pow(a1_mass, 8) -
  4933                  14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  4934                  1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  4935                  6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
  4937                      (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  4939                  1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  4940                  16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  4941                  7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
  4942                  14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  4943                  9. * pow(m_rho, 4) * pow(s, 2) -
  4945                  5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  4946                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  4947                      (2. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) +
  4948                       3. * pow(m_rho, 4) +
  4949                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  4950                       5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  4951                       pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  4952                                          2. * pow(m_rho, 2) + 2. * s)) +
  4954                      (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
  4955                       9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  4956                       pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  4958                      (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  4959                       pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
  4960                       9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
  4963                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  4964                            6. * pow(s, 2))))) *
  4965            atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  4968               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + 4. * pow(a1_mass, 4) +
  4970                pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  4971                2. * pow(m_rho, 2) * s + pow(s, 2) +
  4973                    (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
  4974           0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  4976                    (-4. * pow(a1_mass, 6) +
  4977                     pow(a1_mass, 4) * (12. * pow(
pion_mass, 2) - 6. * s) +
  4979                         (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
  4980                          2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
  4982                         (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  4983                          8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  4986                    (2. * pow(a1_mass, 6) - 2. * pow(
pion_mass, 6) +
  4988                     pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
  4989                     pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
  4990                     pow(a1_mass, 4) * (-6. * pow(
pion_mass, 2) -
  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                          4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
  4997                    (2. * pow(a1_mass, 6) - 2. * pow(
pion_mass, 6) -
  4998                     1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  4999                     pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
  5000                     pow(a1_mass, 4) * (-6. * pow(
pion_mass, 2) +
  5001                                        3. * pow(m_rho, 2) + 3. * s) +
  5003                         (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5004                          pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5005                          2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
  5006               log(abs(-1. * pow(a1_mass, 2) + tmin)) -
  5009                 (-0.5 * pow(a1_mass, 6) - 0.5 * pow(
pion_mass, 6) +
  5010                  0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5012                      (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5014                      (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
  5015             eta1 * (pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) + 0.5 * s) +
  5017                         (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5018                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  5019                          0.5 * pow(m_rho, 2) * s) +
  5021                         (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
  5022                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5023            log(abs(-1. * pow(a1_mass, 2) + tmin))) /
  5024               (pow(a1_mass, 2) - 1. * pow(
pion_mass, 2)) +
  5028                         (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
  5029                     pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
  5030                     pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
  5033                          pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5034                          (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
  5036                     pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) - 0.5 * s) +
  5037                     pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
  5038                     (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
  5039                     pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
  5041                         (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
  5042                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5043            log(abs(-1. * pow(a1_mass, 2) + tmin))) /
  5044               (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  5045                1. * pow(m_rho, 2) + 1. * s) -
  5046           (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5047            (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  5048             0.5 * pow(m_rho, 2) + 0.5 * s) *
  5050                 (-2. * pow(a1_mass, 8) +
  5052                      (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
  5054                      (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
  5055                       4. * pow(m_rho, 2) * s +
  5056                       pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
  5059                       4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5060                       2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
  5061                       2. * pow(m_rho, 2) * pow(s, 2) +
  5063                           (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
  5065                      (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  5066                       pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
  5067                       2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
  5070                           (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
  5075                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5077                                       2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  5078                                       pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5080                      (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5081                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
  5082                       pow(m_rho, 2) * s) +
  5084                      (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
  5085                       5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
  5087                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  5089                           (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
  5093                  2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5095                      (-2. * pow(m_rho, 2) + 2. * s) +
  5097                      (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5099                      (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
  5101                      (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5102                       pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5103                       3. * pow(m_rho, 2) * s) +
  5105                      (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
  5107                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  5109                           (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
  5110                            2. * pow(s, 2))))) *
  5111            log(abs(-1. * pow(a1_mass, 2) + tmin))) /
  5112               (0.25 * pow(
Gammaa1, 2) * pow(a1_mass, 2) + 1. * pow(a1_mass, 4) +
  5114                0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
  5115                0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
  5117                    (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
  5121                      (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5123                      (pow(
pion_mass, 2) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5124                       (-0.25 + 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  5126                      (-1. * 
C4 * pow(m_rho, 6) +
  5128                           (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  5129                       0.25 * 
delta * pow(s, 2) +
  5131                           (-0.25 - 0.375 * 
delta - 1. * 
C4 * s) +
  5132                       pow(m_rho, 4) * (0.75 - 0.125 * 
delta + 2. * 
C4 * s))) +
  5135                      (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  5137                      (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  5139                           (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  5140                       0.25 * 
delta * pow(s, 2) +
  5141                       1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5143                      (pow(
pion_mass, 4) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5144                       s * ((0.25 - 0.125 * 
delta) * pow(m_rho, 2) +
  5145                            (-0.25 + 0.125 * 
delta) * s) +
  5147                           (2. * 
C4 * pow(m_rho, 4) + (0.5 + 0.25 * 
delta) * s +
  5148                            pow(m_rho, 2) * (-1. - 2. * 
C4 * s))))) *
  5149            log(abs(-1. * pow(a1_mass, 2) + tmin))) /
  5152               log(abs(-1. * pow(
pion_mass, 2) + tmin)) +
  5155             8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  5157             5.999999999999999 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 2) *
  5159             1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
  5160            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5161               (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5163            (0. + 
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
  5164             eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
  5165                     pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
  5166            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5167               (pow(a1_mass, 2) - 1. * pow(
pion_mass, 2)) +
  5168           (2. * (-2. + 1. * 
delta) *
  5169            (0. + (-0.25 + 0.125 * 
delta) * pow(m_rho, 2) * s +
  5171                                  pow(m_rho, 2) * (1. + 2. * 
C4 * s))) *
  5172            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5175            (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
  5177                         (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5178                          (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  5180                         (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
  5182                     s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5184             eta2 * (0.5 * pow(m_rho, 6) +
  5185                     pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
  5186                     1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
  5189                         (1.5 * pow(m_rho, 4) - 3. * pow(m_rho, 2) * s +
  5192                         (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  5194                          pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  5195            log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
  5196               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5198                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5199                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5201                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  5203               log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s +
  5208                 (-1. * pow(m_rho, 2) + 1. * s) +
  5210                 (-0.5 * pow(m_rho, 4) +
  5211                  pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5212                  0.5 * pow(m_rho, 2) * s) +
  5214                 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
  5215                  pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  5216            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5217               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  5220            (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  5221                     (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  5222             eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  5224                     pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  5225            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5226               (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  5227                1. * pow(m_rho, 2) + 1. * s) -
  5230             8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  5232             1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 4) * s +
  5234                 (
C4 * (32. - 16. * 
delta) * pow(m_rho, 6) +
  5237                      (-8. + 24. * 
delta - 10. * pow(
delta, 2) + 32. * 
C4 * s -
  5239                  pow(m_rho, 4) * (-16. + 8. * 
delta - 64. * 
C4 * s +
  5241            log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
  5242               (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  5243           0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5245                    (4. * pow(a1_mass, 6) +
  5246                     pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5247                         (-4. * pow(a1_mass, 2) + 4. * pow(
pion_mass, 2) -
  5249                     pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) + 6. * s) +
  5251                         (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
  5252                          2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
  5254                         (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
  5255                          8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
  5258                    (-2. * pow(a1_mass, 6) + 2. * pow(
pion_mass, 6) +
  5259                     3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5261                         (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
  5263                     1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  5264                     1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
  5265                     pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5266                         (2. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) -
  5267                          1. * pow(m_rho, 2) + s) +
  5269                         (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5270                          4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5271                          pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
  5273                    (-2. * pow(a1_mass, 6) + 2. * pow(
pion_mass, 6) -
  5274                     3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5276                         (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
  5279                     pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5280                         (2. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) +
  5281                          pow(m_rho, 2) + s) +
  5283                         (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5284                          2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5285                          pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
  5286               log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  5287                       4. * pow(a1_mass, 2) * pow(
pion_mass, 2) +
  5289                       2. * pow(a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmin) -
  5291                           (-1. * pow(m_rho, 2) + s + tmin) +
  5292                       pow(-1. * pow(m_rho, 2) + s + tmin, 2))) -
  5294            (
eta2 * (pow(
Gammaa1, 2) * pow(a1_mass, 2) * pow(m_rho, 2) *
  5295                         (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5297                         (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  5299                         (pow(
pion_mass, 2) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5300                          (0.25 - 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  5302                         (1. * 
C4 * pow(m_rho, 6) +
  5304                              (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  5305                          0.25 * 
delta * pow(s, 2) +
  5306                          pow(m_rho, 4) * (-0.75 + 0.125 * 
delta - 2. * 
C4 * s) +
  5308                              (0.25 + 0.375 * 
delta + 1. * 
C4 * s))) +
  5309             eta1 * (pow(
Gammaa1, 2) * pow(a1_mass, 2) * pow(m_rho, 2) *
  5310                         (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  5312                         (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5314                         (-0.5 * pow(m_rho, 4) + 1. * 
C4 * pow(m_rho, 6) +
  5316                              (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  5317                          0.25 * 
delta * pow(s, 2) -
  5318                          1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5320                         (pow(
pion_mass, 4) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  5321                          s * ((-0.25 + 0.125 * 
delta) * pow(m_rho, 2) +
  5322                               (0.25 - 0.125 * 
delta) * s) +
  5324                              (-2. * 
C4 * pow(m_rho, 4) +
  5325                               (-0.5 - 0.25 * 
delta) * s +
  5326                               pow(m_rho, 2) * (1. + 2. * 
C4 * s))))) *
  5327            log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) +
  5328                    pow(pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) -
  5329                            1. * pow(m_rho, 2) + s + tmin,
  5334                 (0.5 * pow(
Gammaa1, 4) * pow(a1_mass, 4) -
  5335                  0.5 * pow(a1_mass, 8) + 0.5 * pow(
pion_mass, 8) +
  5336                  0.5 * pow(a1_mass, 4) * pow(
pion_mass, 2) * pow(m_rho, 2) -
  5337                  0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5340                           (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
  5342                       pow(a1_mass, 4) * (-1. * pow(
pion_mass, 2) +
  5343                                          0.5 * pow(m_rho, 2) - 1. * s)) +
  5345                      (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5347                      (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
  5349                 (pow(a1_mass, 6) * (1. * pow(
pion_mass, 2) + 0.5 * s) +
  5350                  pow(a1_mass, 2) * (3. * pow(
pion_mass, 6) +
  5351                                     1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
  5355                       pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  5356                       0.5 * pow(m_rho, 2) * s) +
  5358                      (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5359                       pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
  5360                       0.5 * pow(m_rho, 2) * s) +
  5361                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5363                       pow(a1_mass, 2) * (1. * pow(
pion_mass, 2) + 0.5 * s) -
  5364                       0.5 * pow(m_rho, 2) * s +
  5365                       pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
  5366            log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5368                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5369                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5370                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5371                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5373                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5375               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  5379                     1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5380                     pow(a1_mass, 6) * (1. * pow(
pion_mass, 2) - 0.5 * s) -
  5382                     4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  5383                     1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  5385                     4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  5386                     0.5 * pow(m_rho, 4) * pow(s, 2) -
  5388                     1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  5389                     pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5391                          pow(a1_mass, 2) * (1. * pow(
pion_mass, 2) - 0.5 * s) -
  5392                          0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
  5395                          (1. * pow(m_rho, 2) - 0.5 * s) * s +
  5396                          pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
  5399                          pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
  5400                          0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
  5402                              (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
  5405                 (0.5 * pow(
Gammaa1, 4) * pow(a1_mass, 4) -
  5406                  0.5 * pow(a1_mass, 8) + 0.5 * pow(
pion_mass, 8) +
  5408                      (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
  5410                  pow(a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
  5411                                     (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
  5413                      (-0.5 * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5416                      (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5418                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5419                      (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
  5420                       pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
  5421                       1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
  5422                       pow(a1_mass, 2) * (-1. * pow(
pion_mass, 2) -
  5423                                          1. * pow(m_rho, 2) + 1.5 * s)) +
  5426                       pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
  5428                           (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5430                       s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  5431                            0.5 * pow(s, 2))))) *
  5432            log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5434                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5435                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5436                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5437                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5439                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5441               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5443                pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5444                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5446                    (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  5447           (0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  5449                 (-1. * pow(a1_mass, 10) +
  5451                      (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  5452                  pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  5453                      (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  5454                       0.5 * pow(m_rho, 2) + 0.5 * s) +
  5456                      (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
  5457                       pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
  5458                       6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
  5462                       (1. * pow(m_rho, 2) - 1. * s) * s +
  5463                       pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  5465                      (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  5466                       pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
  5467                       1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
  5469                           (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
  5470                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5471                      (4. * pow(a1_mass, 6) - 4. * pow(
pion_mass, 6) -
  5472                       0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  5473                       1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5474                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  5475                       pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) -
  5476                                          6. * pow(m_rho, 2) + 6. * s) +
  5478                           (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  5481                           (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  5482                            pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  5483                            6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  5485                      (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
  5486                       3.5 * pow(m_rho, 6) * s +
  5487                       4.5 * pow(m_rho, 4) * pow(s, 2) -
  5488                       2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  5489                       pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  5491                           (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
  5494                           (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
  5495                            7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
  5497                 (-1. * pow(a1_mass, 10) +
  5499                      (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  5500                  pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  5501                      (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  5502                       0.5 * pow(m_rho, 2) + 0.5 * s) +
  5504                      (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
  5505                       5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
  5506                       pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  5508                      (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  5509                       pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
  5510                       3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
  5513                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
  5514                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5515                      (4. * pow(a1_mass, 6) - 4. * pow(
pion_mass, 6) -
  5516                       0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  5517                       1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5518                       pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  5519                       pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) -
  5520                                          6. * pow(m_rho, 2) + 6. * s) +
  5522                           (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  5525                           (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  5526                            pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  5527                            6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  5530                       pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
  5532                           (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
  5534                           (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  5537                           (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
  5538                            4.5 * pow(m_rho, 2) * pow(s, 2))) +
  5541                       pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  5543                           (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
  5546                           (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
  5547                            3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
  5548                       s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
  5549                            0.5 * pow(m_rho, 2) * pow(s, 2) +
  5550                            0.5 * pow(s, 3)))) +
  5552                 (2. * pow(a1_mass, 10) +
  5553                  pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  5554                      (-2. * pow(a1_mass, 2) + 2. * pow(
pion_mass, 2) +
  5555                       1. * pow(m_rho, 2) - 1. * s) +
  5557                      (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
  5559                      (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
  5560                       pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
  5561                       8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5563                      (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  5564                       6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
  5565                       pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
  5567                           (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
  5569                      (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
  5570                       3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
  5571                       1. * pow(m_rho, 2) * pow(s, 3) +
  5572                       pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
  5574                           (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
  5576                           (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
  5577                            6. * pow(m_rho, 2) * pow(s, 2))) +
  5579                      (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
  5580                       pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
  5581                       2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
  5584                           (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
  5587                           (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
  5589                  pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5590                      (-8. * pow(a1_mass, 6) + 8. * pow(
pion_mass, 6) +
  5591                       1. * pow(m_rho, 6) +
  5592                       pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
  5593                       pow(a1_mass, 4) * (24. * pow(
pion_mass, 2) +
  5594                                          12. * pow(m_rho, 2) - 12. * s) -
  5595                       3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
  5598                           (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  5601                           (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
  5602                            12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  5604                                (-24. * pow(m_rho, 2) + 24. * s))))) *
  5605            log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5607                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5608                    4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  5609                    pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
  5610                    2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
  5612                        (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  5614               (pow(
Gammaa1, 2) * pow(a1_mass, 2) + 4. * pow(a1_mass, 4) +
  5616                pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  5617                2. * pow(m_rho, 2) * s + pow(s, 2) +
  5619                    (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
  5623          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5625                (-2. * pow(a1_mass, 8) - 2. * pow(
pion_mass, 8) +
  5626                 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5627                 pow(a1_mass, 6) * (8. * pow(
pion_mass, 2) - 4. * s) +
  5629                     (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  5630                      4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  5632                     (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  5633                      8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5636                (1. * pow(a1_mass, 8) + 1. * pow(
pion_mass, 8) -
  5637                 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  5638                 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5640                     (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  5642                     (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5643                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5644                      2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5647                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5648                      pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
  5650                (1. * pow(a1_mass, 8) +
  5652                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5654                     (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  5655                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5656                      4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5659                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5660                      pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  5661                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  5664                      2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5665                      2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  5667                          (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
  5668              (1. * pow(a1_mass, 2) - 1. * tmax) +
  5670           (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  5673           (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
  5674              (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
  5680               0.5 * 
eta2 * pow(m_rho, 2) + 0.5 * 
eta1 * s - 1. * 
eta2 * s) *
  5684                (12. + 1. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  5685                 delta * (-8. + 8. * 
C4 * pow(m_rho, 2))) +
  5686            (-4. - 3. * pow(
delta, 2) - 16. * 
C4 * pow(m_rho, 2) +
  5687             delta * (8. + 8. * 
C4 * pow(m_rho, 2))) *
  5692              (
eta2 * (pow(a1_mass, 2) - 1. * s) +
  5696              (
eta2 * (-1. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * s) +
  5702                       0.5 * pow(m_rho, 2) + 1. * s)) *
  5704          (0.25 * (-2. + 1. * 
delta) *
  5705           (-8. * 
C4 * pow(m_rho, 4) +
  5707            (-2. - 3. * 
delta) * s +
  5708            pow(m_rho, 2) * (2. + 1. * 
delta + 16. * 
C4 * s)) *
  5712           (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
  5713            8 * 
C4 * pow(m_rho, 6) * (-6 + 
delta - 8 * 
C4 * s) +
  5714            2 * 
delta * pow(m_rho, 2) * s * (-6 + 
delta - 8 * 
C4 * s) +
  5715            pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 * 
C4 * (6 + 
delta) * s +
  5716                             32 * pow(
C4, 2) * pow(s, 2))) *
  5721                (0.75 * pow(m_rho, 4) - 0.125 * 
delta * pow(m_rho, 4) -
  5722                 1. * 
C4 * pow(m_rho, 6) +
  5724                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5726                     (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) -
  5727                 0.25 * pow(m_rho, 2) * s - 0.375 * 
delta * pow(m_rho, 2) * s +
  5728                 2. * 
C4 * pow(m_rho, 4) * s + 0.25 * 
delta * pow(s, 2) -
  5729                 1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5730            eta1 * (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  5732                        (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) +
  5734                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5735                    0.25 * 
delta * pow(s, 2) +
  5736                    1. * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5739          0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  5741                   (pow(
Gammaa1, 2) * pow(a1_mass, 2) - 1. * pow(a1_mass, 4) -
  5743                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) +
  5745                        (2. * pow(
pion_mass, 2) + pow(m_rho, 2) - 1. * s) -
  5746                    3. * pow(m_rho, 2) * s + pow(s, 2)) +
  5747               pow(
eta1, 2) * (pow(
Gammaa1, 2) * pow(a1_mass, 2) -
  5748                               1. * pow(a1_mass, 4) - 2. * pow(
pion_mass, 4) -
  5749                               2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  5750                               pow(a1_mass, 2) * (2. * pow(
pion_mass, 2) +
  5751                                                  pow(m_rho, 2) - 1. * s) +
  5752                               pow(m_rho, 2) * s + pow(s, 2)) +
  5754                   (-2. * pow(
Gammaa1, 2) * pow(a1_mass, 2) +
  5755                    2. * pow(a1_mass, 4) + 4. * pow(
pion_mass, 4) +
  5756                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) -
  5758                    pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  5759                                       2. * pow(m_rho, 2) + 2. * s))) *
  5761          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5763                   (-6. * pow(a1_mass, 4) - 12. * pow(
pion_mass, 4) +
  5764                    2. * pow(m_rho, 4) +
  5765                    pow(a1_mass, 2) * (16. * pow(
pion_mass, 2) - 8. * s) +
  5766                    8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5769                   (3. * pow(a1_mass, 4) + 6. * pow(
pion_mass, 4) +
  5771                    pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5772                    4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5773                    pow(a1_mass, 2) * (-8. * pow(
pion_mass, 2) -
  5774                                       4. * pow(m_rho, 2) + 4. * s)) +
  5776                   (3. * pow(a1_mass, 4) + 6. * pow(
pion_mass, 4) +
  5778                    pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5779                    2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5780                    pow(a1_mass, 2) * (-8. * pow(
pion_mass, 2) +
  5781                                       4. * pow(m_rho, 2) + 4. * s))) *
  5783          (0.125 * (-2. + 1. * 
delta) *
  5784           (2. + 1. * 
delta - 8. * 
C4 * pow(m_rho, 2)) * pow(tmax, 2)) /
  5786          0.5 * pow(1. * 
eta1 - 1. * 
eta2, 2) *
  5787              (-0.5 + 1. * 
C4 * pow(m_rho, 2)) * pow(tmax, 2) -
  5789           (0.5 - 0.125 * pow(
delta, 2) - 2. * 
C4 * pow(m_rho, 2) +
  5790            1. * 
C4 * 
delta * pow(m_rho, 2)) *
  5793          0.0625 * pow(1. * 
eta1 - 1. * 
eta2, 4) *
  5794              (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
  5796          0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  5798                       1. * pow(m_rho, 2) - 1. * s) +
  5800                       1. * pow(m_rho, 2) + s)) *
  5802          0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) * pow(tmax, 3) -
  5803          0.020833333333333332 * pow(1. * 
eta1 - 1. * 
eta2, 4) * pow(tmax, 3) +
  5804          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5806                   (2. * pow(
Gammaa1, 2) * pow(a1_mass, 2) -
  5807                    6. * pow(a1_mass, 4) - 4. * pow(
pion_mass, 4) +
  5808                    2. * pow(m_rho, 4) +
  5809                    pow(a1_mass, 2) * (8. * pow(
pion_mass, 2) - 8. * s) +
  5810                    4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
  5812                   (-1. * pow(
Gammaa1, 2) * pow(a1_mass, 2) +
  5813                    3. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) +
  5814                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5815                    4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5816                    pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  5817                                       4. * pow(m_rho, 2) + 4. * s)) +
  5819                   (-1. * pow(
Gammaa1, 2) * pow(a1_mass, 2) +
  5820                    3. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) -
  5821                    2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  5822                    2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  5823                    pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) +
  5824                                       4. * pow(m_rho, 2) + 4. * s))) *
  5825              (-1. * pow(m_rho, 2) + s + tmax) -
  5826          0.03125 * pow(
eta1 - 1. * 
eta2, 3) *
  5827              (
eta2 * (-1. * pow(a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
  5828               eta1 * (pow(a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
  5829              pow(-1. * pow(m_rho, 2) + s + tmax, 2) +
  5830          0.010416666666666666 * pow(
eta1 - 1. * 
eta2, 4) *
  5831              pow(-1. * pow(m_rho, 2) + s + tmax, 3) +
  5833              (-1. + 2. * 
C4 * pow(m_rho, 2)) *
  5834              pow(pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  5839                (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  5840                 0.5 * 
C4 * pow(m_rho, 6) +
  5842                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5844                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5845                 0.125 * pow(m_rho, 2) * s - 0.1875 * 
delta * pow(m_rho, 2) * s +
  5846                 1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  5847                 0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5848            eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  5850                        (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5852                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5853                    0.125 * 
delta * pow(s, 2) +
  5854                    0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5855           (1. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  5856            1. * s + 1. * tmax)) /
  5860                (0.375 * pow(m_rho, 4) - 0.0625 * 
delta * pow(m_rho, 4) -
  5861                 0.5 * 
C4 * pow(m_rho, 6) +
  5863                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5865                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5866                 0.125 * pow(m_rho, 2) * s - 0.1875 * 
delta * pow(m_rho, 2) * s +
  5867                 1. * 
C4 * pow(m_rho, 4) * s + 0.125 * 
delta * pow(s, 2) -
  5868                 0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  5869            eta1 * (0.25 * pow(m_rho, 4) - 0.5 * 
C4 * pow(m_rho, 6) +
  5871                        (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  5873                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) -
  5874                    0.125 * 
delta * pow(s, 2) +
  5875                    0.5 * 
C4 * pow(m_rho, 2) * pow(s, 2))) *
  5876           atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  5882                (-1. * pow(a1_mass, 6) +
  5883                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5884                     (-1. * pow(a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5886                     (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  5887                 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
  5889                     (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
  5892                    pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) + 0.5 * s) +
  5896                        (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5897                         pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
  5898                         1. * pow(m_rho, 2) * s))) *
  5899           atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  5902              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  5906                (-1. * pow(a1_mass, 6) - 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  5907                 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  5909                     (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
  5911                 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5912                 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
  5913                 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  5915                     (-1. * pow(a1_mass, 4) + 0.5 * pow(a1_mass, 2) * s) +
  5917                     (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  5918                      1. * pow(m_rho, 2) * s +
  5919                      pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
  5922                 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  5923                 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
  5924                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5926                 pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) - 0.5 * s) -
  5928                 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  5929                 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
  5930                 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
  5933                      (1. * pow(m_rho, 2) - 1. * s) * s +
  5934                      pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
  5935           atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  5938              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  5940               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  5941               2. * pow(m_rho, 2) * s + pow(s, 2) +
  5943                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
  5944          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  5946                (pow(
Gammaa1, 4) * pow(a1_mass, 4) + pow(a1_mass, 8) +
  5950                     (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
  5952                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5953                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  5954                      2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5957                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5958                      pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
  5959                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5960                     (-6. * pow(a1_mass, 4) - 6. * pow(
pion_mass, 4) -
  5961                      1. * pow(m_rho, 4) +
  5962                      pow(a1_mass, 2) * (12. * pow(
pion_mass, 2) -
  5963                                         6. * pow(m_rho, 2) - 6. * s) +
  5964                      2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  5965                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
  5967                (-2. * pow(
Gammaa1, 4) * pow(a1_mass, 4) - 2. * pow(a1_mass, 8) -
  5969                 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  5970                 pow(a1_mass, 6) * (8. * pow(
pion_mass, 2) - 4. * s) +
  5972                     (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
  5973                      4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
  5975                     (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  5976                      8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  5978                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5979                     (12. * pow(a1_mass, 4) + 12. * pow(
pion_mass, 4) -
  5980                      2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
  5981                      4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
  5982                      pow(a1_mass, 2) * (-24. * pow(
pion_mass, 2) + 12. * s))) +
  5984                (pow(
Gammaa1, 4) * pow(a1_mass, 4) + pow(a1_mass, 8) +
  5986                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  5988                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  5989                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  5990                      4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  5993                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  5994                      pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
  5995                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
  5996                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  5997                     (-6. * pow(a1_mass, 4) - 6. * pow(
pion_mass, 4) -
  5998                      1. * pow(m_rho, 4) +
  5999                      pow(a1_mass, 2) * (12. * pow(
pion_mass, 2) +
  6000                                         6. * pow(m_rho, 2) - 6. * s) +
  6001                      4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  6002                      pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
  6005                      2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  6006                      2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
  6008                          (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
  6009           atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  6015                (-2. * pow(
Gammaa1, 4) * pow(a1_mass, 4) +
  6016                 14. * pow(a1_mass, 8) + 14. * pow(
pion_mass, 8) +
  6017                 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6018                 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6019                 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
  6021                 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6022                 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  6023                 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
  6024                 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6026                 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
  6028                     (-56. * pow(
pion_mass, 2) - 28. * pow(m_rho, 2) + 28. * s) +
  6030                     (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
  6031                      pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
  6032                      36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
  6033                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6034                     (-4. * pow(a1_mass, 4) - 4. * pow(
pion_mass, 4) +
  6035                      pow(a1_mass, 2) * (8. * pow(
pion_mass, 2) +
  6036                                         4. * pow(m_rho, 2) - 4. * s) +
  6037                      (4. * pow(m_rho, 2) - 4. * s) * s +
  6038                      pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
  6040                     (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
  6041                      18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
  6043                      pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
  6045                          (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
  6046                           12. * pow(s, 2)))) +
  6048                (1. * pow(
Gammaa1, 4) * pow(a1_mass, 4) - 7. * pow(a1_mass, 8) -
  6050                 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  6051                 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
  6052                 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
  6054                     (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  6056                 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
  6057                 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  6058                 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
  6059                 1. * pow(m_rho, 4) * pow(s, 2) -
  6061                 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  6062                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6063                     (2. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) +
  6064                      1. * pow(m_rho, 4) +
  6065                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  6066                      1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  6067                      pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  6068                                         2. * pow(m_rho, 2) + 2. * s)) +
  6070                     (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
  6071                      21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6072                      pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  6074                     (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  6075                      pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
  6076                      9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
  6079                          (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
  6082                (1. * pow(
Gammaa1, 4) * pow(a1_mass, 4) - 7. * pow(a1_mass, 8) -
  6084                 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
  6085                 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6086                 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
  6088                     (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
  6090                 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6091                 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
  6092                 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
  6093                 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6094                 9. * pow(m_rho, 4) * pow(s, 2) -
  6096                 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
  6097                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6098                     (2. * pow(a1_mass, 4) + 2. * pow(
pion_mass, 4) +
  6099                      3. * pow(m_rho, 4) +
  6100                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
  6101                      5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
  6102                      pow(a1_mass, 2) * (-4. * pow(
pion_mass, 2) -
  6103                                         2. * pow(m_rho, 2) + 2. * s)) +
  6105                     (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
  6106                      9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6107                      pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
  6109                     (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  6110                      pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
  6111                      9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
  6114                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  6115                           6. * pow(s, 2))))) *
  6116           atan((pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) +
  6119              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + 4. * pow(a1_mass, 4) +
  6121               pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  6122               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6124                   (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
  6125          0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  6127                   (-4. * pow(a1_mass, 6) +
  6128                    pow(a1_mass, 4) * (12. * pow(
pion_mass, 2) - 6. * s) +
  6130                        (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
  6131                         2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
  6133                        (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
  6134                         8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
  6137                   (2. * pow(a1_mass, 6) - 2. * pow(
pion_mass, 6) +
  6139                    pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * 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                         4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
  6148                   (2. * pow(a1_mass, 6) - 2. * pow(
pion_mass, 6) -
  6149                    1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
  6150                    pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
  6152                        (-6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) + 3. * s) +
  6154                        (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  6155                         pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
  6156                         2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
  6157              log(abs(-1. * pow(a1_mass, 2) + tmax)) -
  6160                (-0.5 * pow(a1_mass, 6) - 0.5 * pow(
pion_mass, 6) +
  6161                 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6163                     (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  6165                     (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
  6166            eta1 * (pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) + 0.5 * s) +
  6168                        (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
  6169                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  6170                         0.5 * pow(m_rho, 2) * s) +
  6172                        (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
  6173                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6174           log(abs(-1. * pow(a1_mass, 2) + tmax))) /
  6175              (pow(a1_mass, 2) - 1. * pow(
pion_mass, 2)) +
  6179                        (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
  6180                    pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
  6181                    pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
  6184                         pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6185                         (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
  6187                    pow(a1_mass, 4) * (1. * pow(
pion_mass, 2) - 0.5 * s) +
  6188                    pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
  6189                    (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
  6190                    pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
  6192                        (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
  6193                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6194           log(abs(-1. * pow(a1_mass, 2) + tmax))) /
  6195              (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  6196               1. * pow(m_rho, 2) + 1. * s) -
  6197          (0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  6198           (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) +
  6201                (-2. * pow(a1_mass, 8) +
  6203                     (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
  6205                     (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
  6206                      4. * pow(m_rho, 2) * s +
  6207                      pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
  6210                      4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
  6211                      2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
  6212                      2. * pow(m_rho, 2) * pow(s, 2) +
  6214                          (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
  6216                     (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
  6217                      pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
  6218                      2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
  6221                          (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
  6226                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  6228                                      2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
  6229                                      pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6231                     (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6232                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
  6233                      pow(m_rho, 2) * s) +
  6235                     (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
  6236                      5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
  6238                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  6240                          (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
  6244                 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6246                     (-2. * pow(m_rho, 2) + 2. * s) +
  6248                     (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
  6250                     (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
  6252                     (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
  6253                      pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
  6254                      3. * pow(m_rho, 2) * s) +
  6256                     (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
  6258                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
  6260                          (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
  6261                           2. * pow(s, 2))))) *
  6262           log(abs(-1. * pow(a1_mass, 2) + tmax))) /
  6263              (0.25 * pow(
Gammaa1, 2) * pow(a1_mass, 2) + 1. * pow(a1_mass, 4) +
  6265               0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
  6266               0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
  6268                   (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
  6272                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  6274                     (pow(
pion_mass, 2) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6275                      (-0.25 + 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  6277                     (-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                      pow(m_rho, 2) * s * (-0.25 - 0.375 * 
delta - 1. * 
C4 * s) +
  6282                      pow(m_rho, 4) * (0.75 - 0.125 * 
delta + 2. * 
C4 * s))) +
  6285                     (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  6287                     (0.5 * pow(m_rho, 4) - 1. * 
C4 * pow(m_rho, 6) +
  6289                          (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  6290                      0.25 * 
delta * pow(s, 2) +
  6291                      1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  6293                     (pow(
pion_mass, 4) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6294                      s * ((0.25 - 0.125 * 
delta) * pow(m_rho, 2) +
  6295                           (-0.25 + 0.125 * 
delta) * s) +
  6297                          (2. * 
C4 * pow(m_rho, 4) + (0.5 + 0.25 * 
delta) * s +
  6298                           pow(m_rho, 2) * (-1. - 2. * 
C4 * s))))) *
  6299           log(abs(-1. * pow(a1_mass, 2) + tmax))) /
  6302              log(abs(-1. * pow(
pion_mass, 2) + tmax)) +
  6305            8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  6307            5.999999999999999 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 2) *
  6309            1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
  6310           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6311              (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6313           (0. + 
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
  6314            eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
  6315                    pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
  6316           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6317              (pow(a1_mass, 2) - 1. * pow(
pion_mass, 2)) +
  6318          (2. * (-2. + 1. * 
delta) *
  6319           (0. + (-0.25 + 0.125 * 
delta) * pow(m_rho, 2) * s +
  6321                                 pow(m_rho, 2) * (1. + 2. * 
C4 * s))) *
  6322           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6325           (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
  6327                        (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6328                         (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  6330                        (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
  6332                    s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6335                (0.5 * pow(m_rho, 6) +
  6336                 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
  6337                 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
  6339                 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 4) -
  6340                                      3. * pow(m_rho, 2) * s + 1.5 * pow(s, 2)) +
  6342                     (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  6344                      pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
  6345           log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
  6346              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  6348               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  6349               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6351                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  6353              log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax)) +
  6357                (-1. * pow(m_rho, 2) + 1. * s) +
  6359                (-0.5 * pow(m_rho, 4) +
  6360                 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6361                 0.5 * pow(m_rho, 2) * s) +
  6363                (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
  6364                 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  6365           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6366              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  6369           (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
  6370                    (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
  6371            eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
  6373                    pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
  6374           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6375              (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  6376               1. * pow(m_rho, 2) + 1. * s) -
  6379            8.000000000000002 * pow(2. - 1. * 
delta, 2) * pow(
pion_mass, 4) *
  6381            1. * pow(2. - 1. * 
delta, 2) * pow(m_rho, 4) * s +
  6383                (
C4 * (32. - 16. * 
delta) * pow(m_rho, 6) +
  6386                     (-8. + 24. * 
delta - 10. * pow(
delta, 2) + 32. * 
C4 * s -
  6388                 pow(m_rho, 4) * (-16. + 8. * 
delta - 64. * 
C4 * s +
  6390           log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
  6391              (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
  6392          0.03125 * pow(
eta1 - 1. * 
eta2, 2) *
  6394                   (4. * pow(a1_mass, 6) +
  6395                    pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6396                        (-4. * pow(a1_mass, 2) + 4. * pow(
pion_mass, 2) -
  6398                    pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) + 6. * s) +
  6400                        (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
  6401                         2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
  6403                        (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
  6404                         8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
  6407                   (-2. * pow(a1_mass, 6) + 2. * pow(
pion_mass, 6) +
  6408                    3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6410                        (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
  6412                    1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
  6413                    1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
  6414                    pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6415                        (2. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) -
  6416                         1. * pow(m_rho, 2) + s) +
  6418                        (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6419                         4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  6420                         pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
  6422                   (-2. * pow(a1_mass, 6) + 2. * pow(
pion_mass, 6) -
  6423                    3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
  6425                        (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
  6428                    pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6429                        (2. * pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) +
  6430                         pow(m_rho, 2) + s) +
  6432                        (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6433                         2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
  6434                         pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
  6435              log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  6436                      4. * pow(a1_mass, 2) * pow(
pion_mass, 2) +
  6438                      2. * pow(a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) -
  6439                      4. * pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) +
  6440                      pow(-1. * pow(m_rho, 2) + s + tmax, 2))) -
  6442           (
eta2 * (pow(
Gammaa1, 2) * pow(a1_mass, 2) * pow(m_rho, 2) *
  6443                        (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6445                        (-0.5 * pow(m_rho, 2) + 1. * 
C4 * pow(m_rho, 4)) +
  6447                        (pow(
pion_mass, 2) * (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6448                         (0.25 - 0.125 * 
delta) * (pow(m_rho, 2) + s)) +
  6450                        (1. * 
C4 * pow(m_rho, 6) +
  6452                             (1. * pow(m_rho, 2) - 2. * 
C4 * pow(m_rho, 4)) -
  6453                         0.25 * 
delta * pow(s, 2) +
  6454                         pow(m_rho, 4) * (-0.75 + 0.125 * 
delta - 2. * 
C4 * s) +
  6456                             (0.25 + 0.375 * 
delta + 1. * 
C4 * s))) +
  6458                (pow(
Gammaa1, 2) * pow(a1_mass, 2) * pow(m_rho, 2) *
  6459                     (-0.5 + 1. * 
C4 * pow(m_rho, 2)) +
  6461                     (0.5 * pow(m_rho, 2) - 1. * 
C4 * pow(m_rho, 4)) +
  6463                     (-0.5 * pow(m_rho, 4) + 1. * 
C4 * pow(m_rho, 6) +
  6465                          (-1. * pow(m_rho, 2) + 2. * 
C4 * pow(m_rho, 4)) +
  6466                      0.25 * 
delta * pow(s, 2) -
  6467                      1. * 
C4 * pow(m_rho, 2) * pow(s, 2)) +
  6469                     (pow(
pion_mass, 4) * (0.5 - 1. * 
C4 * pow(m_rho, 2)) +
  6470                      s * ((-0.25 + 0.125 * 
delta) * pow(m_rho, 2) +
  6471                           (0.25 - 0.125 * 
delta) * s) +
  6473                          (-2. * 
C4 * pow(m_rho, 4) + (-0.5 - 0.25 * 
delta) * s +
  6474                           pow(m_rho, 2) * (1. + 2. * 
C4 * s))))) *
  6475           log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) +
  6476                   pow(pow(a1_mass, 2) - 2. * pow(
pion_mass, 2) -
  6477                           1. * pow(m_rho, 2) + s + tmax,
  6482                (0.5 * pow(
Gammaa1, 4) * pow(a1_mass, 4) -
  6483                 0.5 * pow(a1_mass, 8) + 0.5 * pow(
pion_mass, 8) +
  6484                 0.5 * pow(a1_mass, 4) * pow(
pion_mass, 2) * pow(m_rho, 2) -
  6485                 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
  6488                          (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
  6490                      pow(a1_mass, 4) * (-1. * pow(
pion_mass, 2) +
  6491                                         0.5 * pow(m_rho, 2) - 1. * s)) +
  6493                     (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
  6495                     (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
  6496            eta1 * (pow(a1_mass, 6) * (1. * pow(
pion_mass, 2) + 0.5 * s) +
  6497                    pow(a1_mass, 2) * (3. * pow(
pion_mass, 6) +
  6498                                       1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
  6502                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
  6503                         0.5 * pow(m_rho, 2) * s) +
  6505                        (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
  6506                         pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
  6507                         0.5 * pow(m_rho, 2) * s) +
  6508                    pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6510                         pow(a1_mass, 2) * (1. * pow(
pion_mass, 2) + 0.5 * s) -
  6511                         0.5 * pow(m_rho, 2) * s +
  6512                         pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
  6513           log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  6515                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6516                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6517                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6518                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6520                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6522              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) -
  6526                    1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
  6527                    pow(a1_mass, 6) * (1. * pow(
pion_mass, 2) - 0.5 * s) -
  6529                    4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
  6530                    1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
  6532                    4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
  6533                    0.5 * pow(m_rho, 4) * pow(s, 2) -
  6535                    1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  6536                    pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6538                         pow(a1_mass, 2) * (1. * pow(
pion_mass, 2) - 0.5 * s) -
  6539                         0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
  6542                         (1. * pow(m_rho, 2) - 0.5 * s) * s +
  6543                         pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
  6546                         pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
  6547                         0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
  6549                             (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
  6552                (0.5 * pow(
Gammaa1, 4) * pow(a1_mass, 4) -
  6553                 0.5 * pow(a1_mass, 8) + 0.5 * pow(
pion_mass, 8) +
  6555                     (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
  6557                 pow(a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
  6558                                    (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
  6559                 pow(
pion_mass, 4) * (-0.5 * pow(m_rho, 4) +
  6560                                      2. * pow(m_rho, 2) * s - 1.5 * pow(s, 2)) +
  6562                     (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6564                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6565                     (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
  6566                      pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
  6567                      1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
  6568                      pow(a1_mass, 2) * (-1. * pow(
pion_mass, 2) -
  6569                                         1. * pow(m_rho, 2) + 1.5 * s)) +
  6572                      pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
  6574                          (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  6576                      s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
  6577                           0.5 * pow(s, 2))))) *
  6578           log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  6580                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6581                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6582                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6583                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6585                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6587              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  6589               pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
  6590               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6592                   (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
  6593          (0.0625 * pow(
eta1 - 1. * 
eta2, 2) *
  6595                (-1. * pow(a1_mass, 10) +
  6597                     (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  6598                 pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  6599                     (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  6600                      0.5 * pow(m_rho, 2) + 0.5 * s) +
  6602                     (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
  6603                      pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
  6604                      6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
  6608                      (1. * pow(m_rho, 2) - 1. * s) * s +
  6609                      pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  6611                     (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  6612                      pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
  6613                      1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
  6615                          (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
  6616                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6617                     (4. * pow(a1_mass, 6) - 4. * pow(
pion_mass, 6) -
  6618                      0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  6619                      1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  6620                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  6621                      pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) -
  6622                                         6. * pow(m_rho, 2) + 6. * s) +
  6624                          (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  6627                          (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  6628                           pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  6629                           6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  6631                     (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
  6632                      3.5 * pow(m_rho, 6) * s + 4.5 * pow(m_rho, 4) * pow(s, 2) -
  6633                      2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
  6634                      pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  6636                          (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
  6639                          (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
  6640                           7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
  6642                (-1. * pow(a1_mass, 10) +
  6644                     (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
  6645                 pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  6646                     (1. * pow(a1_mass, 2) - 1. * pow(
pion_mass, 2) -
  6647                      0.5 * pow(m_rho, 2) + 0.5 * s) +
  6649                     (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
  6650                      5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
  6651                      pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
  6653                     (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
  6654                      pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
  6655                      3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
  6658                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
  6659                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6660                     (4. * pow(a1_mass, 6) - 4. * pow(
pion_mass, 6) -
  6661                      0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
  6662                      1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
  6663                      pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
  6664                      pow(a1_mass, 4) * (-12. * pow(
pion_mass, 2) -
  6665                                         6. * pow(m_rho, 2) + 6. * s) +
  6667                          (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
  6670                          (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
  6671                           pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
  6672                           6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
  6675                      pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
  6677                          (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
  6679                          (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
  6682                          (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
  6683                           4.5 * pow(m_rho, 2) * pow(s, 2))) +
  6686                      pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
  6688                          (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
  6691                          (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
  6692                           3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
  6693                      s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
  6694                           0.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3)))) +
  6696                (2. * pow(a1_mass, 10) +
  6697                 pow(
Gammaa1, 4) * pow(a1_mass, 4) *
  6698                     (-2. * pow(a1_mass, 2) + 2. * pow(
pion_mass, 2) +
  6699                      1. * pow(m_rho, 2) - 1. * s) +
  6701                     (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
  6703                     (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
  6704                      pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
  6705                      8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
  6707                     (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
  6708                      6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
  6709                      pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
  6711                          (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
  6713                     (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
  6714                      3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
  6715                      1. * pow(m_rho, 2) * pow(s, 3) +
  6716                      pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
  6718                          (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
  6720                          (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
  6721                           6. * pow(m_rho, 2) * pow(s, 2))) +
  6723                     (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
  6724                      pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
  6725                      2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
  6728                          (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
  6731                          (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
  6733                 pow(
Gammaa1, 2) * pow(a1_mass, 2) *
  6734                     (-8. * pow(a1_mass, 6) + 8. * pow(
pion_mass, 6) +
  6735                      1. * pow(m_rho, 6) +
  6736                      pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
  6737                      pow(a1_mass, 4) * (24. * pow(
pion_mass, 2) +
  6738                                         12. * pow(m_rho, 2) - 12. * s) -
  6739                      3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
  6742                          (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
  6745                          (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
  6746                           12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
  6748                               (-24. * pow(m_rho, 2) + 24. * s))))) *
  6749           log(abs(pow(
Gammaa1, 2) * pow(a1_mass, 2) + pow(a1_mass, 4) +
  6751                   4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
  6752                   4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
  6753                   pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
  6754                   2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
  6756                       (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
  6758              (pow(
Gammaa1, 2) * pow(a1_mass, 2) + 4. * pow(a1_mass, 4) +
  6760               pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
  6761               2. * pow(m_rho, 2) * s + pow(s, 2) +
  6763                   (-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