Changeset 1230 for trunk/examples/advanced/brachytherapy
- Timestamp:
- Jan 8, 2010, 3:02:48 PM (16 years ago)
- Location:
- trunk/examples/advanced/brachytherapy
- Files:
-
- 44 edited
-
Brachy.cc (modified) (1 diff)
-
GNUmakefile (modified) (2 diffs)
-
IridiumSourceMacro.mac (modified) (1 diff)
-
LeipzigSourceMacro.mac (modified) (1 diff)
-
include/BrachyAnalysisManager.hh (modified) (2 diffs)
-
include/BrachyDetectorConstruction.hh (modified) (1 diff)
-
include/BrachyDetectorConstructionI.hh (modified) (1 diff)
-
include/BrachyDetectorConstructionIr.hh (modified) (1 diff)
-
include/BrachyDetectorConstructionLeipzig.hh (modified) (1 diff)
-
include/BrachyDetectorMessenger.hh (modified) (1 diff)
-
include/BrachyEventAction.hh (modified) (1 diff)
-
include/BrachyFactory.hh (modified) (1 diff)
-
include/BrachyFactoryI.hh (modified) (1 diff)
-
include/BrachyFactoryIr.hh (modified) (1 diff)
-
include/BrachyFactoryLeipzig.hh (modified) (1 diff)
-
include/BrachyMaterial.hh (modified) (1 diff)
-
include/BrachyPhantomROGeometry.hh (modified) (1 diff)
-
include/BrachyPhantomSD.hh (modified) (1 diff)
-
include/BrachyPhysicsList.hh (modified) (1 diff)
-
include/BrachyPrimaryGeneratorAction.hh (modified) (1 diff)
-
include/BrachyPrimaryGeneratorActionI.hh (modified) (1 diff)
-
include/BrachyPrimaryGeneratorActionIr.hh (modified) (1 diff)
-
include/BrachyPrimaryGeneratorMessenger.hh (modified) (1 diff)
-
include/BrachyRunAction.hh (modified) (1 diff)
-
src/BrachyAnalysisManager.cc (modified) (14 diffs)
-
src/BrachyDetectorConstruction.cc (modified) (1 diff)
-
src/BrachyDetectorConstructionI.cc (modified) (1 diff)
-
src/BrachyDetectorConstructionIr.cc (modified) (1 diff)
-
src/BrachyDetectorConstructionLeipzig.cc (modified) (1 diff)
-
src/BrachyDetectorMessenger.cc (modified) (1 diff)
-
src/BrachyEventAction.cc (modified) (1 diff)
-
src/BrachyFactory.cc (modified) (1 diff)
-
src/BrachyFactoryI.cc (modified) (1 diff)
-
src/BrachyFactoryIr.cc (modified) (1 diff)
-
src/BrachyFactoryLeipzig.cc (modified) (1 diff)
-
src/BrachyMaterial.cc (modified) (1 diff)
-
src/BrachyPhantomROGeometry.cc (modified) (1 diff)
-
src/BrachyPhantomSD.cc (modified) (1 diff)
-
src/BrachyPhysicsList.cc (modified) (4 diffs)
-
src/BrachyPrimaryGeneratorAction.cc (modified) (1 diff)
-
src/BrachyPrimaryGeneratorActionI.cc (modified) (1 diff)
-
src/BrachyPrimaryGeneratorActionIr.cc (modified) (1 diff)
-
src/BrachyPrimaryGeneratorMessenger.cc (modified) (1 diff)
-
src/BrachyRunAction.cc (modified) (1 diff)
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/advanced/brachytherapy/Brachy.cc
r807 r1230 26 26 // 27 27 // $Id: Brachy.cc 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // -------------------------------------------------------------- -
trunk/examples/advanced/brachytherapy/GNUmakefile
r807 r1230 1 1 # -------------------------------------------------------------- 2 # $Id: GNUmakefile,v 1.1 2 2003/05/26 12:36:25 guatelliExp $2 # $Id: GNUmakefile,v 1.13 2008/06/15 18:05:59 cirrone Exp $ 3 3 # -------------------------------------------------------------- 4 4 # GNUmakefile for examples module. Gabriele Cosmo, 06/04/98. … … 19 19 include $(G4INSTALL)/config/binmake.gmk 20 20 21 ifdef G4ANALYSIS_USE22 CPPFLAGS += `aida-config --include`23 LDFLAGS += `aida-config --lib`24 LOADLIBS += `aida-config --lib`25 endif26 21 27 22 -
trunk/examples/advanced/brachytherapy/IridiumSourceMacro.mac
r807 r1230 3 3 /run/verbose 0 4 4 /event/verbose 0 5 /run/beamOn 200005 /run/beamOn 10000 6 6 7 7 -
trunk/examples/advanced/brachytherapy/LeipzigSourceMacro.mac
r807 r1230 3 3 /event/verbose 0 4 4 /source/switch Leipzig 5 /run/beamOn 200005 /run/beamOn 10000 -
trunk/examples/advanced/brachytherapy/include/BrachyAnalysisManager.hh
r807 r1230 24 24 // ******************************************************************** 25 25 // 26 // $Id: BrachyAnalysisManager.hh,v 1.1 3 2007/06/21 14:58:39 gunterExp $27 // GEANT4 tag $Name: $26 // $Id: BrachyAnalysisManager.hh,v 1.14 2008/06/05 13:45:39 cirrone Exp $ 27 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 28 28 29 29 // Code review: MG Pia, 14/05/2207 … … 63 63 64 64 void book(); 65 void FillNtupleWithEnergy(G4double,G4double,G4double,G4 float);65 void FillNtupleWithEnergy(G4double,G4double,G4double,G4double); 66 66 void FillHistogramWithEnergy(G4double,G4double,G4double); 67 67 void PrimaryParticleEnergySpectrum(G4double); -
trunk/examples/advanced/brachytherapy/include/BrachyDetectorConstruction.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyDetectorConstruction.hh,v 1.17 2006/06/29 15:47:11 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // **************************************** -
trunk/examples/advanced/brachytherapy/include/BrachyDetectorConstructionI.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyDetectorConstructionI.hh,v 1.5 2006/06/29 15:47:13 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // **************************************** -
trunk/examples/advanced/brachytherapy/include/BrachyDetectorConstructionIr.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyDetectorConstructionIr.hh,v 1.5 2006/06/29 15:47:16 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // **************************************** -
trunk/examples/advanced/brachytherapy/include/BrachyDetectorConstructionLeipzig.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyDetectorConstructionLeipzig.hh,v 1.5 2006/06/29 15:47:19 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // -
trunk/examples/advanced/brachytherapy/include/BrachyDetectorMessenger.hh
r807 r1230 32 32 // 33 33 // $Id: BrachyDetectorMessenger.hh,v 1.7 2006/06/29 15:47:22 gunter Exp $ 34 // GEANT4 tag $Name: $34 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 35 35 // 36 36 // -
trunk/examples/advanced/brachytherapy/include/BrachyEventAction.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyEventAction.hh,v 1.15 2006/06/29 15:47:28 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // Code developed by: -
trunk/examples/advanced/brachytherapy/include/BrachyFactory.hh
r807 r1230 27 27 // 28 28 // $Id: BrachyFactory.hh,v 1.6 2006/06/29 15:47:30 gunter Exp $ 29 // GEANT4 tag $Name: $29 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 30 30 // 31 31 // ********************************** -
trunk/examples/advanced/brachytherapy/include/BrachyFactoryI.hh
r807 r1230 36 36 // 37 37 // $Id: BrachyFactoryI.hh,v 1.7 2006/06/29 15:47:33 gunter Exp $ 38 // GEANT4 tag $Name: $38 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 39 39 // 40 40 // -
trunk/examples/advanced/brachytherapy/include/BrachyFactoryIr.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyFactoryIr.hh,v 1.6 2006/06/29 15:47:36 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ********************************** -
trunk/examples/advanced/brachytherapy/include/BrachyFactoryLeipzig.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyFactoryLeipzig.hh,v 1.6 2006/06/29 15:47:39 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ********************************** -
trunk/examples/advanced/brachytherapy/include/BrachyMaterial.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyMaterial.hh,v 1.7 2006/06/29 15:47:42 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ********************************** -
trunk/examples/advanced/brachytherapy/include/BrachyPhantomROGeometry.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPhantomROGeometry.hh,v 1.7 2006/06/29 15:47:45 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ************************************ -
trunk/examples/advanced/brachytherapy/include/BrachyPhantomSD.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPhantomSD.hh,v 1.10 2006/06/29 15:47:48 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ******************************** -
trunk/examples/advanced/brachytherapy/include/BrachyPhysicsList.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPhysicsList.hh,v 1.9 2006/11/15 10:02:17 guatelli Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // Code review: MGP, 5 November 2006 (still to be completed) -
trunk/examples/advanced/brachytherapy/include/BrachyPrimaryGeneratorAction.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPrimaryGeneratorAction.hh,v 1.12 2006/06/29 15:47:54 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ******************************************** -
trunk/examples/advanced/brachytherapy/include/BrachyPrimaryGeneratorActionI.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPrimaryGeneratorActionI.hh,v 1.6 2006/06/29 15:47:57 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // -
trunk/examples/advanced/brachytherapy/include/BrachyPrimaryGeneratorActionIr.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyPrimaryGeneratorActionIr.hh,v 1.7 2006/06/29 15:47:59 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // ************************************************* -
trunk/examples/advanced/brachytherapy/include/BrachyPrimaryGeneratorMessenger.hh
r807 r1230 32 32 // 33 33 // $Id: BrachyPrimaryGeneratorMessenger.hh,v 1.2 2006/06/29 15:48:02 gunter Exp $ 34 // GEANT4 tag $Name: $34 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 35 35 // 36 36 // -
trunk/examples/advanced/brachytherapy/include/BrachyRunAction.hh
r807 r1230 26 26 // 27 27 // $Id: BrachyRunAction.hh,v 1.11 2006/06/29 15:48:04 gunter Exp $ 28 // GEANT4 tag $Name: $28 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 29 29 // 30 30 // -
trunk/examples/advanced/brachytherapy/src/BrachyAnalysisManager.cc
r807 r1230 33 33 // ******************************* 34 34 // 35 // $Id: BrachyAnalysisManager.cc,v 1. 19 2007/05/14 09:57:55 pia Exp $36 // GEANT4 tag $Name: $35 // $Id: BrachyAnalysisManager.cc,v 1.22 2009/11/12 10:32:59 pandola Exp $ 36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 #ifdef G4ANALYSIS_USE 39 38 40 #include <stdlib.h> 39 41 #include <fstream> … … 42 44 #include "G4ios.hh" 43 45 44 #ifdef G4ANALYSIS_USE 46 45 47 #include "AIDA/IHistogram1D.h" 46 48 #include "AIDA/IHistogram2D.h" … … 52 54 #include "AIDA/ITree.h" 53 55 #include "AIDA/ITuple.h" 54 #endif 56 55 57 56 58 BrachyAnalysisManager* BrachyAnalysisManager::instance = 0; 57 59 58 BrachyAnalysisManager::BrachyAnalysisManager() 60 BrachyAnalysisManager::BrachyAnalysisManager() : 61 aFact(0), theTree(0), histFact(0), tupFact(0), h1(0), 62 h2(0), h3(0), ntuple(0) 63 59 64 { 60 61 #ifdef G4ANALYSIS_USE 62 aFact = 0; 63 theTree(0); 64 histFact = 0; 65 tupFact = 0; 66 h1 = 0; 67 h2 = 0; 68 ntuple = 0; 69 65 70 66 // Instantiate the factories 71 67 // The factories manage the analysis objects … … 75 71 76 72 // Definition of the output file 77 G4String fileName = "brachytherapy. hbk";78 theTree = treeFact -> create(fileName," hbook",false, true);73 G4String fileName = "brachytherapy.root"; 74 theTree = treeFact -> create(fileName,"ROOT",false, true); 79 75 80 76 delete treeFact; 81 #endif82 77 } 83 78 84 79 BrachyAnalysisManager::~BrachyAnalysisManager() 85 80 { 86 87 #ifdef G4ANALYSIS_USE88 81 delete tupFact; 89 82 tupFact = 0; … … 97 90 delete aFact; 98 91 aFact = 0; 99 #endif100 92 } 101 93 … … 108 100 void BrachyAnalysisManager::book() 109 101 { 110 #ifdef G4ANALYSIS_USE111 102 // Instantiate the histogram and ntuple factories 112 103 histFact = aFact -> createHistogramFactory( *theTree ); … … 129 120 130 121 //defining the ntuple columns' name 131 std::string columnNames = " float energy; float x; float y; float z";122 std::string columnNames = "double energy, x , y , z "; 132 123 std::string options = ""; 133 124 … … 136 127 // check for non-zero ... 137 128 if (ntuple) G4cout<<"The Ntuple is non-zero"<<G4endl; 138 #endif139 129 } 140 130 … … 142 132 G4double yy, 143 133 G4double zz, 144 G4 floaten)134 G4double en) 145 135 { 146 #ifdef G4ANALYSIS_USE147 136 if (ntuple == 0) 148 137 { … … 164 153 165 154 ntuple->addRow(); 166 #endif167 155 } 168 156 … … 171 159 G4double energyDeposit) 172 160 { 173 #ifdef G4ANALYSIS_USE174 161 // 2DHistogram: energy deposit in a voxel which center is fixed in position (x,z) 175 162 h1 -> fill(x,z,energyDeposit); 176 #endif177 163 } 178 164 179 165 void BrachyAnalysisManager::PrimaryParticleEnergySpectrum(G4double primaryParticleEnergy) 180 166 { 181 #ifdef G4ANALYSIS_USE182 167 // 1DHistogram: energy spectrum of primary particles 183 168 h2 -> fill(primaryParticleEnergy); 184 #endif185 169 return; 186 170 } … … 188 172 void BrachyAnalysisManager::DoseDistribution(G4double x,G4double energy) 189 173 { 190 #ifdef G4ANALYSIS_USE191 174 // 1DHistogram: energy spectrum of primary particles 192 175 h3 -> fill(x, energy); 193 #endif194 176 } 195 177 196 178 void BrachyAnalysisManager::finish() 197 179 { 198 #ifdef G4ANALYSIS_USE 180 199 181 // write all histograms to file ... 200 182 theTree -> commit(); … … 202 184 // close (will again commit) ... 203 185 theTree -> close(); 186 } 204 187 #endif 205 }206 188 207 189 … … 215 197 216 198 217 -
trunk/examples/advanced/brachytherapy/src/BrachyDetectorConstruction.cc
r807 r1230 39 39 // 40 40 // $Id: BrachyDetectorConstruction.cc,v 1.29 2006/06/29 15:48:09 gunter Exp $ 41 // GEANT4 tag $Name: $41 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 42 42 // 43 43 #include "G4CSGSolid.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyDetectorConstructionI.cc
r807 r1230 38 38 // 39 39 // $Id: BrachyDetectorConstructionI.cc,v 1.9 2006/06/29 15:48:11 gunter Exp $ 40 // GEANT4 tag $Name: $40 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 41 41 // 42 42 #include "globals.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyDetectorConstructionIr.cc
r807 r1230 39 39 // 40 40 // $Id: BrachyDetectorConstructionIr.cc,v 1.10 2006/06/29 15:48:13 gunter Exp $ 41 // GEANT4 tag $Name: $41 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 42 42 // 43 43 #include "globals.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyDetectorConstructionLeipzig.cc
r807 r1230 40 40 // 41 41 // $Id: BrachyDetectorConstructionLeipzig.cc,v 1.10 2006/06/29 15:48:16 gunter Exp $ 42 // GEANT4 tag $Name: $42 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 43 43 // 44 44 -
trunk/examples/advanced/brachytherapy/src/BrachyDetectorMessenger.cc
r807 r1230 36 36 // 37 37 // $Id: BrachyDetectorMessenger.cc,v 1.12 2006/06/29 15:48:18 gunter Exp $ 38 // GEANT4 tag $Name: $38 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 39 39 // 40 40 // -
trunk/examples/advanced/brachytherapy/src/BrachyEventAction.cc
r807 r1230 34 34 // 35 35 // $Id: BrachyEventAction.cc,v 1.21 2006/06/29 15:48:21 gunter Exp $ 36 // GEANT4 tag $Name: $36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 38 #include "BrachyEventAction.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyFactory.cc
r807 r1230 36 36 // 37 37 // $Id: BrachyFactory.cc,v 1.7 2006/06/29 15:48:23 gunter Exp $ 38 // GEANT4 tag $Name: $38 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 39 39 // 40 40 // Factory of brachytherapic sources -
trunk/examples/advanced/brachytherapy/src/BrachyFactoryI.cc
r807 r1230 34 34 // 35 35 // $Id: BrachyFactoryI.cc,v 1.6 2006/06/29 15:48:26 gunter Exp $ 36 // GEANT4 tag $Name: $36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 38 #include "BrachyFactoryI.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyFactoryIr.cc
r807 r1230 34 34 // 35 35 // $Id: BrachyFactoryIr.cc,v 1.6 2006/06/29 15:48:29 gunter Exp $ 36 // GEANT4 tag $Name: $36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 38 #include "globals.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyFactoryLeipzig.cc
r807 r1230 34 34 // 35 35 // $Id: BrachyFactoryLeipzig.cc,v 1.6 2006/06/29 15:48:31 gunter Exp $ 36 // GEANT4 tag $Name: $36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 38 -
trunk/examples/advanced/brachytherapy/src/BrachyMaterial.cc
r807 r1230 34 34 // 35 35 // $Id: BrachyMaterial.cc,v 1.7 2006/06/29 15:48:33 gunter Exp $ 36 // GEANT4 tag $Name: $36 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 37 37 // 38 38 -
trunk/examples/advanced/brachytherapy/src/BrachyPhantomROGeometry.cc
r807 r1230 35 35 // 36 36 // $Id: BrachyPhantomROGeometry.cc,v 1.11 2006/06/29 15:48:36 gunter Exp $ 37 // GEANT4 tag $Name: $37 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 38 38 // 39 39 #include "BrachyPhantomROGeometry.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyPhantomSD.cc
r807 r1230 34 34 // ******************************** 35 35 // 36 // $Id: BrachyPhantomSD.cc,v 1.1 3 2006/06/29 15:48:39gunter Exp $37 // GEANT4 tag $Name: $36 // $Id: BrachyPhantomSD.cc,v 1.14 2009/02/23 17:34:26 gunter Exp $ 37 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 38 38 // 39 39 #include "BrachyPhantomSD.hh" 40 #ifdef G4ANALYSIS_USE 40 41 #include "BrachyAnalysisManager.hh" 42 #endif 41 43 #include "BrachyDetectorConstruction.hh" 42 44 #include "G4Track.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyPhysicsList.cc
r807 r1230 36 36 // ********************************** 37 37 // 38 // $Id: BrachyPhysicsList.cc,v 1.1 3 2006/11/15 10:02:17 guatelliExp $39 // GEANT4 tag $Name: $38 // $Id: BrachyPhysicsList.cc,v 1.14 2009/11/12 02:50:51 cirrone Exp $ 39 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 40 40 // 41 41 #include "BrachyPhysicsList.hh" … … 47 47 #include "G4UnitsTable.hh" 48 48 #include "G4ios.hh" 49 50 49 #include "G4MultipleScattering.hh" 51 50 // gamma 52 #include "G4LowEnergyRayleigh.hh" 53 #include "G4LowEnergyPhotoElectric.hh" 54 #include "G4LowEnergyCompton.hh" 55 #include "G4LowEnergyGammaConversion.hh" 51 #include "G4PhotoElectricEffect.hh" 52 #include "G4LivermorePhotoElectricModel.hh" 53 54 #include "G4ComptonScattering.hh" 55 #include "G4LivermoreComptonModel.hh" 56 57 #include "G4GammaConversion.hh" 58 #include "G4LivermoreGammaConversionModel.hh" 59 60 #include "G4RayleighScattering.hh" 61 #include "G4LivermoreRayleighModel.hh" 62 56 63 // e- 57 #include "G4LowEnergyIonisation.hh" 58 #include "G4LowEnergyBremsstrahlung.hh" 64 #include "G4eMultipleScattering.hh" 65 #include "G4eIonisation.hh" 66 #include "G4LivermoreIonisationModel.hh" 67 #include "G4eBremsstrahlung.hh" 68 #include "G4LivermoreBremsstrahlungModel.hh" 69 59 70 // e+ 60 71 #include "G4eIonisation.hh" … … 80 91 ConstructBosons(); 81 92 ConstructLeptons(); 93 82 94 } 83 95 … … 114 126 115 127 if (particleName == "gamma") { 128 116 129 // Photon 117 pmanager->AddDiscreteProcess(new G4LowEnergyRayleigh); 118 pmanager->AddDiscreteProcess(new G4LowEnergyPhotoElectric); 119 pmanager->AddDiscreteProcess(new G4LowEnergyCompton); 120 pmanager->AddDiscreteProcess(new G4LowEnergyGammaConversion); 130 G4RayleighScattering* theRayleigh = new G4RayleighScattering(); 131 theRayleigh->SetModel(new G4LivermoreRayleighModel()); //not strictly necessary 132 pmanager->AddDiscreteProcess(theRayleigh); 133 134 G4PhotoElectricEffect* thePhotoElectricEffect = new G4PhotoElectricEffect(); 135 thePhotoElectricEffect->SetModel(new G4LivermorePhotoElectricModel()); 136 pmanager->AddDiscreteProcess(thePhotoElectricEffect); 137 138 G4ComptonScattering* theComptonScattering = new G4ComptonScattering(); 139 theComptonScattering->SetModel(new G4LivermoreComptonModel()); 140 pmanager->AddDiscreteProcess(theComptonScattering); 141 142 G4GammaConversion* theGammaConversion = new G4GammaConversion(); 143 theGammaConversion->SetModel(new G4LivermoreGammaConversionModel()); 144 pmanager->AddDiscreteProcess(theGammaConversion); 121 145 122 146 } else if (particleName == "e-") { 123 147 // Electron 124 G4LowEnergyIonisation* loweIon = new G4LowEnergyIonisation("LowEnergyIoni"); 125 126 G4LowEnergyBremsstrahlung* loweBrem = new G4LowEnergyBremsstrahlung("LowEnBrem"); 127 // Select the Bremsstrahlung angular distribution model (Tsai/2BN/2BS) 128 loweBrem->SetAngularGenerator("tsai"); 129 130 pmanager->AddProcess(new G4MultipleScattering, -1, 1,1); 131 pmanager->AddProcess(loweIon, -1, 2,2); 132 pmanager->AddProcess(loweBrem, -1,-1,3); 148 149 G4eMultipleScattering* msc = new G4eMultipleScattering(); 150 msc->SetStepLimitType(fUseDistanceToBoundary); 151 pmanager->AddProcess(msc,-1, 1, 1); 152 153 // Ionisation 154 G4eIonisation* eIonisation = new G4eIonisation(); 155 eIonisation->SetEmModel(new G4LivermoreIonisationModel()); 156 eIonisation->SetStepFunction(0.2, 100*um); //improved precision in tracking 157 pmanager->AddProcess(eIonisation,-1, 2, 2); 158 159 // Bremsstrahlung 160 G4eBremsstrahlung* eBremsstrahlung = new G4eBremsstrahlung(); 161 eBremsstrahlung->SetEmModel(new G4LivermoreBremsstrahlungModel()); 162 pmanager->AddProcess(eBremsstrahlung, -1,-3, 3); 133 163 134 164 } else if (particleName == "e+") { 135 165 // Positron 136 pmanager->AddProcess(new G4MultipleScattering, -1, 1,1); 137 pmanager->AddProcess(new G4eIonisation, -1, 2,2); 138 pmanager->AddProcess(new G4eBremsstrahlung, -1,-1,3); 139 pmanager->AddProcess(new G4eplusAnnihilation, 0,-1,4); 166 G4eMultipleScattering* msc = new G4eMultipleScattering(); 167 msc->SetStepLimitType(fUseDistanceToBoundary); 168 pmanager->AddProcess(msc,-1, 1, 1); 169 170 // Ionisation 171 G4eIonisation* eIonisation = new G4eIonisation(); 172 eIonisation->SetStepFunction(0.2, 100*um); // 173 pmanager->AddProcess(eIonisation, -1, 2, 2); 174 175 //Bremsstrahlung (use default, no low-energy available) 176 pmanager->AddProcess(new G4eBremsstrahlung(), -1,-1, 3); 177 178 //Annihilation 179 pmanager->AddProcess(new G4eplusAnnihilation(),0,-1, 4); 140 180 141 181 } -
trunk/examples/advanced/brachytherapy/src/BrachyPrimaryGeneratorAction.cc
r807 r1230 38 38 // ******************************************** 39 39 // 40 // $Id: BrachyPrimaryGeneratorAction.cc,v 1.1 8 2006/06/29 15:48:45gunter Exp $41 // GEANT4 tag $Name: $40 // $Id: BrachyPrimaryGeneratorAction.cc,v 1.19 2009/02/23 17:34:26 gunter Exp $ 41 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 42 42 // 43 43 #include "globals.hh" 44 44 #include "BrachyPrimaryGeneratorAction.hh" 45 #ifdef G4ANALYSIS_USE 45 46 #include "BrachyAnalysisManager.hh" 47 #endif 46 48 #include "G4ParticleTable.hh" 47 49 #include "Randomize.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyPrimaryGeneratorActionI.cc
r807 r1230 39 39 // 40 40 // $Id: BrachyPrimaryGeneratorActionI.cc,v 1.10 2006/06/29 15:48:48 gunter Exp $ 41 // GEANT4 tag $Name: $41 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 42 42 // 43 43 #include "BrachyPrimaryGeneratorActionI.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyPrimaryGeneratorActionIr.cc
r807 r1230 39 39 // 40 40 // $Id: BrachyPrimaryGeneratorActionIr.cc,v 1.11 2006/06/29 15:48:51 gunter Exp $ 41 // GEANT4 tag $Name: $41 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 42 42 // 43 43 #include "BrachyPrimaryGeneratorActionIr.hh" -
trunk/examples/advanced/brachytherapy/src/BrachyPrimaryGeneratorMessenger.cc
r807 r1230 36 36 // 37 37 // $Id: BrachyPrimaryGeneratorMessenger.cc,v 1.3 2006/06/29 15:48:54 gunter Exp $ 38 // GEANT4 tag $Name: $38 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 39 39 // 40 40 // -
trunk/examples/advanced/brachytherapy/src/BrachyRunAction.cc
r807 r1230 40 40 // 41 41 // $Id: BrachyRunAction.cc,v 1.18 2006/06/29 15:48:57 gunter Exp $ 42 // GEANT4 tag $Name: $42 // GEANT4 tag $Name: geant4-09-03-cand-01 $ 43 43 // 44 44
Note:
See TracChangeset
for help on using the changeset viewer.
