Estoy realizando algunas de simulación en el que necesita la integración de una sección transversal de la forma
$$\sigma = \int_{-s}^{0}dt \frac{1}{(t-\mu^2)^2} f(\theta,\phi)$$
donde $s$ $t$ son los habituales de mandelstam variables, $\mu^2$ es una constante, y $f$ es no-Lorenz invariante en función de los antecedentes y de la dispersión de los ángulos. Debido a esto $f$ no puedo impulso a la COM marco donde los movimientos se vuelven triviales, por lo tanto quiero, en lugar de integrar
$$ \sigma = \int d\theta d\phi \frac{1}{(t(\theta,\phi)-\mu^2)^2} f(\theta,\phi) $$
lo que me estoy dando cuenta que no trivial es el de determinar el Jacobiano de esta transformación (ya que estamos pasando de un 1-D 2-D integral de la matriz Jacobiana no es cuadrado).
El gran problema es que $t(\theta,\phi)$ es un degenerado función de $\theta$ $\phi$ y por lo tanto la integración de más de $\theta$ $\phi$ doble cuenta mi sección transversal diferencial y por lo tanto el total de la sección transversal I es demasiado grande.
He comprobado esta vía de integración numérica para $f=1$ y en comparación con el resultado analítico $$\sigma = \int_{-s}^{0}dt \frac{1}{(t-\mu^2)^2} = \frac{s}{\mu^2(\mu^2+s)} $$
He intentado un tipo de método de fuerza bruta para encontrar el jacobiano, mediante la extracción de líneas de constante $t$ a partir de una Mathematica contorno de la parcela, y el cálculo de la arclength de estas líneas y, a continuación, dividir la degeneración.
Larga historia corta, no funciona.
tl;dr
Encontrar el jacobiano que convierte a $$ \int_{-s}^{0} dt \rightarrow \int d\theta d\phi$$ en el laboratorio de marco
Alguna idea?
EDITAR: Así que mi ángulos se definen de la siguiente manera, puedo alinear mi eje z con una de las partículas entrantes, especificando la $\theta$ $\phi$ (i.e el vector unitario) en el que se dispersa de forma exclusiva determina su magnitud y, a continuación, trivialmente (por conservación de momento) el otro vector de dispersión. (ya no estamos en COM marco el otro vector entrante no es anti-alineados)
EDIT 2: En un intento de hacer bien planteado puedo hacer el siguiente "cambio"
$$ \int_{-s}^{0} dt = \int_{0}^{2\pi}\int_{-s}^{0} d\psi dt \frac{1}{2\pi} $$
En la COM marco de esta correspondería a la integración sobre el cono de vectores que corresponden a un valor específico de $t$. En principio debería ser capaz de crisis de una Jacobiana ahora, eso si puedo expresar $\psi(\theta,\phi)$ ( que corresponden a las líneas de constante $t$).....hmmmmm....Me voy a quedar con esta idea para un poco.