Geant4 NO3のRunAction.ccの変更：

[1] N03/src/EXN03RunAction.cc　に入って

    G4cout
       << "\n--------------------End of Run------------------------------\n"
  の前に

  #include <iostream> 
  #include <fstream> 
    std::ofstream fout("data.txt", std::ios::app);
    fout << sumEAbs<<"  "<<rmsEAbs<<std::endl;
    fout.close();

  を入れてコンパイルし直す。

[2] また以下のようにマクロ muon_auto.mac をつくる。

# Idle> /control/execute muon_auto.mac
#
/N03/det/setNbOfLayers 1
#
/N03/det/setGapThick 0 cm
/N03/det/setAbsThick 1 cm
/N03/det/setAbsMat Lead
/N03/det/update
#
#
/run/initialize
#
/gun/particle mu-
/gun/energy 0.05 GeV
/run/beamOn 1000
/gun/energy 0.06 GeV
/run/beamOn 1000
/gun/energy 0.07 GeV
/run/beamOn 1000
/gun/energy 0.08 GeV
/run/beamOn 1000
/gun/energy 0.09 GeV
/run/beamOn 1000
/gun/energy 0.1 GeV
/run/beamOn 1000
/gun/energy 0.15 GeV
/run/beamOn 1000
/gun/energy 0.2 GeV
/run/beamOn 1000
/gun/energy 0.3 GeV
/run/beamOn 1000
/gun/energy 0.4 GeV
/run/beamOn 1000
/gun/energy 0.5 GeV
/run/beamOn 1000
/gun/energy 0.6 GeV
/run/beamOn 1000
/gun/energy 0.7 GeV
/run/beamOn 1000
/gun/energy 0.8 GeV
/run/beamOn 1000
/gun/energy 0.9 GeV
/run/beamOn 1000
/gun/energy 1.0 GeV
/run/beamOn 1000
/gun/energy 1.5 GeV
/run/beamOn 1000
/gun/energy 2.0 GeV
/run/beamOn 1000
/gun/energy 3.0 GeV
/run/beamOn 1000
/gun/energy 4.0 GeV
/run/beamOn 1000
/gun/energy 5.0 GeV
/run/beamOn 1000
/gun/energy 7.0 GeV
/run/beamOn 1000
/gun/energy 10. GeV
/run/beamOn 1000
/gun/energy 20. GeV
/run/beamOn 1000
/gun/energy 50. GeV
/run/beamOn 1000
/gun/energy 100. GeV
/run/beamOn 1000

[3] 実行する

/control/execute muon_auto.mac


[4] 計算結果は data.txt に入っている。