《粒子物理及核物理实验中数据分析ppt课件.ppt》由会员分享，可在线阅读，更多相关《粒子物理及核物理实验中数据分析ppt课件.ppt（24页珍藏版）》请在三一办公上搜索。
1、2022/11/19,1,粒子物理与核物理实验中的数据分析,杨振伟清华大学第六讲:Geant4 的探测器模拟介绍(1),2022/11/19,2,本讲要点,Geant4 简介一般介绍和历史基本概念和内核重要的用户类如何安装Geant4如何使用Geant4 (探测器几何与材料),2022/11/19,3,数据分析为什么需要模拟信息？,探测器,数据,重建,物理分析,发现 n,发现 1,刻度,写成 root 或 zbs 等格式,把电学量变为物理量并做各种刻度与修正。,把分立测量联系起来给出四动量等,问题在于我们如何知道所做的工作是正确的？,2022/11/19,4,模拟需要做的主要工作,定义探测器几
2、何(具有有限接受范围的探测器) 指定探测器的物质(有用的/没用的各种材料) 指定物理过程与粒子(什么东西在探测器中干了什么) 产生原始事例(模拟重复实验n次) 数据收集和分析(取出探测器可以获得的信息并处理分析),除了最后的分析，其它部分都可以由Geant4负责处理,2022/11/19,5,其应用领域主要包括：高能物理、核物理、加速器物理以及医学和空间科学的研究等。,什么是 Geant4?,简单地说，Geant4 是模拟粒子穿过物质的工具包。 提供了探测器模拟的完整工具：几何，探测器响应，运行、事例以及径迹管理、图形显示、用户接口等。 提供了极为丰富的可供选择的物理过程 详见http:/ge
3、ant4.cern.ch/support/about.shtml,参考 资料 1) http:/geant4.cern.ch 2)Nuclear Instruments and Methods in Physics Research A 506 (2003) 250-303, and IEEE Transactions on Nuclear Science 53 No. 1 (2006) 270-278. 最新版为9.1版，于2008年2月5日发布,2022/11/19,6,Geant4的特点?,Geant4 发展于 Geant3,充分利用了C+语言中的优势，是粒子与核物理实验模拟最好的工具之
4、一。 免费下载，更新迅速,Geant4 是最早成功使用面向对象环境重新设计的粒子与核物理软件包，并应用于新一代实验。,考虑了各个方面的需求，包括重离子物理、电荷与宇称破缺物理、宇宙线物理、天体物理、空间科学以及医学应用等。,为了满足这些需求，提供了强大的功能和灵活性。,Geant4 的应用远不止于粒子与核物理。,2022/11/19,7,Geant4 的基本概念与内核,Geant4 包含17大类每一大类都有专门独立的研发工作小组负责。各大类之间的界面(例如顶层的设计) 由整体结构工作小组来维护。,Geant4 的内核处理运行号，事例，迹，步，击中，运动轨迹。提供几何表示和物理过程的架构。,20
5、22/11/19,8,Geant4 软件包中的内容,代码1 百万行语句网上可以自由下载还在不断的增加与更新相关文件说明6 份使用工具书网上可以自由下载应用例子与软件包一起附送可以参照工具书自由运行例子程序包含了从简单到较复杂各种模拟程序,平台Linux, MacOSX, SUNWindows/XP需要预安装的软件C+ 编译器 如g+, VC+等CLHEP库GNU Make, shell图像与用户界面X WindowsOpenGL 或 MesaGLDAWN.至少需要其中的一个,2022/11/19,9,安装Geant4(概要),1)建议使用SLC3/4操作系统，支持和兼容性好，在SLC上安装非常
6、简单。以下默认在SLC上工作。2)下载安装CLHEP程序包(这是唯一需要预安装的程序)3)下载Geant4软件包以及相应的数据文件(用于各种物理模型)，按照安装手册进行编译安装,如果系统版本相同(内核版本和g+版本)，把已经编译好的程序直接复制到其它机器上即可使用。比如，对SLC3系统，直接复制training服务器/projects/soft/ext/clhep.tgz和g4.tgz到本地机器，解压缩到相应目录即可。,2022/11/19,10,设置Geant4环境变量,由于Geant4练习过程中有大量图形显示，要求大家在本地机器上进行练习。本地机器Geant4安装在/projects/so
7、ft/ext/g4/geant4.9.0.p01/,运行Geant4之前需要先设置Geant4的环境变量，即运行：source /projects/soft/ext/g4/geant4.9.0.p01/env.sh,Geant4环境变量都以G4开头。比如G4INSTALL : 安装目录G4SYSTEM : 操作系统及c+编译器名称,如Linux-g+G4WORKDIR : 工作目录我们设定的工作目录为/projects/$USER/geant4所以大家需要在/projects/$USER目录下建立geant4子目录。以后所有的geant4工作都在G4WORKDIR下进行。,2022/11/19
8、,11,Geant4自带很多例子，从新手入门、高级、专家级例子： $G4INSTALL/examples/novice目录为入门级7个例子,首先看最简单的例子cd $G4WORKDIRcp -a $G4INSTALL/examples/novice .cd novice/N01gmake (编译)exampleN01 (执行),这是最简单的例子，只是告诉你Geant4如何工作，并且确认你安装的Geant4可以工作。例子中定义了径迹探测器和三明治式量能器，每个事例发射一个虚拟粒子geantino(没有真实的物理作用)。后面我们会具体介绍如何定义探测器几何与材料、如何定义入射粒子、如何添加需要的物
9、理过程、如何记录需要的实验数据。,例一： Novice Example N01,2022/11/19,12,最基本的主函数(exN01),#include G4RunManager.hh#include G4UImanager.hh#include ExN01DetectorConstruction.hh#include ExN01PhysicsList.hh#include ExN01PrimaryGeneratorAction.hhint main() /构造默认的运行管理器(强制) G4RunManager* runManager = new G4RunManager; /用探测器以及物
10、理过程初始化运行管理器(强制) G4VUserDetectorConstruction* detector = new ExN01DetectorConstruction; runManager-SetUserInitialization(detector); G4VUserPhysicsList* physics = new ExN01PhysicsList; runManager-SetUserInitialization(physics); / 设置用户行为(强制) G4VUserPrimaryGeneratorAction* gen_action = new ExN01PrimaryG
11、eneratorAction; runManager-SetUserAction(gen_action); /.接下页,探测器、物理过程以及用户行为的定义随模拟的问题而变，在相应类中指定。,2022/11/19,13,/接上页. / 初始化G4 内核 runManager-Initialize(); / 获取UI管理器的指针，并设置verbosities。 G4UImanager* UI = G4UImanager:GetUIpointer(); UI-ApplyCommand(/run/verbose 1); UI-ApplyCommand(/event/verbose 1); UI-App
12、lyCommand(/tracking/verbose 1); / 开始一次运行 G4int numberOfEvent = 3; runManager-BeamOn(numberOfEvent); / 结束作业 / 释放内存：用户行为、物理过程以及探测器描述属于运行管理器， / 将被运行管理器自动删除，所以不应该在主函数中删除之。 / 只需要删除运行管理器和其它动态指针即可。 delete runManager; return 0;,最基本的主函数(exN01),如果想图形显示探测器以及物理过程的信息，需要在main函数中用到G4VisManager,2022/11/19,14,添加如下代码
13、进行图形显示,G4VisManager* visManager = new G4VisExecutive;visManager-Initialize();UI-ApplyCommand(/control/execute vis.mac);/vis.mac文件中定义了跟图形显示有关的一些命令，比如观看角度、缩放比例等。,为了让事例显示在程序运行后仍然驻留，需要开启一个G4界面(G4UIsession)G4UIsession * session = 0;session = new G4UIterminal();session-SessionStart();,注1：确保已经include了头文件G4
14、VisManager.hh注2：确保已定义了G4VIS_USE环境变量注3：程序结束前，不要忘记删除session和visManager两个动态指针。,2022/11/19,15,如何定义探测器几何和材料,Geant4模拟首先要考虑的是模拟的物理问题主要包括哪些物理过程以及想采用的探测器构造。例N01的探测器是在ExN01DetectorConstruction类中定义的，该类继承于G4VUserDetectorConstruction。探测器构造主要包括：几何、材料、敏感区域、敏感区域的读出机制，本讲只涉及几何与材料。详见该类的定义以及实现，即ExN01DetectorConstructio
15、n.hh(cc),关键： 在Construct()函数中构造必要的物质、必要的几何形状，以及如何将基本单元放置(架构)成探测器。还可能包括磁场、敏感区域定义等。,2022/11/19,16,如何定义探测器材料(1),物质定义：isotopes G4Isotopeelements G4Elementmolecules, compounds, mixtures G4MaterialG4Material属性包括温度、压强、状态、密度等，定义时有时需要明显指出属性。,double density = 1.390*g/cm3;/密度double a = 39.95*g/mole; /原子量double
16、z=18.; /原子序数G4Material* lAr = new G4Material(“liquidArgon”,z,a,dentsity);,G4Element* H = new G4Element(name=Hydrogen,symbol=H , z= 1., a);G4Element* O = new G4Element(name=Oxygen ,symbol=O , z= 8., a);density = 1.000*g/cm3;G4Material* H2O = new G4Material(name=Water, density, ncomponents=2);H2O-AddE
17、lement(H, natoms=2);H2O-AddElement(O, natoms=1); /定义水，给定密度、元素种类数目、添加元素,注：Geant4中的数据如果有单位都要明显给出单位。 比如double length = 1.0*cm = 10.0*mm 二者效果是一样的。方便程序中使用合适的单位，而且看代码的时候比较直观。,2022/11/19,17,/空气的定义 G4Element* N = new G4Element(Nitrogen, N, z=7., a= 14.01*g/mole); G4Element* O = new G4Element(Oxygen , O, z=8
18、., a= 16.00*g/mole); G4Material* Air = new G4Material(Air, density= 1.29*mg/cm3, nel=2); Air-AddElement(N, 70*perCent); Air-AddElement(O, 30*perCent); /指定百分比,如何定义探测器材料(2),/定义真空 G4double a, z; G4double density, temperature, pressure; G4int nel; G4Material* Vacuum = new G4Material(Galactic, z=1., a=1.
19、01*g/mole, density= universe_mean_density, kStateGas, 2.73*kelvin, 3.e-18*pascal);/指定温度、压强等属性,2022/11/19,18,三个概念层面,G4VSolid 形状，尺寸G4LogicalVolume 子物理体积，材料，灵敏区， 磁场，用户限制等等。G4VPhysicalVolume 位置，转动,原则上几何描述的三个层面允许最大限度的信息再利用，以便减小内存空间。,探测器灵敏区应由用户在G4VSensitiveDetector 进行具体的定义并设置到G4LogicalVolume 上。,探测器几何在Gean
20、t4中是由很多“体积”(Volumes)构成的。Volume的概念在几何构造中非常重要。,如何定义探测器几何(1),2022/11/19,19,G4double expHall_x = 3.0*m;G4double expHall_y = 1.0*m;G4double expHall_z = 1.0*m;/Solid，指定几何形状和尺寸G4Box* experimentalHall_box = new G4Box(“expHall_box”,expHall_x,expHall_y,expHall_z); /Logical，指定具体物理特性，如其中物质为Ar气experimentalHall_l
21、og = new G4LogicalVolume(experimentalHall_box, Ar,expHall_log,0,0,0);/Physical，指定放置位置以及旋转角度等experimentalHall_phys = new G4PVPlacement(0,G4ThreeVector(), experimentalHall_log,expHall,0,false,0);,如何定义探测器几何(2),定义一个G4Box，长方体,Geant4中基本的几何形状很多，比如：G4Box, G4Tubs, G4Cons, G4Para, G4Trd, G4Trap, G4Spere, G4To
22、rus, G4Polycone,.,2022/11/19,20,一个 Geant4 体积的调用示意图,2022/11/19,21,作为用户，如何使用Geant4？,首先不要尝试去读懂 Geant4 上百万行语句的含义！,其次不要尝试去读懂所有 Geant4 的例子！,然后根据你的实际情况选择相应的例子，我个人建议大家仔细读懂“Novice”中的第二个例子N02中在 src 与 include 中的每一行语句的含义。,Geant 4 需要用户告诉探测器的几何是什么？材料是什么？研究的粒子是什么？物理过程什么？,用户需要从 Geant4 中得到粒子在反应后或者输运后，可以被探测器记录的击中数，末态
23、粒子数，对应的四动量，能损大小等等。然后把所有物理量都转化为电学量，如光电子数，电子学道数等等。,2022/11/19,22,小结,模拟探测器的必要性Geant4 模拟软件包的基本概念安装Geant4和设置环境变量模拟探测器的简单例子图形显示探测器几何和材料的定义,2022/11/19,23,练习,1. 修改稿例子N01，加上图形显示。 从不同视角查看探测器显示情况。 (hints: 将N02中的vis.mac复制到N01目录，利用其中关于vis的UICommand命令。)2. 修改N01中探测器的物质以及几何形状，在tracker tube的内部中心平行放置一个内径20cm,外径25cm，长
24、度40cm的小圆管，管的材质为硅。3. 编译运行N02例子，查看探测器以及物质构造。4. 阅读Geant4 Users Guide (for Application Developers)第2.1-2.3节，4.1-4.3节。5. 下载Geant4.9.1源码，尝试安装Geant4,2022/11/19,24,参考资料,http:/geant4.cern.chGeant4 Users Guide (for Application Developers)第2.1-2.3节，4.1-4.3节。可在http:/ Instruments and Methods in Physics Research A 506 (2003) 250-303, and IEEE Transactions on Nuclear Science 53 No. 1 (2006) 270-278.,
