4445 const double s_sqrt = sqrt(s);
4446 const double spin_deg_factor = 1.0;
4449 double tmax = mandelstam_t[0];
4450 double tmin = mandelstam_t[1];
4456 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4459 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4462 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4463 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4465 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4466 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4470 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4471 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4473 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4475 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
4476 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4477 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4480 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4481 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
4485 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
4487 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
4488 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4489 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4492 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4493 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4494 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4497 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4498 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4500 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
4501 (1. * pow(
a1_mass, 2) - 1. * tmin) +
4503 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4506 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4507 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
4513 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
4517 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4518 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
4519 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4520 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
4536 0.5 * pow(m_rho, 2) + 1. * s)) *
4538 (0.25 * (-2. + 1. *
delta) *
4539 (-8. *
C4 * pow(m_rho, 4) +
4541 (-2. - 3. *
delta) * s +
4542 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
4546 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
4547 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
4548 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
4549 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
4550 32 * pow(
C4, 2) * pow(s, 2))) *
4555 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
4556 1. *
C4 * pow(m_rho, 6) +
4558 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4560 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
4561 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
4562 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
4563 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4564 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
4566 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
4568 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4569 0.25 *
delta * pow(s, 2) +
4570 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4573 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4576 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4577 2. * pow(m_rho, 4) +
4579 pow(m_rho, 2) - 1. * s) -
4580 3. * pow(m_rho, 2) * s + pow(s, 2)) +
4583 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4585 pow(m_rho, 2) - 1. * s) +
4586 pow(m_rho, 2) * s + pow(s, 2)) +
4590 4. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4591 2. * pow(m_rho, 4) - 2. * pow(s, 2) +
4593 2. * pow(m_rho, 2) + 2. * s))) *
4595 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4598 2. * pow(m_rho, 4) +
4600 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4605 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4606 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4608 4. * pow(m_rho, 2) + 4. * s)) +
4612 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4613 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4615 4. * pow(m_rho, 2) + 4. * s))) *
4617 (0.125 * (-2. + 1. *
delta) *
4618 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmin, 2)) /
4620 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
4621 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmin, 2) -
4623 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
4624 1. *
C4 *
delta * pow(m_rho, 2)) *
4627 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
4628 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
4630 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4632 1. * pow(m_rho, 2) - 1. * s) +
4634 1. * pow(m_rho, 2) + s)) *
4636 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmin, 3) -
4637 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmin, 3) +
4638 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4642 2. * pow(m_rho, 4) +
4644 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
4648 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4649 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4651 4. * pow(m_rho, 2) + 4. * s)) +
4655 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4656 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4658 4. * pow(m_rho, 2) + 4. * s))) *
4659 (-1. * pow(m_rho, 2) + s + tmin) -
4660 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4661 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
4662 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
4663 pow(-1. * pow(m_rho, 2) + s + tmin, 2) +
4664 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
4665 pow(-1. * pow(m_rho, 2) + s + tmin, 3) +
4667 (-1. + 2. *
C4 * pow(m_rho, 2)) *
4669 1. * pow(m_rho, 2) + s + tmin,
4672 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4673 0.5 *
C4 * pow(m_rho, 6) +
4675 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4677 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4678 0.125 * pow(m_rho, 2) * s -
4679 0.1875 *
delta * pow(m_rho, 2) * s +
4680 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4681 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4682 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4684 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4686 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4687 0.125 *
delta * pow(s, 2) +
4688 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4690 1. * s + 1. * tmin)) /
4693 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4694 0.5 *
C4 * pow(m_rho, 6) +
4696 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4698 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4699 0.125 * pow(m_rho, 2) * s -
4700 0.1875 *
delta * pow(m_rho, 2) * s +
4701 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4702 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4703 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4705 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4707 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4708 0.125 *
delta * pow(s, 2) +
4709 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4718 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4720 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4721 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
4723 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
4730 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4731 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
4732 1. * pow(m_rho, 2) * s))) *
4741 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4742 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4744 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
4746 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4747 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
4748 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
4752 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4753 1. * pow(m_rho, 2) * s +
4754 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
4757 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
4758 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4763 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
4764 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
4765 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
4768 (1. * pow(m_rho, 2) - 1. * s) * s +
4769 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
4775 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
4776 2. * pow(m_rho, 2) * s + pow(s, 2) +
4778 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
4779 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4785 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4787 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4788 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4789 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4792 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4793 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
4796 1. * pow(m_rho, 4) +
4798 6. * pow(m_rho, 2) - 6. * s) +
4799 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4800 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
4804 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4807 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4808 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4810 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4811 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4815 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
4816 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
4821 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
4823 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4824 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4825 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4828 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4829 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4830 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4833 1. * pow(m_rho, 4) +
4835 6. * pow(m_rho, 2) - 6. * s) +
4836 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4837 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
4840 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4841 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4843 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
4852 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4853 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4854 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
4856 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4857 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4858 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
4859 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4861 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
4863 28. * pow(m_rho, 2) + 28. * s) +
4865 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
4866 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
4867 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
4871 4. * pow(m_rho, 2) - 4. * s) +
4872 (4. * pow(m_rho, 2) - 4. * s) * s +
4873 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
4875 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
4876 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
4878 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
4880 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
4881 12. * pow(s, 2)))) +
4885 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4886 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
4887 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
4889 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4891 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
4892 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
4893 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
4894 1. * pow(m_rho, 4) * pow(s, 2) -
4896 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4899 1. * pow(m_rho, 4) +
4900 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4901 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4903 2. * pow(m_rho, 2) + 2. * s)) +
4905 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
4906 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4907 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4909 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
4910 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
4911 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
4914 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
4919 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4920 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4921 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
4923 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4925 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4926 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4927 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
4928 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4929 9. * pow(m_rho, 4) * pow(s, 2) -
4931 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4934 3. * pow(m_rho, 4) +
4935 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4936 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4938 2. * pow(m_rho, 2) + 2. * s)) +
4940 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
4941 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4942 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4944 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
4945 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
4946 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
4949 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
4950 6. * pow(s, 2))))) *
4956 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
4957 2. * pow(m_rho, 2) * s + pow(s, 2) +
4959 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
4960 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4965 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
4966 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
4968 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4969 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4974 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
4975 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
4977 3. * pow(m_rho, 2) + 3. * s) +
4979 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4980 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4981 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
4984 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4985 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
4987 3. * pow(m_rho, 2) + 3. * s) +
4989 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4990 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4991 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
4992 log(abs(-1. * pow(
a1_mass, 2) + tmin)) -
4996 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
4998 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5000 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
5003 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5004 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5005 0.5 * pow(m_rho, 2) * s) +
5007 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
5008 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5009 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5014 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
5015 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
5016 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
5019 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5020 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
5023 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
5024 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
5025 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
5027 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
5028 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5029 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5031 1. * pow(m_rho, 2) + 1. * s) -
5032 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5034 0.5 * pow(m_rho, 2) + 0.5 * s) *
5038 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
5040 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
5041 4. * pow(m_rho, 2) * s +
5042 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
5045 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5046 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
5047 2. * pow(m_rho, 2) * pow(s, 2) +
5049 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
5051 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
5052 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
5053 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
5056 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
5061 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5063 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5064 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5066 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5067 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
5068 pow(m_rho, 2) * s) +
5070 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
5071 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
5073 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5075 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
5079 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5081 (-2. * pow(m_rho, 2) + 2. * s) +
5083 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5085 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
5087 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5088 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5089 3. * pow(m_rho, 2) * s) +
5091 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
5093 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5095 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
5096 2. * pow(s, 2))))) *
5097 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5100 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
5101 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
5103 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
5107 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5109 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5110 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5112 (-1. *
C4 * pow(m_rho, 6) +
5114 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
5115 0.25 *
delta * pow(s, 2) +
5117 (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
5118 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
5121 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5123 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5125 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5126 0.25 *
delta * pow(s, 2) +
5127 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5129 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5130 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
5131 (-0.25 + 0.125 *
delta) * s) +
5133 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
5134 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
5135 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5138 log(abs(-1. * pow(
pion_mass, 2) + tmin)) +
5141 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5143 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
5145 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
5146 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5147 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5149 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
5150 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
5151 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
5152 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5154 (2. * (-2. + 1. *
delta) *
5155 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
5157 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
5158 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5161 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
5163 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5164 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5166 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
5168 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5170 eta2 * (0.5 * pow(m_rho, 6) +
5171 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
5172 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
5175 (1.5 * pow(m_rho, 4) - 3. * pow(m_rho, 2) * s +
5178 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5180 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5181 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5184 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5185 2. * pow(m_rho, 2) * s + pow(s, 2) +
5187 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5189 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s +
5194 (-1. * pow(m_rho, 2) + 1. * s) +
5196 (-0.5 * pow(m_rho, 4) +
5197 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5198 0.5 * pow(m_rho, 2) * s) +
5200 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
5201 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5202 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5206 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5207 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5208 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5210 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5211 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5213 1. * pow(m_rho, 2) + 1. * s) -
5216 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5218 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
5220 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
5223 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
5225 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
5227 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5228 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5229 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5237 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
5238 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
5240 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
5241 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
5245 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5247 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
5249 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5250 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
5253 1. * pow(m_rho, 2) + s) +
5255 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5256 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5257 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
5260 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5262 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
5267 pow(m_rho, 2) + s) +
5269 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5270 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5271 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
5275 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmin) -
5277 (-1. * pow(m_rho, 2) + s + tmin) +
5278 pow(-1. * pow(m_rho, 2) + s + tmin, 2))) -
5281 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5283 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5285 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5286 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5288 (1. *
C4 * pow(m_rho, 6) +
5290 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5291 0.25 *
delta * pow(s, 2) +
5292 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
5294 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
5296 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5298 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5300 (-0.5 * pow(m_rho, 4) + 1. *
C4 * pow(m_rho, 6) +
5302 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
5303 0.25 *
delta * pow(s, 2) -
5304 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5306 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5307 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
5308 (0.25 - 0.125 *
delta) * s) +
5310 (-2. *
C4 * pow(m_rho, 4) +
5311 (-0.5 - 0.25 *
delta) * s +
5312 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
5315 1. * pow(m_rho, 2) + s + tmin,
5323 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
5326 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
5329 0.5 * pow(m_rho, 2) - 1. * s)) +
5331 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5333 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
5337 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
5341 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5342 0.5 * pow(m_rho, 2) * s) +
5344 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5345 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
5346 0.5 * pow(m_rho, 2) * s) +
5350 0.5 * pow(m_rho, 2) * s +
5351 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
5354 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5355 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5356 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5357 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5359 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5365 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5368 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
5369 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
5371 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
5372 0.5 * pow(m_rho, 4) * pow(s, 2) -
5374 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5378 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
5381 (1. * pow(m_rho, 2) - 0.5 * s) * s +
5382 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
5385 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
5386 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
5388 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
5394 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
5396 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
5397 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
5399 (-0.5 * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5402 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5405 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
5406 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
5407 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
5409 1. * pow(m_rho, 2) + 1.5 * s)) +
5412 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
5414 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5416 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5417 0.5 * pow(s, 2))))) *
5420 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5421 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5422 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5423 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5425 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5429 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5430 2. * pow(m_rho, 2) * s + pow(s, 2) +
5432 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5433 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5437 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5440 0.5 * pow(m_rho, 2) + 0.5 * s) +
5442 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
5443 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
5444 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
5448 (1. * pow(m_rho, 2) - 1. * s) * s +
5449 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5451 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5452 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
5453 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
5455 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
5458 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5459 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5460 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5462 6. * pow(m_rho, 2) + 6. * s) +
5464 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5467 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5468 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5469 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5471 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
5472 3.5 * pow(m_rho, 6) * s +
5473 4.5 * pow(m_rho, 4) * pow(s, 2) -
5474 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5475 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5477 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
5480 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
5481 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
5485 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5488 0.5 * pow(m_rho, 2) + 0.5 * s) +
5490 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
5491 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
5492 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5494 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5495 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
5496 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
5499 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
5502 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5503 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5504 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5506 6. * pow(m_rho, 2) + 6. * s) +
5508 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5511 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5512 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5513 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5516 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
5518 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
5520 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5523 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
5524 4.5 * pow(m_rho, 2) * pow(s, 2))) +
5527 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5529 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
5532 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
5533 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
5534 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
5535 0.5 * pow(m_rho, 2) * pow(s, 2) +
5536 0.5 * pow(s, 3)))) +
5541 1. * pow(m_rho, 2) - 1. * s) +
5543 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
5545 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
5546 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
5547 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5549 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
5550 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
5551 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
5553 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
5555 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
5556 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
5557 1. * pow(m_rho, 2) * pow(s, 3) +
5558 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
5560 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
5562 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
5563 6. * pow(m_rho, 2) * pow(s, 2))) +
5565 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
5566 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
5567 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
5570 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
5573 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
5577 1. * pow(m_rho, 6) +
5578 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
5580 12. * pow(m_rho, 2) - 12. * s) -
5581 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
5584 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
5587 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
5588 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
5590 (-24. * pow(m_rho, 2) + 24. * s))))) *
5593 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5594 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5595 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5596 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5598 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5602 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
5603 2. * pow(m_rho, 2) * s + pow(s, 2) +
5605 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
5609 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5612 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5615 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5616 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5618 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5619 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5623 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5624 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5626 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5628 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5629 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5630 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5633 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5634 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
5638 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5640 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5641 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5642 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5645 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5646 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5647 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
5650 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5651 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
5653 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
5654 (1. * pow(
a1_mass, 2) - 1. * tmax) +
5656 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5659 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5660 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
5666 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
5670 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5671 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
5672 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5673 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
5688 0.5 * pow(m_rho, 2) + 1. * s)) *
5690 (0.25 * (-2. + 1. *
delta) *
5691 (-8. *
C4 * pow(m_rho, 4) +
5693 (-2. - 3. *
delta) * s +
5694 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
5698 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
5699 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
5700 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
5701 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
5702 32 * pow(
C4, 2) * pow(s, 2))) *
5707 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
5708 1. *
C4 * pow(m_rho, 6) +
5710 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5712 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
5713 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
5714 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
5715 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5716 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5718 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
5720 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5721 0.25 *
delta * pow(s, 2) +
5722 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5725 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5729 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) +
5731 (2. * pow(
pion_mass, 2) + pow(m_rho, 2) - 1. * s) -
5732 3. * pow(m_rho, 2) * s + pow(s, 2)) +
5735 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5737 pow(m_rho, 2) - 1. * s) +
5738 pow(m_rho, 2) * s + pow(s, 2)) +
5742 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) -
5745 2. * pow(m_rho, 2) + 2. * s))) *
5747 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5750 2. * pow(m_rho, 4) +
5752 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5757 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5758 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5760 4. * pow(m_rho, 2) + 4. * s)) +
5764 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5765 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5767 4. * pow(m_rho, 2) + 4. * s))) *
5769 (0.125 * (-2. + 1. *
delta) *
5770 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmax, 2)) /
5772 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
5773 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmax, 2) -
5775 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
5776 1. *
C4 *
delta * pow(m_rho, 2)) *
5779 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
5780 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
5782 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5784 1. * pow(m_rho, 2) - 1. * s) +
5786 1. * pow(m_rho, 2) + s)) *
5788 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmax, 3) -
5789 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmax, 3) +
5790 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5794 2. * pow(m_rho, 4) +
5796 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
5800 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5801 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5803 4. * pow(m_rho, 2) + 4. * s)) +
5807 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5808 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5810 4. * pow(m_rho, 2) + 4. * s))) *
5811 (-1. * pow(m_rho, 2) + s + tmax) -
5812 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5813 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
5814 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
5815 pow(-1. * pow(m_rho, 2) + s + tmax, 2) +
5816 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
5817 pow(-1. * pow(m_rho, 2) + s + tmax, 3) +
5819 (-1. + 2. *
C4 * pow(m_rho, 2)) *
5825 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5826 0.5 *
C4 * pow(m_rho, 6) +
5828 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5830 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5831 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5832 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5833 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5834 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5836 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5838 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5839 0.125 *
delta * pow(s, 2) +
5840 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5842 1. * s + 1. * tmax)) /
5846 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5847 0.5 *
C4 * pow(m_rho, 6) +
5849 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5851 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5852 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5853 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5854 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5855 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5857 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5859 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5860 0.125 *
delta * pow(s, 2) +
5861 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5870 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5872 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5873 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
5875 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
5882 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5883 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
5884 1. * pow(m_rho, 2) * s))) *
5893 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5895 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
5897 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5898 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
5899 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5903 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5904 1. * pow(m_rho, 2) * s +
5905 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
5908 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5909 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5914 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5915 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
5916 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
5919 (1. * pow(m_rho, 2) - 1. * s) * s +
5920 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
5926 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5927 2. * pow(m_rho, 2) * s + pow(s, 2) +
5929 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
5930 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5936 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5938 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5939 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5940 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5943 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5944 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
5947 1. * pow(m_rho, 4) +
5949 6. * pow(m_rho, 2) - 6. * s) +
5950 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5951 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
5955 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5958 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5959 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5961 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5962 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5966 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
5967 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
5972 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5974 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5975 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5976 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5979 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5980 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5981 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
5984 1. * pow(m_rho, 4) +
5986 6. * pow(m_rho, 2) - 6. * s) +
5987 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5988 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
5991 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5992 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
5994 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
6003 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6004 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6005 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
6007 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6008 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6009 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
6010 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6012 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
6014 (-56. * pow(
pion_mass, 2) - 28. * pow(m_rho, 2) + 28. * s) +
6016 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
6017 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
6018 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
6022 4. * pow(m_rho, 2) - 4. * s) +
6023 (4. * pow(m_rho, 2) - 4. * s) * s +
6024 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
6026 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
6027 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
6029 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
6031 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
6032 12. * pow(s, 2)))) +
6036 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6037 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
6038 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
6040 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6042 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
6043 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6044 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
6045 1. * pow(m_rho, 4) * pow(s, 2) -
6047 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6050 1. * pow(m_rho, 4) +
6051 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6052 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6054 2. * pow(m_rho, 2) + 2. * s)) +
6056 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
6057 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6058 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6060 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6061 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
6062 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
6065 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
6070 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6071 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6072 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
6074 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6076 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6077 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6078 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
6079 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6080 9. * pow(m_rho, 4) * pow(s, 2) -
6082 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6085 3. * pow(m_rho, 4) +
6086 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6087 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6089 2. * pow(m_rho, 2) + 2. * s)) +
6091 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
6092 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6093 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6095 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6096 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
6097 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
6100 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6101 6. * pow(s, 2))))) *
6107 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6108 2. * pow(m_rho, 2) * s + pow(s, 2) +
6110 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
6111 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6116 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
6117 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
6119 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
6120 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
6125 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
6126 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
6128 (-6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) + 3. * s) +
6130 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6131 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6132 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
6135 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
6136 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
6138 (-6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) + 3. * s) +
6140 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6141 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
6142 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
6143 log(abs(-1. * pow(
a1_mass, 2) + tmax)) -
6147 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
6149 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6151 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
6154 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
6155 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6156 0.5 * pow(m_rho, 2) * s) +
6158 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
6159 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6160 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6165 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
6166 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
6167 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
6170 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6171 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
6174 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
6175 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
6176 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
6178 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
6179 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6180 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6182 1. * pow(m_rho, 2) + 1. * s) -
6183 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6189 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
6191 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
6192 4. * pow(m_rho, 2) * s +
6193 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
6196 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
6197 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
6198 2. * pow(m_rho, 2) * pow(s, 2) +
6200 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
6202 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6203 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
6204 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
6207 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
6212 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6214 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
6215 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6217 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6218 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
6219 pow(m_rho, 2) * s) +
6221 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
6222 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
6224 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6226 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
6230 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6232 (-2. * pow(m_rho, 2) + 2. * s) +
6234 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6236 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
6238 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6239 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6240 3. * pow(m_rho, 2) * s) +
6242 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
6244 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6246 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
6247 2. * pow(s, 2))))) *
6248 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6251 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
6252 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
6254 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
6258 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6260 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6261 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6263 (-1. *
C4 * pow(m_rho, 6) +
6265 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6266 0.25 *
delta * pow(s, 2) +
6267 pow(m_rho, 2) * s * (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
6268 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
6271 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6273 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
6275 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6276 0.25 *
delta * pow(s, 2) +
6277 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6279 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6280 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
6281 (-0.25 + 0.125 *
delta) * s) +
6283 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
6284 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
6285 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6288 log(abs(-1. * pow(
pion_mass, 2) + tmax)) +
6291 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6293 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
6295 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
6296 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6297 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6299 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
6300 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
6301 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
6302 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6304 (2. * (-2. + 1. *
delta) *
6305 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
6307 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
6308 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6311 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
6313 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6314 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6316 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
6318 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6321 (0.5 * pow(m_rho, 6) +
6322 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
6323 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
6325 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 4) -
6326 3. * pow(m_rho, 2) * s + 1.5 * pow(s, 2)) +
6328 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6330 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6331 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6334 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6335 2. * pow(m_rho, 2) * s + pow(s, 2) +
6337 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6339 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax)) +
6343 (-1. * pow(m_rho, 2) + 1. * s) +
6345 (-0.5 * pow(m_rho, 4) +
6346 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6347 0.5 * pow(m_rho, 2) * s) +
6349 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
6350 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6351 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6355 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6356 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6357 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6359 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6360 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6362 1. * pow(m_rho, 2) + 1. * s) -
6365 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6367 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
6369 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
6372 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
6374 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
6376 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6377 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6378 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6386 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
6387 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
6389 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
6390 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
6394 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6396 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
6398 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
6399 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
6402 1. * pow(m_rho, 2) + s) +
6404 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6405 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6406 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
6409 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6411 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
6416 pow(m_rho, 2) + s) +
6418 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6419 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6420 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
6424 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) -
6425 4. * pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) +
6426 pow(-1. * pow(m_rho, 2) + s + tmax, 2))) -
6429 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6431 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6433 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6434 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6436 (1. *
C4 * pow(m_rho, 6) +
6438 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6439 0.25 *
delta * pow(s, 2) +
6440 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
6442 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
6445 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6447 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6449 (-0.5 * pow(m_rho, 4) + 1. *
C4 * pow(m_rho, 6) +
6451 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6452 0.25 *
delta * pow(s, 2) -
6453 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6455 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6456 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
6457 (0.25 - 0.125 *
delta) * s) +
6459 (-2. *
C4 * pow(m_rho, 4) + (-0.5 - 0.25 *
delta) * s +
6460 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
6463 1. * pow(m_rho, 2) + s + tmax,
6471 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
6474 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
6477 0.5 * pow(m_rho, 2) - 1. * s)) +
6479 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6481 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
6484 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
6488 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6489 0.5 * pow(m_rho, 2) * s) +
6491 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6492 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
6493 0.5 * pow(m_rho, 2) * s) +
6497 0.5 * pow(m_rho, 2) * s +
6498 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
6501 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6502 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6503 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6504 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6506 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6512 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6515 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6516 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6518 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6519 0.5 * pow(m_rho, 4) * pow(s, 2) -
6521 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6525 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
6528 (1. * pow(m_rho, 2) - 0.5 * s) * s +
6529 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
6532 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
6533 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
6535 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
6541 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
6543 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
6544 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
6545 pow(
pion_mass, 4) * (-0.5 * pow(m_rho, 4) +
6546 2. * pow(m_rho, 2) * s - 1.5 * pow(s, 2)) +
6548 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6551 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
6552 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
6553 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
6555 1. * pow(m_rho, 2) + 1.5 * s)) +
6558 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
6560 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6562 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6563 0.5 * pow(s, 2))))) *
6566 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6567 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6568 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6569 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6571 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6575 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6576 2. * pow(m_rho, 2) * s + pow(s, 2) +
6578 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6579 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6583 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6586 0.5 * pow(m_rho, 2) + 0.5 * s) +
6588 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
6589 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
6590 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
6594 (1. * pow(m_rho, 2) - 1. * s) * s +
6595 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6597 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6598 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
6599 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
6601 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
6604 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6605 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6606 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6608 6. * pow(m_rho, 2) + 6. * s) +
6610 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6613 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6614 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6615 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6617 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
6618 3.5 * pow(m_rho, 6) * s + 4.5 * pow(m_rho, 4) * pow(s, 2) -
6619 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6620 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6622 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
6625 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
6626 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
6630 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6633 0.5 * pow(m_rho, 2) + 0.5 * s) +
6635 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
6636 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
6637 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6639 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6640 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
6641 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
6644 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
6647 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6648 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6649 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6651 6. * pow(m_rho, 2) + 6. * s) +
6653 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6656 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6657 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6658 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6661 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
6663 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
6665 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6668 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
6669 4.5 * pow(m_rho, 2) * pow(s, 2))) +
6672 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6674 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
6677 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
6678 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
6679 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
6680 0.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3)))) +
6685 1. * pow(m_rho, 2) - 1. * s) +
6687 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
6689 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
6690 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
6691 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
6693 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6694 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
6695 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
6697 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
6699 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
6700 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
6701 1. * pow(m_rho, 2) * pow(s, 3) +
6702 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
6704 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
6706 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
6707 6. * pow(m_rho, 2) * pow(s, 2))) +
6709 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
6710 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
6711 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
6714 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
6717 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
6721 1. * pow(m_rho, 6) +
6722 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
6724 12. * pow(m_rho, 2) - 12. * s) -
6725 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
6728 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6731 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
6732 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6734 (-24. * pow(m_rho, 2) + 24. * s))))) *
6737 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6738 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6739 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6740 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6742 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6746 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6747 2. * pow(m_rho, 2) * s + pow(s, 2) +
6749 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
constexpr static double C4
Constant in the computation of photon cross sections.
constexpr static double Const
Constant in the computation of photon cross sections.
constexpr static double delta
Constant in the computation of photon cross sections.
constexpr static double Gammaa1
Constant in the computation of photon cross sections.
constexpr static double eta2
Constant in the computation of photon cross sections.
constexpr static double ghat
Constant in the computation of photon cross sections.
constexpr static double Pi
Value of .
constexpr static double eta1
Constant in the computation of photon cross sections.
constexpr double gev2_mb
GeV^-2 <-> mb conversion factor.
constexpr double a1_mass
a1 mass in GeV.
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 pion_mass
Pion mass in GeV.
double cut_off(const double sigma_mb)
Cross section after cut off.