I would like to simulate 10^10 events. Generally, due to a limitation related to the maximal parameter value in BeamOn(…), I could split one Run into a few sub-Runs:
for(int i = 0; i<1000; i++)
but I would like to use G4AnalysisManager in RunAction() and accumulate distributions at EndOfEventAction(…) in 1D histograms.
What is the best way to merge data between runs? Manually load saved histograms at BeginOfRunAction() and save updated at EndOfRunAction()… or there is something else?