コンプトン散乱計算の数値代入について

以前、ここでご教示頂きましたコンプトン散乱の計算式に適当に数値を入れて計算しました。
元の式は、
Tr{(sl[q]+m)γμu(sl[p]+sl[k]+m)γνu(sl[p]+m)γνd( sl[p]+sl[k]+m)γμd}
= Tr[(-2sl[q]+4m)( sl[p]+sl[k]+m)(-2sl[p]+4m)( sl[p]+sl[k]+m)]
です。
これを、mathematicaで計算して

y4 = 64*k0^2*m^2 - 64*k1^2*m^2 - 64*k2^2*m^2 - 64*k3^2*m^2 + 64*m^4 + 64*k0*m^2*p0 - 64*k1*m^2*p1 - 64*k2*m^2*p2 - 64*k3*m^2*p3 -
64*k0*m^2*q0 + 16*k0^2*p0*q0 + 16*k1^2*p0*q0 + 16*k2^2*p0*q0 + 16*k3^2*p0*q0 - 48*m^2*p0*q0 + 32*k0*p0^2*q0 + 16*p0^3*q0 -
32*k0*k1*p1*q0 - 32*k0*p1^2*q0 - 16*p0*p1^2*q0 - 32*k0*k2*p2*q0 - 32*k0*p2^2*q0 - 16*p0*p2^2*q0 - 32*k0*k3*p3*q0 ・・・・・
以下省略


x = 2; y = 1; q0 = x; q1 = y; q2 = y; q3 = y; p0 = x; p1 = y;
p2 = y; p3 = y;k0 = x; k1 = y; k2 = y; k3 = y;
m = Sqrt[q0^2 - q1^2 - q2^2 - q3^2];
Print[N[y4]];

を得て、適当に数値を入れて計算しました。
もちろん適当な数値なので、何の意味もない値が導かれました。そこで質問ですが、この式に意味のある数値を代入して、実験値に近い計算値を導くには、それぞれの変数にどのような値を入れれば良いのでしょうか？

No.1 回答者： noname#108554 回答日時： 2005/09/30 08:03

コンプトン散乱の全断面積を求めたいなら、パラメータは

入射X線のエネルギー（keVオーダー）
電子質量（～0.5MeV）
微細構造定数（～1/137）
しかないのでは？