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