4447 const double s_sqrt = sqrt(s);
4448 const double spin_deg_factor = 1.0;
4451 double tmax = mandelstam_t[0];
4452 double tmin = mandelstam_t[1];
4458 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4461 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4464 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4465 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4467 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4468 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4472 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4473 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4475 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4477 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
4478 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4479 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4482 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4483 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
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 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4494 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4495 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4496 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4499 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4500 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4502 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
4503 (1. * pow(
a1_mass, 2) - 1. * tmin) +
4505 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4508 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
4509 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
4515 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
4519 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4520 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
4521 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
4522 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
4538 0.5 * pow(m_rho, 2) + 1. * s)) *
4540 (0.25 * (-2. + 1. *
delta) *
4541 (-8. *
C4 * pow(m_rho, 4) +
4543 (-2. - 3. *
delta) * s +
4544 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
4548 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
4549 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
4550 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
4551 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
4552 32 * pow(
C4, 2) * pow(s, 2))) *
4557 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
4558 1. *
C4 * pow(m_rho, 6) +
4560 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4562 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
4563 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
4564 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
4565 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4566 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
4568 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
4570 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4571 0.25 *
delta * pow(s, 2) +
4572 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4575 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4578 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4579 2. * pow(m_rho, 4) +
4581 pow(m_rho, 2) - 1. * s) -
4582 3. * pow(m_rho, 2) * s + pow(s, 2)) +
4585 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4587 pow(m_rho, 2) - 1. * s) +
4588 pow(m_rho, 2) * s + pow(s, 2)) +
4592 4. * pow(
pion_mass, 2) * pow(m_rho, 2) +
4593 2. * pow(m_rho, 4) - 2. * pow(s, 2) +
4595 2. * pow(m_rho, 2) + 2. * s))) *
4597 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4600 2. * pow(m_rho, 4) +
4602 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4607 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4608 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4610 4. * pow(m_rho, 2) + 4. * s)) +
4614 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4615 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4617 4. * pow(m_rho, 2) + 4. * s))) *
4619 (0.125 * (-2. + 1. *
delta) *
4620 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmin, 2)) /
4622 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
4623 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmin, 2) -
4625 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
4626 1. *
C4 *
delta * pow(m_rho, 2)) *
4629 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
4630 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
4632 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4634 1. * pow(m_rho, 2) - 1. * s) +
4636 1. * pow(m_rho, 2) + s)) *
4638 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmin, 3) -
4639 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmin, 3) +
4640 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4644 2. * pow(m_rho, 4) +
4646 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
4650 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4651 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4653 4. * pow(m_rho, 2) + 4. * s)) +
4657 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
4658 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4660 4. * pow(m_rho, 2) + 4. * s))) *
4661 (-1. * pow(m_rho, 2) + s + tmin) -
4662 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
4663 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
4664 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
4665 pow(-1. * pow(m_rho, 2) + s + tmin, 2) +
4666 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
4667 pow(-1. * pow(m_rho, 2) + s + tmin, 3) +
4669 (-1. + 2. *
C4 * pow(m_rho, 2)) *
4671 1. * pow(m_rho, 2) + s + tmin,
4674 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4675 0.5 *
C4 * pow(m_rho, 6) +
4677 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4679 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4680 0.125 * pow(m_rho, 2) * s -
4681 0.1875 *
delta * pow(m_rho, 2) * s +
4682 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4683 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4684 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4686 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4688 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4689 0.125 *
delta * pow(s, 2) +
4690 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4692 1. * s + 1. * tmin)) /
4695 (
eta2 * (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
4696 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 * pow(m_rho, 2) * s -
4702 0.1875 *
delta * pow(m_rho, 2) * s +
4703 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
4704 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
4705 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
4707 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
4709 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
4710 0.125 *
delta * pow(s, 2) +
4711 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
4720 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4722 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
4723 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
4725 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
4732 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4733 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
4734 1. * pow(m_rho, 2) * s))) *
4743 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4744 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4746 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
4748 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4749 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
4750 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
4754 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
4755 1. * pow(m_rho, 2) * s +
4756 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
4759 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
4760 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
4765 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
4766 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
4767 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
4770 (1. * pow(m_rho, 2) - 1. * s) * s +
4771 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
4777 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
4778 2. * pow(m_rho, 2) * s + pow(s, 2) +
4780 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
4781 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
4787 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
4789 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4790 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4791 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4794 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4795 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
4798 1. * pow(m_rho, 4) +
4800 6. * pow(m_rho, 2) - 6. * s) +
4801 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4802 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
4806 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4809 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
4810 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
4812 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4813 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4817 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
4818 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
4823 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
4825 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4826 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4827 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
4830 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4831 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
4832 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
4835 1. * pow(m_rho, 4) +
4837 6. * pow(m_rho, 2) - 6. * s) +
4838 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
4839 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
4842 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
4843 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
4845 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
4854 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
4855 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4856 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
4858 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4859 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4860 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
4861 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4863 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
4865 28. * pow(m_rho, 2) + 28. * s) +
4867 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
4868 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
4869 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
4873 4. * pow(m_rho, 2) - 4. * s) +
4874 (4. * pow(m_rho, 2) - 4. * s) * s +
4875 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
4877 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
4878 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
4880 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
4882 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
4883 12. * pow(s, 2)))) +
4887 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4888 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
4889 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
4891 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4893 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
4894 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
4895 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
4896 1. * pow(m_rho, 4) * pow(s, 2) -
4898 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4901 1. * pow(m_rho, 4) +
4902 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4903 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4905 2. * pow(m_rho, 2) + 2. * s)) +
4907 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
4908 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4909 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4911 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
4912 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
4913 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
4916 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
4921 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
4922 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
4923 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
4925 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
4927 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
4928 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
4929 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
4930 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
4931 9. * pow(m_rho, 4) * pow(s, 2) -
4933 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
4936 3. * pow(m_rho, 4) +
4937 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
4938 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
4940 2. * pow(m_rho, 2) + 2. * s)) +
4942 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
4943 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
4944 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
4946 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
4947 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
4948 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
4951 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
4952 6. * pow(s, 2))))) *
4958 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
4959 2. * pow(m_rho, 2) * s + pow(s, 2) +
4961 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
4962 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
4967 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
4968 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
4970 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
4971 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
4976 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
4977 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
4979 3. * pow(m_rho, 2) + 3. * s) +
4981 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4982 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
4983 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
4986 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
4987 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
4989 3. * pow(m_rho, 2) + 3. * s) +
4991 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
4992 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
4993 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
4994 log(abs(-1. * pow(
a1_mass, 2) + tmin)) -
4998 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
5000 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5002 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
5005 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5006 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5007 0.5 * pow(m_rho, 2) * s) +
5009 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
5010 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5011 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5016 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
5017 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
5018 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
5021 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5022 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
5025 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
5026 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
5027 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
5029 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
5030 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5031 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5033 1. * pow(m_rho, 2) + 1. * s) -
5034 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5036 0.5 * pow(m_rho, 2) + 0.5 * s) *
5040 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
5042 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
5043 4. * pow(m_rho, 2) * s +
5044 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
5047 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5048 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
5049 2. * pow(m_rho, 2) * pow(s, 2) +
5051 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
5053 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
5054 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
5055 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
5058 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
5063 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5065 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5066 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5068 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5069 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
5070 pow(m_rho, 2) * s) +
5072 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
5073 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
5075 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5077 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
5081 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5083 (-2. * pow(m_rho, 2) + 2. * s) +
5085 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5087 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
5089 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5090 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5091 3. * pow(m_rho, 2) * s) +
5093 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
5095 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
5097 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
5098 2. * pow(s, 2))))) *
5099 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5102 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
5103 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
5105 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
5109 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5111 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5112 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5114 (-1. *
C4 * pow(m_rho, 6) +
5116 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
5117 0.25 *
delta * pow(s, 2) +
5119 (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
5120 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
5123 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5125 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5127 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5128 0.25 *
delta * pow(s, 2) +
5129 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5131 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5132 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
5133 (-0.25 + 0.125 *
delta) * s) +
5135 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
5136 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
5137 log(abs(-1. * pow(
a1_mass, 2) + tmin))) /
5140 log(abs(-1. * pow(
pion_mass, 2) + tmin)) +
5143 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5145 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
5147 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
5148 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5149 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5151 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
5152 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
5153 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
5154 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5156 (2. * (-2. + 1. *
delta) *
5157 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
5159 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
5160 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5163 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
5165 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5166 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5168 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
5170 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5172 eta2 * (0.5 * pow(m_rho, 6) +
5173 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
5174 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
5177 (1.5 * pow(m_rho, 4) - 3. * pow(m_rho, 2) * s +
5180 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5182 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
5183 log(abs(-1. * pow(
pion_mass, 2) + tmin))) /
5186 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5187 2. * pow(m_rho, 2) * s + pow(s, 2) +
5189 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5191 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s +
5196 (-1. * pow(m_rho, 2) + 1. * s) +
5198 (-0.5 * pow(m_rho, 4) +
5199 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5200 0.5 * pow(m_rho, 2) * s) +
5202 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
5203 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5204 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5208 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
5209 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
5210 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
5212 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
5213 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5215 1. * pow(m_rho, 2) + 1. * s) -
5218 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
5220 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
5222 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
5225 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
5227 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
5229 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmin))) /
5230 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
5231 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5239 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
5240 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
5242 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
5243 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
5247 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5249 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
5251 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5252 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
5255 1. * pow(m_rho, 2) + s) +
5257 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5258 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5259 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
5262 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5264 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
5269 pow(m_rho, 2) + s) +
5271 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5272 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5273 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
5277 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmin) -
5279 (-1. * pow(m_rho, 2) + s + tmin) +
5280 pow(-1. * pow(m_rho, 2) + s + tmin, 2))) -
5283 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5285 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
5287 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5288 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
5290 (1. *
C4 * pow(m_rho, 6) +
5292 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
5293 0.25 *
delta * pow(s, 2) +
5294 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
5296 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
5298 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
5300 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5302 (-0.5 * pow(m_rho, 4) + 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 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5308 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
5309 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
5310 (0.25 - 0.125 *
delta) * s) +
5312 (-2. *
C4 * pow(m_rho, 4) +
5313 (-0.5 - 0.25 *
delta) * s +
5314 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
5317 1. * pow(m_rho, 2) + s + tmin,
5325 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
5328 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
5331 0.5 * pow(m_rho, 2) - 1. * s)) +
5333 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5335 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
5339 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
5343 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
5344 0.5 * pow(m_rho, 2) * s) +
5346 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5347 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
5348 0.5 * pow(m_rho, 2) * s) +
5352 0.5 * pow(m_rho, 2) * s +
5353 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
5356 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5357 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5358 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5359 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5361 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5367 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5370 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
5371 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
5373 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
5374 0.5 * pow(m_rho, 4) * pow(s, 2) -
5376 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5380 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
5383 (1. * pow(m_rho, 2) - 0.5 * s) * s +
5384 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
5387 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
5388 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
5390 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
5396 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
5398 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
5399 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
5401 (-0.5 * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5404 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5407 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
5408 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
5409 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
5411 1. * pow(m_rho, 2) + 1.5 * s)) +
5414 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
5416 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5418 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
5419 0.5 * pow(s, 2))))) *
5422 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5423 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5424 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5425 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5427 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5431 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5432 2. * pow(m_rho, 2) * s + pow(s, 2) +
5434 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
5435 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5439 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5442 0.5 * pow(m_rho, 2) + 0.5 * s) +
5444 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
5445 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
5446 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
5450 (1. * pow(m_rho, 2) - 1. * s) * s +
5451 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5453 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5454 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
5455 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
5457 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
5460 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5461 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5462 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5464 6. * pow(m_rho, 2) + 6. * s) +
5466 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5469 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5470 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5471 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5473 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
5474 3.5 * pow(m_rho, 6) * s +
5475 4.5 * pow(m_rho, 4) * pow(s, 2) -
5476 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
5477 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5479 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
5482 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
5483 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
5487 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
5490 0.5 * pow(m_rho, 2) + 0.5 * s) +
5492 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
5493 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
5494 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
5496 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
5497 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
5498 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
5501 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
5504 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
5505 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5506 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
5508 6. * pow(m_rho, 2) + 6. * s) +
5510 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
5513 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
5514 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
5515 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
5518 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
5520 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
5522 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
5525 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
5526 4.5 * pow(m_rho, 2) * pow(s, 2))) +
5529 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
5531 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
5534 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
5535 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
5536 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
5537 0.5 * pow(m_rho, 2) * pow(s, 2) +
5538 0.5 * pow(s, 3)))) +
5543 1. * pow(m_rho, 2) - 1. * s) +
5545 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
5547 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
5548 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
5549 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5551 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
5552 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
5553 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
5555 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
5557 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
5558 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
5559 1. * pow(m_rho, 2) * pow(s, 3) +
5560 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
5562 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
5564 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
5565 6. * pow(m_rho, 2) * pow(s, 2))) +
5567 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
5568 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
5569 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
5572 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
5575 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
5579 1. * pow(m_rho, 6) +
5580 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
5582 12. * pow(m_rho, 2) - 12. * s) -
5583 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
5586 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
5589 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
5590 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
5592 (-24. * pow(m_rho, 2) + 24. * s))))) *
5595 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5596 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
5597 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmin -
5598 2. * pow(m_rho, 2) * tmin + 2. * s * tmin + pow(tmin, 2) +
5600 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
5604 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
5605 2. * pow(m_rho, 2) * s + pow(s, 2) +
5607 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
5611 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5614 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5617 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5618 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5620 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5621 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5625 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
5626 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5628 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5630 (6. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
5631 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5632 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5635 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5636 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s))) +
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 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5647 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5648 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5649 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
5652 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5653 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
5655 (1. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s))))) /
5656 (1. * pow(
a1_mass, 2) - 1. * tmax) +
5658 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5661 (1. * pow(
pion_mass, 2) - 0.25 * pow(m_rho, 2))) /
5662 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s -
5668 0.5 *
eta2 * pow(m_rho, 2) + 0.5 *
eta1 * s - 1. *
eta2 * s) *
5672 (12. + 1. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5673 delta * (-8. + 8. *
C4 * pow(m_rho, 2))) +
5674 (-4. - 3. * pow(
delta, 2) - 16. *
C4 * pow(m_rho, 2) +
5675 delta * (8. + 8. *
C4 * pow(m_rho, 2))) *
5690 0.5 * pow(m_rho, 2) + 1. * s)) *
5692 (0.25 * (-2. + 1. *
delta) *
5693 (-8. *
C4 * pow(m_rho, 4) +
5695 (-2. - 3. *
delta) * s +
5696 pow(m_rho, 2) * (2. + 1. *
delta + 16. *
C4 * s)) *
5700 (32 * pow(
C4, 2) * pow(m_rho, 8) + 2 * pow(
delta, 2) * pow(s, 2) +
5701 8 *
C4 * pow(m_rho, 6) * (-6 +
delta - 8 *
C4 * s) +
5702 2 *
delta * pow(m_rho, 2) * s * (-6 +
delta - 8 *
C4 * s) +
5703 pow(m_rho, 4) * (12 - pow(
delta, 2) + 8 *
C4 * (6 +
delta) * s +
5704 32 * pow(
C4, 2) * pow(s, 2))) *
5709 (0.75 * pow(m_rho, 4) - 0.125 *
delta * pow(m_rho, 4) -
5710 1. *
C4 * pow(m_rho, 6) +
5712 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5714 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) -
5715 0.25 * pow(m_rho, 2) * s - 0.375 *
delta * pow(m_rho, 2) * s +
5716 2. *
C4 * pow(m_rho, 4) * s + 0.25 *
delta * pow(s, 2) -
5717 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5718 eta1 * (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
5720 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) +
5722 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5723 0.25 *
delta * pow(s, 2) +
5724 1. *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5727 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
5731 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) +
5733 (2. * pow(
pion_mass, 2) + pow(m_rho, 2) - 1. * s) -
5734 3. * pow(m_rho, 2) * s + pow(s, 2)) +
5737 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
5739 pow(m_rho, 2) - 1. * s) +
5740 pow(m_rho, 2) * s + pow(s, 2)) +
5744 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + 2. * pow(m_rho, 4) -
5747 2. * pow(m_rho, 2) + 2. * s))) *
5749 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5752 2. * pow(m_rho, 4) +
5754 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5759 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5760 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5762 4. * pow(m_rho, 2) + 4. * s)) +
5766 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5767 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5769 4. * pow(m_rho, 2) + 4. * s))) *
5771 (0.125 * (-2. + 1. *
delta) *
5772 (2. + 1. *
delta - 8. *
C4 * pow(m_rho, 2)) * pow(tmax, 2)) /
5774 0.5 * pow(1. *
eta1 - 1. *
eta2, 2) *
5775 (-0.5 + 1. *
C4 * pow(m_rho, 2)) * pow(tmax, 2) -
5777 (0.5 - 0.125 * pow(
delta, 2) - 2. *
C4 * pow(m_rho, 2) +
5778 1. *
C4 *
delta * pow(m_rho, 2)) *
5781 0.0625 * pow(1. *
eta1 - 1. *
eta2, 4) *
5782 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 0.5 * s) *
5784 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5786 1. * pow(m_rho, 2) - 1. * s) +
5788 1. * pow(m_rho, 2) + s)) *
5790 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) * pow(tmax, 3) -
5791 0.020833333333333332 * pow(1. *
eta1 - 1. *
eta2, 4) * pow(tmax, 3) +
5792 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5796 2. * pow(m_rho, 4) +
5798 4. * pow(m_rho, 2) * s - 4. * pow(s, 2)) +
5802 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5803 4. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5805 4. * pow(m_rho, 2) + 4. * s)) +
5809 2. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
5810 2. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
5812 4. * pow(m_rho, 2) + 4. * s))) *
5813 (-1. * pow(m_rho, 2) + s + tmax) -
5814 0.03125 * pow(
eta1 - 1. *
eta2, 3) *
5815 (
eta2 * (-1. * pow(
a1_mass, 2) - 1. * pow(m_rho, 2) - 1. * s) +
5816 eta1 * (pow(
a1_mass, 2) - 1. * pow(m_rho, 2) + s)) *
5817 pow(-1. * pow(m_rho, 2) + s + tmax, 2) +
5818 0.010416666666666666 * pow(
eta1 - 1. *
eta2, 4) *
5819 pow(-1. * pow(m_rho, 2) + s + tmax, 3) +
5821 (-1. + 2. *
C4 * pow(m_rho, 2)) *
5827 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5828 0.5 *
C4 * pow(m_rho, 6) +
5830 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5832 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5833 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5834 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5835 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5836 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5838 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5840 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5841 0.125 *
delta * pow(s, 2) +
5842 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5844 1. * s + 1. * tmax)) /
5848 (0.375 * pow(m_rho, 4) - 0.0625 *
delta * pow(m_rho, 4) -
5849 0.5 *
C4 * pow(m_rho, 6) +
5851 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5853 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5854 0.125 * pow(m_rho, 2) * s - 0.1875 *
delta * pow(m_rho, 2) * s +
5855 1. *
C4 * pow(m_rho, 4) * s + 0.125 *
delta * pow(s, 2) -
5856 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2)) +
5857 eta1 * (0.25 * pow(m_rho, 4) - 0.5 *
C4 * pow(m_rho, 6) +
5859 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
5861 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) -
5862 0.125 *
delta * pow(s, 2) +
5863 0.5 *
C4 * pow(m_rho, 2) * pow(s, 2))) *
5872 (-1. * pow(
a1_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5874 (2. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
5875 pow(
pion_mass, 4) * (-1.5 * pow(m_rho, 2) + 1. * s) +
5877 (-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 2. * s)) +
5884 (1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5885 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) +
5886 1. * pow(m_rho, 2) * s))) *
5895 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5897 (2. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) - 1.5 * s) +
5899 3. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5900 0.5 * pow(m_rho, 4) * s - 2. * pow(
pion_mass, 2) * pow(s, 2) -
5901 1. * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
5905 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
5906 1. * pow(m_rho, 2) * s +
5907 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 1. * s))) +
5910 4. * pow(
pion_mass, 4) * pow(m_rho, 2) +
5911 1. * pow(
pion_mass, 2) * pow(m_rho, 4) +
5916 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
5917 0.5 * pow(m_rho, 4) * s + 3. * pow(
pion_mass, 2) * pow(s, 2) +
5918 1. * pow(m_rho, 2) * pow(s, 2) - 0.5 * pow(s, 3) +
5921 (1. * pow(m_rho, 2) - 1. * s) * s +
5922 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 3. * s)))) *
5928 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
5929 2. * pow(m_rho, 2) * s + pow(s, 2) +
5931 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) +
5932 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
5938 (-4. * pow(
pion_mass, 2) + 2. * pow(m_rho, 2) + 2. * s) +
5940 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5941 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
5942 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5945 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5946 pow(
pion_mass, 4) * (6. * pow(m_rho, 2) + 2. * s)) +
5949 1. * pow(m_rho, 4) +
5951 6. * pow(m_rho, 2) - 6. * s) +
5952 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5953 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s))) +
5957 2. * pow(
pion_mass, 4) * pow(m_rho, 4) +
5960 (8. * pow(
pion_mass, 4) - 8. * pow(m_rho, 4) -
5961 4. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s) +
5963 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
5964 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
5968 2. * pow(m_rho, 4) - 8. * pow(
pion_mass, 2) * s -
5969 4. * pow(m_rho, 2) * s + 4. * pow(s, 2) +
5974 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
5976 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
5977 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
5978 4. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
5981 2. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
5982 pow(m_rho, 2) * (2. * pow(m_rho, 2) - 2. * s) * s +
5983 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s)) +
5986 1. * pow(m_rho, 4) +
5988 6. * pow(m_rho, 2) - 6. * s) +
5989 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
5990 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s)) +
5993 2. * pow(
pion_mass, 4) * pow(m_rho, 2) -
5994 2. * pow(m_rho, 6) + 2. * pow(m_rho, 4) * s +
5996 (pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)))) *
6005 28. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6006 20. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6007 10. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) -
6009 16. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6010 12. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6011 4. * pow(m_rho, 6) * s - 4. * pow(
pion_mass, 4) * pow(s, 2) -
6012 6. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6014 4. * pow(m_rho, 2) * pow(s, 3) - 2. * pow(s, 4) +
6016 (-56. * pow(
pion_mass, 2) - 28. * pow(m_rho, 2) + 28. * s) +
6018 (84. * pow(
pion_mass, 4) + 24. * pow(m_rho, 4) +
6019 pow(
pion_mass, 2) * (84. * pow(m_rho, 2) - 72. * s) -
6020 36. * pow(m_rho, 2) * s + 12. * pow(s, 2)) +
6024 4. * pow(m_rho, 2) - 4. * s) +
6025 (4. * pow(m_rho, 2) - 4. * s) * s +
6026 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 8. * s)) +
6028 (-56. * pow(
pion_mass, 6) - 10. * pow(m_rho, 6) +
6029 18. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) -
6031 pow(
pion_mass, 4) * (-84. * pow(m_rho, 2) + 60. * s) +
6033 (-48. * pow(m_rho, 4) + 60. * pow(m_rho, 2) * s -
6034 12. * pow(s, 2)))) +
6038 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6039 7. * pow(
pion_mass, 4) * pow(m_rho, 4) -
6040 2. * pow(
pion_mass, 2) * pow(m_rho, 6) +
6042 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6044 11. * pow(
pion_mass, 4) * pow(m_rho, 2) * s +
6045 6. * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6046 1. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) -
6047 1. * pow(m_rho, 4) * pow(s, 2) -
6049 1. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6052 1. * pow(m_rho, 4) +
6053 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6054 1. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6056 2. * pow(m_rho, 2) + 2. * s)) +
6058 (-42. * pow(
pion_mass, 4) - 9. * pow(m_rho, 4) +
6059 21. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6060 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6062 (28. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6063 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) -
6064 9. * pow(m_rho, 4) * s + 6. * pow(m_rho, 2) * pow(s, 2) +
6067 (18. * pow(m_rho, 4) - 36. * pow(m_rho, 2) * s +
6072 14. * pow(
pion_mass, 6) * pow(m_rho, 2) -
6073 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6074 6. * pow(
pion_mass, 2) * pow(m_rho, 6) + 2. * pow(m_rho, 8) +
6076 (28. * pow(
pion_mass, 2) + 14. * pow(m_rho, 2) - 14. * s) +
6078 1. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6079 16. * pow(
pion_mass, 2) * pow(m_rho, 4) * s -
6080 7. * pow(m_rho, 6) * s + 2. * pow(
pion_mass, 4) * pow(s, 2) +
6081 14. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6082 9. * pow(m_rho, 4) * pow(s, 2) -
6084 5. * pow(m_rho, 2) * pow(s, 3) + 1. * pow(s, 4) +
6087 3. * pow(m_rho, 4) +
6088 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 4. * s) -
6089 5. * pow(m_rho, 2) * s + 2. * pow(s, 2) +
6091 2. * pow(m_rho, 2) + 2. * s)) +
6093 (-42. * pow(
pion_mass, 4) - 3. * pow(m_rho, 4) +
6094 9. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6095 pow(
pion_mass, 2) * (-42. * pow(m_rho, 2) + 36. * s)) +
6097 (28. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6098 pow(
pion_mass, 4) * (42. * pow(m_rho, 2) - 30. * s) +
6099 9. * pow(m_rho, 4) * s - 6. * pow(m_rho, 2) * pow(s, 2) +
6102 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6103 6. * pow(s, 2))))) *
6109 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6110 2. * pow(m_rho, 2) * s + pow(s, 2) +
6112 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)) +
6113 0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6118 (4. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) -
6119 2. * pow(
pion_mass, 2) * s + 2. * pow(m_rho, 2) * s) +
6121 (-12. * pow(
pion_mass, 4) + 2. * pow(m_rho, 4) +
6122 8. * pow(
pion_mass, 2) * s + 4. * pow(m_rho, 2) * s -
6127 pow(m_rho, 2) * (pow(m_rho, 2) - 1. * s) * s +
6128 pow(
pion_mass, 4) * (-3. * pow(m_rho, 2) + s) +
6130 (-6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) + 3. * s) +
6132 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6133 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6134 4. * pow(m_rho, 2) * s + 2. * pow(s, 2))) +
6137 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s +
6138 pow(
pion_mass, 4) * (3. * pow(m_rho, 2) + s) +
6140 (-6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) + 3. * s) +
6142 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6143 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) - 4. * s) -
6144 2. * pow(m_rho, 2) * s + 2. * pow(s, 2)))) *
6145 log(abs(-1. * pow(
a1_mass, 2) + tmax)) -
6149 0.5 * pow(
pion_mass, 4) * pow(m_rho, 2) +
6151 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6153 (0.5 * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 1. * s)) +
6156 (1. * pow(
pion_mass, 4) + 1. * pow(m_rho, 4) +
6157 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6158 0.5 * pow(m_rho, 2) * s) +
6160 (-2. * pow(
pion_mass, 4) - 0.5 * pow(m_rho, 2) * s +
6161 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6162 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6167 (0.5 * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2)) +
6168 pow(
pion_mass, 4) * (0.5 * pow(m_rho, 2) - 1. * s) +
6169 pow(
pion_mass, 2) * (-0.5 * pow(m_rho, 2) + 0.5 * s) * s +
6172 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6173 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s)) +
6176 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 2) - 2. * s) * s +
6177 (-0.5 * pow(m_rho, 2) + 0.5 * s) * pow(s, 2) +
6178 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1.5 * s) +
6180 (-2. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 2) * s +
6181 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6182 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6184 1. * pow(m_rho, 2) + 1. * s) -
6185 (0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6191 (8. * pow(
pion_mass, 2) + 4. * pow(m_rho, 2) - 4. * s) +
6193 (-12. * pow(
pion_mass, 4) - 4. * pow(m_rho, 4) +
6194 4. * pow(m_rho, 2) * s +
6195 pow(
pion_mass, 2) * (-12. * pow(m_rho, 2) + 8. * s)) +
6198 4. * pow(
pion_mass, 4) * pow(m_rho, 2) -
6199 2. * pow(m_rho, 6) + 4. * pow(m_rho, 4) * s -
6200 2. * pow(m_rho, 2) * pow(s, 2) +
6202 (-8. * pow(m_rho, 4) + 8. * pow(m_rho, 2) * s)) +
6204 (8. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) +
6205 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 4. * s) -
6206 2. * pow(m_rho, 4) * s - 2. * pow(m_rho, 2) * pow(s, 2) +
6209 (8. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s -
6214 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6216 2. * pow(
pion_mass, 2) * pow(m_rho, 2) +
6217 pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6219 (6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6220 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) +
6221 pow(m_rho, 2) * s) +
6223 (-4. * pow(
pion_mass, 6) + 2. * pow(m_rho, 6) -
6224 5. * pow(m_rho, 4) * s + 4. * pow(m_rho, 2) * pow(s, 2) -
6226 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6228 (2. * pow(m_rho, 4) - 4. * pow(m_rho, 2) * s +
6232 2. * pow(
pion_mass, 6) * pow(m_rho, 2) +
6234 (-2. * pow(m_rho, 2) + 2. * s) +
6236 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s) +
6238 (3. * pow(m_rho, 4) - 5. * pow(m_rho, 2) * s) +
6240 (6. * pow(
pion_mass, 4) + pow(m_rho, 4) +
6241 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) - 4. * s) -
6242 3. * pow(m_rho, 2) * s) +
6244 (-4. * pow(
pion_mass, 6) + pow(m_rho, 4) * s -
6246 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 2. * s) +
6248 (-2. * pow(m_rho, 4) + 4. * pow(m_rho, 2) * s +
6249 2. * pow(s, 2))))) *
6250 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6253 0.25 * pow(m_rho, 4) - 1. * pow(
pion_mass, 2) * s -
6254 0.5 * pow(m_rho, 2) * s + 0.25 * pow(s, 2) +
6256 (-2. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + 1. * s)) -
6260 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6262 (pow(
pion_mass, 2) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6263 (-0.25 + 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6265 (-1. *
C4 * pow(m_rho, 6) +
6267 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6268 0.25 *
delta * pow(s, 2) +
6269 pow(m_rho, 2) * s * (-0.25 - 0.375 *
delta - 1. *
C4 * s) +
6270 pow(m_rho, 4) * (0.75 - 0.125 *
delta + 2. *
C4 * s))) +
6273 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6275 (0.5 * pow(m_rho, 4) - 1. *
C4 * pow(m_rho, 6) +
6277 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6278 0.25 *
delta * pow(s, 2) +
6279 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6281 (pow(
pion_mass, 4) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6282 s * ((0.25 - 0.125 *
delta) * pow(m_rho, 2) +
6283 (-0.25 + 0.125 *
delta) * s) +
6285 (2. *
C4 * pow(m_rho, 4) + (0.5 + 0.25 *
delta) * s +
6286 pow(m_rho, 2) * (-1. - 2. *
C4 * s))))) *
6287 log(abs(-1. * pow(
a1_mass, 2) + tmax))) /
6290 log(abs(-1. * pow(
pion_mass, 2) + tmax)) +
6293 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6295 5.999999999999999 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 2) *
6297 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 2) * pow(s, 2)) *
6298 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6299 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6301 (0. +
eta2 * pow(
pion_mass, 2) * (4. * pow(m_rho, 2) - 4. * s) +
6302 eta1 * (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s +
6303 pow(
pion_mass, 2) * (-4. * pow(m_rho, 2) + 4. * s))) *
6304 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6306 (2. * (-2. + 1. *
delta) *
6307 (0. + (-0.25 + 0.125 *
delta) * pow(m_rho, 2) * s +
6309 pow(m_rho, 2) * (1. + 2. *
C4 * s))) *
6310 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6313 (
eta1 * (pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 1. * s) +
6315 (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6316 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6318 (-1. * pow(m_rho, 4) + 2.5 * pow(m_rho, 2) * s -
6320 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6323 (0.5 * pow(m_rho, 6) +
6324 pow(
pion_mass, 4) * (1. * pow(m_rho, 2) - 1. * s) -
6325 1.5 * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) -
6327 pow(
pion_mass, 2) * (1.5 * pow(m_rho, 4) -
6328 3. * pow(m_rho, 2) * s + 1.5 * pow(s, 2)) +
6330 (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6332 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s)))) *
6333 log(abs(-1. * pow(
pion_mass, 2) + tmax))) /
6336 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6337 2. * pow(m_rho, 2) * s + pow(s, 2) +
6339 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6341 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax)) +
6345 (-1. * pow(m_rho, 2) + 1. * s) +
6347 (-0.5 * pow(m_rho, 4) +
6348 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6349 0.5 * pow(m_rho, 2) * s) +
6351 (0.5 * pow(m_rho, 4) - 0.5 * pow(m_rho, 2) * s +
6352 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6353 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6357 (
eta1 * (pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 1. * s) +
6358 (-0.5 * pow(m_rho, 2) + 0.5 * s) * s) +
6359 eta2 * (-0.5 * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s -
6361 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1. * s))) *
6362 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6364 1. * pow(m_rho, 2) + 1. * s) -
6367 8.000000000000002 * pow(2. - 1. *
delta, 2) * pow(
pion_mass, 4) *
6369 1. * pow(2. - 1. *
delta, 2) * pow(m_rho, 4) * s +
6371 (
C4 * (32. - 16. *
delta) * pow(m_rho, 6) +
6374 (-8. + 24. *
delta - 10. * pow(
delta, 2) + 32. *
C4 * s -
6376 pow(m_rho, 4) * (-16. + 8. *
delta - 64. *
C4 * s +
6378 log(abs(-1. * pow(
pion_mass, 2) - 1. * pow(m_rho, 2) + s + tmax))) /
6379 (pow(m_rho, 4) - 1. * pow(m_rho, 2) * s) +
6380 0.03125 * pow(
eta1 - 1. *
eta2, 2) *
6388 (-4. * pow(
pion_mass, 4) + 4. * pow(m_rho, 4) +
6389 2. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s) +
6391 (12. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) -
6392 8. * pow(
pion_mass, 2) * s - 4. * pow(m_rho, 2) * s +
6396 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6398 (6. * pow(
pion_mass, 2) + 3. * pow(m_rho, 2) - 3. * s) -
6400 1. * pow(
pion_mass, 2) * pow(m_rho, 2) * s -
6401 1. * pow(m_rho, 4) * s + pow(m_rho, 2) * pow(s, 2) +
6404 1. * pow(m_rho, 2) + s) +
6406 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6407 4. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6408 pow(
pion_mass, 2) * (-6. * pow(m_rho, 2) + 4. * s))) +
6411 3. * pow(
pion_mass, 4) * pow(m_rho, 2) +
6413 (6. * pow(
pion_mass, 2) - 3. * pow(m_rho, 2) - 3. * s) -
6418 pow(m_rho, 2) + s) +
6420 (-6. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6421 2. * pow(m_rho, 2) * s - 2. * pow(s, 2) +
6422 pow(
pion_mass, 2) * (6. * pow(m_rho, 2) + 4. * s)))) *
6426 2. * pow(
a1_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) -
6427 4. * pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + s + tmax) +
6428 pow(-1. * pow(m_rho, 2) + s + tmax, 2))) -
6431 (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6433 (-0.5 * pow(m_rho, 2) + 1. *
C4 * pow(m_rho, 4)) +
6435 (pow(
pion_mass, 2) * (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6436 (0.25 - 0.125 *
delta) * (pow(m_rho, 2) + s)) +
6438 (1. *
C4 * pow(m_rho, 6) +
6440 (1. * pow(m_rho, 2) - 2. *
C4 * pow(m_rho, 4)) -
6441 0.25 *
delta * pow(s, 2) +
6442 pow(m_rho, 4) * (-0.75 + 0.125 *
delta - 2. *
C4 * s) +
6444 (0.25 + 0.375 *
delta + 1. *
C4 * s))) +
6447 (-0.5 + 1. *
C4 * pow(m_rho, 2)) +
6449 (0.5 * pow(m_rho, 2) - 1. *
C4 * pow(m_rho, 4)) +
6451 (-0.5 * pow(m_rho, 4) + 1. *
C4 * pow(m_rho, 6) +
6453 (-1. * pow(m_rho, 2) + 2. *
C4 * pow(m_rho, 4)) +
6454 0.25 *
delta * pow(s, 2) -
6455 1. *
C4 * pow(m_rho, 2) * pow(s, 2)) +
6457 (pow(
pion_mass, 4) * (0.5 - 1. *
C4 * pow(m_rho, 2)) +
6458 s * ((-0.25 + 0.125 *
delta) * pow(m_rho, 2) +
6459 (0.25 - 0.125 *
delta) * s) +
6461 (-2. *
C4 * pow(m_rho, 4) + (-0.5 - 0.25 *
delta) * s +
6462 pow(m_rho, 2) * (1. + 2. *
C4 * s))))) *
6465 1. * pow(m_rho, 2) + s + tmax,
6473 0.5 * pow(
pion_mass, 6) * pow(m_rho, 2) +
6476 (1. * pow(
pion_mass, 2) + 1.5 * pow(m_rho, 2) -
6479 0.5 * pow(m_rho, 2) - 1. * s)) +
6481 (1. * pow(
pion_mass, 2) + 0.5 * pow(m_rho, 2) - 1. * s) +
6483 (-1. * pow(
pion_mass, 2) - 0.5 * pow(m_rho, 2) + 1. * s)) +
6486 1. * pow(
pion_mass, 2) * pow(m_rho, 4) -
6490 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 0.5 * s) -
6491 0.5 * pow(m_rho, 2) * s) +
6493 (-1. * pow(
pion_mass, 4) - 1. * pow(m_rho, 4) +
6494 pow(
pion_mass, 2) * (1. * pow(m_rho, 2) - 0.5 * s) +
6495 0.5 * pow(m_rho, 2) * s) +
6499 0.5 * pow(m_rho, 2) * s +
6500 pow(
pion_mass, 2) * (-1. * pow(m_rho, 2) + 1.5 * s)))) *
6503 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6504 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6505 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6506 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6508 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6514 1. * pow(
pion_mass, 4) * pow(m_rho, 4) +
6517 4. * pow(
pion_mass, 4) * pow(m_rho, 2) * s -
6518 1.5 * pow(
pion_mass, 2) * pow(m_rho, 4) * s +
6520 4. * pow(
pion_mass, 2) * pow(m_rho, 2) * pow(s, 2) +
6521 0.5 * pow(m_rho, 4) * pow(s, 2) -
6523 1. * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6527 0.5 * pow(
pion_mass, 2) * s + 0.5 * pow(s, 2)) +
6530 (1. * pow(m_rho, 2) - 0.5 * s) * s +
6531 pow(
pion_mass, 2) * (-2. * pow(m_rho, 2) + 2.5 * s)) +
6534 pow(
pion_mass, 4) * (2. * pow(m_rho, 2) - 1.5 * s) -
6535 0.5 * pow(m_rho, 4) * s + 0.5 * pow(s, 3) +
6537 (1. * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s -
6543 (1. * pow(
pion_mass, 2) + 1. * pow(m_rho, 2) - 0.5 * s) +
6545 pow(
a1_mass, 4) * (-0.5 * pow(m_rho, 4) +
6546 (-0.5 * pow(
pion_mass, 2) + 0.5 * s) * s) +
6547 pow(
pion_mass, 4) * (-0.5 * pow(m_rho, 4) +
6548 2. * pow(m_rho, 2) * s - 1.5 * pow(s, 2)) +
6550 (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6553 (1. * pow(
pion_mass, 4) + 0.5 * pow(m_rho, 4) +
6554 pow(
pion_mass, 2) * (2. * pow(m_rho, 2) - 1.5 * s) -
6555 1. * pow(m_rho, 2) * s + 0.5 * pow(s, 2) +
6557 1. * pow(m_rho, 2) + 1.5 * s)) +
6560 pow(
pion_mass, 4) * (-1. * pow(m_rho, 2) + 0.5 * s) +
6562 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6564 s * (0.5 * pow(m_rho, 4) - 1. * pow(m_rho, 2) * s +
6565 0.5 * pow(s, 2))))) *
6568 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6569 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6570 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6571 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6573 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6577 pow(m_rho, 4) - 2. * pow(
pion_mass, 2) * s -
6578 2. * pow(m_rho, 2) * s + pow(s, 2) +
6580 (-2. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s)) -
6581 (0.0625 * pow(
eta1 - 1. *
eta2, 2) *
6585 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6588 0.5 * pow(m_rho, 2) + 0.5 * s) +
6590 (10. * pow(
pion_mass, 6) - 2.5 * pow(m_rho, 6) +
6591 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) +
6592 6. * pow(m_rho, 4) * s - 4.5 * pow(m_rho, 2) * pow(s, 2) +
6596 (1. * pow(m_rho, 2) - 1. * s) * s +
6597 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6599 (1. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6600 pow(
pion_mass, 4) * (2.5 * pow(m_rho, 2) - 0.5 * s) -
6601 1. * pow(m_rho, 4) * s + 0.5 * pow(m_rho, 2) * pow(s, 2) +
6603 (2. * pow(m_rho, 4) - 2. * pow(m_rho, 2) * s)) +
6606 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6607 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6608 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6610 6. * pow(m_rho, 2) + 6. * s) +
6612 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6615 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6616 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6617 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6619 (-5. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) -
6620 3.5 * pow(m_rho, 6) * s + 4.5 * pow(m_rho, 4) * pow(s, 2) -
6621 2.5 * pow(m_rho, 2) * pow(s, 3) + 0.5 * pow(s, 4) +
6622 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6624 (-2. * pow(m_rho, 4) + 1. * pow(m_rho, 2) * s +
6627 (3. * pow(m_rho, 6) - 8. * pow(m_rho, 4) * s +
6628 7. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)))) +
6632 (5. * pow(
pion_mass, 2) + 2.5 * pow(m_rho, 2) - 2.5 * s) +
6635 0.5 * pow(m_rho, 2) + 0.5 * s) +
6637 (-10. * pow(
pion_mass, 4) - 2. * pow(m_rho, 4) +
6638 5. * pow(m_rho, 2) * s - 1. * pow(s, 2) +
6639 pow(
pion_mass, 2) * (-10. * pow(m_rho, 2) + 8. * s)) +
6641 (10. * pow(
pion_mass, 6) + 0.5 * pow(m_rho, 6) +
6642 pow(
pion_mass, 4) * (15. * pow(m_rho, 2) - 9. * s) -
6643 3. * pow(m_rho, 4) * s + 1.5 * pow(m_rho, 2) * pow(s, 2) +
6646 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s)) +
6649 0.5 * pow(m_rho, 6) + 1.5 * pow(m_rho, 4) * s -
6650 1.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3) +
6651 pow(
pion_mass, 4) * (-6. * pow(m_rho, 2) + 6. * s) +
6653 6. * pow(m_rho, 2) + 6. * s) +
6655 (-3. * pow(m_rho, 4) + 6. * pow(m_rho, 2) * s -
6658 (12. * pow(
pion_mass, 4) + 3. * pow(m_rho, 4) +
6659 pow(
pion_mass, 2) * (12. * pow(m_rho, 2) - 12. * s) -
6660 6. * pow(m_rho, 2) * s + 3. * pow(s, 2))) +
6663 pow(
pion_mass, 6) * (2.5 * pow(m_rho, 2) - 0.5 * s) +
6665 (4. * pow(m_rho, 4) - 6. * pow(m_rho, 2) * s) +
6667 (-1. * pow(m_rho, 4) + 2. * pow(m_rho, 2) * s -
6670 (1.5 * pow(m_rho, 6) - 6. * pow(m_rho, 4) * s +
6671 4.5 * pow(m_rho, 2) * pow(s, 2))) +
6674 pow(
pion_mass, 6) * (-10. * pow(m_rho, 2) + 4. * s) +
6676 (-8. * pow(m_rho, 4) + 13. * pow(m_rho, 2) * s +
6679 (-1. * pow(m_rho, 6) + 6. * pow(m_rho, 4) * s -
6680 3. * pow(m_rho, 2) * pow(s, 2) - 2. * pow(s, 3)) +
6681 s * (0.5 * pow(m_rho, 6) - 0.5 * pow(m_rho, 4) * s -
6682 0.5 * pow(m_rho, 2) * pow(s, 2) + 0.5 * pow(s, 3)))) +
6687 1. * pow(m_rho, 2) - 1. * s) +
6689 (-10. * pow(
pion_mass, 2) - 5. * pow(m_rho, 2) + 5. * s) +
6691 (20. * pow(
pion_mass, 4) + 6. * pow(m_rho, 4) +
6692 pow(
pion_mass, 2) * (20. * pow(m_rho, 2) - 16. * s) -
6693 8. * pow(m_rho, 2) * s + 2. * pow(s, 2)) +
6695 (-20. * pow(
pion_mass, 6) - 4. * pow(m_rho, 6) +
6696 6. * pow(m_rho, 4) * s - 2. * pow(s, 3) +
6697 pow(
pion_mass, 4) * (-30. * pow(m_rho, 2) + 18. * s) +
6699 (-18. * pow(m_rho, 4) + 18. * pow(m_rho, 2) * s)) +
6701 (-2. * pow(
pion_mass, 8) - 1. * pow(m_rho, 8) +
6702 3. * pow(m_rho, 6) * s - 3. * pow(m_rho, 4) * pow(s, 2) +
6703 1. * pow(m_rho, 2) * pow(s, 3) +
6704 pow(
pion_mass, 6) * (-5. * pow(m_rho, 2) + 1. * s) +
6706 (-10. * pow(m_rho, 4) + 10. * pow(m_rho, 2) * s) +
6708 (-6. * pow(m_rho, 6) + 12. * pow(m_rho, 4) * s -
6709 6. * pow(m_rho, 2) * pow(s, 2))) +
6711 (10. * pow(
pion_mass, 8) + 1. * pow(m_rho, 8) +
6712 pow(
pion_mass, 6) * (20. * pow(m_rho, 2) - 8. * s) -
6713 2. * pow(m_rho, 6) * s + 2. * pow(m_rho, 2) * pow(s, 3) -
6716 (22. * pow(m_rho, 4) - 20. * pow(m_rho, 2) * s -
6719 (8. * pow(m_rho, 6) - 12. * pow(m_rho, 4) * s +
6723 1. * pow(m_rho, 6) +
6724 pow(
pion_mass, 4) * (12. * pow(m_rho, 2) - 12. * s) +
6726 12. * pow(m_rho, 2) - 12. * s) -
6727 3. * pow(m_rho, 4) * s + 3. * pow(m_rho, 2) * pow(s, 2) -
6730 (6. * pow(m_rho, 4) - 12. * pow(m_rho, 2) * s +
6733 (-24. * pow(
pion_mass, 4) - 6. * pow(m_rho, 4) +
6734 12. * pow(m_rho, 2) * s - 6. * pow(s, 2) +
6736 (-24. * pow(m_rho, 2) + 24. * s))))) *
6739 4. * pow(
pion_mass, 2) * pow(m_rho, 2) + pow(m_rho, 4) -
6740 4. * pow(
pion_mass, 2) * s - 2. * pow(m_rho, 2) * s +
6741 pow(s, 2) - 4. * pow(
pion_mass, 2) * tmax -
6742 2. * pow(m_rho, 2) * tmax + 2. * s * tmax + pow(tmax, 2) +
6744 (-4. * pow(
pion_mass, 2) - 2. * pow(m_rho, 2) + 2. * s +
6748 pow(m_rho, 4) - 4. * pow(
pion_mass, 2) * s -
6749 2. * pow(m_rho, 2) * s + pow(s, 2) +
6751 (-8. * pow(
pion_mass, 2) - 4. * pow(m_rho, 2) + 4. * s)))) /
constexpr static double C4
Constant in the computation of photon cross sections.
constexpr static double Const
Constant in the computation of photon cross sections.
constexpr static double delta
Constant in the computation of photon cross sections.
constexpr static double Gammaa1
Constant in the computation of photon cross sections.
constexpr static double eta2
Constant in the computation of photon cross sections.
constexpr static double ghat
Constant in the computation of photon cross sections.
constexpr static double Pi
Value of .
constexpr static double a1_mass
Value of a1 mass in GeV.
constexpr static double eta1
Constant in the computation of photon cross sections.
constexpr double gev2_mb
GeV^-2 <-> mb conversion factor.
std::array< T, 2 > get_t_range(const T sqrts, const T m1, const T m2, const T m3, const T m4)
Get the range of Mandelstam-t values allowed in a particular 2->2 process, see PDG 2014 booklet,...
constexpr double pion_mass
Pion mass in GeV.
double cut_off(const double sigma_mb)
Cross section after cut off.