PHENIX muon tracker no tracking algorithm no rikai kentaro watanabe
Muon Tracking System
phnx_z South: 〜 St1 ( 1.800m ) 〜 St2 ( 3.000m ) 〜 St1 ( 4.600m ) North: 〜 St1 ( 1.800m ) 〜 St2 ( 3.470m ) 〜 St1 ( 6.125m )
At North Arm, the position of octant and half octant ① ② ③ ④ ⑤ ⑥ ⑦ ⑧
Muon Magnet によってできる磁場の向きは動径 (radial) 方向を向いている ため ( 図 ) 、粒子は 円周 (azimuthal) 方向に曲げられる。よって、この方向に 良い位置分解能を持つように設計されている。
図に示すように、 1 つのギャップを構成する cathode 面は、 1 つは動径方向に張られて いる (non-stereo plane と呼ぶ ) が、もう 1 つは動径方向から少し角度をつけて張ってあ る (stereo plane と呼ぶ ) 。これに よって 2 次元座標が得られるようになっている。更に、 stereo plane の角度は各ギャップによって変化をつけてあり ( 傾きの具体的な値について は表参照 ) 、多重度の高いイベントにおける不定性を減らしている。 100μm の位置分解 能は non-stereo plane に対する要求であり、 stereo plane のストリップの角度は各 ギャップによって異なるため、位置分解能は 300 μm である。
cluster の Mathieson Fit カソードストリップに誘起される電荷は 1 ヒットに対して通常 2~3 ストリップ に広がっている。これを Mathieson Fit([12] 参照 ) することで正確な粒子の通過位 置情報を得ることが出来、この方法で 100 μm の位置分解能が得られることはテ ストベンチにおいて確認されている 3 。このとき重要になるのが、電荷が 誘起 された strip の中で最も多く電荷が誘起された strip ( ピークストリップ ) に対する ノイズレベルである。 これが 1 % 以下でなければ 100 μm の位置分解能を得る ことが出来ない [11] 。誘起される電荷の典型的な 値は 100 fC であるため、ノイ ズは 1 fC 以下である必要があるらしい
hit cluster coord stub Track Track Fit TMutHit TMutClus TMutCoord TMutStub TMutGapCoord TMuTrk TMuTrkRes TMutTrkPar
radial azimuthal radial
Sagitta analysis on the Zero field run motivation
θ sagitta_x sagitta_y sagitta_w sagitta_r track point extarnal point
external point error track_z[0] track_z[1] track_z[2] track_x[0] track_x[1] external_x track_x[2]
Sagitta analysis on the Zero field run tracking parameter study
parameter study Entries
Run9 March Zero field run track entries par half arm March half octant = 1 half octant = 2 South Octant 7-h1 is little extremely!!
track distribution phnx X:Y satation 1
track distribution phnx X:Y satation 2
track distribution phnx X:Y satation 3 South Octant 7-h1 is little extremely!! that reason is here. And octant 3’s out side acceptance is dead. It is important information.
Run9 March Zero field run track entries par half arm March North Octant 2-h1,2 is little extremely!! half octant = 1 half octant = 2
track distribution phnx X:Y satation 1
track distribution phnx X:Y satation 2
track distribution phnx X:Y satation 3
Nouth Oct 0 half 1 track_half_octant==0 track_half_octant==1 track_half_octant==0 track_half_octant==1 Run9 march Run9 may
parameter study Muon ID hit
衝突点で発生した粒子はまず NeosCone( 銅 ) と Central Magnet( 鉄 ) を通る。これらの 吸収材による反応 長は ̃5 程度でこの時点でハドロンの数は 1/100 になる。また、 5 層の MuID には 1 層ごとに鉄が挟まっ ており、 μ 粒子が MuID の 5 層目まで到達す るには 2.5GeV/C の運動量が必要となる。 (← ハドロン吸収材を Run11 から導入した から変わった? ) 2.5GeV/c ← first layer 0 origin
parameter study vertex position
# TMutTrkPar *trk_par_vtx = mut_trk → get_trk_par_vtx( ); float vtx_x; → TMutTrkPar*trk_par_vtx.get_x(); arm 0 (South) arm 1 (North)
float vtx_y; → TMutTrkPar*trk_par_vtx.get_y(); arm 0 (South) arm 1 (Nouth)
float vtx_z; → TMutTrkPar*trk_par_vtx.get_z(); arm 0 (South) arm 1 (North)
parameter study clstr size
ハイ!!犯人は south arm station3 gap1 octant 2,3,5 っす!! Conclusion !! cluster size 16 strange peak from South Arm Station2, octant5,gap0. South Arm Station3, octant2, 3, 5,gap1. parameter study 2 cluster size
parameter study sagitta_s3_w
sagitta distribution south
sagitta distribution North
Sagitta analysis on the Zero field run sagitta error study
多重散乱 ( multiple scattering ) 多重散乱による平均散乱角 物質中で粒子は原子核のクーロン場によって多数回の散乱を受ける。 散乱微分断面積 ( ラザフォード公式 ) ・放射長 荷電粒子が制動放射などをして減速する際、エネルギー が 1/e ( e は自然対数の底)に減少するまでに通過する平 均距離。
collision external point stub point st3 sagitta stub point st2 stub point st1 stub point st3 external point つまり、この間での多重散乱 の効果を考慮すればよい。 この間での多重散乱の効果は sagitta に影響しない。 station3 での sagitta を考慮する際の multiple scattering の効 果
Station2 での放射長は元々のデザイン状は 0.1% であるが、その後ノイズ対策のためチェンバー表面にアルミホイルを 貼っ ている。この追加により現在ではトータルの放射長は 0.2% となっている。従って多重散乱によるト ラック分布の広がりは station3 での sagitta を考慮する際の multiple scattering の効 果 station2 内での散乱 と求まる。 North と South arm では station2-3 間の距離が違うので、この角度の広がりを元にそれぞれ Station-3 へ投影し た場合の分布の広がりをそれぞれ計算しその結果を以下のテーブルに示す 更に今回は station2-3 間での多重散乱を考慮する。 next
空気の放射長は 37g cm−2 である。従って √x/X0 を各 station 毎に計算すると North South
Sagitta analysis on the Zero field run fitting study
従来は sagitta に対して2 gausiam で fit を行っていたが、どうも fit が上手く決まら ないという問題があった。そこで今回は 2gaus の他に gaus+pol2 と 2gaus+pol2 の二つ を新たに試した。 以下にその結果を示す。
2gausiam
gausiam + pol2
2gausiam + pol2
総じて
south octant8 half2 south octant3 half2 south octant1 half2 sample of high improvement south octant8 half2 south octant3 half2 south octant1 half2
what is the composition gaus + gaus + pol2 sample: south octant8 half2
Run9 South March mean position direction ?
Run9 North March mean position direction ?
Sagitta analysis on the Zero field run compere OASys parameter
IKEDA san の Run9 Zero field Run analysis に対して segment の不足分を補った事による各 half octant の number of track の増加量 (March) ※ octant ID number octant 1.0 = octant1 half oct 1 octant 1.5 = octant1 half oct 2 South Run9 March 平均して各 half octant 当たり 2.9 倍、統計量 にして 30,000 event 増 えた。 IDΔincreaceincreace IDΔincreaceincreace North Run9 March 平均して各 half octant 当たり 2.8 倍、統計量 にして 35,000 event 増 えた。
IDΔincreaceincreace South Run9 May 平均して各 half octant 当たり 4.2 倍、統計量 にして 5,500 event 増 えた。 South Run9 May 平均して各 half octant 当たり 4.2 倍、統計量 にして 8,000 event 増 えた。 IDΔincreaceincreace
統計量が増えた事で、どのような変化が生じたか? sample is Run9 March South octant 0 half oct 0 entries Rchis st gaus mean [cm] mean_er ±0.0021[cm] sigma 0.19[cm] entries Rchis st gaus mean [cm] mean_er ± [cm] sigma 0.14[cm]
improvement parameters 当然ながら統計が増えた事で中心値決定精度は向上した、特に North arm に関し ては ±50μm のオーダーに大分が入ったので、これは嬉しい情報である。
increase event from ex-analysis South Run9 March to MayNorth Run9 March to May
armΔnew(μm)Δold(μm)Δnew_er(μm)Δold_er(μm) south
armΔnew(μm)Δold(μm)Δnew_er(μm)Δold_er(μm) north
ここまで segment を全て使う事で統計量が増えてい る事はすでに知っていた。その増加した イベントは特定の octant に偏らない事を確 認した。 これよりは全ての segment を用いたデータ を解析に用いる。