Skip to content

Commit

Permalink
update readme, add some results form server, clean up whitespace
Browse files Browse the repository at this point in the history
  • Loading branch information
@studouglas
studouglas committed Feb 23, 2016
1 parent bc9b855 commit 2031f45
Show file tree
Hide file tree
Showing 5 changed files with 543 additions and 53 deletions.
35 changes: 35 additions & 0 deletions Examples/Hadr04/src/run_50events.mac
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,35 @@
# $Id$
#
# Macro file for "Hadr04.cc"
# (can be run in batch, without graphic)
#
# neutron 2 MeV; all processes
#
/control/verbose 2
/run/verbose 2
#
/testhadr/det/setMat Water_ts
/testhadr/det/setSize 1 m
#
/testhadr/phys/thermalScattering true
/random/setSeeds 9876 54321
#
/run/initialize
#
/process/list
#
/gun/particle neutron
/gun/energy 2 MeV
#
/analysis/setFileName Water_thermal_44
/analysis/h1/set 1 100 0. 100. none #nb colli >1eV
/analysis/h1/set 2 100 0. 100. cm #track len >1eV
/analysis/h1/set 3 100 0. 5. mus #time of flight >1eV
/analysis/h1/set 4 100 0. 1000. none #nb colli <1eV
/analysis/h1/set 5 100 0. 500. cm #track len <1eV
/analysis/h1/set 6 100 0. 1000. mus #time of flight <1eV
/analysis/h1/set 7 100 0. 500. meV #energy dist <1eV
#
/run/printProgress 10
#
/run/beamOn 50
3 changes: 2 additions & 1 deletion README.md
View file
Original file line number Diff line number Diff line change
Expand Up @@ -57,11 +57,12 @@ extra `\)` before the `;\n";`. The version in the repo is fixed, but if you
downloaded marshalgen from another source, you need to fix this)

**Install GEANT-4**<br>
If installing on McMaster's servers, add `. /opt/rh/devtoolset-3/enable` to your bash_profile to use the newer version of gcc.
1. (Fedora only) `yum install expat-devel`<br>
2. `mkdir /path/to/GEANT4-GPU/geant4.10.00.p02-build /path/to/GEANT4-GPU/
geant4.10.00.p02-install`<br>
3. `cd /path/to/GEANT4-GPU/geant4.10.00.p02-build`<br>
4. `cmake -DGEANT4_INSTALL_DATA=ON -DGEANT4_ENABLE_CUDA=ON -DUSE_SYSTEM_EXPAT=OFF -DCMAKE_INSTALL_PREFIX=/path/to/GEANT4-GPU/geant4.10.00.p02-install /path/to/GEANT4-GPU/geant4.10.00.p02`<br>
4. `cmake -DGEANT4_INSTALL_DATA=ON -DGEANT4_ENABLE_CUDA=ON -DGEANT4_USE_SYSTEM_EXPAT=OFF -DCMAKE_INSTALL_PREFIX=/path/to/GEANT4-GPU/geant4.10.00.p02-install /path/to/GEANT4-GPU/geant4.10.00.p02`<br>
IF installing on McMaster's server, you must add flag `-DCUDA_HOST_COMPILER=/usr/bin/g++`<br>
5. `make -jN` where `N` is the number of processors on your computer<br>
6. `make install`
Expand Down
224 changes: 224 additions & 0 deletions Results/Hadr04-Feb23_12pm-gpu1-CPU.txt
View file
Original file line number Diff line number Diff line change
@@ -0,0 +1,224 @@
*************************************************************
Geant4 version Name: geant4-10-02 (4-December-2015)
Copyright : Geant4 Collaboration
Reference : NIM A 506 (2003), 250-303
WWW : http://cern.ch/geant4
*************************************************************

/run/verbose 2
#
/testhadr/det/setMat Water_ts
/testhadr/det/setSize 1 m
/run/reinitializeGeometry
#
/testhadr/phys/thermalScattering true
/random/setSeeds 9876 54321
#
/run/initialize
userDetector->Construct() start.

The Box is 1 m of Water_ts

Material: Water_ts density: 1.000 g/cm3 RadL: 36.084 cm Nucl.Int.Length: 75.376 cm
Imean: 78.000 eV

---> Element: TS_H_of_Water (H) Z = 1.0 N = 1 A = 1.008 g/mole
---> Isotope: H1 Z = 1 N = 1 A = 1.01 g/mole abundance: 99.989 %
---> Isotope: H2 Z = 1 N = 2 A = 2.01 g/mole abundance: 0.011 %
ElmMassFraction: 11.19 % ElmAbundance 66.67 %

---> Element: Oxygen (O) Z = 8.0 N = 16 A = 16.000 g/mole
---> Isotope: O16 Z = 8 N = 16 A = 15.99 g/mole abundance: 99.757 %
---> Isotope: O17 Z = 8 N = 17 A = 17.00 g/mole abundance: 0.038 %
---> Isotope: O18 Z = 8 N = 18 A = 18.00 g/mole abundance: 0.205 %
ElmMassFraction: 88.81 % ElmAbundance 33.33 %

Water_ts is registered to the default region.
physicsList->Construct() start.
@@@ G4ParticleHPInelasticData instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5
NeutronHP: /Capture file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/8_17_Oxygen
NeutronHP: /Elastic file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/8_17_Oxygen
NeutronHP: /Inelastic file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic/CrossSection/8_17_Oxygen
NeutronHP: /Capture file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/6_nat_Carbon
NeutronHP: /Elastic file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/6_nat_Carbon
NeutronHP: /Inelastic file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic/CrossSection/6_nat_Carbon
/u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5
@@@ G4ParticleHPInelastic instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic
physicsList->CheckParticleList() start.
physicsList->setCut() start.
#
/process/list
Transportation, hadElastic, neutronInelastic, nCapture
nFission
#
/gun/particle neutron
/gun/energy 2 MeV
#
/analysis/setFileName Water_thermal_44
Set file name: Water_thermal_44
/analysis/h1/set 1 100 0. 100. none
/analysis/h1/set 2 100 0. 100. cm
/analysis/h1/set 3 100 0. 5. mus
/analysis/h1/set 4 100 0. 1000. none
/analysis/h1/set 5 100 0. 500. cm
/analysis/h1/set 6 100 0. 1000. mus
/analysis/h1/set 7 100 0. 500. meV
#
/run/printProgress 100
#
/run/beamOn 500
NeutronHP: /Capture file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/8_17_Oxygen
NeutronHP: /Elastic file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/8_17_Oxygen
NeutronHP: /Inelastic file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic/CrossSection/8_17_Oxygen
NeutronHP: /Capture file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/6_nat_Carbon
NeutronHP: /Elastic file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/6_nat_Carbon
NeutronHP: /Inelastic file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic/CrossSection/6_nat_Carbon

Neutron HP Thermal Scattering Data: Following material-element pairs and/or elements are registered.
Element TS_H_of_Water, internal thermal scattering id 0
Element TS_C_of_Graphite, internal thermal scattering id 1

NeutronHP: /Elastic file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/8_17_Oxygen
NeutronHP: /Elastic file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Elastic/CrossSection/6_nat_Carbon

Neutron HP Thermal Scattering: Following material-element pairs or elements are registered.
Element TS_H_of_Water, internal thermal scattering id 0
Element TS_C_of_Graphite, internal thermal scattering id 1


Neutron HP Thermal Scattering Data: Following material-element pairs and/or elements are registered.
Element TS_H_of_Water, internal thermal scattering id 0
Element TS_C_of_Graphite, internal thermal scattering id 1


====================================================================
HADRONIC PROCESSES SUMMARY (verbose level 1)

---------------------------------------------------
Hadronic Processes for neutron

Process: hadElastic
Model: NeutronHPElastic: 4 eV ---> 20 MeV
Model: NeutronHPThermalScattering: 0 meV ---> 4 eV
Cr_sctns: NeutronHPThermalScatteringData: 0 meV ---> 4 eV
Cr_sctns: NeutronHPElasticXS: 0 meV ---> 20 MeV
Cr_sctns: GheishaElastic: 0 meV ---> 100 TeV

Process: neutronInelastic
Model: NeutronHPInelastic: 0 meV ---> 20 MeV
Cr_sctns: NeutronHPInelasticXS: 0 meV ---> 20 MeV
Cr_sctns: GheishaInelastic: 0 meV ---> 100 TeV

Process: nCapture
Model: NeutronHPCapture: 0 meV ---> 20 MeV
Cr_sctns: NeutronHPCaptureXS: 0 meV ---> 20 MeV
Cr_sctns: GheishaCaptureXS: 0 meV ---> 100 TeV

Process: nFission
Model: NeutronHPFission: 0 meV ---> 20 MeV
Cr_sctns: NeutronHPFissionXS: 0 meV ---> 20 MeV
Cr_sctns: GheishaFissionXS: 0 meV ---> 100 TeV

================================================================
/u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5
@@@ G4ParticleHPInelastic instantiated for particle neutron data directory variable is G4NEUTRONHPDATA pointing to /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Inelastic
NeutronHP: /Capture file for Z = 8, A = 18 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/8_17_Oxygen
NeutronHP: /Capture file for Z = 6, A = 12 is not found and NeutronHP will use /u50/dougls2/GEANT4-GPU/geant4.10.02-install/share/Geant4-10.2.0/data/G4NDL4.5/Capture/CrossSection/6_nat_Carbon

Region <DefaultRegionForTheWorld> -- -- appears in <Water_ts> world volume
This region is in the mass world.
Root logical volume(s) : Water_ts
Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
Materials : Water_ts
Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm

Region <DefaultRegionForParallelWorld> -- -- is not associated to any world.
Root logical volume(s) :
Pointers : G4VUserRegionInformation[0], G4UserLimits[0], G4FastSimulationManager[0], G4UserSteppingAction[0]
Materials :
Production cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm

========= Table of registered couples ==============================

Index : 0 used in the geometry : Yes
Material : Water_ts
Range cuts : gamma 1 mm e- 1 mm e+ 1 mm proton 0 fm
Energy thresholds : gamma 2.94056 keV e- 351.877 keV e+ 342.545 keV proton 0 meV
Region(s) which use this couple :
DefaultRegionForTheWorld

====================================================================

Start closing geometry.
G4GeometryManager::ReportVoxelStats -- Voxel Statistics

Total memory consumed for geometry optimisation: 0 kByte
Total CPU time elapsed for geometry optimisation: 0 seconds
### Run 0 starts.

--------- Ranecu engine status ---------
Initial seed (index) = 0
Current couple of seeds = 9876, 54321
----------------------------------------
... open Root analysis file : Water_thermal_44.root - done
--> Event 0 starts.
--> Event 100 starts.
--> Event 200 starts.
--> Event 300 starts.
--> Event 400 starts.
Run terminated.
Run Summary
Number of events processed : 500
User=14.19s Real=14.2s Sys=0s

The run is 500 neutron of 2 MeV through 50 cm of Water_ts (density: 1 g/cm3 )

Process calls frequency :
hadElastic= 108818 nCapture= 500

Parcours of incident neutron:
nb of collisions E>1*eV= 16.63 E<1*eV= 202.01 total= 218.64
track length E>1*eV= 19.976 cm E<1*eV= 75.695 cm total= 95.671 cm
time of flight E>1*eV= 753.28 ns E<1*eV= 207.82 mus total= 208.57 mus

List of generated particles:
O16: 1512 Emean = 43.921 keV ( 2.1992 meV --> 446.59 keV)
O17: 2 Emean = 202.95 eV ( 187.15 eV --> 218.76 eV )
O18: 1 Emean = 9.3256 eV ( 9.3256 eV --> 9.3256 eV )
deuteron: 498 Emean = 1.3195 keV ( 1.2775 keV --> 1.7025 keV)
gamma: 502 Emean = 2.2188 MeV ( 870.8 keV --> 2.2244 MeV)
proton: 11474 Emean = 81.412 keV ( 0.14188 meV --> 1.9993 MeV)
... write Root file : Water_thermal_44.root - done

--------- Ranecu engine status ---------
Initial seed (index) = 0
Current couple of seeds = 1392330422, 48169388
----------------------------------------
G4 kernel has come to Quit state.
UserDetectorConstruction deleted.
UserPhysicsList deleted.
UserActionInitialization deleted.
UserRunAction deleted.
UserPrimaryGenerator deleted.
RunManager is deleting RunManagerKernel.
EventManager deleted.
Units table cleared.
Total navigation history collections cleaned: 4
================== Deleting memory pools ===================
Pool ID '20G4NavigationLevelRep', size : 0.00385 MB
Pool ID '24G4ReferenceCountedHandleIvE', size : 0.000961 MB
Pool ID '7G4Event', size : 0.000961 MB
Pool ID '15G4PrimaryVertex', size : 0.000961 MB
Pool ID '17G4PrimaryParticle', size : 0.000961 MB
Pool ID '17G4DynamicParticle', size : 0.0154 MB
Pool ID '7G4Track', size : 0.0308 MB
Pool ID '18G4TouchableHistory', size : 0.000961 MB
Pool ID '15G4CountedObjectIvE', size : 0.000961 MB
Pool ID '17G4ReactionProduct', size : 0.000961 MB
Number of memory pools allocated: 10; of which, static: 0
Dynamic pools deleted: 10 / Total memory freed: 0.057 MB
============================================================
G4Allocator objects are deleted.
UImanager deleted.
StateManager deleted.
RunManagerKernel is deleted. Good bye :)
Loading

0 comments on commit 2031f45

Please sign in to comment.