Changeset 1337 for trunk/examples/extended/medical/GammaTherapy
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/medical/GammaTherapy
- Files:
-
- 23 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/medical/GammaTherapy/GNUmakefile
r1230 r1337 1 # $Id: GNUmakefile,v 1. 3 2008/06/11 15:17:57 vnivanchExp $1 # $Id: GNUmakefile,v 1.4 2010/01/11 16:43:05 gcosmo Exp $ 2 2 # -------------------------------------------------------------- 3 3 # GNUmakefile for examples module. Gabriele Cosmo, 06/04/98. … … 9 9 10 10 ifndef G4INSTALL 11 G4INSTALL = ../ geant411 G4INSTALL = ../../../.. 12 12 endif 13 13 -
trunk/examples/extended/medical/GammaTherapy/GammaTherapy.cc
r807 r1337 91 91 if (argc==1) // Define UI terminal for interactive mode 92 92 { 93 visManager = new G4VisExecutive( );93 visManager = new G4VisExecutive("Quiet"); 94 94 visManager->Initialize(); 95 95 G4UIsession * session; -
trunk/examples/extended/medical/GammaTherapy/GammaTherapy.in
r807 r1337 13 13 /testem/phys/addPhysics emstandard 14 14 #/testem/phys/addPhysics emstandard_opt2 15 #/testem/phys/addPhysics penelope16 #/testem/phys/addPhysics livermore15 #/testem/phys/addPhysics empenelope 16 #/testem/phys/addPhysics emlivermore 17 17 #/testem/phys/addPhysics step_limit 18 18 # -
trunk/examples/extended/medical/GammaTherapy/GammaTherapy.out
r807 r1337 1 2 ############################################ 3 !!! WARNING - FPE detection is activated !!! 4 ############################################ 1 5 2 6 ************************************************************* 3 Geant4 version Name: g lobal-V09-00-03 (9-May-2008)7 Geant4 version Name: geant4-09-03-ref-06 (25-June-2010) 4 8 Copyright : Geant4 Collaboration 5 9 Reference : NIM A 506 (2003), 250-303 … … 10 14 ***** Table : Nb of materials = 6 ***** 11 15 12 Material: G4_Be density: 1.848 g/cm3 RadL: 35.276 cm Imean: 63.700 eV16 Material: G4_Be density: 1.848 g/cm3 RadL: 35.276 cm Nucl.Int.Length: 39.449 cm Imean: 63.700 eV 13 17 ---> Element: Be (Be) Z = 4.0 N = 9.0 A = 9.01 g/mole 14 ---> Isotope: BeZ = 4 N = 9 A = 9.01 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %15 16 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Imean: 727.000 eV18 ---> Isotope: Be9 Z = 4 N = 9 A = 9.01 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 19 20 Material: G4_W density: 19.300 g/cm3 RadL: 3.504 mm Nucl.Int.Length: 10.310 cm Imean: 727.000 eV 17 21 ---> Element: W (W) Z = 74.0 N = 183.9 A = 183.84 g/mole 18 ---> Isotope: WZ = 74 N = 180 A = 179.95 g/mole abundance: 0.12 %19 ---> Isotope: WZ = 74 N = 182 A = 181.95 g/mole abundance: 26.50 %20 ---> Isotope: WZ = 74 N = 183 A = 182.95 g/mole abundance: 14.31 %21 ---> Isotope: WZ = 74 N = 184 A = 183.95 g/mole abundance: 30.64 %22 ---> Isotope: WZ = 74 N = 186 A = 185.95 g/mole abundance: 28.43 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %23 24 Material: G4_He density: 0.166 mg/cm3 RadL: 5.671km Imean: 41.800 eV temperature: 273.15 K pressure: 1.00 atm22 ---> Isotope: W180 Z = 74 N = 180 A = 179.95 g/mole abundance: 0.12 % 23 ---> Isotope: W182 Z = 74 N = 182 A = 181.95 g/mole abundance: 26.50 % 24 ---> Isotope: W183 Z = 74 N = 183 A = 182.95 g/mole abundance: 14.31 % 25 ---> Isotope: W184 Z = 74 N = 184 A = 183.95 g/mole abundance: 30.64 % 26 ---> Isotope: W186 Z = 74 N = 186 A = 185.95 g/mole abundance: 28.43 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 27 28 Material: G4_He density: 0.166 kg/m3 RadL: 5.671 km Nucl.Int.Length: 3.343 km Imean: 41.800 eV temperature: 273.15 K pressure: 1.00 atm 25 29 ---> Element: He (He) Z = 2.0 N = 4.0 A = 4.00 g/mole 26 ---> Isotope: HeZ = 2 N = 3 A = 3.02 g/mole abundance: 0.00 %27 ---> Isotope: HeZ = 2 N = 4 A = 4.00 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 %28 29 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Imean: 75.000 eV30 ---> Isotope: He3 Z = 2 N = 3 A = 3.02 g/mole abundance: 0.00 % 31 ---> Isotope: He4 Z = 2 N = 4 A = 4.00 g/mole abundance: 100.00 % ElmMassFraction: 100.00 % ElmAbundance 100.00 % 32 33 Material: G4_WATER H_2O density: 1.000 g/cm3 RadL: 36.083 cm Nucl.Int.Length: 75.505 cm Imean: 78.000 eV 30 34 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 31 ---> Isotope: HZ = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %32 ---> Isotope: HZ = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 11.19 % ElmAbundance 66.67 %35 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 36 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 11.19 % ElmAbundance 66.67 % 33 37 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 34 ---> Isotope: OZ = 8 N = 16 A = 15.99 g/mole abundance: 99.76 %35 ---> Isotope: OZ = 8 N = 17 A = 17.00 g/mole abundance: 0.04 %36 ---> Isotope: OZ = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 %37 38 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm38 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 39 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 40 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 88.81 % ElmAbundance 33.33 % 41 42 Material: G4_AIR density: 1.205 mg/cm3 RadL: 303.921 m Nucl.Int.Length: 710.137 m Imean: 85.700 eV temperature: 273.15 K pressure: 1.00 atm 39 43 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 40 ---> Isotope: CZ = 6 N = 12 A = 12.00 g/mole abundance: 98.93 %41 ---> Isotope: CZ = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 0.01 % ElmAbundance 0.02 %44 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 45 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 0.01 % ElmAbundance 0.02 % 42 46 ---> Element: N (N) Z = 7.0 N = 14.0 A = 14.01 g/mole 43 ---> Isotope: NZ = 7 N = 14 A = 14.00 g/mole abundance: 99.63 %44 ---> Isotope: NZ = 7 N = 15 A = 15.00 g/mole abundance: 0.37 % ElmMassFraction: 75.53 % ElmAbundance 78.44 %47 ---> Isotope: N14 Z = 7 N = 14 A = 14.00 g/mole abundance: 99.63 % 48 ---> Isotope: N15 Z = 7 N = 15 A = 15.00 g/mole abundance: 0.37 % ElmMassFraction: 75.53 % ElmAbundance 78.44 % 45 49 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 46 ---> Isotope: OZ = 8 N = 16 A = 15.99 g/mole abundance: 99.76 %47 ---> Isotope: OZ = 8 N = 17 A = 17.00 g/mole abundance: 0.04 %48 ---> Isotope: OZ = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 23.18 % ElmAbundance 21.07 %50 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 51 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 52 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 23.18 % ElmAbundance 21.07 % 49 53 ---> Element: Ar (Ar) Z = 18.0 N = 40.0 A = 39.95 g/mole 50 ---> Isotope: ArZ = 18 N = 36 A = 35.97 g/mole abundance: 0.34 %51 ---> Isotope: ArZ = 18 N = 38 A = 37.96 g/mole abundance: 0.06 %52 ---> Isotope: ArZ = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 1.28 % ElmAbundance 0.47 %53 54 Material: G4_MYLAR density: 1.400 g/cm3 RadL: 28.536 cm Imean: 78.700 eV54 ---> Isotope: Ar36 Z = 18 N = 36 A = 35.97 g/mole abundance: 0.34 % 55 ---> Isotope: Ar38 Z = 18 N = 38 A = 37.96 g/mole abundance: 0.06 % 56 ---> Isotope: Ar40 Z = 18 N = 40 A = 39.96 g/mole abundance: 99.60 % ElmMassFraction: 1.28 % ElmAbundance 0.47 % 57 58 Material: G4_MYLAR density: 1.400 g/cm3 RadL: 28.536 cm Nucl.Int.Length: 55.961 cm Imean: 78.700 eV 55 59 ---> Element: H (H) Z = 1.0 N = 1.0 A = 1.01 g/mole 56 ---> Isotope: HZ = 1 N = 1 A = 1.01 g/mole abundance: 99.99 %57 ---> Isotope: HZ = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 4.20 % ElmAbundance 36.36 %60 ---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.99 % 61 ---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.01 % ElmMassFraction: 4.20 % ElmAbundance 36.36 % 58 62 ---> Element: C (C) Z = 6.0 N = 12.0 A = 12.01 g/mole 59 ---> Isotope: CZ = 6 N = 12 A = 12.00 g/mole abundance: 98.93 %60 ---> Isotope: CZ = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 62.50 % ElmAbundance 45.46 %63 ---> Isotope: C12 Z = 6 N = 12 A = 12.00 g/mole abundance: 98.93 % 64 ---> Isotope: C13 Z = 6 N = 13 A = 13.00 g/mole abundance: 1.07 % ElmMassFraction: 62.50 % ElmAbundance 45.46 % 61 65 ---> Element: O (O) Z = 8.0 N = 16.0 A = 16.00 g/mole 62 ---> Isotope: OZ = 8 N = 16 A = 15.99 g/mole abundance: 99.76 %63 ---> Isotope: OZ = 8 N = 17 A = 17.00 g/mole abundance: 0.04 %64 ---> Isotope: OZ = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 33.30 % ElmAbundance 18.18 %66 ---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.76 % 67 ---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.04 % 68 ---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.20 % ElmMassFraction: 33.30 % ElmAbundance 18.18 % 65 69 66 70 … … 78 82 ### PhysicsList Add Physics <emstandard> emBuilderIsRegisted= 0 79 83 PhysicsList::AddPhysicsList <emstandard> 80 #/testem/phys/addPhysics emstandard_opt281 #/testem/phys/addPhysics penelope82 #/testem/phys/addPhysics livermore83 #/testem/phys/addPhysics step_limit84 #/testem/phys/addPhysics emstandard_opt2 85 #/testem/phys/addPhysics empenelope 86 #/testem/phys/addPhysics emlivermore 87 #/testem/phys/addPhysics step_limit 84 88 # 85 89 /testem/target1Material G4_Be … … 90 94 /testem/target2Z 6 mm 91 95 # 92 #/testem/histoName be_w_st.paw93 # 94 #/testem/stepMax 1.0mm96 #/testem/histoName be_w_st.paw 97 # 98 #/testem/stepMax 1.0 mm 95 99 # 96 100 /gun/particle e- … … 126 130 Index : 0 used in the geometry : Yes recalculation needed : No 127 131 Material : G4_AIR 128 Range cuts : gamma 100 um e- 100 um e+ 100 um129 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV132 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 133 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV 130 134 Region(s) which use this couple : 131 135 DefaultRegionForTheWorld … … 133 137 Index : 1 used in the geometry : Yes recalculation needed : No 134 138 Material : G4_Be 135 Range cuts : gamma 100 um e- 100 um e+ 100 um136 Energy thresholds : gamma 1.02934 keV e- 107.795 keV e+ 106.476 keV139 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 140 Energy thresholds : gamma 1.05172 keV e- 109.289 keV e+ 107.83 keV proton 0 eV 137 141 Region(s) which use this couple : 138 142 DefaultRegionForTheWorld … … 140 144 Index : 2 used in the geometry : Yes recalculation needed : No 141 145 Material : G4_He 142 Range cuts : gamma 100 um e- 100 um e+ 100 um143 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV146 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 147 Energy thresholds : gamma 990 eV e- 990 eV e+ 990 eV proton 0 eV 144 148 Region(s) which use this couple : 145 149 DefaultRegionForTheWorld … … 147 151 Index : 3 used in the geometry : Yes recalculation needed : No 148 152 Material : G4_W 149 Range cuts : gamma 100 um e- 100 um e+ 100 um150 Energy thresholds : gamma 36.0479 keV e- 351.438 keV e+ 334.551 keV153 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 154 Energy thresholds : gamma 36.0599 keV e- 351.877 keV e+ 337.972 keV proton 0 eV 151 155 Region(s) which use this couple : 152 156 DefaultRegionForTheWorld … … 154 158 Index : 4 used in the geometry : Yes recalculation needed : No 155 159 Material : G4_MYLAR 156 Range cuts : gamma 100 um e- 100 um e+ 100 um157 Energy thresholds : gamma 1.19365 keV e- 98.2879 keV e+ 96.4896 keV160 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 161 Energy thresholds : gamma 1.21113 keV e- 99.4755 keV e+ 98.1475 keV proton 0 eV 158 162 Region(s) which use this couple : 159 163 DefaultRegionForTheWorld … … 161 165 Index : 5 used in the geometry : Yes recalculation needed : No 162 166 Material : G4_WATER 163 Range cuts : gamma 100 um e- 100 um e+ 100 um164 Energy thresholds : gamma 1.09571 keV e- 84.2696 keV e+ 83.2385 keV167 Range cuts : gamma 100 um e- 100 um e+ 100 um proton 0 fm 168 Energy thresholds : gamma 1.10981 keV e- 84.6589 keV e+ 83.5287 keV proton 0 eV 165 169 Region(s) which use this couple : 166 170 DefaultRegionForTheWorld … … 178 182 Run Summary 179 183 Number of events processed : 10000 180 User= 14.88s Real=14.89s Sys=0.01s184 User=8.44s Real=10.22s Sys=0s 181 185 RunAction: End of run actions are started 182 186 Histo: End of run actions are started 183 187 ======================================================== 184 188 Number of events 10000 185 Average number of e- 12.2 6186 Average number of gamma 21. 05187 Average number of e+ 0.7 483189 Average number of e- 12.24 190 Average number of gamma 21.28 191 Average number of e+ 0.7963 188 192 Average number of neutrons 0 189 Average number of steps in absorber 3 3.78190 Average number of e- steps in target 18. 76191 Average number of g produced in the target 20. 18192 Average number of e- produced in the target 7. 089193 Average number of g produced in the phantom 0.5 544194 Average number of e- produced in the phantom 2. 927195 Total gamma fluence in front of phantom 9.022MeV193 Average number of steps in absorber 31.91 194 Average number of e- steps in target 18.67 195 Average number of g produced in the target 20.44 196 Average number of e- produced in the target 7.238 197 Average number of g produced in the phantom 0.5203 198 Average number of e- produced in the phantom 2.795 199 Total gamma fluence in front of phantom 8.671 MeV 196 200 ======================================================== 197 201 198 202 199 203 G4 kernel has come to Quit state. 204 -
trunk/examples/extended/medical/GammaTherapy/History
r807 r1337 1 $Id: History,v 1.2 5 2007/07/04 11:14:33 vnivanchExp $1 $Id: History,v 1.29 2010/06/07 05:40:46 perl Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 15 15 * Reverse chronological order (last date on top), please * 16 16 ---------------------------------------------------------- 17 18 06-06-10 J.Perl (gtherapy-V09-03-01) 19 - Remove unused variable in EventAction 20 21 04-06-10 J.Perl (gtherapy-V09-03-00) 22 - Updated vis usage 23 24 21-11-09 V.Ivant (gtherapy-V09-02-00) 25 - update EM physics according to 9.3 26 27 11-06-08 V.Ivant (gtherapy-V09-01-00) 28 - Remove AIDA from GNUmakefile 29 - Added opt3 variant of EM physics 17 30 18 31 04-07-07 V.Ivant (gtherapy-V09-00-00) -
trunk/examples/extended/medical/GammaTherapy/README
r807 r1337 42 42 Cu W target). 43 43 44 Following Physics Lists can be activated via UI commands: 44 Only Physics Lists from physics_list kernal library can be activated. 45 Corresponding UI commands are following 45 46 46 /testem/phys/addPhysics standard47 /testem/phys/addPhysics g4v5248 /testem/phys/addPhysics penelope49 /testem/phys/addPhysics lowenergy47 /testem/phys/addPhysics emstandard 48 /testem/phys/addPhysics emstandard_opt3 49 /testem/phys/addPhysics empenelope 50 /testem/phys/addPhysics emlivermore 50 51 51 52 For interactive mode G4 visualization options and variables should be -
trunk/examples/extended/medical/GammaTherapy/be.in
r807 r1337 11 11 /testem/phys/addPhysics emstandard 12 12 #/testem/phys/addPhysics emstandard_opt2 13 #/testem/phys/addPhysics livermore13 #/testem/phys/addPhysics emlivermore 14 14 #/testem/phys/addPhysics elastic 15 15 #/testem/phys/addPhysics binary -
trunk/examples/extended/medical/GammaTherapy/be_w.in
r807 r1337 13 13 /testem/phys/addPhysics emstandard 14 14 #/testem/phys/addPhysics emstandard_opt2 15 #/testem/phys/addPhysics livermore15 #/testem/phys/addPhysics emlivermore 16 16 #/testem/phys/addPhysics elastic 17 17 #/testem/phys/addPhysics binary … … 31 31 /gun/energy 50.0 MeV 32 32 # 33 /run/initialize34 /run/beamOn 100000033 #/run/initialize 34 #/run/beamOn 1000000 35 35 -
trunk/examples/extended/medical/GammaTherapy/cu_w.in
r807 r1337 13 13 /testem/phys/addPhysics emstandard 14 14 #/testem/phys/addPhysics emstandard_opt2 15 #/testem/phys/addPhysics lowenergy16 15 # 17 16 /testem/target1Material G4_W -
trunk/examples/extended/medical/GammaTherapy/include/EventActionMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.hh,v 1.2 2006/06/29 17:25:54 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/include/G4StepLimiterBuilder.hh
r1230 r1337 26 26 // 27 27 // $Id: G4StepLimiterBuilder.hh,v 1.2 2006/06/29 17:26:25 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/include/G4StepLimiterMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: G4StepLimiterMessenger.hh,v 1.2 2006/06/29 17:26:28 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/include/G4StepLimiterPerRegion.hh
r1230 r1337 25 25 // 26 26 // $Id: G4StepLimiterPerRegion.hh,v 1.4 2007/05/16 16:27:53 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/include/PhysicsList.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsList.hh,v 1.6 2006/11/16 14:42:20 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 // Modified: -
trunk/examples/extended/medical/GammaTherapy/include/PhysicsListMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.hh,v 1.3 2006/06/29 17:26:41 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/src/EventAction.cc
r807 r1337 41 41 42 42 #include "G4UImanager.hh" 43 #include "G4TrajectoryContainer.hh"44 #include "G4Trajectory.hh"45 #include "G4VVisManager.hh"46 43 #include "G4ios.hh" 47 44 … … 118 115 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo.... 119 116 120 void EventAction::EndOfEventAction(const G4Event* evt)117 void EventAction::EndOfEventAction(const G4Event*) 121 118 { 122 119 (Histo::GetPointer())->SaveEvent(); 123 G4VVisManager* pVVisManager = G4VVisManager::GetConcreteInstance();124 125 if(pVVisManager) {126 G4TrajectoryContainer* trjc = evt->GetTrajectoryContainer();127 G4int n_trajectories = 0;128 if (trjc) n_trajectories = trjc->entries();129 130 for(G4int i=0; i<n_trajectories; i++) {131 G4Trajectory* t = (G4Trajectory*)((*(evt->GetTrajectoryContainer()))[i]);132 if (drawFlag == "all") t->DrawTrajectory(1000);133 else if ((drawFlag == "charged")&&(t->GetCharge() != 0.))134 t->DrawTrajectory(1000);135 else if ((drawFlag == "neutral")&&(t->GetCharge() == 0.))136 t->DrawTrajectory(1000);137 else if ((drawFlag == "charged+n")&&((t->GetCharge() != 0.)||138 (t->GetCharge()==0.&&t->GetParticleName()=="neutron")))139 t->DrawTrajectory(1000);140 }141 }142 120 143 121 if(verbose > 0) { -
trunk/examples/extended/medical/GammaTherapy/src/EventActionMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.cc,v 1.2 2006/06/29 17:27:19 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/src/G4StepLimiterBuilder.cc
r1230 r1337 26 26 // 27 27 // $Id: G4StepLimiterBuilder.cc,v 1.3 2006/06/29 17:27:48 gunter Exp $ 28 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $28 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 29 29 // 30 30 //--------------------------------------------------------------------------- -
trunk/examples/extended/medical/GammaTherapy/src/G4StepLimiterMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: G4StepLimiterMessenger.cc,v 1.2 2006/06/29 17:27:50 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/src/G4StepLimiterPerRegion.cc
r1230 r1337 25 25 // 26 26 // $Id: G4StepLimiterPerRegion.cc,v 1.4 2007/05/16 16:27:53 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/src/PhysicsList.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.cc,v 1.1 6 2008/06/11 15:17:57 vnivanch Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: PhysicsList.cc,v 1.17 2009/11/21 16:47:07 vnivanch Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //--------------------------------------------------------------------------- … … 46 46 #include "PhysicsListMessenger.hh" 47 47 48 #include "ParticlesBuilder.hh"49 48 #include "G4EmStandardPhysics.hh" 50 49 #include "G4EmStandardPhysics_option1.hh" 51 50 #include "G4EmStandardPhysics_option2.hh" 52 51 #include "G4EmStandardPhysics_option3.hh" 53 #include " PhysListEmLivermore.hh"54 #include " PhysListEmPenelope.hh"52 #include "G4EmLivermorePhysics.hh" 53 #include "G4EmPenelopePhysics.hh" 55 54 #include "G4StepLimiterBuilder.hh" 56 55 #include "G4DecayPhysics.hh" … … 60 59 #include "G4EmExtraPhysics.hh" 61 60 #include "G4QStoppingPhysics.hh" 62 //#include "PhysListEmModelPai.hh"63 61 64 62 #include "G4UnitsTable.hh" … … 89 87 90 88 // Add Physics builders 91 RegisterPhysics(new ParticlesBuilder());89 RegisterPhysics(new G4DecayPhysics()); 92 90 } 93 91 … … 156 154 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 157 155 158 } else if (name == "livermore" && !emBuilderIsRegisted) { 159 RegisterPhysics(new PhysListEmLivermore()); 160 emBuilderIsRegisted = true; 161 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 162 163 } else if (name == "penelope" && !emBuilderIsRegisted) { 164 RegisterPhysics(new PhysListEmPenelope()); 165 emBuilderIsRegisted = true; 166 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 167 /* 168 } else if (name == "pai" && !emBuilderIsRegisted) { 169 RegisterPhysics(new PhysListEmModelPai()); 170 emBuilderIsRegisted = true; 171 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 172 */ 173 174 } else if (name == "decay" && !decayIsRegisted && emBuilderIsRegisted) { 175 RegisterPhysics(new G4DecayPhysics()); 176 decayIsRegisted = true; 156 } else if (name == "emlivermore" && !emBuilderIsRegisted) { 157 RegisterPhysics(new G4EmLivermorePhysics()); 158 emBuilderIsRegisted = true; 159 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 160 161 } else if (name == "empenelope" && !emBuilderIsRegisted) { 162 RegisterPhysics(new G4EmPenelopePhysics()); 163 emBuilderIsRegisted = true; 177 164 G4cout << "PhysicsList::AddPhysicsList <" << name << ">" << G4endl; 178 165 -
trunk/examples/extended/medical/GammaTherapy/src/PhysicsListMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.cc,v 1.3 2006/06/29 17:28:04 gunter Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... -
trunk/examples/extended/medical/GammaTherapy/vis.mac
r1230 r1337 27 27 /gun/energy 50.0 MeV 28 28 # 29 # Create empty scene ("world" is default)30 /vis/ scene/create29 # Use this open statement to create an OpenGL view: 30 /vis/open OGL 600x600-0+0 31 31 # 32 # Create a scene handler for a specific graphics system 33 # Edit the next line(s) to choose another graphic system 32 # Use this open statement to create a .prim file suitable for 33 # viewing in DAWN: 34 #/vis/open DAWNFILE 34 35 # 35 ####/vis/open DAWNFILE 36 /vis/open OGLIX 36 # Use this open statement to create a .heprep file suitable for 37 # viewing in HepRApp: 38 #/vis/open HepRepFile 37 39 # 38 # Draw scene 39 /vis/viewer/set/viewpointThetaPhi -90 0 deg 40 # Use this open statement to create a .wrl file suitable for 41 # viewing in a VRML viewer: 42 #/vis/open VRML2FILE 43 # 44 # Disable auto refresh and quieten vis messages whilst scene and 45 # trajectories are established: 46 /vis/viewer/set/autoRefresh false 47 /vis/verbose errors 48 # 49 # Draw geometry: 50 /vis/drawVolume 51 # 52 # Specify view angle: 53 /vis/viewer/set/viewpointThetaPhi -90. 0. 54 # 55 # Specify zoom value: 40 56 /vis/viewer/zoom 1.4 41 /vis/viewer/flush42 57 # 43 # for drawing the tracks 44 # if too many tracks cause core dump => storeTrajectory 0 45 /tracking/storeTrajectory 1 58 # Specify style (surface or wireframe): 59 #/vis/viewer/set/style wireframe 60 # 61 # Draw coordinate axes: 62 #/vis/scene/add/axes 0 0 0 1 m 63 # 64 # Draw smooth trajectories at end of event, showing trajectory points 65 # as markers 2 pixels wide: 66 /vis/scene/add/trajectories smooth 67 /vis/modeling/trajectories/create/drawByCharge 68 /vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true 69 /vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 2 70 # (if too many tracks cause core dump => /tracking/storeTrajectory 0) 71 # 72 # Draw hits at end of event: 73 #/vis/scene/add/hits 74 # 75 # To draw only gammas: 76 #/vis/filtering/trajectories/create/particleFilter 77 #/vis/filtering/trajectories/particleFilter-0/add gamma 78 # 79 # To invert the above, drawing all particles except gammas, 80 # keep the above two lines but also add: 81 #/vis/filtering/trajectories/particleFilter-0/invert true 82 # 83 # Many other options are available with /vis/modeling and /vis/filtering. 84 # For example, to select colour by particle ID: 85 #/vis/modeling/trajectories/create/drawByParticleID 86 #/vis/modeling/trajectories/drawByParticleID-0/set e- blue 87 # 88 # To superimpose all of the events from a given run: 46 89 /vis/scene/endOfEventAction accumulate 47 90 # 48 /testem/event/drawTracks all 91 # Re-establish auto refreshing and verbosity: 92 /vis/viewer/set/autoRefresh true 93 /vis/verbose warnings 94 # 95 # For file-based drivers, use this to create an empty detector view: 96 #/vis/viewer/flush
Note: See TracChangeset
for help on using the changeset viewer.