《dm粒子物理与核物理实验中的数据分析》由会员分享，可在线阅读，更多相关《dm粒子物理与核物理实验中的数据分析（24页珍藏版）》请在金锄头文库上搜索。
1、2019/5/24,1,粒子物理与核物理实验中的数据分析,杨振伟 清华大学 第六讲:Geant4 的探测器模拟介绍(1),2019/5/24,2,本讲要点,Geant4 简介 一般介绍和历史 基本概念和内核 重要的用户类 如何安装Geant4 如何使用Geant4 (探测器几何与材料),2019/5/24,3,数据分析为什么需要模拟信息？,探测器,数据,重建,物理 分析,发现 n,发现 1,刻度,写成 root 或 zbs 等格式,把电学量变为 物理量并做各种刻度与修正。,把分立测量 联系起来给 出四动量等,问题在于我们如何知道所做的工作是正确的？,2019/5/24,4,模拟需要做的主要工作
2、,定义探测器几何(具有有限接受范围的探测器) 指定探测器的物质(有用的/没用的各种材料) 指定物理过程与粒子(什么东西在探测器中干了什么) 产生原始事例(模拟重复实验n次) 数据收集和分析(取出探测器可以获得的信息并处理分析),除了最后的分析，其它部分都可以由Geant4负责处理,2019/5/24,5,其应用领域主要包括：高能物理、核物理、加速器物理以及医学和空间科学的研究等。,什么是 Geant4?,简单地说，Geant4 是模拟粒子穿过物质的工具包。 提供了探测器模拟的完整工具：几何，探测器响应，运行、事例以及径迹管理、图形显示、用户接口等。 提供了极为丰富的可供选择的物理过程 详见ht
3、tp:/geant4.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日发布,2019/5/24,6,Geant4的特点?,Geant4 发展于 Geant3,充分利用了C+语言中的优势，是粒子与核物理实验模拟最
4、好的工具之一。 免费下载，更新迅速,Geant4 是最早成功使用面向对象环境重新设计的粒子与核物理软件包，并应用于新一代实验。,考虑了各个方面的需求，包括重离子物理、电荷与宇称破缺物理、宇宙线物理、天体物理、空间科学以及医学应用等。,为了满足这些需求，提供了强大的功能和灵活性。,Geant4 的应用远不止于粒子与核物理。,2019/5/24,7,Geant4 的基本概念与内核,Geant4 包含17大类 每一大类都有专门独立的研发工作小组负责。 各大类之间的界面(例如顶层的设计) 由整体结构工作小组来维护。,Geant4 的内核 处理运行号，事例，迹，步，击中，运动轨迹。 提供几何表示和物理过
5、程的架构。,2019/5/24,8,Geant4 软件包中的内容,代码 1 百万行语句 网上可以自由下载 还在不断的增加与更新 相关文件说明 6 份使用工具书 网上可以自由下载 应用例子 与软件包一起附送 可以参照工具书自由运行例子程序 包含了从简单到较复杂各种模拟程序,平台 Linux, MacOSX, SUN Windows/XP 需要预安装的软件 C+ 编译器 如g+, VC+等 CLHEP库 GNU Make, shell 图像与用户界面 X Windows OpenGL 或 MesaGL DAWN . 至少需要其中的一个,2019/5/24,9,安装Geant4(概要),1)建议使用
6、SLC3/4操作系统，支持和兼容性好，在SLC上安装非常简单。以下默认在SLC上工作。 2)下载安装CLHEP程序包(这是唯一需要预安装的程序) 3)下载Geant4软件包以及相应的数据文件(用于各种物理模型)，按照安装手册进行编译安装,如果系统版本相同(内核版本和g+版本)，把已经编译好的程序直接复制到其它机器上即可使用。 比如，对SLC3系统，直接复制training服务器/projects/soft/ext/clhep.tgz和g4.tgz到本地机器，解压缩到相应目录即可。,2019/5/24,10,设置Geant4环境变量,由于Geant4练习过程中有大量图形显示，要求大家在本地机器上
7、进行练习。本地机器Geant4安装在/projects/soft/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子目录。 以后
8、所有的geant4工作都在G4WORKDIR下进行。,2019/5/24,11,Geant4自带很多例子，从新手入门、高级、专家级例子： $G4INSTALL/examples/novice目录为入门级7个例子,首先看最简单的例子 cd $G4WORKDIR cp -a $G4INSTALL/examples/novice . cd novice/N01 gmake (编译) exampleN01 (执行),这是最简单的例子，只是告诉你Geant4如何工作，并且确认你安装的Geant4可以工作。 例子中定义了径迹探测器和三明治式量能器，每个事例发射一个虚拟粒子geantino(没有真实的物理作
9、用)。 后面我们会具体介绍如何定义探测器几何与材料、如何定义入射粒子、如何添加需要的物理过程、如何记录需要的实验数据。,例一： Novice Example N01,2019/5/24,12,最基本的主函数(exN01),#include “G4RunManager.hh“ #include “G4UImanager.hh“ #include “ExN01DetectorConstruction.hh“ #include “ExN01PhysicsList.hh“ #include “ExN01PrimaryGeneratorAction.hh“ int main() /构造默认的运行管理器(
10、强制) G4RunManager* runManager = new G4RunManager; /用探测器以及物理过程初始化运行管理器(强制) G4VUserDetectorConstruction* detector = new ExN01DetectorConstruction; runManager-SetUserInitialization(detector); G4VUserPhysicsList* physics = new ExN01PhysicsList; runManager-SetUserInitialization(physics); / 设置用户行为(强制) G4VU
11、serPrimaryGeneratorAction* gen_action = new ExN01PrimaryGeneratorAction; runManager-SetUserAction(gen_action); /.接下页,探测器、物理过程以及用户行为的定义随模拟的问题而变，在相应类中指定。,2019/5/24,13,/接上页. / 初始化G4 内核 runManager-Initialize(); / 获取UI管理器的指针，并设置verbosities。 G4UImanager* UI = G4UImanager:GetUIpointer(); UI-ApplyCommand(“/
12、run/verbose 1“); UI-ApplyCommand(“/event/verbose 1“); UI-ApplyCommand(“/tracking/verbose 1“); / 开始一次运行 G4int numberOfEvent = 3; runManager-BeamOn(numberOfEvent); / 结束作业 / 释放内存：用户行为、物理过程以及探测器描述属于运行管理器， / 将被运行管理器自动删除，所以不应该在主函数中删除之。 / 只需要删除运行管理器和其它动态指针即可。 delete runManager; return 0; ,最基本的主函数(exN01),如果
13、想图形显示探测器以及物理过程的信息，需要在main函数中用到G4VisManager,2019/5/24,14,添加如下代码进行图形显示,G4VisManager* visManager = new G4VisExecutive; visManager-Initialize(); UI-ApplyCommand(“/control/execute vis.mac“); /vis.mac文件中定义了跟图形显示有关的一些命令，比如观看角度、缩放比例等。,为了让事例显示在程序运行后仍然驻留，需要开启一个G4界面(G4UIsession) G4UIsession * session = 0; sess
14、ion = new G4UIterminal(); session-SessionStart();,注1：确保已经include了头文件G4VisManager.hh 注2：确保已定义了G4VIS_USE环境变量 注3：程序结束前，不要忘记删除session和visManager两个动态指针。,2019/5/24,15,如何定义探测器几何和材料,Geant4模拟首先要考虑的是模拟的物理问题主要包括哪些物理过程以及想采用的探测器构造。例N01的探测器是在ExN01DetectorConstruction类中定义的，该类继承于G4VUserDetectorConstruction。 探测器构造主要
15、包括：几何、材料、敏感区域、敏感区域的读出机制，本讲只涉及几何与材料。详见该类的定义以及实现，即ExN01DetectorConstruction.hh(cc),关键： 在Construct()函数中构造必要的物质、必要的几何形状，以及如何将基本单元放置(架构)成探测器。还可能包括磁场、敏感区域定义等。,2019/5/24,16,如何定义探测器材料(1),物质定义： isotopes G4Isotope elements G4Element molecules, compounds, mixtures G4Material G4Material属性包括温度、压强、状态、密度等，定义时有时需要明
16、显指出属性。,double density = 1.390*g/cm3;/密度 double a = 39.95*g/mole; /原子量 double 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-Ad
