A method for luminosity determination based on real-time hit reconstruction with the LHCb silicon pixel detector

LHCb 実験のアップグレード版において、2024 年の物理ランから実運用されている、検出器読み出しレベルでのリアルタイムヒット再構成に基づく新しいルミノシティ測定手法は、統計的分解能が 1% 未満であり、非線形性が無視できるレベルで、LHCb の Run 3 運転条件を含む広範囲のルミノシティ領域において高い精度と応答性を示すものである。

要約

素粒子の「雨」を数える新しい方法：LHCb 実験の革新的な輝度測定

この論文は、スイスにある巨大な粒子加速器「LHC（大型ハドロン衝突型加速器）」で行われている実験、LHCb について書かれています。特に、「加速器の中でどれくらい頻繁に粒子が衝突しているか（これを『輝度』と呼びます）」を、これまでとは全く違う方法で、リアルタイムかつ高精度に測る技術を紹介しています。

まるで、激しい雨の中を歩いているときに、「どれくらい雨粒が当たっているか」を数える方法を工夫したような話です。

---

1. 背景：なぜ「輝度」を測る必要があるの？

LHCb 実験は、素粒子の不思議な振る舞いを調べるために、プロトン（水素の原子核）同士を激しくぶつけ合っています。
実験を成功させるためには、「今、衝突がどれくらい起きているか」を常に正確に把握し、それを一定に保つ（レベリングする）ことが不可欠です。

• 従来の方法： 以前は、衝突によって発生する「熱」や「光」を別のセンサーで測って推定していました。
• 今回の課題： 実験がさらに高性能化（アップグレード）したため、より速く、より正確に、かつ直接「衝突そのもの」を数える新しい方法が必要になりました。

2. 新技術の核心：「VELO」という超高速カメラ

LHCb には、衝突点のすぐそばに**「VELO（ベロ）」**という特殊な検出器があります。これは、シリコンのピクセル（画素）でできた、非常に解像度の高い「カメラ」のようなものです。

• 従来の壁： 以前は、このカメラが捉えたデータをすべてコンピューターに送り、ソフトウェアで解析してから「衝突があった」と判断していました。しかし、衝突が 1 秒間に 4000 万回（40MHz）も起きていると、ソフトウェアが追いつきません。
• 今回の breakthrough（ブレイクスルー）： 研究チームは、「カメラの読み出しボードそのもの（FPGA というチップ）」に、データ解析機能を組み込んでしまいました。
  • イメージ： 従来の方法は、雨粒が地面に落ちるのをすべて写真に撮り、後でコンピューターで「雨粒の数」を数えていたのに対し、今回の方法は**「カメラ自体が『今、雨粒が 100 個落ちたよ』と瞬時に計算して報告する」**ようなものです。

3. 具体的な仕組み：雨粒の「集まり」を数える

この新しい方法は、2 つのステップで輝度を測ります。

ステップ 1：ピクセルを「クラスター（集まり）」にする

衝突すると、検出器のピクセルに信号が走ります。隣り合ったピクセルが同時に反応したら、それは「1 つの粒子（雨粒）」が通った証拠だと判断します。これを**「クラスター」**と呼びます。
この「クラスター化」を、データが流れてくる瞬間に、ハードウェア（チップ）の中で即座に行います。

ステップ 2：「空っぽの部屋」を数える（ゼロカウント法）

ここが最も面白い部分です。輝度を測るために、検出器の特定のエリア（内側と外側）を「部屋」のように区切ります。

• 平均カウント法： 「この部屋に、1 秒間に何個のクラスターが入ってきたか」を単純に数えます。
• ゼロカウント法（log0 法）： 「この部屋に、全く何も入ってこなかった瞬間（空っぽの部屋）」を数えます。

なぜ「空っぽ」を数えるの？
もし衝突が激しすぎて、部屋が常に満員（クラスターが重なり合っている）だと、正確な数が分かりません（飽和現象）。しかし、「全く何も入ってこなかった瞬間」の割合を統計的に分析すれば、「実際にはどれくらいの雨粒が降っていたか」を逆算して、正確に推測できるという数学的な裏付けがあります。

• アナロジー： 激しい雨の中で、傘を差している人が「雨に濡れていない瞬間」を数えれば、その人がどれだけの雨にさらされていたかを推測できる、といった感じです。

4. 結果：どれくらいすごいのか？

この新技術は、2024 年の実験で実用化され、驚異的な性能を示しました。

• 超高速： 100 ミリ秒（0.1 秒）未満の時間で結果が出ます。まるで、雨の降り方をリアルタイムで追っているようです。
• 高精度： 統計的な誤差は 1% 未満。これは、100 個の雨粒を数えて、1 個も間違えないレベルです。
• 頑丈さ： 衝突が激しすぎてデータが混雑しても（鉛イオン衝突など）、正確に測れます。また、検出器の一部が壊れても、他の「部屋」のデータを使って補正するため、全体として安定しています。
• 柔軟性： 検出器の位置が少しずれても、ソフトウェアで補正して正確な値を出せます。

5. まとめ：なぜこれが重要なのか？

この論文で紹介された方法は、**「ハードウェア（チップ）の力で、データ解析を前倒しして行う」**という画期的なアプローチです。

• 従来： データを全部集めてから「どれくらい衝突したか」を計算していた。
• 今回： データが流れてくる瞬間に「衝突数」を計算し、それを即座に実験の制御にフィードバックしている。

これにより、LHCb 実験は、より安全に、より効率的に、そしてより多くの「素粒子の謎」を解明できるようになりました。まるで、嵐の中で航海する船が、最新のレーダーと自動操縦システムを備えて、波の揺れを瞬時に検知し、進路を微調整できるようになったようなものです。

この技術は、将来の他の実験や、より複雑なデータ処理が必要な分野でも応用が期待される、非常に重要な進歩です。

技術概要

CERN LHCb 実験における VELO ピクセル検出器を用いたリアルタイム輝度測定手法の技術的概要

本論文は、CERN の LHCb 実験（アップグレード版）において、ビームクロスイングレート（40 MHz）で Vertex Locator（VELO）ピクセル検出器の全ヒットをリアルタイムに再構成し、それに基づいて即座に輝度（Luminosity）を決定する新しい手法を提案・実装したものです。以下に、問題提起、手法、主要な貢献、結果、および意義について詳細をまとめます。

1. 背景と課題

LHCb 実験は、重陽子ハドロン崩壊における CP 対称性の破れの精密測定を主目的としていますが、その物理プログラムは量子色力学や重イオン衝突など多岐にわたります。LHCb は他の LHC 実験（ATLAS, CMS）に比べて低い輝度で運転されるように設計されており、これを一定に保つ「輝度レベリング（Luminosity Levelling）」が不可欠です。

• 課題: Run 3（2022 年以降）では、従来の L0 トリガー段階が廃止され、トリガーレス読み出しアーキテクチャへ移行しました。これにより、従来のカロリメータベースのリアルタイム輝度測定が困難になり、より高精度（5% 以上の精度）かつ高速なリアルタイム輝度モニタリングが求められました。
• 既存手法の限界: 従来の手法は、特定の検出器（PLUME など）や HLT1 での高レベル再構成に依存しており、検出器レベルでの低レベル情報（ヒット）を直接利用した高速な手法は不足していました。

2. 提案手法と実装

本研究では、VELO 検出器の読み出しボード（FPGA）上で実行される「オンザフライ・クラスターリング（On-the-fly Clustering）」技術を活用し、再構成されたヒット（クラスター）の数を直接カウントする手法を開発しました。

• ハードウェア実装:
  • VELO の読み出しボード（TELL40）に搭載された Intel Arria-10 FPGA 上で、40 MHz のビームクロスイングレートで全ヒットの 2 次元クラスター再構成を実行しています。
  • これにより、イベントビルディングやソフトウェアトリガー（HLT1/HLT2）を待たずに、検出器レベルで約 $10^{11}$ ヒット/秒のデータを処理可能です。
• カウンター設計:
  • VELO の 26 層それぞれに対して、ビーム軸から約 14 mm（Inner）と 26 mm（Outer）の位置に 4 つずつの蓄積領域（Accumulation regions）を定義し、合計 208 個のカウンターを FPGA 内に実装しました。
  • 平均カウンター法（Average Method）: 各ビームクロスイングタイプ（bb, be, eb, ee）ごとに再構成されたクラスター数を累積し、平均値を算出します。これが主たる輝度推定器として採用されました。
  • Log0 法（Zero-counting）: クラスターが検出されないイベントの割合から輝度を推定する手法です。飽和効果に強い特性を持ち、平均法の線形性の検証用クロスチェックとして利用されました。
• 統計処理:
  • 読み出しボード上で「移動平均（Running Mean）」を計算し、統計的揺らぎを低減します。パラメータは $\lambda=0.875$, $N=1024 \times N_{bb}$ として設定され、約 90ms ごとに更新されます。
  • 最終的なグローバル推定器として、外れ値を除去した「トリムド平均（Trimmed Mean）」を採用し、特定のカウンターノイズや故障に対する頑健性を確保しています。

3. 主要な成果と結果

2024 年の LHCb 物理運転期間（pp 衝突および PbPb 衝突）において、この手法は実運用され、以下の性能が確認されました。

• 統計的分解能:
  • 輝度測定値の統計的分解能は1% 未満（実際には 0.27% 程度）であり、これはビームレベリングに必要な 5% の精度を十分に満たしています。
  • 時間分解能は 100ms 未満（分析用には 3 秒にサンプリング）であり、極めて高速な応答性を示しています。
• 線形性:
  • 平均カウンター法は、LHCb の運転点（Run 3）を含む広範囲の輝度において、非線形性がゼロとみなせるレベル（$<0.5\%$）であることを確認しました。
  • 飽和効果（ヒットの重なり）は、適切な距離（ビーム軸から 11mm 以上）の領域を選択することで無視できるレベルに抑えられています。
• 安定性と校正:
  • van der Meer（vdM）スキャンを用いた校正により、各カウンター間の相対的なばらつきは pp 衝突で 0.8%、PbPb 衝突で 0.9% 以内と非常に良好でした。
  • 2 週間の運転データにおいて、PLUME 検出器によるオンライン輝度との比は**0.3%**の安定性を示しました（ビームスポット位置の補正後）。
• 重イオン衝突（PbPb）への適用:
  • PbPb 衝突は pp 衝突に比べて約 70 倍のクラスター密度を持ちますが、FPGA アルゴリズムの最適化と平均法の特性により、この手法も高占有率環境で有効に機能することが確認されました。
• 可用性:
  • システムの可用性は安定ビーム時間の約 93% に達しており、ダウンタイムは極めて少なくなっています。

4. 技術的意義と貢献

• 検出器レベルでの独立した輝度測定: ソフトウェア再構成や外部トリガーに依存せず、検出器読み出しレベル（FPGA ファームウェア）で完結する独立した輝度モニタリングを実現しました。これにより、レベリングフィードバックループの信頼性が向上し、クロス校正が可能になりました。
• 柔軟性と拡張性: FPGA 実装であるため、蓄積領域の形状、サイズ、位置、積分パラメータを柔軟に調整可能です。将来的には、ビームスポット位置のリアルタイム推定など、他の物理量への応用も可能です。
• リアルタイム性能: 従来のオフライン解析や高レベルトリガーを待たずに、ミリ秒単位で輝度情報を提供できるため、実験運転の最適化や異常検知に極めて有効です。
• 非線形性の抑制: 従来のヒットカウント法で見られた飽和や非線形性を、適切な領域選択と統計手法（Log0 法との比較、トリムド平均）によって効果的に制御・評価することに成功しました。

結論

本論文で提案された VELO ピクセルヒットに基づくリアルタイム輝度測定手法は、LHCb Run 3 以降の運転において、高精度・高安定性・高可用性を兼ね備えた主要な輝度モニタリング手段として確立されました。この手法は、FPGA 技術を活用したデータ取得システムの高度化を示す成功例であり、将来の高輝度実験におけるリアルタイム制御システムの重要な基盤技術となっています。