北京谱仪的R测量及QCD检验

1强子产生子LUARLW参数协调及辐射修正计算中值得注意的几个问题胡海明R&QCD研讨会黄山学院2015.07.27-30北京谱仪上的R值测量和QCD实验检验R值的定义辐射修正问题费曼图方案用于产生器的优点结构函数方案用于产生器的不自然真空极化处理对狭窄共振态line-shape测量的影响LUARLW简介LUARLW参数协调及效率数值的困惑初始强子及短寿命粒子衰变由产生器模拟+探测器模拟初始强子由产生器模拟+短寿命粒子由EvnGen模拟+探测器模拟效率的差别R值误差估计事例选择误差与效率误差综合在一起估计未测量零叉事例的误差主要内容2R值测量原理实验测量公式：实验测量的基本工作是：数据获取(BEPC/BES运行取数,触发表和触发效率,原始数据的刻度和重建)数据分析(本底排除,强子事例选择,积分亮度测量)物理分析(强子化模型,强子效率,辐射修正因子,本底计算)误差分析(分析测量公式中各量的不确定性给出R值总误差)理论定义：3关于辐射修正辐射修正（基于Feynman图）R值的定义:(a)树图~对应R值；(b)韧致辐射；(c)电子子能；(d)顶角和真空极化。5物理上感兴趣的、反映过程基本物理性质的强子截面是Feynman图中树图所表示的相对QED过程最低阶截面(但包括QCD过程的任意阶高阶截面)。但实验数据分析中并不能把最低阶过程的事例与高阶过程的事例区分开来，只能选择所有的强子事例。从理论上计算出高阶过程产生事例所占的比例，即辐射修正因子(1+):用于计算辐射修正因子（基于Feynman图）实验上观测到的强子过程总截面：物理总截面：零阶截面(R值)：除以辐射修正因子，即扣除高阶费曼图的贡献按Feynman图，初态辐射修正来自如下各项的贡献：其中，虚辐射修正项：6Bose-Einstein关联用于确定蒙特卡罗模拟辐射事例辐射硬光子的角度和动量分布为零阶强子截面考虑了连续态和共振态的贡献.注：对于辐射事例,强子产生在其自身质心系中模拟，然后再把所有强子能-动量boost到总质心系(实验室系)。7在MC模拟中，与费曼图相对应，所产生的事例分为两类：①无辐射事例(零阶过程，软辐射，虚辐射。强子系统的能量等于初始正负电子的质心能量)权重为②有辐射事例(硬辐射强子系统的能量小于初始正负电子的质心能量)权重为（~45%）（~55%）辐射修正中的真空极化因子真空极化效应改变了电磁相互作用的强度，即把电磁耦合常数α变为跑动耦合常数α(s).共振态Born截面可表示为Breit-Wigner形式，在轻子宽度仍采用理论值来表示时，共振态本征截面为右图为用PDG所给的共振参数计算的从强子产生阈到5GeV的Breit-Wigner截面在狭窄共振态峰附近的电磁跑动耦合常数α(s)，左：J/Ψ峰附近右：Ψ(3686)峰附近在狭窄共振态峰附近的真空极化因子Π(s)实部(红)ReΠ(s)虚部Im(蓝色)左：J/Ψ峰附近右：Ψ(3686)附近经真空极化修正Breit-Wigner截面(把Breit-Wigner公式中轻子宽度的固定精细结构常数换成能量相关的跑动耦合常数时的BreitWigner截面)红色线为原来的Breit-Wigner截面黑色线为真空极化修正因子蓝色线为考虑了真空极化修正后的Breit-Wigner截面，左图：J/Ψ峰附近；右图：Ψ(3686)峰附近真空极化因子对Breit-Wigner截面影响LUARLW中辐射修正的改进计算一般文献认为，辐射修正项是相对于树图项是小量，考虑了辐射修正的总截面可以表示为δvert顶点修正，δvac真空极化，称为虚辐射修正，δγ轫致辐射。但e+e-质心能量在狭窄共振态附近，真空极化项将随能量迅速变化，其绝对值与1相比不仅不小，甚至会远大于1，因此上述近似假设将不成立。因此，在还原为原始的精确表达式后，考虑真空极化的作用下，把精细结构常数变为跑动耦合常数后，得到初态辐射修正因子为：12把真空极化近似展开还原为精确表示采用PDG数值，分别用原来方案(红色)和改进方案(蓝色)的狭窄共振态的总截面左图J/Ψ峰附近，右图Ψ(3686)峰附近辐射修正改进计算对line-shape的影响14辐射修正改进计算的(1+)2.000000049.32868355.3075081.1212038322.099999942.73715150.3610591.1783906632.150000140.77360248.2792241.1840804442.175000039.84110647.2681961.1864177752.200000038.93959046.2814351.1885444962.232399937.81520745.0405821.1910706172.309400135.32769042.2653251.1963795482.386400033.07452939.7248051.2010693892.395999932.80862139.4236081.20162342102.400000132.69875439.2990821.20185258112.500000030.12042836.3644621.20730233122.644399926.89991132.6687221.21445464132.700000025.79575831.3943471.21703526142.800000023.97412629.2836711.22146983152.900000122.34100927.3820091.22563884162.950000021.58946826.5030151.22758997172.980999921.14582825.9820911.22871004183.000000020.88384325.6730021.22932364………………………………………..……………Ecm0(nb)tot(nb)(1+)15BESIII重粲结构更精细扫描特别精细扫描区域的步长：1~2MeV；如有新结构，所获取的样本能提供足够的信息以确定新粒子态的全部参数。干涉的处理衰变道fc:相干性相加:振幅:非相干性相加:共振截面:棘手的理论问题：考虑干涉时的幺正性强子宽度的能量相关性辐射修正精度问题ISR的误差应该考虑上述两方面的总误差的精度会被的误差所淹没方案本身的精度在Feynman图的微扰计算中，精确至3阶的计算精度为~1%。零阶截面的误差QED过程：零阶截面是准确的(零误差)。强子过程：测量值误差，模型误差，参数值误差等。17关于产生子和参数协调19正负电子对撞反应过程强子事例的例子也可以通过中间机制产生(先产生中间态,再衰变为末态),例如:20正负电子湮没强子产生过程强子态可以直接产生,例如:强子末态21e+e-探测器：原始数据：事例重建：数据分析：物理结果：蒙特卡罗产生器模拟事例产生探测器模拟（粒子与物质相互作用，吸收，产生，散射，输运，噪音及本底真实化）蒙特卡罗模拟数据加速器粒子对撞击中，相互作用，产生，衰变，热噪声，本底，信号触发，数字化，数据各种信息MCMC探测效率事例重建数据、本底分析刻度，寻迹本底，信号产生事例高能物理实验流程系统不确定性(误差)测量值误差真实过程模拟过程探测效率问题从粒子对撞到获得物理结果经历如下过程：对撞产生Ngen个强子事例（有物理意义的量，未知）经实验测量得到Nobs个强子事例(测量值)真实的探测效率(实验分析中丢失事例的比例)是未知的！问题：如何知道实验中探测效率的数值？解决办法:MC！物理截面表示在对撞顶点强子产生的概率。但在实验测量中由于各种原因(事例重建，事例选择条件等)，所观测到的强子事例数少于产生的事例数。在数据分析中把所观测到的强子事例数与对撞产生的事例数之比称为探测效率22蒙特卡罗模拟两个基本要素蒙特卡罗包括两个方面：①事例产生模拟：物理理论和唯象模型要正确,或者是很好的近似。有些过程的产生器是现成可用的(例如e+e-,+等QED过程)；有些是标准程序库里没有的,需要实验者编写和调试(例如强子过程产生器)。②探测器模拟：必须反映真实的探测器结构、物质分布、性能、状态等,可靠地模拟粒子进入探测器后的散射、传播、与物质相互作用、衰变等。通常是已经编写、调试、检查好放入软件库中，实验者可以现成调用。蒙特卡罗技术可以对实验全过程进行模拟，从而可以得到探测效率.当粒子产生模型是正确合理的，探测器模拟是真实可靠的，则可以用蒙特卡罗的效率代替真实的探测效率。但必须给出蒙特卡罗效率的系统误差。23Lund面积定律弦碎裂图象量子真空涨落碎裂强子O轻夸克(u,d,s)质量近似于0，以光速沿光锥运动。但对重夸克(c,b)，质量很大，它们不沿光锥运动，因此面积需要修正。24LUND强子产生模型LUARLW强子产生过程：“矩阵元”：其中，Lund模型横动量分布(高斯形式)：，纵向动量分布(Lund面积定律)：，25弦碎裂为n强子的概率(类似于量子场论中的微分截面)：其中，碎裂强子的光锥动量分数完成对相空间变量的积分后得到弦碎裂为不同末态的分布。HuHaiming将选择的调节参数关于inclusive多重数:关于4-动量:关于单粒子产额比:PARJ(1-17)取其它实验(OPAL)的调节值，或考虑个别作微调LUARLW的模拟功能LUARLW可以按单举(部分过程按遍举)的方式模拟BEPC能区连续态和JPC=1共振态的产生和衰变过程。28参产生子初始强子态@4.6GeV........................统计到13体29参产生子衰变末态@4.6GeV........................统计到13体目前参数协调中的困惑方法：把真实数据与模拟数据经过相同的过程后进行比较TriggertrgRawdataGeneratorBESsimulationEventselectionNobshadNgenhadNgenMCNobsMCBESrawdateTuneparameters如果：1.强子产生模型正确2.探测器模拟真实可靠hadNobshadNgenhad=NobsMCNgenMC则：所有的蒙特卡罗模拟与真实数据在探测器层次的分布应很好符合好参数组判据。强子产生器参数协调3132参数协调的基本思路LUARLW/JETSET中包含众多的唯象参数，它们虽然有明确的物理意义，但其数值是未知的，需要通过与实验数据的比较来确定，此过程就是参数协调。从原则上说，强子事例的任何末态分布(不管是直接测量量，如能动量、顶点空间位置，角度、粒子产额比等，还是组合物理量，如快度、球度、冲度，扁度、非平面度等)都可参与协调，所用的分布越多，符合得越好，得到的参数越合理，效率越可靠。目前，参数协调目的是将LUARLW用于确定R值测量的强子效率，并具有较小的系统误差。鉴于误差分析是在确定的事例选择条件下比较数据与与MC的差别，因此，最实用的办法就是参数协调时选用全部与事例选择条件有关的分布，再加上尽可能多的其它分布，多多益善。33参数协调的具体方法桥梁式在所有能量点上用一套参数(或粲阈以下、以上各一套)统一调好，虽不可能做到在所有点上都非常好，但能保证在所有能量点上整体都好，效率随能量变化是一致的。跷跷板式选择在某一个能量点上调好(可以调得非常好)，然后运用于其它能量点，但不能保证都仍然好，很可能是很不好。34参产生子与探测器模拟的界面BESII实验方案LUARLW/JETSET能给出粒子的5-动量及径迹顶点的时空坐标，还有一套描述束流管的几何参数。可以做到产生子负责模拟初始强子在顶点的产生，以及寿命小于束流管尺度的粒子衰变。在LUARLW产生与探测器模拟SIMBES界面处写入的是寿命大于束流管尺寸的稳定粒子与长寿命粒子。SIMBES模拟的过程与真实实验中BESII所接收到的粒子一致，所得效率合理。35参产生子与探测器模拟的界面BESIII模仿BESIILUARLW/JETSET