[Gate-users] The question about the speed of calculation

[hao qiushi]

Dear Gate Expert,
I have simulated the radiotherapy of carbon ion beam to hit the mouse model,but the simulation speed is only a few particles per second,and my project needs to simulate 10E9 particles,I do not know if I'm the only one that's so slow,could someone help me to check my macro file? If possible ,please help me to  verify the calculation speed ,thanks in advance, the macro file is attached below


/control/execute mac/verbose.mac
#/control/execute mac/visu.mac

/gate/geometry/setMaterialDatabase /home/gate/Downloads/vgate/dosimetry/data/GateMaterials.db

# World=======================

/gate/world/geometry/setXLength 2 m
/gate/world/geometry/setYLength 2 m
/gate/world/geometry/setZLength 2 m
/gate/world/setMaterial Vacuum

# digimouse ==================

/gate/HounsfieldMaterialGenerator/SetMaterialTable                  /home/gate/Downloads/vgate/dosimetry/data/Schneider2000MaterialsTable.txt
/gate/HounsfieldMaterialGenerator/SetDensityTable                   /home/gate/Downloads/vgate/dosimetry/data/Schneider2000DensitiesTable.txt
/gate/HounsfieldMaterialGenerator/SetDensityTolerance               0.1 g/cm3
/gate/HounsfieldMaterialGenerator/SetOutputMaterialDatabaseFilename /home/gate/Downloads/vgate/dosimetry/data/patient-HUmaterials.db
/gate/HounsfieldMaterialGenerator/SetOutputHUMaterialFilename       /home/gate/Downloads/vgate/dosimetry/data/patient-HU2mat.txt
/gate/HounsfieldMaterialGenerator/Generate

/gate/world/daughters/name                        pmma
/gate/world/daughters/insert                      ImageNestedParametrisedVolume
/gate/geometry/setMaterialDatabase                /home/gate/Downloads/vgate/dosimetry/data/patient-HUmaterials.db
/gate/pmma/geometry/setHUToMaterialFile           /home/gate/Downloads/vgate/dosimetry/data/patient-HU2mat.dat
/gate/pmma/geometry/setImage                      /home/gate/Downloads/vgate/dosimetry/data/ct_data/ct_380x992x208.hdr

/gate/pmma/placement/setTranslation               0 0 500 mm
/gate/pmma/placement/setRotationAxis              1 0 0
/gate/pmma/placement/setRotationAngle             -90 deg

# PHYSICS=========

/gate/physics/addPhysicsList QGSP_BERT_EMV

/gate/physics/addPhysicsList QGSP_BERT_HP

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

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

# output===================

/gate/actor/addActor                               DoseActor  doseDistribution
/gate/actor/doseDistribution/save                  /home/gate/Downloads/vgate/dosimetry/output_txt/output-C100.mhd
/gate/actor/doseDistribution/attachTo         pmma
#/gate/actor/doseDistribution/setSize              38 99.2 20.8 mm
/gate/actor/doseDistribution/stepHitType           random
/gate/actor/doseDistribution/setPosition           0 0 0 cm
/gate/actor/doseDistribution/setResolution         380 992 208
/gate/actor/doseDistribution/enableEdep            true
/gate/actor/doseDistribution/enableDose            true
/gate/actor/doseDistribution/enableNumberOfHits    true

/gate/actor/addActor                               DoseActor  doseDistribution_txt
/gate/actor/doseDistribution_txt/save                  /home/gate/Downloads/vgate/dosimetry/output_txt/output-C100.txt
/gate/actor/doseDistribution_txt/attachTo         pmma
#/gate/actor/doseDistribution_txt/setSize               38 99.2 20.8 mm
/gate/actor/doseDistribution_txt/stepHitType           random
/gate/actor/doseDistribution_txt/setPosition           0 0 0 cm
/gate/actor/doseDistribution_txt/setResolution         1 992 1
/gate/actor/doseDistribution_txt/enableEdep            true
/gate/actor/doseDistribution_txt/enableDose            true
/gate/actor/doseDistribution_txt/enableNumberOfHits    true



# ProductionAndStoppingActor==================

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_11
/gate/actor/prdandstp_11/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop11.mhd
/gate/actor/prdandstp_11/attachTo           pmma
/gate/actor/prdandstp_11/setResolution      380 992 208
/gate/actor/prdandstp_11/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_10
/gate/actor/prdandstp_10/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop10.mhd
/gate/actor/prdandstp_10/attachTo           pmma
/gate/actor/prdandstp_10/setResolution      380 992 208
/gate/actor/prdandstp_10/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_15
/gate/actor/prdandstp_15/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop15.mhd
/gate/actor/prdandstp_15/attachTo           pmma
/gate/actor/prdandstp_15/setResolution      380 992 208
/gate/actor/prdandstp_15/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_11_txt
/gate/actor/prdandstp_11_txt/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop11.txt
/gate/actor/prdandstp_11_txt/attachTo           pmma
/gate/actor/prdandstp_11_txt/setResolution      1 992 1
/gate/actor/prdandstp_11_txt/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_10_txt
/gate/actor/prdandstp_10_txt/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop10.txt
/gate/actor/prdandstp_10_txt/attachTo           pmma
/gate/actor/prdandstp_10_txt/setResolution      1 992 1
/gate/actor/prdandstp_10_txt/stepHitType        random

/gate/actor/addActor                     ProductionAndStoppingActor    prdandstp_15_txt
/gate/actor/prdandstp_15_txt/save               /home/gate/Downloads/vgate/dosimetry/output_txt/product_stop15.txt
/gate/actor/prdandstp_15_txt/attachTo           pmma
/gate/actor/prdandstp_15_txt/setResolution      1 992 1
/gate/actor/prdandstp_15_txt/stepHitType        random

/gate/actor/addActor               SimulationStatisticActor stat
/gate/actor/stat/save              /home/gate/Downloads/vgate/dosimetry/output_txt/stat-C100.txt
/gate/actor/stat/saveEveryNSeconds 10

# filtter====================

/gate/actor/prdandstp_11/addFilter                     particleFilter
/gate/actor/prdandstp_11/particleFilter/addParticleZ      6
/gate/actor/prdandstp_11/particleFilter/addParticleA      11

/gate/actor/prdandstp_10/addFilter                     particleFilter
/gate/actor/prdandstp_10/particleFilter/addParticleZ      6
/gate/actor/prdandstp_10/particleFilter/addParticleA      10

/gate/actor/prdandstp_15/addFilter                     particleFilter
/gate/actor/prdandstp_15/particleFilter/addParticleZ      8
/gate/actor/prdandstp_15/particleFilter/addParticleA      15

/gate/actor/prdandstp_11_txt/addFilter                     particleFilter
/gate/actor/prdandstp_11_txt/particleFilter/addParticleZ      6
/gate/actor/prdandstp_11_txt/particleFilter/addParticleA      11

/gate/actor/prdandstp_10_txt/addFilter                     particleFilter
/gate/actor/prdandstp_10_txt/particleFilter/addParticleZ      6
/gate/actor/prdandstp_10_txt/particleFilter/addParticleA      10

/gate/actor/prdandstp_15_txt/addFilter                     particleFilter
/gate/actor/prdandstp_15_txt/particleFilter/addParticleZ      8
/gate/actor/prdandstp_15_txt/particleFilter/addParticleA      15
# INITIALISATION============

/gate/run/initialize
#/gate/physics/processList Available
#/gate/physics/processList Enabled


# BEAMS==================

/gate/source/addSource     mybeam  PencilBeam
/gate/source/mybeam/setParticleType GenericIon
/gate/source/mybeam/setIonProperties 6 12 6 0
/gate/source/mybeam/setEnergy 2400 MeV
/gate/source/mybeam/setSigmaEnergy 2 MeV
/gate/source/mybeam/setPosition 0 0 0 mm
/gate/source/mybeam/setSigmaX 2 mm
/gate/source/mybeam/setSigmaY 2 mm
/gate/source/mybeam/setSigmaTheta  0.01 mrad
/gate/source/mybeam/setSigmaPhi    0.01 mrad
/gate/source/mybeam/setEllipseXThetaEmittance 0.01 mm*mrad
/gate/source/mybeam/setEllipseYPhiEmittance   0.01 mm*mrad
#/gate/source/mybeam/setIntensity 100000000


# JamesRandom Ranlux64 MersenneTwister
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto

/gate/application/noGlobalOutput
#/gate/application/readTimeSlicesIn               data/timeslice.dat
#/gate/application/readNumberOfPrimariesInAFile   data/primary.dat
/gate/application/setTotalNumberOfPrimaries 500000
/gate/application/start
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.opengatecollaboration.org/pipermail/gate-users/attachments/20210822/cebcb32f/attachment-0001.html>