LHEReader: Simplified Conversion from Les Houches Event Files to ROOT Format

إليك ملخص تقني مفصل للورقة البحثية "LHEReader: تبسيط التحويل من ملفات أحداث Les Houches إلى تنسيق ROOT" بقلم أمان ديساي.

1. بيان المشكلة

تنتج محاكاة فيزياء الطاقات العالية (HEP) عادةً سجلات الأحداث على مستوى البارتونات باستخدام مولدات أحداث مونت كارلو (MC) مثل MadGraph5_aMC@NLO. تُخرج هذه المولدات البيانات بتنسيق ملفات Les Houches Event (LHE)، وهو تنسيق نصي قياسي يحتوي على هويات الجسيمات، رموز الحالة، الزخم الرباعي، اللف المغزلي، معلومات اللون، وأوزان الحدث.

بينما يعد تنسيق LHE معيارياً لعملية التوليد، فإن مراحل التحليل اللاحقة (مثل عملية التذرية أو "parton showering"، أو الهدرنة عبر Pythia8، أو محاكاة الكاشف عبر MadAnalysis5) غالباً ما تتطلب البيانات بتنسيق أكثر كفاءة، ثنائي وثنائي الأبعاد، ومنظم ككائنات (object-oriented) مناسب للتحليل الإحصائي عالي الأداء. ويعد إطار عمل ROOT (الذي طوره CERN) هو المعيار الصناعي لهذا التحليل. ومع ذلك، فإن تحويل ملفات LHE إلى ملفات ROOT يتطلب غالباً كتابة سكربتات مخصصة ومعقدة أو الاعتماد على مكتبات خارجية قد لا تكون محسنة لسير عمل محدد. كانت هناك حاجة لأداة خفيفة الوزن، سهلة الوصول، وفعالة لسد هذه الفجوة، مما يسمح للباحثين بتحويل ملفات LHE مباشرة إلى هياكل TTree بتنسيق ROOT للتحليل الفوري.

2. المنهجية

تقدم الورقة البحثية LHEReader، وهو برنامج قائم على لغة بايثون مصمم لتحليل ملفات LHE وتحويلها إلى تنسيق ROOT.

• بنية البرمجيات:

  • اللغة: بايثون (تم اختباره مع الإصدار 3.10).
  • التبعيات: يتطلب CERN ROOT و PyROOT.
  • المدخلات: يدعم ملفات .lhe غير المضغوطة وملفات .lhe.gz المضغوطة.
  • المخرجات: يُنشئ ملف .root يحتوي على كائن TTree باسم events.

• تحويل بنية البيانات:
  يقوم البرنامج برسم خرائط لهيكل LHE الشبيه بـ XML إلى TTree مهيكل مع فروع (branches) محددة:

  • فروع على مستوى الحدث (Event-level branches): numParticles (عدد الجسيمات)، eventweight (وزن الحدث)، scale (المقياس)، alpha_qed (ثابت البنية الدقيقة الكهرومغناطيسية)، alpha_qcd (ثابت البنية الدقيقة القوية).
  • فروع على مستوى الجسيم (Particle-level branches - Vectors): pid (معرف الجسيم)، status (الحالة)، mother1/2 (الأم 1/2)، color1/2 (اللون 1/2)، px/py/pz (الزخم)، energy (الطاقة)، mass (الكتلة)، tau (الزمن النسبي)، و spin (اللف المغزلي).

• سير عمل التحليل (EventLoop):
  يوفر المؤلف نموذجاً بلغة ++C يسمى EventLoop لتوضيح كيفية معالجة ملف ROOT الناتج. تقوم هذه الحلقة بـ:

  1. التكرار عبر أحداث الـ TTree.
  2. التكرار عبر الجسيمات داخل كل حدث.
  3. تصفية الجسيمات بناءً على رموز الحالة (على سبيل المثال، status == 1 للجسيمات في الحالة النهائية).
  4. تحديد الجسيمات باستخدام معرفات PDG (على سبيل المثال، abs(pid) == 5 للكواركات القاعية b-quarks، و 11/13 للبتونات).
  5. حساب المتغيرات الكينماتيكية، بما في ذلك طاقة المفقود المستعرضة (MET) بناءً على تحديد النيوترينو.

3. المساهمات الرئيسية

• أداة LHEReader: أداة مفتوحة المصدر ومتاحة للجمهور (عبر GitHub) تبسط عملية التحويل من LHE إلى ROOT، مما يلغي الحاجة إلى محللات (parsers) مخصصة ومعقدة.
• اختبار الأداء: تقدم الورقة بيانات تجريبية على سرعات التحويل لعمليتين فيزيائيتين مختلفتين:
  • $pp \to jj$ (إنتاج كواركات ثنائية QCD).
  • $e^+e^- \to ZH$ (إنتاج هيجز عبر إلكترون-بوزيترون).
• عرض فيزيائي شامل (End-to-End): لا تكتفي الورقة بعرض تحويل الملفات فحسب؛ بل تستعرض سلسلة تحليل كاملة:
  1. توليد أحداث $pp \to Z(\to \ell^+\ell^-)H(\to b\bar{b})$ باستخدام MadGraph5_aMC@NLO.
  2. عمليات التذرية والهدرنة عبر Pythia8.
  3. تجميع النفاثات (Jet clustering) عبر FastJet (ضمن MadAnalysis5).
  4. التحويل إلى ROOT عبر LHEReader.
  5. تحليل التوزيعات الكينماتيكية باستخدام EventLoop المقدم.
• مورد تعليمي: إن تضمين مقتطفات برمجية لكل من سكربت التحويل وحلقة تحليل ++C يعمل كدليل تعليمي للطلاب والباحثين حول كيفية التعامل مع بيانات أحداث MC في بيئة ROOT.

4. النتائج

• أداء التحويل: قيمت الدراسة الوقت المطلوب لتحويل الأحداث إلى تنسيق ROOT عبر أحجام عينات متغيرة ($10^4$ إلى $10^6$ حدث).

  • $pp \to jj$: انخفض متوسط وقت التحويل لكل حدث من 0.27 مللي ثانية (عند $10^4$ حدث) إلى 0.18 مللي ثانية (عند $10^6$ حدث).
  • $e^+e^- \to ZH$: انخفض متوسط وقت التحويل لكل حدث من 0.14 مللي ثانية (عند $10^4$ حدث) إلى 0.03 مللي ثانية (عند $10^6$ حدث).
  • ملاحظة: الأداة تتوسع بكفاءة، حيث تظهر الدفعات الأكبر تحسناً في أوقات المعالجة لكل حدث، ويرجع ذلك على الأرجح إلى تقليل تأثير عبء الإدخال والإخراج (I/O overhead).

• التحقق من التحليل الفيزيائي: تم تطبيق البرنامج لمحاكاة عملية $pp \to ZH$ عند طاقة مركز الكتلة $\sqrt{s} = 13.6$ تريليون إلكترون فولت (ظروف تشغيل LHC Run 3).

  • قيود الاختيار (Selection Cuts): تم تصفية الأحداث للحصول على ليبتونين متضادين في الشحنة ومتماثلين في النوع ومعزولين ($p_T > 25$ GeV، $|\eta| < 2.5$) مع كتلة ثابتة في نافذة بوزون Z [80, 100] GeV.
  • إعادة البناء: تم تحديد الـ b-jets (على مستوى الحقيقة للتحليل البارتوني، وبناءً على PID للتحليل المعاد بناؤه) واستخدامها لإعادة بناء مرشح هيجز.
  • التوزيعات: تعرض الورقة توزيعات معيرة لـ:
    • الزخوم المستعرضة ($p_T$) للميونات، الإلكترونات، والـ b-jets.
    • الكتلة الثابتة لأزواج الإلكترون-إلكترون والميون-ميون (التي تبلغ ذروتها عند كتلة Z).
    • الكتلة الثابتة لزوج الـ b-jet (التي تبلغ ذروتها عند كتلة هيجز).
  • أكدت النتائج أن ملفات ROOT المحولة حافظت بدقة على الخصائص الكينماتيكية لأحداث LHE الأصلية، حيث أظهرت التوزيعات المعاد بناؤها الاتساع المتوقع نتيجة لتأثيرات تجميع النفاثات ومحاكاة الكاشف.

5. الأهمية

• كفاءة سير العمل: تقلل LHEReader بشكل كبير من حواجز الدخول لتحليل بيانات أحداث MC. من خلال أتمتة التحويل إلى ROOT، تتيح للفيزيائيين الاستفادة من كامل قوة إطار عمل تحليل ROOT (الرسوم البيانية، الملاءمة، الأدوات الإحصائية) دون الحاجة لكتابة محللات XML منخفضة المستوى.
• تعدد الاستخدامات: تدعم الأداة ملفات LHE القياسية والمضغوطة، وهي متوافقة مع مخرجات مختلف المولدات (MadGraph, MadAnalysis5)، مما يجعلها أداة متعددة الاستخدامات لمختلف مشاريع أبحاث فيزياء الطاقات العالية (بما في ذلك دراسات LHC و FCCee المستقبلية).
• قابلية التكرار: من خلال توفير الكود المصدري وسكربتات التحليل المحددة (EventLoop)، تعزز الورقة البحثية قابلية تكرار التحليلات في فيزياء الطاقات العالية، مما يسمح للباحثين الآخرين باختبار التوزيعات الكينماتيكية بسرعة أو التحقق من نماذج الفيزياء الجديدة مقابل خلفيات النموذج القياسي.
• القيمة التعليمية: تعمل الورقة كدليل عملي للانتقال من توليد الأحداث إلى تحليل البيانات، وهي مهارة حاسمة لطلاب الدراسات العليا والباحثين في بداية مسيرتهم في فيزياء الجسيمات.