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