geant4二次粒子信息获取

在 Geant4 中，可以通过实现 G4SteppingVerbose 类来获取二次粒子的信息。这个类是 Geant4 中用于输出每一步模拟信息的基类。你可以继承 G4SteppingVerbose 类并覆盖其成员函数，以输出你需要的信息。

以下是一个简单的例子，演示了如何获取二次粒子的信息：

#include "G4SteppingVerbose.hh"
#include "G4Track.hh"

class MySteppingVerbose: public G4SteppingVerbose
{
public:
    void TrackingStarted() override
    {
        G4SteppingVerbose::TrackingStarted();
        G4Track* track = fpSteppingManager->GetTrack();
        G4cout << "Primary particle: " << track->GetDefinition()->GetParticleName() << G4endl;
    }

    void StepInfo() override
    {
        G4SteppingVerbose::StepInfo();
        const G4Step* step = fpSteppingManager->GetStep();
        G4Track* track = fpSteppingManager->GetTrack();
        if (track->GetCurrentStepNumber() == 1) {
            G4cout << "Secondary particle: " << step->GetSecondary()->at(0)->GetDefinition()->GetParticleName() << G4endl;
        }
    }
};

int main()
{
    // ...
    MySteppingVerbose* mySteppingVerbose = new MySteppingVerbose();
    runManager->SetUserVerbose(mySteppingVerbose);
    // ...
    return 0;
}

在上面的例子中，我们继承了 G4SteppingVerbose 类，并覆盖了其中的 TrackingStarted() 和 StepInfo() 函数。TrackingStarted() 函数在开始跟踪一个粒子时调用，我们在这里输出了主粒子的信息。StepInfo() 函数在每一步结束后调用，我们在这里判断当前步骤是否为第一步，如果是，则输出第一个二次粒子的信息。

注意，为了使用自定义的 G4SteppingVerbose 类，我们需要在创建 G4RunManager 对象之后，调用 SetUserVerbose() 函数，将其设置为用户定义的详细模式。