抽取信号的χ2分布拟合方法

BESⅡ上测量ψ(2S)辐射衰变时, 鉴于本底道的污染非常严重, 提出了一种抽取信号事例数的χ²拟合方法. 大量的蒙特卡罗模拟表明, 信号道与各个本底道的χ²分布形状各不相同. 这些特点提供了从数据的χ²分布中抽出信号道事例数的χ²拟合方法. 从蒙特卡罗模拟的结果来看, χ²拟合方法的误差小于2%.     

北京谱仪主漂移室的径迹重建

本文阐述了北京谱仪主漂移室的径迹寻找和径迹拟合程序的原理及使用情况.从五百万北京谱仪J/ψ事例以及相应的J/ψ衰变道的蒙特卡罗数据分析结果,证明了程序的可靠性.     

北京谱仪Ⅲ电磁量能器触发系统的模拟研究

利用北京谱仪Ⅲ探测器模拟和触发系统模拟程序,研究了北京谱仪Ⅲ电磁量能器的触发方案.根据物理目标设计了中性事例触发,巴巴事例判选和带电事例触发的方案,优化了各触发条件的参数,并仔细研究了对本底的排斥能力.研究了一批典型的物理道,给出了它们的触发效率,同时给出了预期的本底触发事例率.     

R值测量中束流相关本底的扣除方法

利用北京谱仪R值扫描数据,研究了束流相关本底的特征及其扣除方法.采用f因子方法和拟合强子事例顶点z分布的方法来扣除残余束流相关本底,并讨论其对R值测量误差的贡献.结果表明,两种方法对R值测量误差的影响约为0.3%至2.3%.     

纯中性事例事例顶点z坐标的最小二乘拟合

对于相互作用顶点z坐标未知的纯中性事例, 发展了一种运动学拟合的方法. 经蒙特卡罗样本检验, 拟合所用的程序是合理可靠的. 把该方法用于ψ′→γπ⁰→γγγ蒙特卡罗样本, 得到了更高的选择效率. 该方法可普遍应用于正负电子对撞实验中纯中性事例的分析.     

双光束干涉仪中椭圆拟合估算的参数精度研究

双光束干涉仪的相位调制解调通常会用到椭圆拟合的方法来计算相应的解调参数,提出了一种用于评价这类参数的计算精度的方法。根据自适应的蒙特卡罗法,计算不同噪声水平下椭圆拟合估算的参数的最大偏差;根据实验测得的干涉信号噪声本底的大小,判断拟合的参数的精度。以基于3×3耦合器的光纤水听器的两路输出方案为例,该方案为提高水听器的解调性能,需研究参数的数值随器件性能的变化规律,因此首先要判定实验中通过椭圆拟合所求得的参数的准确性。将该精度评定方法运用在水听器测量实验中,测量的干涉信号在20 Hz~4 k Hz处的噪声本底接近-70 d B·V2/Hz,根据自适应的蒙特卡罗法,得到参数的偏差不超过0.06%和0.08°。     

BESⅢ触发系统MDC子系统的物理设计

介绍北京谱仪BESⅢ触发系统MDC子系统的物理设计. 该设计根据BESⅢ的物理研究目标,用软件模拟的方法研究了不同触发方案的可行性,并优化了触发参数. 对触发方案的效率进行了仔细研究,给出一些重要反应道的触发效率和本底事例率.     

神经元网络法重建LHC上H→ZZ(WW)→llvv的不变质量

在LHCpp对撞机上通过H→ZZ(WW)→llvv道寻找重Higgs粒子的主要困难在于,该反应道具有大的丢失能量,用传统方法无法重建它的不变质量,其存在只能由p_T或m_T分布的宽Jacobi峰显现。这里作为基于蒙特卡罗模拟的唯象性研究,设计了一个前馈式神经元网络,用以重建Higgs粒子的不变质量,不仅得到了正确的M_H质量峰值及比较好的宽度,并同时具有抗本底事例干扰能力,可应用于实验的共振态新粒子寻找及粒子质量的精确测量。     

BEPCⅡ束流相关本底的蒙特卡罗模拟研究(Ⅱ)—束流-气体相互作用

介绍了用蒙特卡罗方法研究正常对撞模式下束流本底中的束流-气体相互作用的原理与效果.使用通用的模拟工具,研究了BEPCⅡ的束流-气体本底对设计中的BESⅢ的影响.结果显示,在储存环上设置挡板后,探测器的本底在允许的范围之内.模拟结果对BEPCⅡ-BESⅢ的设计建造具有指导意义.     

高温热辐射输运模拟的蒙特卡罗全局降方差方法

针对隐式蒙特卡罗方法模拟高温热辐射输运问题时存在的“辐射强度计算误差时间空间分布严重不均匀”现象,通过理论分析和数值模拟手段,找到决定误差大小的主要因素为“网格内的辐射径迹长度记录数”.据此提出“隐式蒙特卡罗全局降方差方法”并推导相应的计算公式.该方法主要包含如下3个关键技术:1)针对辐射输运蒙特卡罗模拟的自适应动态分配源粒子方法; 2)与自适应动态分配源粒子相匹配的动态权窗设计技术及粒子权无偏估计算法; 3)辐射强度的解析估计降方差方法.针对这3个关键技术,设计蒙特卡罗数值模拟方案,编写相应的数值模拟程序模块.典型辐射输运问题模拟结果显示:隐式蒙特卡罗全局降方差方法能够使网格辐射强度计算误差在整个时空范围内分布相对比较均匀,最大误差可控,计算效率提升10倍左右.新方法在激光惯性约束聚变的黑腔辐射输运模拟应用中取得了显著效果.     

硅纳米薄膜中声子弹道扩散导热的蒙特卡罗模拟

通过建立声子散射概率函数描述声子在输运过程中的散射,提出了一种模拟声子弹道扩散导热的蒙特卡罗方法,并将其应用于硅纳米薄膜中的稳态和瞬态弹道扩散导热过程的研究. 提出的蒙特卡罗方法对边界发射的声子束进行跟踪,根据散射概率函数模拟声子束在传播区域内经历的散射过程,并通过统计声子束的分布得到温度分布. 稳态导热过程的模拟发现,尺寸效应会引起边界温度跳跃,其值随着Knudsen数的增大而增大;计算的硅纳米薄膜的热导率随着厚度的增大而增大,与文献中的实验数据和理论模型相符. 通过瞬态导热过程的模拟得到了纳米薄膜内的温度分布随时间的变化,发现瞬态导热过程中的热波现象与空间尺度相关,材料尺寸越小,弹道输运越强,薄膜中的热波现象也越显著. 关键词： 纳米薄膜 弹道扩散导热 蒙特卡罗模拟 尺寸效应     

中子照相CCD成像系统反射镜引入本底分析

在基于CCD相机的中子照相系统中,反射镜距离闪烁屏太近就会将部分荧光反射回闪烁屏,将闪烁屏照亮,入射中子束的少部分还会被反向散射回闪烁屏,形成图像本底叠加在图像上,对定量分析和CT重建结果产生影响。为此建立了反射镜所引入反射分量的计算方法,可根据闪烁屏受照分布和反射镜参数计算反射分量的分布,并对中子反向散射进行了蒙特卡罗模拟,计算结果与实验测量的本底分布规律相符。     

混浊物质中光子传输的快速模拟与干涉特性分析

建立了混浊物质中低相干干涉系统的蒙特卡罗模型,结合相干面光强的高斯分布和低相干光源的相干条件,对干涉信号进行了数值分析.提出在蒙特卡罗计算的基线数据基础上快速模拟其他参数条件下干涉信号的方法.快速模拟结果与蒙特卡罗计算结果进行了比较,结果显示在满足前向散射条件下,两种方法所得模拟信号强度基本吻合.当各向异性因子g=0....     

噪声对漂移室性能影响的模拟研究

利用蒙特卡罗方法模拟了漂移室的随机噪声和对其性能的影响,分别考虑了击中丝和非击中丝的时间道和电荷道的噪声,给出了径迹效率和动量分辨随噪声水平和噪声分布的变化情况.     

蒙特卡罗方法模拟薄膜电致发光器件中碰撞离化的作用

基于经验赝势法得到的能带结构数据，采用分段多项式拟合获得ZnS能带结构的解析表达式 ，建立解析能带模型.使用建立的模型计算得到各能谷的态密度和总的散射速率,并与文献 的计算结果进行了对比，验证该解析能带模型既具有非抛物型多能谷能带模型运算速度快、 使用方便的优势,又具有与采用全导带模型相近的计算精度.进一步利用该模型进行蒙特卡罗 模拟，得到第一导带和第二导带中电子数随电场强度的变化、不同电场中能量分布函数以及 包含与不包含碰撞离化情况下电子能量随时间变化的曲线.讨论在外加电场下，电子在导带 内各个能谷间和 关键词： 蒙特卡罗模拟 解析能带模型 多项式拟合 碰撞离化     

反应堆中微子实验的基线优化

介绍了在χ2假设检验中引入本底和系统误差的方法，并通过最大化混合角θ13的测量敏感度对反应堆中微子实验的基线长度进行优化. 考察了对基线有影响的几个因素，包括振荡效应，统计误差，系统误差，本底等. 尽管实验条件可能很不相同，最优基线仍可被确定在一个较小的范围内.     

中子照相CCD成像系统反射镜引入本底分析

 在基于CCD相机的中子照相系统中,反射镜距离闪烁屏太近就会将部分荧光反射回闪烁屏,将闪烁屏照亮,入射中子束的少部分还会被反向散射回闪烁屏,形成图像本底叠加在图像上,对定量分析和CT重建结果产生影响。为此建立了反射镜所引入反射分量的计算方法,可根据闪烁屏受照分布和反射镜参数计算反射分量的分布,并对中子反向散射进行了蒙特卡罗模拟,计算结果与实验测量的本底分布规律相符。     

蒙特卡罗粒子输运问题的全局降方差方法

精确获得中子通量的时间-空间-能量精细分布对于一类非定常粒子输运蒙特卡罗模拟至关重要.以零方差理论为基础,采取化整为零的策略,把蒙特卡罗方法与离散纵标方法相耦合,用近似重要函数指导蒙特卡罗模拟,给出一种实现蒙特卡罗全局降方差的计算方案.数值算例表明,该计算方案获得了全局降方差的效果.     

蒙特卡罗模拟元素特征X射线注量的展宽技术

蒙特卡罗模拟技术可以有效提高核辐射探测分析效率、降低分析成本,在促进X荧光分析技术发展过程中也发挥了重要作用。在应用蒙特卡罗模拟能量色散X荧光分析中,为使蒙特卡罗模拟元素特征X射线注量谱与实验能谱之间的物理差异更小,提出了一种模拟注量展宽算法。首先依据离散统计物理理论,建立了Si-PIN半导体探测器对X射线的响应函数模型,并利用加权非线性最小二乘法拟合得到各系数,再对系数进行归一化处理,得到归一化后的响应函数R(E, E′)。然后利用R(E, E′)对模拟注量进行展宽转换,得到脉冲幅度谱(pulse height spectrum, PHS),使模拟谱更趋近于实验谱,这样就可以通过蒙特卡罗模拟得到与实测能谱一致性较高的模拟能谱。研究中最后利用建立的探测器响应函数模型对Ti, V, Cr, Mn, Fe, Ni, Cu, Zn, As, Sn元素K系全能峰拟合χ2_r值在1.04～1.18范围内,对铜合金样品中六种元素特征X射线模拟注量转换结果与实验值一致性较好,并实现了对多元素混合样品复杂重迭谱的有效模拟和展宽。     

蒙特卡罗剂量分布去噪中三维高斯和Savitzky-Golay滤波器的改进与混合

对模拟粒子轨迹数较少模拟时间较短的蒙特卡罗粗糙剂量分布进行三维滤波,可以加速其收敛速度.结合蒙特卡罗剂量分布特征,改进三维高斯和Savitzky-Golay滤波器,建立三维混合滤波方法,并比较并联和级联两种基本混合方式.根据卷积性质,提出用等效卷积核简化混合滤波器结构的方法.结果表明,改进后的高斯和Savitzky-Golay滤波器的整体去噪效果得以增强,混合滤波器进一步降低滤波结果的局部误差,两种混合滤波器都能够大幅度抑制MC粗糙剂量分布中的噪声,级联混合滤波器降噪效果略优于并联混合滤波器.