粒子物理与核物理但实验中数据分析.pptVIP

粒子物理与核物理实验中的数据分析 杨振伟 清华大学 第八讲:Geant4 的探测器模拟介绍(3) 上讲回顾 本讲要点 产生主事例 G4HEPEvtInterface 敏感探测器 取出敏感探测器的数据，并存入ROOT格式文件 事例产生子接口 G4HEPEvtInterface 很多时候，事例产生子已经存在，而且是Fortran语言。Geant4并不直接链接这些Fortran程序，而是提供了一个接口： G4HEPEvtInterface读取事例产生子生成的ASCII文件中的信息，重新生成G4PrimaryParticle对象，并关联到对应的G4PrimaryVertex 也就是说，G4HEPEvtInterface将/HEPEVT/公共块的信息转换为一个O-O数据结构。这个公共块在高能物理中被广泛使用。 用/HEPEVT/公共块生成ASCII文件 以HEPEVT格式输出的ASCII文件 使用HEPEvtInterface的例子 敏感探测器(Sensitive Detector) SD的首要任务是通过粒子“迹”(track)上的“步”(step)的信息，构造“击中”(hit)。 这些击中经过数字化，被读出几何读出的信息是真正的模拟结果。(当然在模拟中我们也可以忽略数字化而直接读出hit的信息) 用户敏感探测器继承于抽象基类G4VSensitiveDetector，用户需要完成3个主要函数： ProcessHits(G4Step* aStep, G4TouchableHistory*) 构造“击中”，被G4SteppingManager调用 Initialize(G4HCofThisEvent* HCE) 初始化，事例开始时调用，指定构造的“集中”与当前事例关联起来 EndOfEvent(G4HCofThisEvent*) 事例结束时调用 定义和添加敏感探测器 读取敏感探测器的信息 将模拟结果写入root文件 小结 练习 在例子N02的基础上，将模拟的信息存储到root文件中。这些信息包括：粒子的PDG号、质量、能量沉积、径迹长度。生成root文件后画出这些信息的直方图，并进行分析 修改探测器物质和入射粒子，重新运行，得到新的root文件，并画出储存信息的直方图。 在N03的基础上，加入敏感探测器。 参考资料 Geant4应用开发手册3.6节 Geant4应用开发手册4.4节 Geant4例子novice/N02,N04 * * 粒子定义 G4ParticleDefinition 6大类粒子：G4LeptonConstructor G4BosonConstructor G4MesonConstructor G4BaryonConstructor G4IonConstructor G4ShortlivedConstructor 产生主事例：G4ParticleGun 物理过程：电磁、强作用、衰变、光轻子-强子作用、光学、参数化、输运(必要过程) common block 将以下量写入文件中 第一行：NHEP，当前事例粒子数(包括中间态) 随后的NHEP行：每个粒子的ISTHEP,IDHEP,JDAHEP,PHEP信息 ISTHEP：粒子状态；IDHEP：粒子PDG号；JDAHEP：粒子衰变产物 位置的指针；PHEP(1-3,5)：粒子x,y,z动量，能量，质量 102 3 11 0 0 0+00 0+00 0+03 0-03 3 -11 0 0 0+00 0+00 -0+03 0-03 3 11 0 0 0+00 0+00 0+03 0+00 3 -11 0 0 0+00 0+00 -0+03 0+00 3 11 0 0 0.3