<p>
        Dear Gate Users:
</p>
<p>
        I previously used Pencil Beam Source to simulate single-energy pencil beam for carbon-ion therapy. The features of the pencil beam source such as spot size, divergence and beam emittance were defined as follows: 3 mm spot size in both x and y directions, 4 mrad divergence and 15 mm·mrad beam emittance. The energy spread was set to 0.1% of the respective energy. Except for dipole scanning magnets,vacuum pipe and window, main and 
sub-dose monitors, a multi-strip ionization chamber (IC) and a mini 
ridge filter were modeled in my MC simulations. The source position was located before the vacuum window and a water box was located at isocenter.
</p>
<p>
        #################Pencil Beam Source########################################
</p>
<p>
        # SOURCE PBS<br>
/gate/source/addSource C12 PencilBeam<br>
<br>
#particle type<br>
/gate/source/C12/setParticleType GenericIon<br>
/gate/source/C12/setIonProperties 6 12 6 0<br>
<br>
#energy spectrum<br>
/gate/source/C12/setEnergy 2687.52 MeV<br>
/gate/source/C12/setSigmaEnergy 2.68752 MeV<br>
<br>
#position<br>
/gate/source/C12/setPosition 0 0 -1.4059 m<br>
<br>
# spot size<br>
/gate/source/C12/setSigmaX 3 mm<br>
/gate/source/C12/setSigmaY 3 mm<br>
<br>
#divergence<br>
/gate/source/C12/setSigmaTheta 4 mrad<br>
/gate/source/C12/setSigmaPhi 4 mrad<br>
<br>
#emittance<br>
/gate/source/C12/setEllipseXThetaEmittance 15 mm*mrad<br>
/gate/source/C12/setEllipseYPhiEmittance 15 mm*mrad<br>
<br>
#positive-> converging beam; negative diverging beam<br>
/gate/source/C12/setEllipseXThetaRotationNorm positive<br>
/gate/source/C12/setEllipseYPhiRotationNorm positive
</p>
<p>
        <br>
</p>
<p>
        Now, I want to use TPS Pencil Beam Source to simulate realistic treatment plan for carbon ion therapy. In Source-properties.txt file, it contains 8 polynomial equations, each equation being a function of the energy at the nozzle entrance.
</p>
<p>
        #################TPS Pencil Beam Source-Source-properties########################################
</p>
<p>
        # Source properties<br>
<br>
# Nozzle exit to Isocenter distance<br>
1405.9<br>
# SMX to Isocenter distance<br>
4000<br>
# SMY to Isocenter distance<br>
4000<br>
<br>
# ENERGY<br>
# mean energy<br>
# polynomial order<br>
1<br>
# polynomial parameters (highest to lowest)<br>
1<br>
0<br>
<br>
# energy spread<br>
# polynomial order<br>
percent<br>
1<br>
# polynomial parameters (highest to lowest)<br>
0.001<br>
0<br>
<br>
# SPOT x y theta phi<br>
# x<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
3<br>
# Theta<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
0.004<br>
# y<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
3<br>
# Phi<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
0.004<br>
<br>
# SPOT emmittances<br>
# x-theta emmittance<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
0.015<br>
# y-phi emmittance<br>
# polynomial order<br>
1<br>
# polynomial parameters<br>
0<br>
0.015
</p>
<p>
        <br>
</p>
<p>
        #################TPS Pencil Beam Source-PlanDescriptionToGate########################################
</p>
<p>
        #TREATMENT-PLAN-DESCRIPTION<br>
#PlanName<br>
IMPT<br>
#NumberOfFractions            ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
1<br>
##FractionID                ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
1<br>
##NumberOfFields<br>
1<br>
###FieldsID<br>
1<br>
<br>
#TotalMetersetWeightOfAllFields        ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
662421<br>
<br>
#FIELD-DESCRIPTION<br>
###FieldID<br>
1<br>
###FinalCumulativeMeterSetWeight    ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
336732<br>
###GantryAngle (in degrees)<br>
0<br>
###PatientSupportAngle            ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
0<br>
###IsocenterPosition            ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
0 0 0<br>
###NumberOfControlPoints<br>
1<br>
<br>
#SPOTS-DESCRIPTION<br>
####ControlPointIndex<br>
0<br>
####SpotTunnedID            ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
0<br>
####CumulativeMetersetWeight        ***    CURRENTLY THIS VALUE IS NOT USED    ***<br>
0<br>
####Energy (MeV)<br>
2687.52<br>
####NbOfScannedSpots<br>
1<br>
####X Y Weight (spot position at isocenter in mm, with weight in MU (default) or number of ions "setSpotIntensityAsNbIons true")<br>
0 0 100<br>
</p>
<p>
        I simulated the same energy beam in my PlanDescriptionToGate.txt file and the spot position is at isocenter. However, I find that the dose distributions with these two sources in water box were different. I'm puzzled about the parameter settings in TPS Pencil Beam Source if I want to make them consistent with those in Pencil Beam Source settings. How can I do?
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        Best regards,
</p>
<p>
        Hui Zhang
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<p>
        <br>
</p>
<span class="spnEditorSign">
<hr class="signature-separator" style="margin:0.5em 0;width:10em;height:1px;background-color:#999;border:none;" align="left">
</span>