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