Orz日記 by Akio Morita

_ [SAD]Trackingの謎

SADで RFSWとRADをONにしたTrackingを 実行するとTraclkingの中で実現する平衡分布が z方向へずれるという謎現象があるのですが、ヒントを得ました。

どんなヒントかといえば、トラッキングのために KEKBOpticsをロード(自動的にCALCが走る)した後、 事前にファイルに保存してある粒子分布を読み出し TrackParticles[]へ突っ込んだ場合(1)と、 その場で粒子分布を生成して TrackParticles[]へ突っ込んだ場合(2)で、 粒子分布の重心運動が異なるというものです。 (1)で事前に用意した粒子分布自体は、(2)で生成した 分布をファイルに書き出したものなわけで... File I/Oが破滅的にイカレてない限り粒子分布は 同じものになります。念のため、TrackParticles[]直前で 擬似乱数列も初期化しているので、プログラム的には まったく同じ出力を返すべきところですが、z方向に 30mm謎のドリフトをしてくれます。

注目すべきは、Opticsロードから TrackParticles[]開始までの 粒子分布生成の方法しか違わない点で、分布生成に Emittance[]の吐き出す座標変換行列を使っているので、 TrackParticles[]の直前に Emittance[]を差し込んだ所、 謎のドリフトが消失...トラッキングルーチンで Emittance[]が書き出す大域変数を参照している ようです。そんな仕様、聞いてないよ Orz

NORADの場合

NORADにした場合でも、謎のドリフトは発生する。