CodyMilne G4Proj - Revision history

Milncody at 19:35, 28 May 2016

2016-05-28T19:35:14Z

Milncody: /* Running Types */

2016-05-28T19:34:34Z

Running Types

Milncody: /* 500 keV */

2016-05-28T19:34:13Z

500 keV

Milncody at 19:33, 28 May 2016

2016-05-28T19:33:50Z

Milncody: /* Photon Energy deposition in a pure Ge crystal */

2016-05-28T19:31:14Z

Photon Energy deposition in a pure Ge crystal

Milncody: /* Runing Types */

2016-05-28T19:28:45Z

Runing Types

Milncody: /* 10 MeV */

2016-05-28T19:21:22Z

10 MeV

Milncody: /* Simulation Set Up - Detector Construction */

2016-05-28T19:19:43Z

Simulation Set Up - Detector Construction

Milncody at 19:15, 28 May 2016

2016-05-28T19:15:01Z

Milncody at 19:09, 28 May 2016

2016-05-28T19:09:13Z

@@ Line 29: / Line 29: @@
 ==100 keV==
-Start with<math> E_{gamma}</math> = 100 keV
+Start with <math> E_{gamma}</math> = 100 keV
@@ Line 40: / Line 40: @@
 ==500 keV==
-Start with<math> E_{gamma}</math> = 500 keV
+Start with <math> E_{gamma}</math> = 500 keV
 Change units on Y-axis to Energy (eV?)
@@ Line 61: / Line 61: @@
 ==1 MeV==
-Start with<math> E_{gamma}</math> = 1 MeV
+Start with <math> E_{gamma}</math> = 1 MeV
@@ Line 71: / Line 71: @@
 ==3 MeV==
-Start with<math> E_{gamma}</math> = 3 MeV
+Start with <math> E_{gamma}</math> = 3 MeV

← Older revision		Revision as of 19:34, 28 May 2016
Line 90:		Line 90:

	To submit batch job to minerve:		To submit batch job to minerve:
−	nohup ../exampleN02 run1.mac > /dev/null &	+	nohup ../exampleN02 run1.mac > /dev/null &

← Older revision		Revision as of 19:34, 28 May 2016
Line 45:		Line 45:

	"No Backscattering" Do another plot with momentum cut Pz>0		"No Backscattering" Do another plot with momentum cut Pz>0
−
−	~~Add 9000 to Z-position to start Ge barrier at zero~~

	[[File:CM_500keVPOSZvsEVENTS_2-7-2016.png \| 200 px]]		[[File:CM_500keVPOSZvsEVENTS_2-7-2016.png \| 200 px]]

@@ Line 23: / Line 23: @@
 Photons are fired from a cylindrical gps from within the inner radius of the lead cask. This cylinder has a height of 40.25 cm and a radius of 0.5 cm. The goal of this gun is to simulate decay from a fuel rod on the inside of the lead cylinder. Angular cuts are made on the gun to limit wasteful gammas from being calculated. There is a phi (x direction) cut of 22 degrees, centered on the hole leading to the Germanium block, and a theta cut of about 45 degrees. The phi coincides with the 1mm side of the hole, and the theta coincides with the 1cm side of the hole.
 =Energy Deposition=
-Send photons from inside an infinite crystal at 100, 500, 1000,3000 keV energy.  Add up all the energy loss and see if you get delta functions or a delta + tail because GEANT4 stopped taking the photon and thus the energy is not accounted for.
+Send photons from inside an infinite crystal at 100, 500, 1000, 3000 keV energy.  Add up all the energy loss and see if you get delta functions or a delta + tail because GEANT4 stopped taking the photon and thus the energy is not accounted for.
 ==100 keV==

@@ Line 9: / Line 9: @@
 Perform a fit to the 2nd hit distribution using an exponential fitting function.  Do this for 100keV, 500keV, 1MeV, 3Mev and for this test, please run some 10 MeV gammas.
 —the result of the fit should be equal to the mass attenuation coefficient multiplied by the density of Ge.  Therefore, dividing the fit result by 5.323 yields the mass attenuation coefficient at each energy and we can compare them directly with the NIST database values for validation of the GEANT4 total cross sections and sampling.
 Using S2WN.Fit("expo");
-I got for 3MeV constant: -1.88934 and slope: -1.84926e-02 +/- 1.63709e-02. Research shows this says these values are on a "log scale", not sure how exactly to read this.
+   I got for
-MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
+MeV constant: -1.88934 and slope: -1.84926e-02 +/- 1.63709e-02.
-k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
+MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
+k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
-) for the same events (2nd hits), make a plot of the number of times that interaction was photoelectric, compton or pair production. Make the same 3-bin plot for each of the incident gamma energies. The relative rates of the first interactions evolve as a function of incident photon energy (see the plot of linear absorption coefficient on the wiki).
+   Research shows this says these values are on a "log scale", not sure how exactly to read this.
 =Simulation Set Up - Detector Construction =

@@ Line 120: / Line 120: @@
 Confirmed, GEANT keeps track of all energy deposited in the crystal, with little error that scales with initial photon energy. Could be minor leakage of secondary photons.
-==Runing Types==
+==Running Types==
 nohup ../exampleN02 run1.mac > /dev/null &
-Divide canvas into 2 x 5
+For parallel processing on plexi:
+   Run
 [[User_talk:Foretony#CodyMilne_G4Proj]]

@@ Line 117: / Line 117: @@
 [[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
 [[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]
 ==Runing Types==

← Older revision		Revision as of 19:19, 28 May 2016
Line 23:		Line 23:

	[[File:CM_Detectorsetup.png \| 200 px]]		[[File:CM_Detectorsetup.png \| 200 px]]
		+
		+	Photons are fired from a cylindrical gps from within the inner radius of the lead cask. This cylinder has a height of 40.25 cm and a radius of 0.5 cm. The goal of this gun is to simulate decay from a fuel rod on the inside of the lead cylinder. Angular cuts are made on the gun to limit wasteful gammas from being calculated. There is a phi (x direction) cut of 22 degrees, centered on the hole leading to the Germanium block, and a theta cut of about 45 degrees. The phi coincides with the 1mm side of the hole, and the theta coincides with the 1cm side of the hole.

	=Results of energy loss in Germanium=		=Results of energy loss in Germanium=

@@ Line 22: / Line 22: @@
 A pure block of 50x50x50cm Germanium is placed in front of a lead cylinder, inner radius 0.75cm, outer radius 50cm of height 55cm. There is a 1mmx1cm hole (made of Argon) from the center of the lead cask to the outside, pointed toward the middle of the Germanium crystal.
+[[File:CM_Detectorsetup.png | 200 px]]
 =Results of energy loss in Germanium=

@@ Line 18: / Line 18: @@
 ) for the same events (2nd hits), make a plot of the number of times that interaction was photoelectric, compton or pair production. Make the same 3-bin plot for each of the incident gamma energies. The relative rates of the first interactions evolve as a function of incident photon energy (see the plot of linear absorption coefficient on the wiki).
 —The ratios of event types can be compared with published values.  I think these "spot checks" at those five energies are all we need to complete our initial validation of the GEANT4 sub-sampling.
 =Results of energy loss in Germanium=
 ==10 eV==