[thread display] [new arrival display] [word search] [past log] [管理用]
| Subject | ： TrackParticles[]の進行方向の計算 |
| Date | ： 2007/11/30(Fri) 02:38:42 |
| Contributor | ： Hiroyuki Harada < > |
広島大学の原田寛之です。
粒子の初期分布を
a = {{0},{0},{0},{0},{0},{0.01},{1}};
とし、いくつかの計算より、進行方向の変位ｚを求め、その差を比較しました。
しかしながら、TrackParticlesによる結果を理解することができません。
ソースコードの方を眺めてみましたが、理解するには時間がかかりそうです。
そこで、今回、TrackParticlesによる進行方向の計算に関して、お聞きしたいことがあり、
ご連絡いたしました。
計算手法として以下の（１）〜（３）を行いました。
手法（１）：TrackParticlesを用いた計算
TrackParticlesを用いて、初期分布を1回周回させると、運動量のずれの分、同期位相Z=0から
ずれが生じます。
現在計算されるスリッページファクターは、η＜０であるため、zは増加します。
ハーモニックナンバーを１としているため、位相は２πz/Lとなります。Lは周長です。
＜実際の計算＞
b = TrackParticles[{1,a},1];
z = b[5,1];
手法（２）：スリッページファクターηからの計算
ηとΔp/pから進行方向の変位zは、
z = -η Δp/p L
となりますので、以下のように計算しました。
＜実際の計算＞
slip = MomentumCompaction /. Emittance[Emittance->Flase][6];
dp = 0.01;
L = 348.333;
z = -slip * dp * L;
手法（３）：６次元行列による計算
SADから6次元行列を取り出し、行列計算で変位zを求めました。
＜実際の計算＞
m = TransferMatrices /. Last[Emittance[Matrix->True]];
m = Last[m]; !!! 一周の行列の導出
b = m.Take[a,{1,6}];
z = b[5];
上記の（１）、（２）、（３）を比較すると、（２）と（３）は一致しますので、zは進行方向の
位相だと考えております。しかし、（１）つまりTrackParticlesにおけるzが何を意味し、計算を
行っているのか、わかりません。
何か特別なおまじないが抜けているのか、最初からTrackParticlesにおけるzの定義が異なるのか
教えていただきたいです。
お手数ですが、宜しく御願い致します。