粒子物理与核物理实验中数.pdfVIP

粒子物理与核物理实验中的数 据分析 杨振伟 清华大学 第八讲:Geant4 的探测 器模拟介绍(3) 2009-4-22 1 上讲回顾 粒子定义 G4ParticleDefinition 6大类粒子：G4LeptonConstructor G4BosonConstructor G4MesonConstructor G4BaryonConstructor G4IonConstructor G4ShortlivedConstructor 产生事例：G4ParticleGun 物理过程：电磁、强作用、衰变、光轻子-强 子作用、光学、参数化、输运(必要过程) 2009-4-22 2 本讲要点 产生主事例 G4HEPEvtInterface 灵敏探测器 取出灵敏探测器的数据，并存入ROOT 格式文件 2009-4-22 3 事例产生子接口 G4HEPEvtInterface 很多时候，事例产生子已经存在，而且是Fortran 语言。Geant4并不直接链接这些Fortran程序， 而是提供了一个接口： G4HEPEvtInterface读取事例产生子生成的 ASCII文件中的信息，重新生成 G4PrimaryParticle对象，并关联到对应的 G4PrimaryVertex 也就是说，G4HEPEvtInterface将/HEPEVT/公共 块的信息转换为一个O-O数据结构。这个公共块 在高能物理中被广泛使用。 2009-4-22 4 用/HEPEVT/公共块生成ASCII文件 common block 将以下量写入文件中 第一行：NHEP，当前事例粒子数(包括中间态) 随后的NHEP行：每个粒子的ISTHEP,IDHEP,JDAHEP,PHEP信息 ISTHEP：粒子状态；IDHEP：粒子PDG号；JDAHEP：粒子衰变产物 位置的指针；PHEP(1-3,5)：粒子x,y,z动量，能量，质量 2009-4-22 5 以HEPEVT格式输出的ASCII文件 比如：下面这个事例表示该事例共102个粒子(包括中间态)， 随后的102行分别为这102个粒子的具体信息： 第一列为粒子状态(3：对撞入射粒子或其它；2：衰变了；1： 存在的粒子；0：空)， 第2列为粒子PDG号， 最后4列分别为粒子的x，y，z方向动量和质量。 102 3 11 0 0 0+00 0+00 0+03 0-03 3 -11 0 0 0.