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

金录嘉 jinlujia at pku.edu.cn
Tue Aug 18 15:48:00 CEST 2020


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
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20200818/1402f756/attachment-0001.html>


More information about the Gate-users mailing list