Afterpulse prediction for SUBMET experiment

إليك ملخص تقني مفصل لورقة "التنبؤ بالنبضات اللاحقة لتجربة SUBMET" من قبل Campagnari وآخرون.

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

تهدف تجربة SUB-Millicharge ExperimenT (SUBMET) في J-PARC إلى الكشف عن الجسيمات ذات الشحنة المليونية ($\chi$) التي تقل كتلتها عن $1.6 \text{ GeV}/c^2$ وشحنتها $Q_\chi < 10^{-3}e$. يستخدم الكاشف 160 وحدة، تتكون كل منها من ومض بلاستيكي (Eljen-200) متصل بأنبوب مضاعف ضوئي (PMT) من طراز Hamamatsu R7725.

يتمثل التحدي الحرج في هذه التجربة في وجود النبضات اللاحقة (afterpulses) للأنبوب المضاعف الضوئي. وهي عبارة عن إشارات متأخرة ناتجة عندما تؤين الإلكترونات الضوئية شوائب الغاز بين الكاثود الضوئي والداينودات، مما يؤدي إلى إنشاء إلكترونات ضوئية ثانوية بعد مئات النانو ثواني.

• المشكلة: لا يمكن تمييز النبضات اللاحقة عن الإشارة المستهدفة (نبضات الفوتوإلكترون المنفرد أو SPE).
• السياق: يتكون شعاع J-PARC من ثماني حزم تفصل بينها فترة زمنية قدرها $\sim 60$ نانو ثانية. تولد النبضات الكبيرة (الارتفاع $> 140$ ملي فولت) الناتجة عن الجسيمات عالية الطاقة في حزمة واحدة نبضات لاحقة يمكن أن تتداخل مع إشارات SPE في الحزم اللاحقة، مما يزيد بشكل كبير من معدل الخلفية.
• الهدف: تطوير نموذج تنبؤي لمعدلات النبضات اللاحقة بناءً على معاملات قابلة للقياس لتقدير مساهمات الخلفية بدقة وتحسين موثوقية تحليل بيانات التجربة.

2. المنهجية

قام المؤلفون بتطوير نموذج تنبؤي يعتمد على بيانات الشعاع التي جُمعت في يونيو 2024، باستخدام تأخير تحفيز (trigger delay) مخصص للتركيز على النافذة الزمنية التي تلي آخر حزمة شعاع. تكونت مجموعة البيانات من 50,000 حدث، تم تقسيمها بالتساوي بين ملاءمة النموذج والتحقق من الأداء.

الملاحظات الرئيسية والافتراضات

يفترض النموذج أن معدل النبضات اللاحقة ($dn/dt$) يعتمد على:

1. مساحة النبضة الكبيرة ($A$): الشحنة المتكاملة للنبضة الأولية التي تولد النبضات اللاحقة.
2. الثابت الزمني ($\tau$): معدل اضمحلال احتمالية النبضة اللاحقة بمرور الوقت.
3. النافذة الزمنية: يركز التحليل على النافذة $t \in [720 \text{ ns}, 3475.2 \text{ ns}]$ بعد النبضة الكبيرة. يتم استبعاد أول 720 نانو ثانية (المنطقة 1) بسبب تقلبات خط الأساس (baseline).

بناء النموذج

يتضمن إطار التنبؤ مكونين رئيسيين:

1. الاعتماد على مساحة النبضة ($A$):
   وُجد أن العدد الإجمالي للنبضات اللاحقة ($n$) في نافذة زمنية ثابتة يرتبط إيجابياً بمساحة النبضة الكبيرة ($A$). تم اختبار شكلين وظيفيين:

   • الملاءمة الخطية (Linear Fit): $n_{pred}^{lin}(A) = p_0^{lin} + p_1^{lin} A$
   • الملاءمة الأسية (Exponential Fit): $n_{pred}^{expo}(A) = \exp(p_0^{expo} + p_1^{expo} A)$

2. البنية الزمنية:
   يتبع توزيع الفرق الزمني ($\Delta t$) بين النبضة الكبيرة والنبضة اللاحقة اضمحلالاً أسياً: $N(\Delta t) \propto \exp(-\Delta t / \tau)$. وُجد أن الثابت الزمني $\tau$ مستقل عن مساحة النبضة $A$ ويُقاس بشكل فردي لكل وحدة من الوحدات الـ 160.

3. صيغة التنبؤ:
   من خلال دمج اعتماد المساحة مع البنية الزمنية، يتم حساب عدد النبضات اللاحقة المتوقع في نافذة زمنية تعسفية $[t_i, t_f]$ كالتالي:
   $$n_{pred}(A, t_i, t_f) = \frac{1}{\tau} \int_{t_i}^{t_f} n_0(A) e^{-t'/\tau} dt'$$
   حيث يتم اشتقاق $n_0(A)$ عن طريق تطبيع (normalization) الملاءمة المعتمدة على المساحة ($n_{pred}(A)$) عبر نطاق التكامل $[t_0, t_1]$.

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

• إطار العمل مزدوج النماذج: تقدم الورقة نموذجين تنبؤيين متميزين (خطي وأسي) يربطان مساحة النبضة الكبيرة بعدد النبضات اللاحقة، وتم التحقق من صحتهما مقابل بيانات الشعاع الحقيقية.
• المعايرة الخاصة بكل وحدة: إدراكاً بأن خصائص الأنبوب المضاعف الضوئي تختلف، قام المؤلفون بقياس الثابت الزمني $\tau$ ومعاملات الملاءمة ($p_0, p_1$) بشكل مستقل لجميع وحدات الكاشف الـ 160.
• استراتيجية تخفيف الخلفية: تسمح هذه الطريقة للتجربة بتضمين الأحداث التي تحتوي على "نبضات كبيرة" (والتي كانت سابقاً عرضة للاستبعاد أو التعامل معها كخلفية عالية) في التحليل من خلال طرح مساهمة النبضات اللاحقة المتوقعة بدقة.

4. النتائج

• نطاقات المعاملات:
  • معاملات النموذج الخطي: $p_0 < 4$، $p_1 \in [3, 14] \times 10^{-5}$.
  • معاملات النموذج الأسي: $p_0 < 2$، $p_1 \in [1, 4] \times 10^{-5}$.
  • متوسط الثابت الزمني $\tau$: $\approx 1450$ نانو ثانية (الانحراف المعياري: 140 نانو ثانية).
• أداء النموذج:
  • أعاد كل من النموذجين الخطي والأسي إنتاج معدلات النبضات اللاحقة المرصودة بدقة تقارب 20%.
  • كانت توزيعات قيمة "كا تربيع المختزلة" ($\chi^2/ndf$) لكلا النموذجين متقاربة، حيث بلغت القيم المتوسطة 0.731 (للخطي) و 0.750 (للأسي).
  • تم اختيار النموذج الخطي لمزيد من التحليل نظراً لبساطته.
• التحقق من الصحة:
  • يظهر الشكل 8 (في الورقة) إجمالي النبضات اللاحقة المتوقلة مقابل المرصودة لجميع الوحدات الـ 160، مما يظهر توافقاً قوياً ضمن حدود عدم اليقين الإحصائي.
  • نجح النموذج في التنبؤ بالتوزيع الزمني للنبضات اللاحقة، حيث وقعت الأعداد المتوقعة ضمن نطاقات عدم اليقين $\pm 1\sigma$ للبيانات المرصودة.

5. الأهمية

• استغلال البيانات: بما أن أكثر من 70% من أحداث التصادم تحتوي على نبضة كبيرة واحدة على الأقل، فإن القدرة على التنبؤ بالنبضات اللاحقة وطرحها بدقة أمر بالغ الأهمية. بدون هذا النموذج، ستكون غالبية بيانات التصادم غير قابلة للاستخدام بسبب ارتفاع عدم اليقين في الخلفية.
• موثوقية الخلفية: تعمل هذه الطريقة على تحسين موثوقية توقعات الخلفية للبحث عن الجسيمات ذات الشحنة المليونية بشكل كبير، مما يضمن عدم حجب الإشارات المحتملة بواسطة النبضات اللاحقة التي قد تُخطئ في تحديد هويتها.
• القابلية للتطبيق العام: يوفر هذا النهج إطاراً قوياً لتوصيف وتصحيح تأثيرات النبضات اللاحقة للأنبوب المضاعف الضوئي في تجارب فيزياء الطاقة العالية الأخرى التي تستخدم أنظمة قراءة مشابهة للمومضات والأنبوب المضاعف الضوئي.

في الختام، تؤسس هذه الورقة لطريقة تنبؤ معتمدة وخاصة بكل وحدة، مما يقلل من عدم اليقين في خلفية النبضات اللاحقة إلى حوالي 20%، مما يمكّن تجربة SUBMET من تعظيم حساسيتها للجسيمات ذات الشحنة المليونية في مساحة بارامترية لم تُستكشف من قبل.