Changeset 1337 for trunk/examples/extended/medical/fanoCavity2
- Timestamp:
- Sep 30, 2010, 2:47:17 PM (14 years ago)
- Location:
- trunk/examples/extended/medical/fanoCavity2
- Files:
-
- 41 edited
Legend:
- Unmodified
- Added
- Removed
-
trunk/examples/extended/medical/fanoCavity2/History
r807 r1337 1 $Id: History,v 1. 4 2007/11/12 18:19:30 maireExp $1 $Id: History,v 1.18 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 (fano2-V09-03-03) 19 - Remove unused variable in EventAction 20 21 04-06-10 J.Perl (fano2-V09-03-02) 22 - Updated vis usage 23 24 20-01-10 mma (fano2-V09-03-00) 25 - introduce PhysListEmStandard_option2 : UrbanMscModel93 without options 26 PhysListEmStandard_WVI : WentzelVIModel 27 - put in run01.mac all physics lists choice 28 29 16-11-09 mma (fano2-V09-02-01) 30 - PhysListEmStandard_option3 : use UrbanMscModel93 31 32 31-10-09 mma (fano2-V09-02-00) 33 - splite PhysicsList in 2, as others TestEm examples. 34 - local copy of PhysListEmStandard_option0, PhysListEmStandard_option3, 35 PhysListEmStandard_SS 36 - new PhysListEmStandard_GS for GoudsmitSaundersonMsc model 37 - PhysicsList: add cut for proton recoil in case of Single Scattering 38 39 16-11-08 mma (fano2-V09-01-06) 40 - Em options: 20bins/decade LinLossLimit=default=0.01 41 42 01-11-08 mma (fano2-V09-01-05) 43 - PhysicsList : add comments for Em options 44 45 26-10-08 mma (fano2-V09-01-04) 46 - PhysicsList : modified Em options 47 48 28-09-08 mma (fano2-V09-01-03) 49 - HistoManager : root format by default 50 51 11-06-08 V.Ivant (fano2-V09-01-02) 52 - Remove AIDA from GNUmakefile 53 54 07-05-08 mma (fano2-V09-01-01) 55 - README : update Aida informations 56 57 11-04-08 mma 58 - HistoManager::Scale() : suppress warning 59 60 04-04-08 mma(fano2-V09-01-00) 61 - PhysicsList : replace G4MultipleScattering by new G4eMultipleScattering 62 --> allows to use MscUrbanModel2 17 63 18 64 12-11-07 mma (fano2-V09-00-01) -
trunk/examples/extended/medical/fanoCavity2/README
r807 r1337 1 $Id: README,v 1. 8 2007/11/12 18:19:30maire Exp $1 $Id: README,v 1.11 2009/10/31 18:05:56 maire Exp $ 2 2 ------------------------------------------------------------------- 3 3 … … 82 82 - Bremsstrahlung : Fano conditions imply no energy transfer via 83 83 bremsstrahlung radiation. Therefore this process is not registered in the 84 physics list. However, it is always possible to include it via an UI85 command. See PhysicsListMessenger class.84 physics list. However, it is always possible to include it. 85 See PhysListEm classes. 86 86 87 87 - Ionization : In order to have same stopping power in wall and cavity, one … … 125 125 126 126 One can control the name of the histograms file with the command: 127 /testem/histo/setFileName name (default fano Cavity)127 /testem/histo/setFileName name (default fanocavity2) 128 128 129 129 It is possible to choose the format of the histogram file (hbook, root, XML) 130 with the command /testem/histo/setFileType ( hbookby default)130 with the command /testem/histo/setFileType (root by default) 131 131 132 132 It is also possible to print selected histograms on an ascii file: … … 156 156 157 157 To use histograms, at least one of the AIDA implementations should be 158 available (see http://aida.freehep.org). 159 160 8a - PI 161 162 A package including AIDA and extended interfaces also using Python is PI, 163 available from: http://cern.ch/pi 164 165 Once installed PI or PI-Lite in a specified local area $MYPY, it is required 166 to add the installation path to $PATH, i.e. for example, for release 1.2.1 of 167 PI: 168 setenv PATH ${PATH}:$MYPI/1.2.1/app/releases/PI/PI_1_2_1/rh73_gcc32/bin 169 170 CERN users can use the PATH to the LCG area on AFS. 171 Before running the example the command should be issued: 172 eval `aida-config --runtime csh` 173 174 8b - OpenScientist 175 176 OpenScientist is available at http://OpenScientist.lal.in2p3.fr. 177 178 You have to "setup" the OpenScientist AIDA implementation before compiling 179 (then with G4ANALYSIS_USE set) and running your Geant4 application. 180 181 On UNIX you setup, with a csh flavoured shell : 182 csh> source <<OpenScientist install path>/aida-setup.csh 183 or with a sh flavoured shell : 184 sh> . <<OpenScientist install path>/aida-setup.sh 185 On Windows : 186 DOS> call <<OpenScientist install path>/aida-setup.bat 187 188 You can use various file formats for writing (AIDA-XML, hbook, root). 189 These formats are readable by the Lab onx interactive program 190 or the OpenPAW application. See the web pages. 191 192 193 With OpenPAW, on a run.hbook file, one can view the histograms 194 with something like : 195 OS> opaw 196 opaw> h/file 1 run.hbook ( or opaw> h/file 1 run.aida or run.root) 197 opaw> zone 2 2 198 opaw> h/plot 1 199 opaw> h/plot 2 158 available. See InstallAida.txt -
trunk/examples/extended/medical/fanoCavity2/essai.mac
r1230 r1337 1 # $Id: essai.mac,v 1. 5 2009/10/31 18:05:56maire Exp $1 # $Id: essai.mac,v 1.6 2010/01/20 15:41:14 maire Exp $ 2 2 # 3 3 /control/verbose 2 4 4 /run/verbose 2 5 5 # 6 /testem/phys/addPhysics standard_ opt36 /testem/phys/addPhysics standard_WVI 7 7 # 8 8 /run/initialize … … 10 10 /gun/energy 1 MeV 11 11 # 12 /testem/event/printModulo 10000 012 /testem/event/printModulo 10000 13 13 # 14 /run/beamOn 100000 014 /run/beamOn 100000 -
trunk/examples/extended/medical/fanoCavity2/fanoCavity2.cc
r1230 r1337 25 25 // 26 26 // $Id: fanoCavity2.cc,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/fanoCavity2.in
r807 r1337 1 # $Id: fanoCavity2.in,v 1. 2 2007/11/05 13:19:16 maire Exp $1 # $Id: fanoCavity2.in,v 1.3 2009/10/31 18:05:56 maire Exp $ 2 2 # 3 3 /control/verbose 2 4 4 /run/verbose 2 5 # 6 /testem/phys/addPhysics standard_opt3 5 7 # 6 8 /run/initialize -
trunk/examples/extended/medical/fanoCavity2/fanoCavity2.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: Water density: 1.000 g/cm3 RadL: 36.092 cm Imean: 75.000 eV16 Material: Water density: 1.000 g/cm3 RadL: 36.092 cm Nucl.Int.Length: 75.416 cm Imean: 75.000 eV 13 17 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 14 18 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 15 19 16 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.92 5m Imean: 75.000 eV temperature: 273.15 K pressure: 1.00 atm20 Material: Water_gas density: 1.000 mg/cm3 RadL: 360.924 m Nucl.Int.Length: 754.163 m Imean: 75.000 eV temperature: 273.15 K pressure: 1.00 atm 17 21 ---> Element: Hydrogen (H) Z = 1.0 N = 1.0 A = 1.01 g/mole ElmMassFraction: 11.21 % ElmAbundance 66.67 % 18 22 ---> Element: Oxygen (O) Z = 8.0 N = 16.0 A = 16.00 g/mole ElmMassFraction: 88.79 % ElmAbundance 33.33 % 19 23 20 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Imean: 78.000 eV24 Material: Graphite density: 2.265 g/cm3 RadL: 18.850 cm Nucl.Int.Length: 35.387 cm Imean: 81.000 eV 21 25 ---> Element: Graphite ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 22 26 23 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.496 m Imean: 78.000 eV temperature: 273.15 K pressure: 1.00 atm27 Material: Graphite_gas density: 2.265 mg/cm3 RadL: 188.496 m Nucl.Int.Length: 353.873 m Imean: 81.000 eV temperature: 273.15 K pressure: 1.00 atm 24 28 ---> Element: Graphite_gas ( ) Z = 6.0 N = 12.0 A = 12.01 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 25 29 26 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Imean: 166.400 eV30 Material: Aluminium density: 2.700 g/cm3 RadL: 8.893 cm Nucl.Int.Length: 38.879 cm Imean: 166.000 eV 27 31 ---> Element: Aluminium ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 28 32 29 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.925 m Imean: 166.400 eV temperature: 273.15 K pressure: 1.00 atm33 Material: Aluminium_gas density: 2.700 mg/cm3 RadL: 88.925 m Nucl.Int.Length: 388.793 m Imean: 166.000 eV temperature: 273.15 K pressure: 1.00 atm 30 34 ---> Element: Aluminium_gas ( ) Z = 13.0 N = 27.0 A = 26.98 g/mole ElmMassFraction: 100.00 % ElmAbundance 100.00 % 31 35 32 36 33 37 /run/verbose 2 38 # 39 /testem/phys/addPhysics standard_opt3 40 PhysicsList::AddPhysicsList: <standard_opt3> 34 41 # 35 42 /run/initialize … … 44 51 physicsList->Construct() start. 45 52 physicsList->Construct() start. 53 physicsList->CheckParticleList() start. 46 54 physicsList->setCut() start. 47 55 PhysicsList::SetCuts:CutLength : 10 km … … 52 60 # 53 61 /run/beamOn 40000 54 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for gamma 55 The cut in range [1e+07 (mm)] is too big for material idx=0 56 The cut in energy is set1.79769e+305GeV 57 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e- 58 The cut in range [1e+07 (mm)] is too big for material idx=0 59 The cut in energy is set1.79769e+305GeV 60 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for e+ 61 The cut in range [1e+07 (mm)] is too big for material idx=0 62 The cut in energy is set1.79769e+305GeV 63 G4VRangeToEnergyConverter::ConvertCutToKineticEnergy for gamma 64 The cut in range [1e+07 (mm)] is too big for material idx=1 65 The cut in energy is set1.79769e+305GeV 66 67 phot: Total cross sections from Sandia parametrisation. 68 Sampling according PhotoElectric model 69 70 compt: Total cross sections has a good parametrisation from 10 KeV to (100/Z) GeV 71 Sampling according Klein-Nishina model 72 tables are built for gamma 73 Lambda tables from 100 eV to 10 GeV in 800 bins. 74 75 conv: Total cross sections has a good parametrisation from 1.5 MeV to 100 GeV for all Z; 76 sampling secondary e+e- according Bethe-Heitler model 77 tables are built for gamma 78 Lambda tables from 1.022 MeV to 10 GeV in 800 bins. 79 80 msc: Model variant of multiple scattering for e- 81 Lambda tables from 100 eV to 10 GeV in 800 bins. 82 LateralDisplacementFlag= 1 Skin= 2 83 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2 84 85 eIoni: tables are built for e- 86 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 87 Lambda tables from threshold to 10 GeV in 800 bins. 88 Delta cross sections and sampling from myMollerBhabha model 89 Good description from 1 KeV to 100 GeV. 90 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 91 CSDA range table up to 10 GeV in 800 bins. 92 93 eIoni: tables are built for e+ 94 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 95 Lambda tables from threshold to 10 GeV in 800 bins. 96 Delta cross sections and sampling from myMollerBhabha model 97 Good description from 1 KeV to 100 GeV. 98 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 99 CSDA range table up to 10 GeV in 800 bins. 100 101 msc: Model variant of multiple scattering for proton 102 Lambda tables from 100 eV to 10 GeV in 800 bins. 103 LateralDisplacementFlag= 1 Skin= 2 104 Boundary/stepping algorithm is active with RangeFactor= 0.02 Step limit type 2 105 106 hIoni: tables are built for proton 107 dE/dx and range tables from 100 eV to 10 GeV in 800 bins. 108 Lambda tables from threshold to 10 GeV in 800 bins. 109 Scaling relation is used from proton dE/dx and range. 110 Delta cross sections and sampling from BetheBloch model for scaled energy > 2 MeV 111 Parametrisation from Bragg for protons below. NuclearStopping= 1 112 Step function: finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1 113 CSDA range table up to 10 GeV in 800 bins. 62 63 phot: for gamma SubType= 12 64 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 65 PhotoElectric : Emin= 0 eV Emax= 10 GeV 66 67 compt: for gamma SubType= 13 68 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 69 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 70 Klein-Nishina : Emin= 0 eV Emax= 10 GeV 71 72 conv: for gamma SubType= 14 73 Lambda tables from 1.022 MeV to 10 GeV in 160 bins, spline: 1 74 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 75 Bethe-Heitler : Emin= 0 eV Emax= 10 GeV 76 77 msc: for e- SubType= 10 78 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 79 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 80 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 81 UrbanMsc93 : Emin= 0 eV Emax= 10 GeV 82 83 eIoni: for e- SubType= 2 84 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 85 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 86 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 87 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 88 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 89 CSDA range table up to 10 GeV in 160 bins 90 91 eIoni: for e+ SubType= 2 92 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 93 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 94 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 95 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 96 myMollerBhabha : Emin= 0 eV Emax= 10 GeV 97 CSDA range table up to 10 GeV in 160 bins 98 99 annihil: for e+ SubType= 5 100 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 101 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 102 eplus2gg : Emin= 0 eV Emax= 10 GeV 103 104 msc: for proton SubType= 10 105 Lambda tables from 100 eV to 10 GeV in 160 bins, spline: 1 106 RangeFactor= 0.02, stepLimitType: 2, latDisplacement: 1, skin= 3, geomFactor= 2.5 107 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 108 UrbanMsc90 : Emin= 0 eV Emax= 10 GeV 109 110 hIoni: for proton SubType= 2 111 dE/dx and range tables from 100 eV to 10 GeV in 160 bins 112 Lambda tables from threshold to 10 GeV in 160 bins, spline: 1 113 finalRange(mm)= 0.01, dRoverRange= 0.2, integral: 1, fluct: 1, linLossLimit= 0.01 114 ===== EM models for the G4Region DefaultRegionForTheWorld ====== 115 Bragg : Emin= 0 eV Emax= 2 MeV 116 BetheBloch : Emin= 2 MeV Emax= 10 GeV 117 CSDA range table up to 10 GeV in 160 bins 114 118 115 119 Region <DefaultRegionForTheWorld> -- appears in <Wall> world volume … … 117 121 Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0] 118 122 Materials : Water Water_gas 119 Production cuts : gamma 10 km e- 10 km e+10 km123 Production cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 120 124 121 125 ========= Table of registered couples ============================== … … 123 127 Index : 0 used in the geometry : Yes recalculation needed : No 124 128 Material : Water 125 Range cuts : gamma 10 km e- 10 km e+10 km126 Energy thresholds : gamma 2.88022e+295 J e- 2.88022e+295 J e+ 2.88022e+295 J129 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 130 Energy thresholds : gamma 10 GeV e- 10 GeV e+ 10 GeV proton 1 TeV 127 131 Region(s) which use this couple : 128 132 DefaultRegionForTheWorld … … 130 134 Index : 1 used in the geometry : Yes recalculation needed : No 131 135 Material : Water_gas 132 Range cuts : gamma 10 km e- 10 km e+10 km133 Energy thresholds : gamma 2.88022e+295 J e- 18.4792 GeV e+ 17.5913 GeV136 Range cuts : gamma 10 km e- 10 km e+ 10 km proton 10 km 137 Energy thresholds : gamma 10 GeV e- 10 GeV e+ 10 GeV proton 1 TeV 134 138 Region(s) which use this couple : 135 139 DefaultRegionForTheWorld … … 163 167 Start Run processing. 164 168 165 --->evntNb= 10000 Nwall= 9997 Ncav= 4 Ic/Iw= 2.09883 Ne-_cav= 2823 doseCavity/Ebeam= 1.13192 (100*(ratio-1) = 13.19193%)166 167 --->evntNb= 20000 Nwall= 1999 6 Ncav= 5 Ic/Iw= 1.31164 Ne-_cav= 5677 doseCavity/Ebeam= 1.23974 (100*(ratio-1) = 23.97404%)168 169 --->evntNb= 30000 Nwall= 2999 1 Ncav= 10 Ic/Iw= 1.74903 Ne-_cav= 8523 doseCavity/Ebeam= 1.29183 (100*(ratio-1) = 29.18310 %)169 --->evntNb= 10000 Nwall= 10000 Ncav= 1 Ic/Iw= 0.52456 Ne-_cav= 2787 doseCavity/Ebeam= 1.12338 (100*(ratio-1) = 12.33820 %) 170 171 --->evntNb= 20000 Nwall= 19999 Ncav= 2 Ic/Iw= 0.52458 Ne-_cav= 5507 doseCavity/Ebeam= 1.00206 (100*(ratio-1) = 0.20619 %) 172 173 --->evntNb= 30000 Nwall= 29998 Ncav= 3 Ic/Iw= 0.52459 Ne-_cav= 8236 doseCavity/Ebeam= 0.98084 (100*(ratio-1) = -1.91610 %) 170 174 Run terminated. 171 175 Run Summary 172 176 Number of events processed : 40000 173 User= 13s Real=13s Sys=0.01s174 175 Process calls frequency ---> msc= 351630 eIoni= 635349 Transportation= 22715177 User=9.6s Real=12s Sys=0s 178 179 Process calls frequency ---> msc= 507867 eIoni= 620396 Transportation= 22034 176 180 177 181 Charged particle flow in cavity : 178 Enter --> nbParticles = 11 358 Energy = 7.238GeV179 Exit --> nbParticles = 11 357 Energy = 7.238GeV180 181 beamFluence in wall = 3999 0 in cavity = 10 Icav/Iwall = 1.31170 energyFluence = 38.12063MeV*cm2/mg182 183 Total edep in cavity = 9.95706 MeV +- 6.35210%184 Total dose in cavity = 49.78531 MeV*cm2/mg +- 6.35210%185 186 DoseCavity/EnergyFluence = 1. 30599 +- 0.08296187 188 Total charged trackLength in cavity = 40.0185 m (mean value = 3.5234mm )189 190 StepSize of ch. tracks in wall = 1 83.6241 um +- 185.5614 um (nbSteps/track = 23.8988)191 StepSize of ch. tracks in cavity = 7 41.3298 um +- 3.0344 mm (nbSteps/track = 4.7528)182 Enter --> nbParticles = 11019 Energy = 7.036 GeV 183 Exit --> nbParticles = 11015 Energy = 7.034 GeV 184 185 beamFluence in wall = 39994 in cavity = 6 Icav/Iwall = 0.78695 energyFluence = 38.12014 MeV*cm2/mg 186 187 Total edep in cavity = 7.76475 MeV +- 5.14221 % 188 Total dose in cavity = 38.82375 MeV*cm2/mg +- 5.14221 % 189 190 DoseCavity/EnergyFluence = 1.01846 +- 0.05237 191 192 Total charged trackLength in cavity = 33.1218 m (mean value = 3.0059 mm ) 193 194 StepSize of ch. tracks in wall = 158.8831 um +- 182.9748 um (nbSteps/track = 27.6011) 195 StepSize of ch. tracks in cavity = 713.5388 um +- 1.1439 mm (nbSteps/track = 4.2126) 192 196 193 197 --------- Ranecu engine status --------- 194 198 Initial seed (index) = 0 195 Current couple of seeds = 580686433, 475796940199 Current couple of seeds = 991846787, 685450694 196 200 ---------------------------------------- 197 201 UserDetectorConstruction deleted. … … 201 205 G4 kernel has come to Quit state. 202 206 EventManager deleted. 203 Default detector region deleted.204 207 UImanager deleted. 205 208 Units table cleared. … … 207 210 RunManagerKernel is deleted. 208 211 RunManager is deleting. 212 -
trunk/examples/extended/medical/fanoCavity2/include/DetectorConstruction.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/DetectorMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/EventAction.hh
r1230 r1337 25 25 // 26 26 // $Id: EventAction.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/EventActionMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/HistoManager.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.hh,v 1.3 2007/11/12 18:19:30 maire 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/fanoCavity2/include/HistoMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.hh,v 1.2 2007/11/12 18:19:30 maire 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/fanoCavity2/include/MyMollerBhabhaModel.hh
r1230 r1337 25 25 // 26 26 // $Id: MyMollerBhabhaModel.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/PhysicsList.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsList.hh,v 1.2 2009/10/31 18:05:56 maire 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/fanoCavity2/include/PhysicsListMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.hh,v 1.2 2009/10/31 18:05:56 maire 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/fanoCavity2/include/PrimaryGeneratorAction.hh
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/ProcessesCount.hh
r1230 r1337 26 26 // 27 27 // $Id: ProcessesCount.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/RunAction.hh
r1230 r1337 25 25 // 26 26 // $Id: RunAction.hh,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/include/StepMax.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMax.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/StepMaxMessenger.hh
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/SteppingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.hh,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/include/SteppingVerbose.hh
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.hh,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/include/TrackingAction.hh
r1230 r1337 25 25 // 26 26 // $Id: TrackingAction.hh,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/run01.mac
r1230 r1337 1 # $Id: run01.mac,v 1. 3 2009/10/31 18:05:56maire Exp $1 # $Id: run01.mac,v 1.4 2010/01/20 15:41:14 maire Exp $ 2 2 # 3 3 /control/verbose 2 4 4 /run/verbose 2 5 5 # 6 ###/testem/phys/addPhysics standard_opt0 7 ###/testem/phys/addPhysics standard_opt2 6 8 /testem/phys/addPhysics standard_opt3 9 ###/testem/phys/addPhysics standard_GS 10 ###/testem/phys/addPhysics standard_WVI 11 ###/testem/phys/addPhysics standard_SS 7 12 # 8 13 /run/initialize -
trunk/examples/extended/medical/fanoCavity2/src/DetectorConstruction.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorConstruction.cc,v 1.2 2007/11/05 13:19:16 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 29 29 // -
trunk/examples/extended/medical/fanoCavity2/src/DetectorMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: DetectorMessenger.cc,v 1.2 2007/11/05 13:19:16 maire 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/fanoCavity2/src/EventAction.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: EventAction.cc,v 1. 1 2007/10/15 16:20:23 maireExp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: EventAction.cc,v 1.3 2010/06/07 05:40:46 perl Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 37 37 38 38 #include "G4Event.hh" 39 #include "G4TrajectoryContainer.hh"40 #include "G4Trajectory.hh"41 #include "G4VVisManager.hh"42 39 43 40 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 69 66 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... 70 67 71 void EventAction::EndOfEventAction(const G4Event* evt)68 void EventAction::EndOfEventAction(const G4Event*) 72 69 { 73 //trajectories74 //75 if (G4VVisManager::GetConcreteInstance())76 {77 G4TrajectoryContainer* trajectoryContainer = evt->GetTrajectoryContainer();78 G4int n_trajectories = 0;79 if (trajectoryContainer) n_trajectories = trajectoryContainer->entries();80 for(G4int i=0; i<n_trajectories; i++)81 { G4Trajectory* trj = (G4Trajectory*)82 ((*(evt->GetTrajectoryContainer()))[i]);83 if (drawFlag == "all") trj->DrawTrajectory(1000);84 else if ((drawFlag == "charged")&&(trj->GetCharge() != 0.))85 trj->DrawTrajectory();86 }87 }88 70 } 89 71 -
trunk/examples/extended/medical/fanoCavity2/src/EventActionMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: EventActionMessenger.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/HistoManager.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoManager.cc,v 1.5 2008/09/28 16:15:39 maire 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/fanoCavity2/src/HistoMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: HistoMessenger.cc,v 1.2 2007/11/12 18:19:30 maire 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/fanoCavity2/src/MyMollerBhabhaModel.cc
r1230 r1337 25 25 // 26 26 // $Id: MyMollerBhabhaModel.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/PhysicsList.cc
r1230 r1337 24 24 // ******************************************************************** 25 25 // 26 // $Id: PhysicsList.cc,v 1. 7 2009/10/31 18:05:01maire Exp $27 // GEANT4 tag $Name: geant4-09-0 3-cand-01 $26 // $Id: PhysicsList.cc,v 1.8 2010/01/20 15:41:14 maire Exp $ 27 // GEANT4 tag $Name: geant4-09-04-beta-01 $ 28 28 // 29 29 //....oooOO0OOooo........oooOO0OOooo........oooOO0OOooo........oooOO0OOooo...... … … 34 34 35 35 #include "PhysListEmStandard_option0.hh" 36 #include "PhysListEmStandard_option2.hh" 36 37 #include "PhysListEmStandard_option3.hh" 37 38 #include "PhysListEmStandard_GS.hh" 39 #include "PhysListEmStandard_WVI.hh" 38 40 #include "PhysListEmStandard_SS.hh" 39 41 … … 149 151 emPhysicsList = new PhysListEmStandard_option0(name,detector); 150 152 153 } else if (name == "standard_opt2") { 154 155 emName = name; 156 delete emPhysicsList; 157 emPhysicsList = new PhysListEmStandard_option2(name,detector); 158 151 159 } else if (name == "standard_opt3") { 152 160 … … 154 162 delete emPhysicsList; 155 163 emPhysicsList = new PhysListEmStandard_option3(name,detector); 156 164 157 165 } else if (name == "standard_GS") { 158 166 … … 160 168 delete emPhysicsList; 161 169 emPhysicsList = new PhysListEmStandard_GS(name,detector); 162 170 171 } else if (name == "standard_WVI") { 172 173 emName = name; 174 delete emPhysicsList; 175 emPhysicsList = new PhysListEmStandard_WVI(name,detector); 176 163 177 } else if (name == "standard_SS") { 164 178 165 179 emName = name; 166 180 delete emPhysicsList; 167 emPhysicsList = new PhysListEmStandard_SS(name,detector); 181 emPhysicsList = new PhysListEmStandard_SS(name,detector); 182 168 183 } else { 169 184 -
trunk/examples/extended/medical/fanoCavity2/src/PhysicsListMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: PhysicsListMessenger.cc,v 1.2 2009/10/31 18:05:01 maire 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/fanoCavity2/src/PrimaryGeneratorAction.cc
r1230 r1337 25 25 // 26 26 // $Id: PrimaryGeneratorAction.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/RunAction.cc
r1230 r1337 25 25 // 26 26 // $Id: RunAction.cc,v 1.3 2007/11/05 13:19:16 maire 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/fanoCavity2/src/StepMax.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMax.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/StepMaxMessenger.cc
r1230 r1337 25 25 // 26 26 // $Id: StepMaxMessenger.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/SteppingAction.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingAction.cc,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/src/SteppingVerbose.cc
r1230 r1337 25 25 // 26 26 // $Id: SteppingVerbose.cc,v 1.1 2007/10/15 16:20:23 maire 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/fanoCavity2/src/TrackingAction.cc
r1230 r1337 25 25 // 26 26 // $Id: TrackingAction.cc,v 1.2 2007/10/31 16:16:20 maire 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/fanoCavity2/vis.mac
r807 r1337 10 10 /run/initialize 11 11 # 12 # Create empty scene ("world" is default)13 /vis/ scene/create12 # Use this open statement to create an OpenGL view: 13 /vis/open OGL 600x600-0+0 14 14 # 15 # Create a scene handler for a specific graphics system 16 # Edit the next line(s) to choose another graphic system 17 # 18 /vis/open OGLIX 19 # 15 # Use this open statement to create a .prim file suitable for 16 # viewing in DAWN: 20 17 #/vis/open DAWNFILE 21 18 # 22 # Draw scene 23 /vis/viewer/set/viewpointThetaPhi 90 180 deg 19 # Use this open statement to create a .heprep file suitable for 20 # viewing in HepRApp: 21 #/vis/open HepRepFile 22 # 23 # Use this open statement to create a .wrl file suitable for 24 # viewing in a VRML viewer: 25 #/vis/open VRML2FILE 26 # 27 # Disable auto refresh and quieten vis messages whilst scene and 28 # trajectories are established: 29 /vis/viewer/set/autoRefresh false 30 /vis/verbose errors 31 # 32 # Draw geometry: 33 /vis/drawVolume 34 # 35 # Specify view angle: 36 /vis/viewer/set/viewpointThetaPhi 90. 180. 37 # 38 # Specify zoom value: 24 39 /vis/viewer/zoom 1.4 25 /vis/viewer/flush26 40 # 27 # for drawing the tracks 28 # if too many tracks cause core dump => storeTrajectory 0 29 /tracking/storeTrajectory 1 41 # Specify style (surface or wireframe): 42 #/vis/viewer/set/style wireframe 43 # 44 # Draw coordinate axes: 45 #/vis/scene/add/axes 0 0 0 1 m 46 # 47 # Draw smooth trajectories at end of event, showing trajectory points 48 # as markers 2 pixels wide: 49 /vis/scene/add/trajectories smooth 50 /vis/modeling/trajectories/create/drawByCharge 51 /vis/modeling/trajectories/drawByCharge-0/default/setDrawStepPts true 52 /vis/modeling/trajectories/drawByCharge-0/default/setStepPtsSize 1 53 # (if too many tracks cause core dump => /tracking/storeTrajectory 0) 54 # 55 # Draw hits at end of event: 56 #/vis/scene/add/hits 57 # 58 # To draw only gammas: 59 #/vis/filtering/trajectories/create/particleFilter 60 #/vis/filtering/trajectories/particleFilter-0/add gamma 61 # 62 # To invert the above, drawing all particles except gammas, 63 # keep the above two lines but also add: 64 #/vis/filtering/trajectories/particleFilter-0/invert true 65 # 66 # Many other options are available with /vis/modeling and /vis/filtering. 67 # For example, to select colour by particle ID: 68 #/vis/modeling/trajectories/create/drawByParticleID 69 #/vis/modeling/trajectories/drawByParticleID-0/set e- blue 70 # 71 # To superimpose all of the events from a given run: 30 72 /vis/scene/endOfEventAction accumulate 31 73 # 32 /testem/event/drawTracks all 33 74 # Re-establish auto refreshing and verbosity: 75 /vis/viewer/set/autoRefresh true 76 /vis/verbose warnings 77 # 78 # For file-based drivers, use this to create an empty detector view: 79 #/vis/viewer/flush
Note: See TracChangeset
for help on using the changeset viewer.