Dear experts
Recently i am design an carbon pencil beam scanning model in GATE tools. However, an cross-section error happens . Here is the detail
The aim of this project is to establish an carbon pencil beam scanning model in gate simulation, so i designed an pencil beam first. The code runs well at first when i used GATE singel pencil beam commond setting, while when i change the beam setting in plan description&source property file the code will crash and remind "No isotope cross section found for C12 off Element Nitrogen ".
To avoid the code structure error, i change the ion typr from genericion(C-ion) to proton, the code runs well again. Then i use a lower energy in carbon(10MeV/U).Code also runs without error information. It realy confuse me, would u mind help me fix this problem? Thanks a lot for your kind help.Here is the code for check
main.mac
#=====================================================
VERBOSITY
#=====================================================
#/control/execute mac/verbose.mac
#=====================================================
VISUALISATION
#=====================================================
#/control/execute mac/visu.mac
=====================================================
GEOMETRY
=====================================================
/gate/geometry/setMaterialDatabase /home/vgate/Desktop/3002/myGateMaterials.db
WORLD
/gate/world/setMaterial Vacuum
/gate/world/geometry/setXLength 2 m
/gate/world/geometry/setYLength 20 m
/gate/world/geometry/setZLength 2 m
/gate/world/vis/setVisible 0
LINAC HEAD
World origin is in the center of the beam e- source (not the real
isocenter of the gantry room)
/gate/world/daughters/name waterbox
/gate/world/daughters/insert box
/gate/waterbox/geometry/setXLength 100. mm
/gate/waterbox/geometry/setYLength 300 mm
/gate/waterbox/geometry/setZLength 100. mm
/gate/waterbox/placement/setTranslation 0 -332 0 mm
/gate/waterbox/setMaterial Water
/gate/waterbox/vis/setVisible 1
/gate/waterbox/vis/setColor blue
=====================================================
PHYSICS
=====================================================
/gate/physics/addPhysicsList QGSP_BIC_EMY
/gate/physics/Gamma/SetCutInRegion world 10 m
/gate/physics/Electron/SetCutInRegion world 10 m
/gate/physics/Positron/SetCutInRegion world 10 m
/gate/physics/SetMaxStepSizeInRegion waterbox 1 mm
/gate/physics/ActivateStepLimiter proton
/gate/physics/ActivateStepLimiter deuteron
/gate/physics/ActivateStepLimiter triton
/gate/physics/ActivateStepLimiter alpha
/gate/physics/ActivateStepLimiter GenericIon
/gate/physics/displayCuts
#=====================================================
DETECTORS
#=====================================================
----------------------------------------------------
/gate/actor/addActor DoseActor doseprofile
/gate/actor/doseprofile/save output/output.txt
/gate/actor/doseprofile/attachTo waterbox
/gate/actor/doseprofile/stepHitType random
/gate/actor/doseprofile/setPosition 0 0 0 cm
/gate/actor/doseprofile/setResolution 1 3000 1
/gate/actor/doseprofile/setSize 10 30 10 cm
/gate/actor/doseprofile/saveEveryNSeconds 60
/gate/actor/doseprofile/enableEdep true
/gate/actor/doseprofile/enableDose true
#=====================================================
INITIALISATION
#=====================================================
/gate/run/initialize
#=====================================================
BEAMS
#=====================================================
#/gate/source/addSource PBS PencilBeam
#/gate/source/PBS/setParticleType GenericIon
#/gate/source/PBS/setIonProperties 6 12 6 0
#/gate/source/PBS/setEnergy 1443.12 MeV
#/gate/source/PBS/setSigmaEnergy 39.17652707 MeV
#/gate/source/PBS/setPosition 0 -7000 0 mm
#/gate/source/PBS/setSigmaX 1.2 mm
#/gate/source/PBS/setSigmaY 1.2 mm
#/gate/source/PBS/setSigmaTheta 1 mrad
#/gate/source/PBS/setSigmaPhi 1 mrad
#/gate/source/PBS/setEllipseXThetaEmittance 0.15 mm*mrad
#/gate/source/PBS/setEllipseXThetaRotationNorm negative
#/gate/source/PBS/setEllipseYPhiEmittance 0.15 mm*mrad
#/gate/source/PBS/setEllipseYPhiRotationNorm negative
#/gate/source/PBS/setRotationAxis 1 0 0
#/gate/source/PBS/setRotationAngle 270 deg
/gate/source/addSource PBS TPSPencilBeam
#/gate/source/PBS/setParticleType GenericIon
/gate/source/PBS/setIonProperties 6 12 6 0
/gate/source/PBS/setTestFlag false
/gate/source/PBS/setPlan data/PlanDescriptionToGATE.txt
#/gate/source/PBS/setNotAllowedFieldID 1
/gate/source/PBS/setFlatGenerationFlag false
/gate/source/PBS/setSpotIntensityAsNbIons true
/gate/source/PBS/setSortedSpotGenerationFlag false
/gate/source/PBS/setSigmaEnergyInMeVFlag false
/gate/source/PBS/setSourceDescriptionFile data/Source-Properties.txt
#=====================================================
START BEAMS
#=====================================================
/gate/random/setEngineName MersenneTwister
/gate/random/setEngineSeed auto
/gate/application/noGlobalOutput
/gate/application/setTotalNumberOfPrimaries 100000
/gate/application/start
PlanDescriptionToGATE.txt
#TREATMENT-PLAN-DESCRIPTION
#PlanName
IMCT
#NumberOfFractions *** CURRENTLY THIS VALUE IS NOT USED ***
1
##FractionID *** CURRENTLY THIS VALUE IS NOT USED ***
1
##NumberOfFields
1
###FieldsID
1
#TotalMetersetWeightOfAllFields *** CURRENTLY THIS VALUE IS NOT USED ***
662421
#FIELD-DESCRIPTION
###FieldID
1
###FinalCumulativeMeterSetWeight *** CURRENTLY THIS VALUE IS NOT USED ***
336732
###GantryAngle (in degrees)
0
###PatientSupportAngle *** CURRENTLY THIS VALUE IS NOT USED ***
0
###IsocenterPosition *** CURRENTLY THIS VALUE IS NOT USED ***
0 -638 0
###NumberOfControlPoints
1
#SPOTS-DESCRIPTION
####ControlPointIndex
0
####SpotTunnedID *** CURRENTLY THIS VALUE IS NOT USED ***
0
####CumulativeMetersetWeight *** CURRENTLY THIS VALUE IS NOT USED ***
0
####Energy (MeV)
4000
####NbOfScannedSpots
1
####X Y Weight (spot position at isocenter in mm, with weight in MU (default) or number of ions “setSpotIntensityAsNbIons true”)
0 0 10000
Source property.txt
Source properties
Nozzle exit to Isocenter distance
638
SMX to Isocenter distance
2000
SMY to Isocenter distance
2000
ENERGY
mean energy
polynomial order
1
polynomial parameters (highest to lowest)
1
0
energy spread
polynomial order
percent
1
polynomial parameters (highest to lowest)
0.028
0
SPOT x y theta phi
x
polynomial order
1
polynomial parameters
0
1.2
Theta
polynomial order
1
polynomial parameters
0
0.001
y
polynomial order
1
polynomial parameters
0
1.2
Phi
polynomial order
1
polynomial parameters
0
0.001
SPOT emmittances
x-theta emmittance
polynomial order
1
polynomial parameters
0
0.00015
y-phi emmittance
polynomial order
1
polynomial parameters
0
0.00015