1. GEANT4模拟在MATHUSLA探测器中的核心作用粒子物理实验的设计与优化离不开精确的蒙特卡洛模拟。GEANT4作为高能物理领域事实标准的模拟工具包通过跟踪粒子在物质中的输运过程为探测器性能评估提供了可靠的计算框架。在MATHUSLA探测器的设计中GEANT4模拟主要承担三个关键任务背景量化精确计算宇宙射线(CR)、LHC对撞产生的μ子、大气中微子等背景源在探测器内的相互作用率信号响应建模模拟长寿命粒子(LLP)衰变产物在探测器中的行为包括径迹形成、能量沉积等物理过程几何优化通过参数化扫描确定探测器各组件闪烁体层、支撑结构、否决层等的最佳布局方案提示GEANT4模拟中需要特别注意材料边界的二次散射效应。在MATHUSLA的模拟中所有结构材料如铝制闪烁体外壳、钢支撑柱都需精确建模其厚度误差控制在1mm以内。2. MATHUSLA探测器几何建模详解2.1 基本结构参数MATHUSLA采用分层式设计主要包含以下功能单元天花板跟踪模块6层闪烁体阵列层间距80cm距地面高度11.8m开始布置单层覆盖面积42.8m×42.8m闪烁体条尺寸1cm(厚)×4.5cm(宽)×450cm(长)后墙跟踪模块6层垂直放置的闪烁体阵列距地面1.85m起始向上延伸9m采用与天花板相同的闪烁体规格否决层系统地面否决层2层全封闭式闪烁体阵列前墙否决层2层全封闭式设计层间距80cm确保对穿透粒子的高效探测2.2 材料与支撑结构探测器材料预算的精确模拟对背景抑制至关重要闪烁体封装每层包含0.3cm铝制外壳模拟实际的光导和机械支撑结构钢支撑柱垂直支柱14英寸方钢壁厚5/8英寸水平横梁8英寸方钢壁厚1/2英寸地面结构85.47m厚的混凝土层模拟CMS实验大厅上方地层2.3 几何简化与验证为平衡计算效率与精度模拟中进行了合理简化省略了前墙闪烁体的交错排列设计支撑柱与地面探测器间保留1cm气隙地面探测器置于地表而非地下数米通过对比5层/6层配置的模拟结果验证了几何简化对物理结论的影响小于5%。关键验证指标包括几何接受度(Geometric acceptance)顶点重建效率(Vertex reconstruction efficiency)背景否决率(Background rejection power)3. 背景模拟与抑制策略3.1 主要背景源建模3.1.1 宇宙射线背景使用CRY/PARMA生成器模拟初级宇宙射线质子通量1.04×10¹¹个/60m×60m平面对应5年HL-LHC运行时间的1/32中子通量8×10¹¹个/60m×60m平面对应5年HL-LHC运行时间的1/27能谱覆盖10MeV-1TeV符合511m海拔日内瓦地区的预期分布3.1.2 LHC μ子背景通过MadGraph5Pythia8链模拟产生机制s-channel γ*/Z*/W*交换截面积分pT10GeV|η|2.5K因子1.3NLO修正总样本量6.1×10⁹个μ子相当于HL-LHC总μ子通量的1/103.1.3 大气中微子采用GENIEFLUKA3D组合通量模型三维大气中微子角分布靶材料空气、钢、土壤、探测器材料预期事件率6次/年整个探测器体积内3.2 背景否决技术通过七级筛选实现背景压制空间否决边界收缩0.3m的基准体积钢支柱周围1m²区域排除时间否决顶点前100ns时间窗内允许≤1个hit有效抑制宇宙射线二次相互作用速度筛选剔除β5/6的非相对论性径迹排除表观超光速的假径迹方向一致性顶点锥角与CMS对撞点方向偏差0.3rad利用LLP衰变产物的前冲特性径迹多重性要求顶点包含≥4条径迹可排除80%的CR质子/中子背景向下径迹比例事件中向下径迹占比阈值控制平衡信号效率与背景抑制非关联hit比例剔除与顶点时空关联性差的hit聚集针对中子散射特有的hit分布特征4. 信号重建关键技术4.1 径迹重建算法采用基于Kalman滤波的递推算法种子形成要求至少2个不同层的hit构成双联点空间间隔ΔR15cm时间间隔Δt10ns径迹延伸逐层搜索符合预测的hit考虑多次散射效应等效0.5GeV μ子参数提取最终拟合包含≥4个hit角度分辨率1°对10GeV/c μ子性能指标效率90%p0.5GeV/c误建率0.1%6层配置4.2 顶点重建流程顶点种子由径迹对构成要求最小接近距离(DCA)10cm到达时间差5ns顶点扩展逐步添加符合χ²阈值的径迹最大包含径迹数20条四维拟合最小二乘法确定最佳(x,y,z,t)分辨率10cm(横向)11cm(纵向)4.3 电子学读出模型为模拟真实探测器响应空间分辨率沿条方向2.25cm高斯展宽横向13cm取最近条中心时间分辨率1ns高斯分布能量阈值单条600keV20ns积分探测效率基准值95%范围90-100%噪声模拟基准300kHz/ch等效宇宙射线本底极端情况3MHz/ch10倍噪声5. 物理性能评估5.1 信号效率分析对h→XXX→bb过程mX15-55GeV质量(GeV)几何接受度重建效率总效率1598.2%53%52%2597.8%49%48%3597.5%42%41%4597.1%37%36%5596.7%32%31%效率下降主因高质量LLP的boost因子降低衰变产物角分布更各向同性5.2 背景抑制能力通过ABCD方法估计全期背景背景源原始事件率剩余事件(5年)宇宙射线质子6.4×10⁷0.2宇宙射线中子2.2×10⁷0.1LHC μ子2.4×10⁴0.2大气中微子240关键发现径迹多重性筛选贡献最大否决功率方向一致性对LHC μ子特别有效四层配置会使背景增加10倍5.3 物理灵敏度展望对40m×40m基准设计可探测分支比Br(h→XX)~10⁻⁶寿命覆盖cτ~10-1000m相比LHC主探测器提升1-2个数量级在暗胶球搜索方面可探测暗QCD能标Λ~10-100GeV对Fraternal Twin Higgs模型敏感区伴顶夸克质量1TeV暗胶球质量20-60GeV6. 经验总结与优化方向在实际模拟中发现几个关键经验点材料效应处理钢支撑柱导致的多次散射会恶化径迹分辨率解决方案在支柱周围增设包裹式探测器层时间关联优化600ns读出窗口可平衡本底与效率过短会截断长寿命次级粒子径迹几何布局改进后墙跟踪层提升接受度15-20%地面下挖设计可进一步抑制宇宙射线未来优化方向引入机器学习辅助顶点重建开发快速模拟框架(FastSim)研究塑料闪烁体替代方案降低成本