[Gate-users] Run Aborted after *** events processed

金录嘉 jinlujia at pku.edu.cn
Wed Aug 19 11:52:40 CEST 2020


Hello Gate Users,
I am the questioner and I am here again. At someone's prompt, I turned all vis function commands commented, and I set: 


/gate/verbose Error     2
/gate/verbose Warning   2
/run/verbose  2
/event/verbose 1
/tracking/verbose 0


...to see more details of the error during run. 
Everything ran well until after the 9996th event was processed. The system prompted "0 primaries are passed from G4EventTransformer", followed by "Run aborted after 9996 events processed". The complete log file of the running process is shown in the attachment.

I tried the ecat system and the cylindricalPET system, and also tried different scanners and phantoms setups, but the results were all the same. I tried running it on vGate v9.0, vGate v8.0 and the Linux server with Gate9 installed, but the problem still exists. I think this does not seem to be a problem of operating environment or memory overflow.
I need help urgently. If anyone can make valuable suggestions, I would be very grateful!

Lujia Jin,

Peking University




-----原始邮件-----
发件人:"金录嘉" <jinlujia at pku.edu.cn>
发送时间:2020-08-18 21:48:00 (星期二)
收件人: gate-users at lists.opengatecollaboration.org
抄送:
主题: Run Aborted after *** events processed



Hello Gate Users,

 I built an Ecat system for PET simulation and I successfully got the output file in .root format. But there seemed to be some errors when the program ended:




====================================================================
G4VisManager: Using G4TrajectoryDrawByCharge as fallback trajectory model.
See commands in /vis/modeling/trajectories/ for other options.
### Run 0 starts.
[Core-1]  Begin Of Run 0

### ===  Deexcitation model UAtomDeexcitation is activated for 5 regions:
          DefaultRegionForTheWorld  1  1  1
          ecat  1  1  1
          block  1  1  1
          crystal  1  1  1
          phantom  1  1  1
### ===  Ignore cuts flag:   0
WARNING: G4VisManager::EndOfEvent: Automatic event keeping suspended.
  The number of events exceeds the maximum, 100, that may be kept by
  the vis manager.
********************* WARNING! ********************
*  Display list limit reached in OpenGL.
*  Continuing drawing WITHOUT STORING. Scene only partially refreshable.
*  Current limit: 50000 primitives.  Change with "/vis/ogl/set/displayListLimit".


***************************************************
 Run terminated.
Run Summary
  Run Aborted after 9996 events processed.
  User=5.880000s Real=6.027382s Sys=0.120000s
[Core-1]  End Of Run 0
[Core-1]  Run 0  ---  Elapsed time = 6.1
100 events have been kept for refreshing and/or reviewing.

  "/vis/reviewKeptEvents" to review them one by one.
  "/vis/enable", then "/vis/viewer/flush" or "/vis/viewer/rebuild" to see them accumulated.
WARNING: G4VisManager::EndOfRun: Automatic event keeping was suspended.
  The number of events in the run exceeded the maximum, 100, that may be
  kept by the vis manager.
  The number of events kept by the vis manager can be changed with
  "/vis/scene/endOfEventAction accumulate <N>", where N is the
  maximum number you wish to allow.  N < 0 means "unlimited".
[Core-0] End of macro PET_Ecat_System.mac
Graphics systems deleted.
Visualization Manager deleting...
G4 kernel has come to Quit state.
================== Deleting memory pools ===================
Number of memory pools allocated: 22; of which, static: 6
Dynamic pools deleted: 16 / Total memory freed: 3.2 MB
============================================================
RunManagerKernel is deleted. Good bye :)




 I think "Run Aborted after *** events processed" means that the program has exited abnormally. Changing the settings of scanner or phantom would not eliminate this problem. However, the cylindrical PET system was able to operate successfully. So maybe someone has encountered a similar problem or knows a solution?


 The content of my entire mac file is as follows:

#=====================================================
# VISUALISATION
#=====================================================
/vis/open                   OGLSX
/vis/viewer/set/viewpointThetaPhi   60 60
/vis/viewer/zoom               6.0
/vis/drawVolume
/tracking/storeTrajectory        1
/vis/scene/add/trajectories 
/vis/scene/endOfEventAction       accumulate

#=====================================================
# VERBOSITY
#=====================================================
/gate/verbose Physic    1
/gate/verbose Cuts      1
/gate/verbose SD        1
/gate/verbose Actions   1
/gate/verbose Actor     1
/gate/verbose Step      1
/gate/verbose Error     1
/gate/verbose Warning   1
/gate/verbose Output    1
/gate/verbose Beam      1
/gate/verbose Volume    1
/gate/verbose Image     1
/gate/verbose Geometry  1
/gate/verbose Core      1

/run/verbose 1
/event/verbose 0
/tracking/verbose 0      

#=====================================================
# GEOMETRY
#=====================================================
/gate/geometry/setMaterialDatabase    GateMaterials.db

#
#     W O R L D
#
/gate/world/geometry/setXLength       400. cm
/gate/world/geometry/setYLength       400. cm
/gate/world/geometry/setZLength       400. cm

#     E C A T
/gate/world/daughters/name ecat
/gate/world/daughters/insert cylinder
/gate/ecat/setMaterial Air
/gate/ecat/geometry/setRmax 44.2 cm
/gate/ecat/geometry/setRmin 41.2 cm
/gate/ecat/geometry/setHeight 15.52 cm
/gate/ecat/vis/forceWireframe

#     B L O C K
/gate/ecat/daughters/name block
/gate/ecat/daughters/insert box
/gate/block/placement/setTranslation 427.0 0.0 0.0 mm
/gate/block/geometry/setXLength 30.0 mm
/gate/block/geometry/setYLength 35.8594 mm
/gate/block/geometry/setZLength 38.7 mm
/gate/block/setMaterial Air
/gate/block/vis/forceWireframe

#C R Y S T A L
/gate/block/daughters/name crystal
/gate/block/daughters/insert box
/gate/crystal/geometry/setXLength 30.0 mm
/gate/crystal/geometry/setYLength 4.4 mm
/gate/crystal/geometry/setZLength 4.75 mm
/gate/crystal/setMaterial BGO
/gate/crystal/vis/setColor yellow

#R E P E A T    C R Y S T A L
/gate/crystal/repeaters/insert cubicArray
/gate/crystal/cubicArray/setRepeatNumberX 1
/gate/crystal/cubicArray/setRepeatNumberY 8
/gate/crystal/cubicArray/setRepeatNumberZ 8
/gate/crystal/cubicArray/setRepeatVector 0. 4.4942 4.85 mm

#R E P E A T    BLOCK
/gate/block/repeaters/insert linear
/gate/block/linear/setRepeatNumber 4
/gate/block/linear/setRepeatVector 0. 0. 38.8 mm
/gate/block/repeaters/insert ring
/gate/block/ring/setRepeatNumber 72

# A T T A C H    S Y S T E M 
/gate/systems/ecat/block/attach block
/gate/systems/ecat/crystal/attach crystal

#A T T A C H    C R Y S T A L  SD
/gate/crystal/attachCrystalSD

#=====================================================
# PHANTOM
#=====================================================
/gate/world/daughters/name phantom
/gate/world/daughters/insert cylinder
/gate/phantom/setMaterial Water
/gate/phantom/vis/forceSolid
/gate/phantom/vis/setColor blue
/gate/phantom/geometry/setRmax   10.0 cm
/gate/phantom/geometry/setHeight 20. cm
/gate/phantom/attachPhantomSD

#=====================================================
# PHYSICS
#=====================================================
/gate/physics/addPhysicsList QGSP_BERT_EMV
/gate/physics/addProcess PhotoElectric
/gate/physics/processes/PhotoElectric/setModel StandardModel

/gate/physics/addProcess Compton
/gate/physics/processes/Compton/setModel StandardModel

/gate/physics/addProcess RayleighScattering
/gate/physics/processes/RayleighScattering/setModel PenelopeModel

/gate/physics/addProcess ElectronIonisation
/gate/physics/processes/ElectronIonisation/setModel StandardModel e-
/gate/physics/processes/ElectronIonisation/setModel StandardModel e+

/gate/physics/addProcess Bremsstrahlung
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e-
/gate/physics/processes/Bremsstrahlung/setModel StandardModel e+

/gate/physics/addProcess PositronAnnihilation

/gate/physics/addProcess MultipleScattering e+
/gate/physics/addProcess MultipleScattering e-

/gate/physics/processList Enabled
/gate/physics/processList Initialized

#=====================================================
# CUTS
#=====================================================
/gate/physics/Gamma/SetCutInRegion      crystal 1.0 cm
/gate/physics/Electron/SetCutInRegion   crystal 1.0 cm
/gate/physics/Positron/SetCutInRegion   crystal 1.0 cm

/gate/physics/Gamma/SetCutInRegion      phantom 0.1 mm
/gate/physics/Electron/SetCutInRegion   phantom 0.1 mm
/gate/physics/Positron/SetCutInRegion   phantom 0.1 mm

/gate/physics/SetMaxStepSizeInRegion    phantom 0.01 mm

#=====================================================
# INITIALISATION
#=====================================================
/gate/run/initialize

#=====================================================
# DIGITIZER
#=====================================================
#A D D E R
/gate/digitizer/Singles/insert adder

#       R E A D O U T 
/gate/digitizer/Singles/insert readout
/gate/digitizer/Singles/readout/setDepth 1

#      E N E R G Y   B L U R R I N G
/gate/digitizer/Singles/insert blurring
/gate/digitizer/Singles/blurring/setResolution 0.26
/gate/digitizer/Singles/blurring/setEnergyOfReference 511. keV

#       E N E R G Y   C U T
/gate/digitizer/Singles/insert thresholder
/gate/digitizer/Singles/thresholder/setThreshold 250. keV
/gate/digitizer/Singles/insert upholder
/gate/digitizer/Singles/upholder/setUphold 750. keV

#       D E A D    T I M E 
/gate/digitizer/Singles/insert deadtime
/gate/digitizer/Singles/deadtime/setDeadTime 3000000. ps
/gate/digitizer/Singles/deadtime/setMode paralysable
/gate/digitizer/Singles/deadtime/chooseDTVolume block

#C O I N C I    S O R T E R 
/gate/digitizer/Coincidences/setWindow 10. ns
/gate/digitizer/Coincidences/setOffset 0. ns
/gate/digitizer/Coincidences/describe

/gate/digitizer/name delay
/gate/digitizer/insert coincidenceSorter
/gate/digitizer/delay/setWindow 10. ns
/gate/digitizer/delay/setOffset 500. ns
/gate/digitizer/delay/describe

/gate/digitizer/name finalCoinc
/gate/digitizer/insert coincidenceChain
/gate/digitizer/finalCoinc/addInputName delay
/gate/digitizer/finalCoinc/addInputName Coincidences
/gate/digitizer/finalCoinc/usePriority true
/gate/digitizer/finalCoinc/describe

#  R A N D O M
/gate/random/setEngineName JamesRandom
/gate/random/setEngineSeed default
/gate/random/verbose 1

#=====================================================
# OUTPUT
#=====================================================
# S E T U P  - R O O T  F I L E
/gate/output/root/enable
/gate/output/root/setFileName YourPETSimu
/gate/output/root/setRootHitFlag 1
/gate/output/root/setRootSinglesFlag 1
/gate/output/root/setRootCoincidencesFlag 1

#=====================================================
# SOURCE
#=====================================================
# F A S T   P R O D U C T I O N   M O D U L E   F O R   F - 1 8   I S O T O P E S
/gate/source/addSource F18CylinderSource
/gate/source/F18CylinderSource/setActivity 10000. becquerel
/gate/source/F18CylinderSource/gps/particle e+
/gate/source/F18CylinderSource/setForcedUnstableFlag true
/gate/source/F18CylinderSource/setForcedHalfLife 6586 s
/gate/source/F18CylinderSource/gps/energytype Fluor18
/gate/source/F18CylinderSource/gps/type Volume
/gate/source/F18CylinderSource/gps/shape Cylinder
/gate/source/F18CylinderSource/gps/radius 0.5 mm
/gate/source/F18CylinderSource/gps/halfz 0.5 mm
/gate/source/F18CylinderSource/gps/angtype iso
/gate/source/F18CylinderSource/gps/centre 0.0 0.0 0. cm

/gate/source/list

#
# EXPERIMENT
#
/gate/application/setTimeSlice     1.  s
/gate/application/setTimeStart     0.  s
/gate/application/setTimeStop      1.  s

#
# LET'S RUN THE SIMULATION!
#
/gate/application/startDAQ




 Thank you very much for taking time out of your busy schedule to read my mail. I would be grateful if anyone can provide some ideas to solve the problem!




Lujia Jin

Peking University






Lujia Jin

Department of Biomedical Engineering

College of Engineering

Peking University

Beijing, China, 100871

TEL +86-18813118588

 

金录嘉

北京大学工学院

生物医学工程系

中国 北京 100871

电话 18813118588

 

 

 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200819/3220e9a0/attachment-0001.html>
-------------- next part --------------
An embedded and charset-unspecified text was scrubbed...
Name: log.txt
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200819/3220e9a0/attachment-0001.txt>


More information about the Gate-users mailing list