geant4用G4Scintillation获取随机粒子打在已知闪烁体上的波形，并将输出的波形数据保存在txt文件中，用代码演示一下

以下是一个简单的示例代码，展示了如何使用Geant4中的G4Scintillation类获取随机粒子在闪烁体中产生的光子，并将其波形数据保存到txt文件中。

#include "G4Scintillation.hh"
#include "G4Scintillation.hh"
#include "G4OpticalPhoton.hh"
#include "G4ParticleDefinition.hh"
#include "G4ParticleTable.hh"
#include "G4Event.hh"
#include "G4RunManager.hh"
#include "G4SystemOfUnits.hh"
#include <iostream>
#include <fstream>

int main()
{
    // 定义材料和粒子
    G4Material* scintillator = G4Material::GetMaterial("G4_PLASTIC_SC_VINYLTOLUENE");
    G4ParticleDefinition* particle = G4OpticalPhoton::Definition();

    // 定义闪烁过程
    G4Scintillation* scintillation = new G4Scintillation();
    scintillation->SetScintillationYieldFactor(1.0);
    scintillation->SetTrackSecondariesFirst(true);

    // 初始化随机数生成器
    CLHEP::HepRandom::setTheEngine(new CLHEP::MTwistEngine);
    G4int seed = time(NULL);
    CLHEP::HepRandom::setTheSeed(seed);

    // 定义输出文件
    std::ofstream outfile("waveform.txt");

    // 生成一个事件
    G4Event event;
    G4RunManager::GetRunManager()->BeamOn(1, &event);

    // 获取所有在闪烁体中产生的光子
    std::vector<G4DynamicParticle*>* secondaries = event.GetTrajectoryContainer()[0]->GetSecondaries();
    for (auto secondary : *secondaries)
    {
        // 如果是光子
        if (secondary->GetDefinition() == particle)
        {
            // 获取该光子在闪烁体中产生的波形
            std::vector<G4double> waveform = scintillation->GetPhotonWaveform(secondary, scintillator);

            // 将波形数据保存到文件中
            for (G4double data : waveform)
            {
                outfile << data << " ";
            }
            outfile << std::endl;
        }
    }

    outfile.close();

    return 0;
}

这个示例代码假设你已经定义了一个包含G4_PLASTIC_SC_VINYLTOLUENE材料的几何体，以及在该几何体中发射了一些粒子。它将获取每个在闪烁体中产生的光子的波形，并将其写入名为“waveform.txt”的文件中。请注意，这个示例代码只生成了一个事件，因此它只会输出一个光子的波形数据。如果你需要获取更多光子的波形数据，需要修改代码以生成更多事件。