基于蒙特卡洛模拟的真空康普顿探测器性能研究与优化

基于蒙特卡洛模拟的真空康普顿探测器性能研究与优化一、引言1.1研究背景与意义在现代物理学和核技术应用领域，射线探测技术一直占据着至关重要的地位，从基础科学研究到工业检测、医学诊断与治疗、安全检查等诸多实际应用场景，都离不开射线探测技术的支持。例如，在医学影像诊断中，X射线、γ射线等被用于生成人体内部结构的图像，帮助医生检测疾病；在工业无损检测中，射线用于检测材料内部的缺陷，确保产品质量。而康普顿探测器作为射线探测的重要工具之一，凭借其独特的工作原理和性能优势，在强流脉冲γ辐射场探测等方面发挥着关键作用。真空康普顿探测器（VacuumComptonDetector，VCD）是基于康普顿效应设计的一种探测器。康普顿效应指的是当光子与物质中的电子相互作用时，光子将一部分能量传递给电子，自身能量降低且散射方向改变。真空康普顿探测器利用这一效应，通过探测散射后的光子或电子来获取入射射线的相关信息，如能量、方向等。相比于其他类型的探测器，真空康普顿探测器具有时间响应快的特点，能够快速捕捉射线信号，适用于探测短脉冲辐射；同时，它属于电荷收集型探测器，工作在电流模式下，在强流脉冲γ辐射场探测中表现出色，自20世纪60年代起，就在美国强流γ射线测量中得到了应用。然而，随着科学技术的不断发展和应用需求的日益增长，对真空康普顿探测器的性能提出了更高的要求。在宽能谱的脉冲γ辐射场强度测量中，为了准确获取更为精确的强度信息，需要探测器在宽能谱范围内具有能量响应平坦的特点。虽然基于康普顿效应的各种探测器，如VCD145、VCD-I～VCD-Ⅳ、ICI探测器、复合金属发射极真空康普顿探测器等的研制使γ测量技术取得了重要进展，但在能量响应方面，由于以往的应用需求不同，多数文献仅给出了灵敏度能量响应特性，并未开展以实现能量响应平坦特性为目的的深入研究。2005年出现的能量响应相对平坦的介质康普顿探测器，设计了前后金属盖、散射体和吸收体厚度可调结构，在0.5-5.0MeVγ能区内灵敏度变化波动小于17%，但其能量响应平坦性评价结果不是按照本征灵敏度对比给出，结果可信度有待进一步验证。2008年发展的基于散射电子和闪烁薄膜的探测技术，设计的背散射法探测器系统，在0.4-5.0MeV能量区间，灵敏度变化小于15%，虽然是探测器能量响应方面研究取得的重要进展，但适用的能谱范围仍不够宽。因此，开展宽能谱范围的能量响应平坦的真空康普顿探测器的研究具有重要的现实意义。蒙特卡洛模拟方法作为一种基于概率统计的数值计算方法，在解决复杂物理问题中具有独特的优势。它通过对大量随机事件的模拟来近似求解问题，能够处理包含多种相互作用和复杂几何结构的系统。在真空康普顿探测器的研究中，蒙特卡洛模拟可以用于建立探测器的3D动态模型，对探测器射线的能量及其转换、信号检测等各方面的参数进行模拟计算。同时，还可以对射线入射角度、位置等参数进行参数化，以便于后续的参数优化。通过蒙特卡洛模拟，能够在实际制造探测器之前，深入了解探测器的性能，预测不同参数条件下探测器的响应，为探测器的优化设计提供理论依据，从而节省实验成本和时间，提高研究效率。因此，运用蒙特卡洛模拟对真空康普顿探测器进行研究，对于提升探测器性能、拓展其应用领域具有重要的科学意义和实用价值。1.2国内外研究现状真空康普顿探测器的研究一直是国内外核探测领域的重要课题。国外对真空康普顿探测器的研究起步较早，自20世纪60年代起，美国就将其应用于强流γ射线测量中。早期的研究主要集中在探测器的基本原理验证和性能初步测试方面。随着技术的不断发展，研究人员开始关注如何提高探测器的性能，如能量分辨率、探测效率等。在蒙特卡洛模拟应用于真空康普顿探测器研究方面，国外学者开展了大量工作。1989年，HsuH.-H.等人进行了蒙特卡洛模拟，计算真空康普顿探测器（VCD）对单能伽马射线和电子的静态灵敏度，并且使用^{60}Co光谱作为输入的计算结果与测量值吻合良好，同时还计算了由单能电子束产生的轫致辐射光谱的探测器灵敏度，并与实验数据进行比较，得出该方法可准确可靠地计算VCD各种灵敏度的结论。这一研究为后续利用蒙特卡洛模拟优化探测器性能奠定了基础。此后，国外研究团队不断利用蒙特卡洛模拟对探测器的结构、材料等参数进行优化，以提升探测器在不同应用场景下的性能。例如，在一些核物理实验中，通过蒙特卡洛模拟精确设计探测器，使其能够更好地适应复杂的辐射环境，准确测量射线参数。国内对于真空康普顿探测器的研究虽然起步相对较晚，但近年来发展迅速。研究人员在借鉴国外先进技术的基础上，结合国内实际需求，开展了一系列有针对性的研究工作。在探测器的设计方面，尝试采用新的结构和材料组合，以实现更好的性能。2012年，西北核技术研究所的翁秀峰等人以开展能量响应平坦的真空康普顿探测器为目的，采用厚薄材料叠加补偿设计思想，使用蒙特卡罗方法，对真空康普顿探测器的灵敏度能量响应进行了优化设计。优化后的探测器系统在0.4-7.0MeVγ能区内，探测器本征灵敏度极值变化小于10.7％，能量响应平坦性优于当时各种常见的γ探测器。这一成果展示了国内在利用蒙特卡洛模拟改进真空康普顿探测器性能方面的重要进展。此外，国内研究人员还利用蒙特卡洛模拟研究探测器在不同辐射场下的响应特性，为探测器在工业检测、医学放疗剂量监测等实际应用提供理论支持。在工业无损检测中，通过模拟不同缺陷情况下射线与探测器的相互作用，优化探测器的设置，提高对缺陷的检测能力；在医学放疗剂量监测中，模拟射线在人体组织和探测器中的传输过程，确保探测器能够准确测量放疗剂量，保障患者的治疗安全。然而，目前国内外的研究仍存在一些不足之处。在能量响应平坦性方面，虽然已经取得了一定的进展，但仍有提升空间，现有探测器在宽能谱范围内的能量响应平坦度还不能完全满足一些高精度测量的需求。在探测器的小型化和便携化研究方面，虽然蒙特卡洛模拟可以在理论上进行优化设计，但在实际制造过程中，还面临着工艺和成本等方面的挑战。未来，需要进一步深入研究蒙特卡洛模拟方法，结合新型材料和制造工艺，不断改进真空康普顿探测器的性能，以满足日益增长的科学研究和实际应用需求。1.3研究内容与方法本文围绕真空康普顿探测器的蒙特卡洛模拟展开深入研究，主要研究内容包括以下几个方面：真空康普顿探测器原理深入剖析：全面系统地阐述真空康普顿探测器的工作原理，包括康普顿效应在探测器中的具体作用机制，光子与探测器内物质相互作用的过程，如散射、吸收等，以及这些相互作用如何产生可检测的信号。详细分析探测器的结构组成，各部件的功能及对探测器性能的影响，为后续的蒙特卡洛模拟和探测器优化设计奠定坚实的理论基础。蒙特卡洛模型构建：运用蒙特卡洛模拟方法，精心构建真空康普顿探测器的3D动态模型。在模型构建过程中，充分考虑探测器的实际几何结构、材料特性以及射线与物质相互作用的各种物理过程，如康普顿散射、光电效应等。对探测器射线的能量及其转换、信号检测等各方面的参数进行精确模拟计算，确保模型能够准确反映探测器的实际工作情况。同时，对射线入射角度、位置等参数进行参数化处理，为后续的参数优化提供便利条件。参数优化及结果分析：依据建立的蒙特卡洛模型，对探测器的相关参数进行系统的调整和优化。研究电场强度、放射性样品位置、探测器各部件的材料和尺寸等参数变化对探测器性能的影响规律，通过比较不同参数设置下探测器的信号强度、能量分辨率、探测效率等性能参数，确定探测器的最优参数组合。深入分析模拟结果，挖掘数据背后的物理意义，为探测器的实际设计和改进提供有价值的参考依据。模拟结果可视化展示：借助三维动态展示、图形自动生成等先进的技术手段，将模拟结果以直观、可视化的形式呈现出来。生成各种参数变化对结果影响的散点图、直方图、容积图等，使复杂的模拟数据更加易于理解和分析。通过可视化展示，能够更清晰地观察探测器在不同条件下的性能表现，及时发现问题并进行优化调整。在研究方法上，本文主要采用以下几种方法：理论研究法：通过查阅大量的国内外相关文献资料，深入学习和研究真空康普顿探测器的工作原理、蒙特卡洛模拟方法的基本理论以及相关的物理学知识。对康普顿效应、射线与物质相互作用的理论进行系统梳理，为后续的研究提供坚实的理论支撑。数值模拟法：运用蒙特卡洛模拟软件，如Geant4等，进行真空康普顿探测器的模拟研究。在模拟过程中，严格按照探测器的实际结构和物理过程进行建模，设置合理的模拟参数，确保模拟结果的准确性和可靠性。通过多次模拟计算，获取不同条件下探测器的性能数据，为参数优化和结果分析提供数据基础。对比分析法：在参数优化和结果分析阶段，采用对比分析的方法。对比不同参数设置下探测器的性能差异，找出影响探测器性能的关键因素。同时，将模拟结果与已有的实验数据或其他研究成果进行对比，验证模拟方法和模型的正确性，进一步优化模拟结果。可视化方法：利用专业的可视化软件，将模拟结果转化为直观的图形和图像，如探测器内部射线的传输轨迹、能量沉积分布、信号强度分布等。通过可视化展示，能够更直观地了解探测器的工作过程和性能特点，为研究人员提供更清晰的研究思路和方向。本文通过以上研究内容和方法，旨在深入探究真空康普顿探测器的性能特点，利用蒙特卡洛模拟实现探测器的优化设计，为提高探测器的性能和拓展其应用领域提供理论支持和技术指导。二、真空康普顿探测器原理2.1康普顿效应基础康普顿效应是真空康普顿探测器工作的核心物理基础，它深刻地揭示了光子与电子之间的相互作用机制。1923年，美国物理学家康普顿在研究X射线通过实物物质发生散射的实验时，首次发现了这一重要的物理现象。实验中，康普顿观察到散射光中除了存在与原波长相同的成分外，还出现了波长变长的成分，且波长的增量随散射角的不同而变化，这一现象无法用经典电磁理论进行解释。从微观角度来看，康普顿效应的过程可描述为：当高能光子（如X射线或γ射线光子）与物质中的自由电子或准自由电子相互作用时，光子将其部分能量和动量传递给电子，自身能量降低，波长变长，同时电子获得动能并向其他方向移动，这一过程类似于两个弹性小球的碰撞。在这个过程中，能量和动量守恒定律起着关键作用。从能量守恒的角度出发，设入射光子的能量为E_{γ0}，散射光子的能量为E_{γ}，电子的静止能量为m_0c^2，碰撞后电子的动能为E_{k}，则能量守恒方程可表示为：E_{γ0}=E_{γ}+E_{k}+m_0c^2由于光子能量E=hν=\frac{hc}{\lambda}（其中h为普朗克常数，c为光速，ν为频率，\lambda为波长），将其代入上式可得：\frac{hc}{\lambda_{0}}=\frac{hc}{\lambda}+E_{k}+m_0c^2从动量守恒的角度，建立直角坐标系，将静止的电子定位在原点。假设入射光子具有初始动量\vec{p}_{γ0}，碰撞后光子的动量为\vec{p}_{γ}，并以相对于x轴的角度\theta的方向离开；而电子的动量是\vec{p}_{e}，并以相对于x轴的角度\varphi的方向离开。根据动量守恒定律，在x和y方向上的分量分别满足：p_{γ0x}=p_{γx}+p_{ex}p_{γ0y}=p_{γy}+p_{ey}将上述动量分量式平方并相加，经过一系列推导，并结合能量守恒方程以及爱因斯坦的能量-动量关系E^2=p^2c^2+m_0^2c^4（对于光子，m_0=0，则E=pc），最终可以得到著名的康普顿散射公式，用于描述散射前后光子波长的关系：\Delta\lambda=\lambda-\lambda_{0}=\frac{h}{m_0c}(1-\cos\theta)其中，\Delta\lambda为波长的改变量，即康普顿位移；m_0为电子的静止质量；\theta为散射角，即散射光子与入射光子方向之间的夹角。该公式清晰地表明，波长的改变量\Delta\lambda只与散射角\theta有关，与入射光子的波长\lambda_{0}无关。当散射角\theta增大时，\cos\theta减小，(1-\cos\theta)增大，从而\Delta\lambda增大，即散射光的波长随散射角的增大而变长。例如，当\theta=0时，\cos\theta=1，\Delta\lambda=0，散射光波长与入射光波长相同；当\theta=180^{\circ}时，\cos\theta=-1，\Delta\lambda达到最大值，此时散射光波长最长。康普顿效应的发现具有极其重要的科学意义，它不仅为光的粒子性提供了有力的实验证据，进一步证实了爱因斯坦的光子假说，而且在量子力学的发展历程中起到了关键的推动作用。通过康普顿效应，人们更加深入地认识到光子不仅具有能量，还具有动量，光子与电子的相互作用是量子化的过程，这为后续研究微观世界的物理规律奠定了坚实的基础。2.2真空康普顿探测器工作机制真空康普顿探测器的工作机制基于康普顿效应，其核心在于利用射线与探测器内物质相互作用产生的物理现象来实现射线探测。当准直后的窄束伽马射线入射到真空康普顿探测器中时，会经历一系列复杂的物理过程。首先，伽马射线光子与探测器内的金属发射极发生相互作用，这种相互作用主要包括光电效应、康普顿效应和电子对效应等。在这些效应中，康普顿效应起着关键作用。当高能伽马射线光子与金属发射极中的自由电子或准自由电子碰撞时，光子将一部分能量和动量传递给电子，自身能量降低，波长变长，成为散射光子，而电子则获得动能成为反冲电子。这一过程满足能量守恒和动量守恒定律，其能量和动量的变化关系遵循前文所述的康普顿效应相关公式。以能量为E_{γ0}的伽马射线光子为例，假设其与静止的电子发生康普顿散射，散射后光子的能量变为E_{γ}，电子获得动能E_{k}。根据能量守恒定律E_{γ0}=E_{γ}+E_{k}+m_0c^2，其中m_0c^2为电子的静止能量。在动量守恒方面，建立直角坐标系，设入射光子动量为\vec{p}_{γ0}，散射光子动量为\vec{p}_{γ}，反冲电子动量为\vec{p}_{e}，在x和y方向上分别满足p_{γ0x}=p_{γx}+p_{ex}和p_{γ0y}=p_{γy}+p_{ey}。这些产生的反冲电子从金属发射极表面发射出来，进而引起金属发射极的电势变化。由于电子的发射，金属发射极失去部分负电荷，其电势会发生改变。这种电势变化与入射伽马射线的强度、能量等信息密切相关。例如，入射伽马射线强度越高，产生的反冲电子数量可能越多，金属发射极的电势变化也就越明显；而入射伽马射线能量的不同，会导致反冲电子能量分布的差异，从而也会影响金属发射极电势变化的具体情况。为了测量这种电势变化，探测器通常配备有信号检测和处理电路。信号检测部分会将金属发射极的电势变化转化为电信号，例如电压信号或电流信号。这些电信号随后被传输到信号处理电路中，在信号处理电路中，电信号会经过放大、滤波等一系列处理步骤。放大过程可以增强信号的强度，使其更容易被后续的分析和测量；滤波则可以去除信号中的噪声和干扰，提高信号的质量。经过处理后的电信号，最终被转换为数字信号，以便于数据采集系统进行采集和分析。通过对采集到的数据进行处理和分析，就可以获取入射伽马射线的相关信息，如射线的强度、能量分布、时间特性等。在整个工作过程中，探测器的各个部分协同工作，从射线入射到电信号输出，每一个环节都对探测器的性能有着重要影响。例如，金属发射极的材料和结构会影响射线与电子相互作用的概率和效率，进而影响探测器的灵敏度；信号检测和处理电路的性能则直接决定了能否准确、快速地获取和处理电信号，影响探测器的时间响应和测量精度。2.3探测器结构组成真空康普顿探测器主要由真空容器、电场装置、探测器元件和信号处理电路等部分组成，各部分相互协作，共同实现对射线的有效探测。真空容器是探测器的基础结构，通常采用金属材料制成，如不锈钢等。其主要功能是为探测器内部的物理过程提供一个真空环境，减少气体分子对射线和电子的散射与吸收，从而提高探测器的探测效率和精度。例如，在真空中，射线与物质相互作用产生的电子能够更自由地运动，减少了因气体碰撞而导致的能量损失和散射，使得探测器能够更准确地捕捉到电子信号。同时，真空容器还起到电磁屏蔽的作用，防止外界电磁干扰影响探测器内部的信号传输和处理。电场装置在探测器中起着关键作用，一般由电源和电极结构组成。电源为电极提供稳定的电压，从而在探测器内部建立起特定强度和方向的电场。电场的主要作用是对射线与物质相互作用产生的电子进行加速和引导，使其能够按照预定的路径运动到探测器元件上，便于信号的检测和收集。例如，在一些真空康普顿探测器中，通过在金属发射极和收集电极之间施加一定的电压，形成电场，使得从金属发射极发射出的反冲电子在电场力的作用下快速向收集电极运动，提高了信号的产生效率和响应速度。探测器元件是直接与射线相互作用并产生信号的关键部分，主要包括金属发射极等。金属发射极通常采用高原子序数的金属材料，如金（Au）、铂（Pt）等。当射线入射到金属发射极时，会与其中的电子发生康普顿效应等相互作用，产生反冲电子。这些反冲电子从金属发射极表面发射出来，成为探测器检测的信号源。金属发射极的材料特性和几何形状对探测器的性能有着重要影响，不同的材料和形状会导致射线与电子相互作用的概率和产生的反冲电子能量分布不同，进而影响探测器的灵敏度和能量分辨率。信号处理电路用于对探测器元件产生的电信号进行处理和分析，通常包括前置放大器、滤波器、主放大器、模数转换器（ADC）等部分。前置放大器的作用是在信号传输的初始阶段，对探测器输出的微弱电信号进行初步放大，提高信号的强度，以便后续的处理；滤波器则用于去除信号中的噪声和干扰，通过设置合适的滤波参数，如截止频率、带宽等，可以有效滤除与信号无关的高频或低频噪声，提高信号的质量；主放大器进一步对信号进行放大，使其达到适合模数转换的电平范围；模数转换器将模拟信号转换为数字信号，便于计算机进行数据采集和处理。例如，在实际应用中，经过前置放大器和滤波器处理后的信号，进入主放大器进行放大，然后由模数转换器将其转换为数字信号，传输到计算机中进行分析和存储，通过对数字信号的处理和分析，可以获取射线的能量、强度、时间等信息。三、蒙特卡洛模拟方法3.1基本原理蒙特卡洛模拟方法作为一种强大的数值计算手段，其核心思想建立在随机抽样和统计计算的基础之上，旨在通过模拟大量的随机事件来近似求解复杂系统或过程中的各种问题。该方法的理论基石源于概率论中的大数定律，即随着试验次数的不断增加，事件发生的频率会逐渐趋近于其理论概率。以估算不规则图形的面积为例，假设存在一个形状不规则的图形，难以通过常规的几何公式直接计算其面积。我们可以将这个不规则图形放置在一个已知面积的正方形区域内，通过在正方形区域内随机生成大量的点，统计落在不规则图形内的点的数量与总点数的比例。根据大数定律，当随机生成的点足够多时，这个比例会趋近于不规则图形面积与正方形面积的比值，从而可以通过已知的正方形面积来估算不规则图形的面积。在数学计算领域，蒙特卡洛模拟可用于求解复杂的积分问题。例如，对于一个难以通过解析方法求解的积分\int_{a}^{b}f(x)dx，可以通过在区间[a,b]上随机生成大量的点x_i，计算函数值f(x_i)，然后根据这些随机样本计算积分的近似值，即\hat{I}=\frac{1}{N}\sum_{i=1}^{N}f(x_i)，其中N为随机样本的数量，\hat{I}为积分的蒙特卡洛估计。在物理领域，蒙特卡洛模拟被广泛应用于模拟粒子在物质中的输运过程。例如，在研究中子在核反应堆中的运动时，由于中子与反应堆内物质的相互作用十分复杂，包含散射、吸收等多种过程，难以通过精确的数学模型进行描述。蒙特卡洛模拟通过随机抽样的方式，模拟中子的初始位置、速度和方向等参数，以及中子与物质相互作用的类型和结果。通过大量的模拟，统计中子在不同位置的分布、能量变化等信息，从而深入了解核反应堆内的物理过程，为反应堆的设计和优化提供重要依据。在金融领域，蒙特卡洛模拟可用于评估投资组合的风险和收益。假设投资者持有多种不同的资产，每种资产的收益率都受到多种不确定因素的影响，如市场波动、宏观经济状况等。通过蒙特卡洛模拟，可以随机生成各种可能的市场情景，计算在不同情景下投资组合的收益率，进而统计投资组合收益率的概率分布，评估投资组合面临的风险，为投资决策提供参考。蒙特卡洛模拟方法的优势在于其通用性强，能够处理各种复杂的、难以用传统解析方法解决的问题，不受问题复杂度和维度的限制；同时，该方法简单易实现，无需复杂的数学推导和解析公式，只需定义问题的概率模型和随机抽样方式即可进行模拟。然而，蒙特卡洛模拟也存在一些局限性，其结果的准确性依赖于大量的随机样本，计算量较大，且由于结果具有随机性，每次模拟的结果可能存在一定的差异。3.2实现步骤蒙特卡洛模拟在真空康普顿探测器研究中的实现是一个系统且严谨的过程，主要涵盖问题定义、随机数生成、模拟过程设计和结果统计分析等关键步骤，每个步骤都紧密相连，对准确模拟探测器性能起着不可或缺的作用。在问题定义阶段，需精准确定研究目标，即深入了解真空康普顿探测器在不同条件下的性能表现，包括但不限于其对不同能量射线的探测效率、能量分辨率以及信号响应特性等。同时，明确模拟过程中涉及的各种物理量和参数，例如射线的能量范围、入射角度分布，探测器各部件的材料特性（如原子序数、密度等）、几何尺寸（如金属发射极的厚度、真空容器的形状和尺寸等），以及电场强度和方向等。这些参数将直接影响探测器与射线相互作用的物理过程，进而决定探测器的性能。以射线能量范围为例，不同能量的射线与探测器内物质相互作用的方式和概率不同，低能射线可能主要发生光电效应，而高能射线则更易发生康普顿效应和电子对效应，明确射线能量范围有助于准确模拟探测器对不同能量射线的响应。随机数生成是蒙特卡洛模拟的基础环节，其质量直接关系到模拟结果的准确性和可靠性。在模拟过程中，大量的物理过程具有随机性，如射线与物质相互作用的位置、散射角度、能量损失等，这些都需要通过随机数来模拟。通常使用伪随机数生成器来产生满足特定概率分布的随机数，常见的概率分布包括均匀分布、正态分布等。在模拟射线与物质相互作用的位置时，可根据探测器的几何形状和尺寸，利用均匀分布的随机数来确定相互作用点在探测器内的坐标；而在模拟散射角度时，可能需要根据康普顿散射的理论概率分布，通过特定的算法将均匀分布的随机数转换为符合散射角度分布的随机数。为了确保随机数的独立性和均匀性，需要对随机数生成器进行严格的测试和验证，采用一些标准的测试方法，如频率测试、游程测试等，以保证生成的随机数能够准确模拟物理过程的随机性。模拟过程设计是蒙特卡洛模拟的核心部分，需要建立详细的物理模型来描述射线与真空康普顿探测器相互作用的全过程。这个物理模型应全面考虑各种物理过程，如康普顿散射、光电效应、电子对效应等，以及探测器内的电场分布、电子输运等因素。以康普顿散射过程为例，根据康普顿散射公式，结合随机数生成的散射角度和能量损失，计算散射后光子和电子的能量和方向。在模拟电子输运时，考虑电子在电场中的加速、与气体分子或探测器材料原子的碰撞等过程，确定电子的运动轨迹和最终到达探测器元件的位置。同时，利用计算机编程实现模拟过程，选择合适的编程语言和模拟软件，如C++、Python结合Geant4等专业模拟工具。在编程过程中，按照物理模型的逻辑，依次实现射线源的定义、粒子的发射、相互作用的模拟、信号的产生和收集等功能模块，确保模拟过程的准确性和高效性。结果统计分析是蒙特卡洛模拟的重要环节，通过对大量模拟结果的统计分析，获取有价值的物理信息和探测器性能参数。统计分析的内容包括但不限于探测效率、能量分辨率、信号强度分布等。以探测效率为例，统计在一定模拟次数内，探测器成功探测到射线的次数与总发射射线次数的比值，得到探测效率；对于能量分辨率，分析探测器输出信号的能量分布，计算能量分辨率的相关指标，如半高宽等。为了更直观地展示模拟结果，采用多种数据分析和可视化方法，生成散点图展示射线能量与探测效率的关系，绘制直方图呈现信号强度的分布情况，利用容积图展示探测器内不同位置的能量沉积分布等。通过这些可视化手段，能够更清晰地观察探测器性能随参数的变化规律，为探测器的优化设计提供有力依据。同时，对模拟结果进行不确定性分析，评估由于随机数生成、物理模型简化等因素导致的结果不确定性，确定模拟结果的置信区间，提高模拟结果的可靠性和可信度。3.3优势与局限性蒙特卡洛模拟在真空康普顿探测器研究中展现出多方面的显著优势，为深入探究探测器性能提供了有力支持。其强大的处理复杂问题能力是一大突出优势，在真空康普顿探测器中，射线与探测器内物质的相互作用涉及多种复杂物理过程，如康普顿散射、光电效应、电子对效应等，同时探测器的结构和材料特性也较为复杂。蒙特卡洛模拟能够全面考虑这些因素，通过建立详细的物理模型，对各种相互作用过程进行精确模拟，从而准确预测探测器在不同条件下的性能。例如，在模拟射线在探测器内的传输过程时，蒙特卡洛模拟可以根据探测器各部件的材料和几何形状，准确计算射线与物质相互作用的位置、散射角度和能量损失等参数，这是传统解析方法难以实现的。蒙特卡洛模拟无需复杂的数学推导，这使得研究人员能够更专注于探测器的物理过程和性能优化。在传统的探测器研究中，对于复杂的物理问题，往往需要进行大量的数学推导和近似处理，这不仅耗时费力，还可能引入误差。而蒙特卡洛模拟通过随机抽样和统计计算的方式，直接对物理过程进行模拟，避免了复杂的数学推导过程，降低了研究的难度和工作量。例如，在计算探测器对不同能量射线的探测效率时，蒙特卡洛模拟只需根据设定的物理模型和随机数生成规则，进行大量的模拟计算，然后通过统计分析即可得到探测效率的结果，无需进行复杂的数学积分和推导。此外，蒙特卡洛模拟具有高度的灵活性，能够方便地对探测器的各种参数进行调整和优化。研究人员可以根据实际需求，在模拟过程中改变射线的能量、入射角度、探测器各部件的材料和尺寸等参数，快速评估不同参数组合对探测器性能的影响，从而找到最优的设计方案。这种灵活性为探测器的设计和改进提供了极大的便利，能够大大缩短探测器的研发周期，降低研发成本。例如，在研究探测器的能量分辨率时，可以通过蒙特卡洛模拟，快速改变探测器的电场强度、电极间距等参数，观察这些参数变化对能量分辨率的影响，从而确定最佳的参数设置。然而，蒙特卡洛模拟也存在一些局限性。计算效率较低是其主要不足之一，为了获得较为准确的结果，蒙特卡洛模拟通常需要进行大量的随机抽样和模拟计算，这导致计算量巨大，计算时间较长。尤其是在模拟复杂的探测器系统或进行高精度的性能分析时，计算时间可能会达到数小时甚至数天，这对于一些需要快速得到结果的研究场景来说，是一个较大的限制。例如，在对探测器进行全面的性能评估时，需要模拟不同能量、不同入射角度的大量射线与探测器的相互作用，计算量非常庞大，会耗费大量的计算资源和时间。蒙特卡洛模拟的结果存在一定的随机性，由于模拟过程依赖于随机数生成，每次模拟的结果可能会略有不同。虽然可以通过增加模拟次数来减小这种随机性的影响，但无法完全消除。这就需要研究人员在分析模拟结果时，考虑到结果的不确定性，进行多次模拟并对结果进行统计分析，以获得更可靠的结论。例如，在模拟探测器的探测效率时，每次模拟得到的探测效率可能会在一定范围内波动，需要进行多次模拟，取平均值作为最终结果，并评估结果的误差范围。四、真空康普顿探测器蒙特卡洛模型建立4.1模型假设与简化在构建真空康普顿探测器的蒙特卡洛模型时，为了使模型既能够精确反映探测器的关键特性，又能在合理的计算资源和时间范围内进行有效分析，需要对探测器的结构和物理过程做出一系列合理的假设与简化。从探测器结构方面来看，首先假设真空容器为理想的封闭结构，忽略其可能存在的微小缝隙或缺陷。实际上，在制造工艺中，真空容器难以做到绝对完美，可能存在一些细微的缝隙，但这些缝隙对探测器内部的物理过程影响极小，在模型中可忽略不计。这样的假设能够简化对真空容器内部电场分布和粒子输运过程的模拟，使得模拟过程更加高效且易于实现。对于电场装置，假定其产生的电场是均匀且稳定的，不考虑电场在边缘处可能出现的畸变以及随时间的微小波动。在实际情况中，由于电极的形状、材料以及电源的稳定性等因素，电场在边缘区域可能会出现一定程度的畸变，并且在长时间运行过程中可能会有微小的波动。然而，通过合理设计电场装置和选择高质量的电源，可以使这些影响控制在较小范围内，在模型中进行简化处理不会对模拟结果产生显著影响。例如，在一些高精度的探测器中，通过优化电极形状和采用高精度的稳压电源，电场的均匀性和稳定性可以得到很好的保证。探测器元件中的金属发射极，假设其材料均匀且表面光滑，不考虑材料内部可能存在的杂质以及表面的微观粗糙度。尽管在实际的金属发射极中，不可避免地会存在一定的杂质和微观粗糙度，但这些因素对射线与电子相互作用过程的影响相对较小。在模型中进行这样的简化，能够减少模拟过程中的参数数量，提高计算效率，同时也不会对探测器的主要性能模拟产生较大偏差。在物理过程方面，主要考虑康普顿散射、光电效应和电子对效应这三种射线与物质相互作用的主要过程，而忽略其他一些相对次要的相互作用，如瑞利散射等。这是因为在真空康普顿探测器所涉及的能量范围内，康普顿散射、光电效应和电子对效应是主要的相互作用方式，瑞利散射等过程发生的概率极低，对探测器的整体性能影响可以忽略不计。以能量为1MeV的γ射线为例，在与常见的金属发射极材料相互作用时，康普顿散射和光电效应是主要的作用过程，瑞利散射的概率非常小，在模型中不考虑瑞利散射不会影响对探测器性能的准确模拟。对于电子在探测器内的输运过程，假设电子在真空中的运动是自由的，不考虑电子与残留气体分子的碰撞以及电子之间的相互作用。虽然在实际的真空环境中，仍然会存在极少量的残留气体分子，电子在运动过程中可能会与这些分子发生碰撞，同时电子之间也可能存在相互作用。但是，通过良好的真空系统和高真空度的维持，可以使残留气体分子的数量极少，电子与残留气体分子的碰撞概率非常低，电子之间的相互作用也相对较弱。在模型中忽略这些因素，能够简化电子输运过程的模拟，提高计算效率，并且对模拟结果的准确性影响较小。通过以上对探测器结构和物理过程的假设与简化，建立的蒙特卡洛模型既能抓住真空康普顿探测器的核心特性，又能在可接受的计算资源和时间成本下进行高效模拟，为后续深入研究探测器性能和参数优化提供了坚实的基础。4.2参数设定在真空康普顿探测器的蒙特卡洛模型中，关键参数的设定对于准确模拟探测器性能至关重要，这些参数涵盖射线相关参数以及探测器材料和结构参数等多个方面，它们的取值依据充分考虑了探测器的实际应用场景、相关物理理论以及前人的研究成果。射线能量是影响探测器性能的关键因素之一，根据真空康普顿探测器在实际应用中常涉及的射线能量范围，如在工业无损检测、医学放疗剂量监测等领域，射线能量通常在keV到MeV量级。为了全面研究探测器在不同能量下的性能，设定射线能量范围为0.1-10MeV，在此范围内选取多个离散的能量点进行模拟计算，如0.1MeV、0.5MeV、1MeV、2MeV、5MeV、10MeV等。这样的能量范围设定能够覆盖探测器在大多数实际应用中的射线能量情况，同时选取多个能量点进行模拟，可以更详细地分析探测器性能随射线能量的变化规律。例如，在医学放疗中，常用的γ射线能量约为1-5MeV，工业无损检测中使用的射线能量也多在这个范围内，通过模拟这些能量点下探测器的性能，能够为实际应用提供有针对性的参考。射线入射角度对探测器的响应也有显著影响。在实际情况中，射线可能以不同的角度入射到探测器上。为了研究这种影响，设定射线入射角度范围为0°-90°，以10°为间隔进行模拟。选择这样的角度范围和间隔，一方面是因为在实际应用中，射线入射角度通常不会超过90°，另一方面，10°的间隔能够在保证计算精度的前提下，减少计算量。例如，当射线垂直入射（0°）时，探测器的响应与斜入射时会有所不同，通过模拟不同入射角度下探测器的信号强度、能量分辨率等性能参数，可以了解探测器对不同入射角度射线的适应性。探测器材料参数方面，金属发射极作为探测器中与射线直接相互作用的关键部件，其材料的选择对探测器性能起着决定性作用。常见的金属发射极材料有金（Au）、铂（Pt）等，这些材料具有高原子序数和良好的导电性，能够提高射线与电子相互作用的概率，增强探测器的灵敏度。在模型中，选取金作为金属发射极材料，其密度为19.32g/cm³，原子序数为79。这些参数是根据金的实际物理性质确定的，准确的材料参数设定能够保证模型对探测器物理过程的模拟准确性。例如，金的高原子序数使得射线与其中电子发生康普顿效应和光电效应的概率相对较高，从而能够产生更多的可检测信号。真空容器材料通常采用不锈钢，其密度约为7.9g/cm³，原子序数根据不锈钢的主要成分确定。不锈钢具有良好的机械性能和真空密封性，能够为探测器内部提供稳定的真空环境。在模型中，准确设定不锈钢的材料参数，有助于模拟射线在真空容器内的传输过程以及探测器整体的性能。电场强度是探测器中的重要参数，它影响着电子在探测器内的运动轨迹和能量分布。根据探测器的设计要求和实际应用场景，设定电场强度范围为100-1000V/cm。在这个范围内调整电场强度，研究其对探测器性能的影响。例如，增大电场强度可以加速电子的运动，提高探测器的时间响应速度，但同时也可能增加电子与探测器部件的碰撞概率，影响探测器的能量分辨率。通过模拟不同电场强度下探测器的性能，可以找到最佳的电场强度设置，以满足探测器在不同应用场景下的需求。4.3模拟程序实现在实现真空康普顿探测器的蒙特卡洛模拟时，选用了功能强大的Geant4软件结合C++编程语言。Geant4是一款基于C++开发的开源蒙特卡洛模拟工具包，广泛应用于粒子物理、医学物理、材料科学等领域的粒子输运模拟，具有丰富的物理模型库和灵活的几何建模能力，能够准确模拟射线与物质的相互作用过程。而C++语言以其高效的执行效率和对底层硬件资源的良好控制能力，为模拟程序的快速运行和复杂算法的实现提供了有力支持。模拟程序首先定义射线源，利用Geant4提供的相关类，设置射线的类型（如γ射线）、能量范围（根据前文设定的0.1-10MeV）以及入射角度范围（0°-90°）。例如，通过以下代码片段实现射线源的定义：G4ParticleGun*particleGun=newG4ParticleGun(1);G4ParticleTable*particleTable=G4ParticleTable::GetParticleTable();G4ParticleDefinition*particle=particleTable->FindParticle("gamma");particleGun->SetParticleDefinition(particle);G4doubleenergyMin=0.1*CLHEP::MeV;G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4ParticleTable*particleTable=G4ParticleTable::GetParticleTable();G4ParticleDefinition*particle=particleTable->FindParticle("gamma");particleGun->SetParticleDefinition(particle);G4doubleenergyMin=0.1*CLHEP::MeV;G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4ParticleDefinition*particle=particleTable->FindParticle("gamma");particleGun->SetParticleDefinition(particle);G4doubleenergyMin=0.1*CLHEP::MeV;G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);particleGun->SetParticleDefinition(particle);G4doubleenergyMin=0.1*CLHEP::MeV;G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4doubleenergyMin=0.1*CLHEP::MeV;G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4doubleenergyMax=10*CLHEP::MeV;G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4doubleangleMin=0*CLHEP::degree;G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4doubleangleMax=90*CLHEP::degree;//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);//设置能量和角度的随机分布G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4RandomDirection*direction=newG4RandomDirection();particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);particleGun->SetParticleMomentumDirection(direction);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);G4RandomEnergy*energy=newG4RandomEnergy(energyMin,energyMax);particleGun->SetParticleEnergy(energy);particleGun->SetParticleEnergy(energy);上述代码中，首先创建了一个粒子枪G4ParticleGun，并指定每次发射一个粒子。然后获取粒子表G4ParticleTable，从中找到γ射线粒子定义G4ParticleDefinition并设置到粒子枪中。接着定义了射线的能量范围和入射角度范围，并通过G4RandomEnergy和G4RandomDirection类实现能量和角度的随机分布设置。在定义探测器几何结构时，充分利用Geant4的几何建模功能。对于真空容器，使用相应的几何类创建其形状（如圆柱体），并设置其尺寸参数（半径、高度等）以及材料（不锈钢）。对于金属发射极，同样创建其几何形状（如平板），设置尺寸和材料（金），并将其放置在真空容器内合适的位置。例如：G4Material*stainlessSteel=newG4Material("StainlessSteel",7.9*g/cm3,2);//添加不锈钢的主要成分及比例G4NistManager*nist=G4NistManager::Instance();stainlessSteel->AddElement(nist->FindOrBuildElement("Fe"),0.7);stainlessSteel->AddElement(nist->FindOrBuildElement("Cr"),0.18);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);//添加不锈钢的主要成分及比例G4NistManager*nist=G4NistManager::Instance();stainlessSteel->AddElement(nist->FindOrBuildElement("Fe"),0.7);stainlessSteel->AddElement(nist->FindOrBuildElement("Cr"),0.18);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4NistManager*nist=G4NistManager::Instance();stainlessSteel->AddElement(nist->FindOrBuildElement("Fe"),0.7);stainlessSteel->AddElement(nist->FindOrBuildElement("Cr"),0.18);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);stainlessSteel->AddElement(nist->FindOrBuildElement("Fe"),0.7);stainlessSteel->AddElement(nist->FindOrBuildElement("Cr"),0.18);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);stainlessSteel->AddElement(nist->FindOrBuildElement("Cr"),0.18);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4Box*vacuumVesselBox=newG4Box("VacuumVesselBox",5*cm,5*cm,10*cm);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4LogicalVolume*vacuumVesselLogical=newG4LogicalVolume(vacuumVesselBox,stainlessSteel,"VacuumVesselLogical");G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4VPhysicalVolume*vacuumVesselPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,0),vacuumVesselLogical,"VacuumVesselPhysical",0,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4Material*gold=newG4Material("Gold",19.32*g/cm3,1);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);gold->AddElement(nist->FindOrBuildElement("Au"),1.0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4Box*emitterBox=newG4Box("EmitterBox",4*cm,4*cm,0.1*cm);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4LogicalVolume*emitterLogical=newG4LogicalVolume(emitterBox,gold,"EmitterLogical");G4VPhysicalVolume*emitterPhysical=newG4PVPlacement(0,G4ThreeVector(0,0,5*cm),emitterLogical,"EmitterPhysical",vacuumVesselLogical,false,0);G4VPhysicalVolume