class ACB5_2012_6X(TestTRS398Photon, TestCase):
temp = 21.7
press = trs398.mmHg2kPa(757.2)
nd_w = 5.446
m_reference = 25.27
m_opposite = -25.19
m_reduced = 25.17
clinical_pdd_zref = 66.8
tpr2010 = trs398.tpr2010_from_pdd2010(pdd2010=38.4/66.8)
tissue_correction = 0.99
dose_mu_zref = 0.679
dose_mu_zmax = 1.0159
class IMMCTB_15X(TestTRS398Photon, TestCase):
energy = 15
temp = 22.4
press = trs398.mmHg2kPa(748.1)
nd_w = 5.394
m_reference = 14.307
m_opposite = -14.323
m_reduced = 14.220
clinical_pdd_zref = 76.7
tpr2010 = trs398.tpr2010_from_pdd2010(pdd2010=(49.9/76.9)) * (76.79/76.9)
mu = 100
dose_mu_zref = 0.770
dose_mu_zmax = 1.004
class ACB5_2012_18X(TestTRS398Photon, TestCase):
energy = 18
temp = 21.7
press = trs398.mmHg2kPa(757.2)
tpr2010 = trs398.tpr2010_from_pdd2010(pdd2010=52.5/79.4)
nd_w = 5.446
m_reference = 30.67
m_opposite = -30.65
m_reduced = 30.50
clinical_pdd_zref = 79.7
tissue_correction = 0.99
dose_mu_zref = 0.807
dose_mu_zmax = 1.0125
class MDA_TB1_2015_10x(TestTRS398Photon, TestCase):
energy = 10
temp = 21
press = trs398.mmHg2kPa(763)
nd_w = 5.393
k_elec = 1.003
m_reference = 27.727
m_opposite = 27.784
m_reduced = 27.635
clinical_pdd_zref = 73.5
dose_mu_zref = 0.734
dose_mu_zmax = 0.998
tpr2010 = (73.42/73.7) * trs398.tpr2010_from_pdd2010(pdd2010=46.3/73.7)
class MDA_TB2_2015_15x(TestTRS398Photon, TestCase):
energy = 15
temp = 20.5
press = trs398.mmHg2kPa(760)
nd_w = 5.444
k_elec = 1.002
m_reference = 29.28
m_opposite = -29.33
m_reduced = 29.10
dose_mu_zref = 0.779
dose_mu_zmax = 1.007
clinical_pdd_zref = 77.4
tpr2010 = 0.762
class IMMCTB_10FFF(TestTRS398Photon, TestCase):
energy = 10
fff = True
temp = 22.4
press = trs398.mmHg2kPa(748.1)
nd_w = 5.394
m_reference = 13.00067
m_opposite = -13.013
m_reduced = 12.867
tpr2010 = trs398.tpr2010_from_pdd2010(pdd2010=(43/71.2))
clinical_pdd_zref = 71.1
mu = 100
dose_mu_zref = 0.712
dose_mu_zmax = 1.0005
class IMMC_TB_20E(TestTRS398Electron, TestCase):
energy = 20
cone = '15x15'
mu = 100
temp = 22.1
press = trs398.mmHg2kPa(748.2)
k_elec = 0.999
nd_w = 5.394
m_reference = 19.670 * 0.99354 # * Pgradient because TRS398 is done at dref+0.5cm
m_opposite = 19.707 * 0.99354
m_reduced = 19.437 * 0.99354
i_50 = 8.22
clinical_pdd_zref = 96.8
dose_mu_zref = 0.972
dose_mu_zmax = 1.004
class IMMCTB_6FFF(TestTRS398Photon, TestCase):
energy = 6
fff = True
temp = 22.5
press = trs398.mmHg2kPa(749)
tpr2010 = (64.16 / 63.6) * trs398.tpr2010_from_pdd2010(pdd2010=34.5 / 63.6)
nd_w = 5.394
m_reference = 11.610
m_opposite = -11.613
m_reduced = 11.533
clinical_pdd_zref = 63.5
mu = 100
dose_mu_zref = 0.638
dose_mu_zmax = 1.005
print_data = True
class ACB5_2011_6x(TestTRS398Photon, TestCase):
temp = 22
press = trs398.mmHg2kPa(751.2)
nd_w = 5.450
tpr2010 = trs398.tpr2010_from_pdd2010(pdd2010=38.4/66.8)
m_reference = 24.82
m_opposite = -24.83
m_reduced = 24.79
clinical_pdd_zref = 66.8
tissue_correction = 0.99
dose_mu_zref = 0.673
dose_mu_zmax = 1.007
def test_zmax_adjusted(self):
self.trs398.m_reference_adjusted = 24.65
self.assertAlmostEqual(self.trs398.dose_mu_zmax_adjusted, 1.000, delta=0.0005)
def test_zref_adjusted(self):
self.trs398.m_reference_adjusted = 24.65
self.assertAlmostEqual(self.trs398.dose_mu_zref_adjusted, 0.668, delta=0.0005)
"""A script to calculate TRS-398 dose using pylinac classes and following the TRS-398 photon form"""
from pylinac.calibration import trs398
ENERGY = 6
TEMP = 22.1
PRESS = trs398.mmHg2kPa(755.0)
CHAMBER = '30013' # PTW
K_ELEC = 1.000
ND_w = 5.443 # Gy/nC
MU = 200
CLINICAL_PDD = 66.5
trs398_6x = trs398.TRS398Photon(
unit='TrueBeam1',
setup='SSD',
chamber=CHAMBER,
temp=TEMP, press=PRESS,
n_dw=ND_w,
clinical_pdd_zref=CLINICAL_PDD,
tpr2010=(38.2/66.6),
energy=ENERGY,
fff=False,
k_elec=K_ELEC,
voltage_reference=-300, voltage_reduced=-150,
m_reference=(25.65, 25.66, 25.65),
m_opposite=(25.64, 25.65, 25.65),
m_reduced=(25.64, 25.63, 25.63),
mu=MU, tissue_correction=1.0
)
"""A script to calculate TRS-398 dose using pylinac functions and following the TRS-398 photon form"""
from pylinac.calibration import trs398
TEMP = 22.1
PRESS = trs398.mmHg2kPa(755.0)
CHAMBER = '30013' # PTW
K_ELEC = 1.000
ND_w = 5.443 # Gy/nC
MU = 200
# Section 3 (dosimeter corrections)
k_tp = trs398.k_tp(temp=TEMP, press=PRESS)
k_pol = trs398.k_pol(m_reference=(25.66, 25.67, 25.66),
m_opposite=(25.65, 25.66, 25.66))
k_s = trs398.k_s(voltage_reference=300,
voltage_reduced=150,
m_reference=(25.66, 25.67, 25.66),
m_reduced=(25.63, 25.65, 25.64))
m_corrected = trs398.m_corrected(m_reference=(25.66, 25.67, 25.66),
k_tp=k_tp, k_elec=K_ELEC, k_pol=k_pol, k_s=k_s) \
/ MU
# Section 4 (kQ + dose at zref)
kq = trs398.kq_photon(chamber=CHAMBER, tpr=(39.2 / 68.1))
dose_mu_zref = m_corrected * ND_w * kq
# Section 5 (Dose at zmax)
# SSD setup
CLINICAL_PDD = 66.5
dose_mu_zmax = dose_mu_zref * 100 / CLINICAL_PDD