New IAC Wiki - User contributions [en]

CodyMilne G4Proj

2016-05-28T19:35:14Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
Research shows this says these values are on a "log scale", not sure how exactly to read this.

=Simulation Set Up - Detector Construction =
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]]

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.

==100 keV==
Start with <math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with <math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"No Backscattering" Do another plot with momentum cut Pz>0

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

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with <math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with <math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:34:34Z

Milncody: /* Running Types */

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
Research shows this says these values are on a "log scale", not sure how exactly to read this.

=Simulation Set Up - Detector Construction =
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]]

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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"No Backscattering" Do another plot with momentum cut Pz>0

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

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:34:13Z

Milncody: /* 500 keV */

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
Research shows this says these values are on a "log scale", not sure how exactly to read this.

=Simulation Set Up - Detector Construction =
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]]

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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"No Backscattering" Do another plot with momentum cut Pz>0

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

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:33:50Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
Research shows this says these values are on a "log scale", not sure how exactly to read this.

=Simulation Set Up - Detector Construction =
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]]

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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:31:14Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01
Research shows this says these values are on a "log scale", not sure how exactly to read this.

=Simulation Set Up - Detector Construction =
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]]

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=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:28:45Z

Milncody: /* Runing Types */

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

=Simulation Set Up - Detector Construction =
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]]

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=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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.

==Running Types==

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

For parallel processing on plexi:
Run
module load geant/geant4-10.2.0
when you first log in to set up the environment for Geant4
Use the command
G4make
when you configure your application (instead of 'cmake') to provide the necessary default options for Plexi.
To run:
salloc -n <number_of_processes> mpirun /absolute/path/to/your/application

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:21:22Z

Milncody: /* 10 MeV */

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

=Simulation Set Up - Detector Construction =
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]]

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=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

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==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:19:43Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

=Simulation Set Up - Detector Construction =
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]]

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=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

File:CM Detectorsetup.png

2016-05-28T19:15:12Z

Milncody:

CodyMilne G4Proj

2016-05-28T19:15:01Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

=Simulation Set Up - Detector Construction =
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=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:09:13Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

=Simulation Set Up - Detector Construction =
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.

=Results of energy loss in Germanium=

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-05-28T19:05:00Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium ------ [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium ------------- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-04-01T03:46:32Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hourp
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|-
|3/18/16||1.5 hours
|-
|3/30/16||2 hours
|-
|3/31/16||1 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-30T22:00:05Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hourp
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|-
|3/18/16||1.5 hours
|-
|3/30/16||2 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-23T07:31:27Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hourp
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|-
|3/18/16||1.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-19T21:42:05Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|-
|3/18/16||1.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-18T19:47:58Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|
|3/18/16||1.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-11T20:38:57Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|-
|3/11/16||2.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

File:CM 10Menergydepositionzoom 3-11-2016.png

2016-03-11T20:38:14Z

Milncody:

File:CM 10Menergydeposition 3-11-2016.png

2016-03-11T20:38:03Z

Milncody:

CodyMilne G4Proj

2016-03-11T20:37:49Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==10 MeV==
This is looking at total energy deposition from the child particles in Germanium, from a incident gamma shot from within the crystal.

[[File:CM_10Menergydeposition_3-11-2016.png | 200 px]]
[[File:CM_10Menergydepositionzoom_3-11-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-11T05:32:35Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|-
|3/10/16||1.5 hours
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-07T19:41:20Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.
1MeV: -2.852121 slope: -3.04127e-02 +/- 4.63910e-02
500k: -3.66350 slope: -4.25739e-02 +/- 1.10890e-01

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-07T03:22:25Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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. Research shows this says these values are on a "log scale", not sure how exactly to read this.
1MeV: -2.852121 slope: -3.04127e-02
500k: -3.66350 slope: -4.25739e-02

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|-
|3/06/16||1 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-07T03:21:22Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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. Research shows this says these values are on a "log scale", not sure how exactly to read this.
1MeV: -2.852121 slope: -3.04127e-02
500k: -3.66350 slope: -4.25739e-02

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-03-07T03:13:22Z

Milncody:

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|-
|3/02/16||1 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-21T02:04:32Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

1)Determine total energy scoring in the Germanium block.

2) 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.

3) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-21T02:04:17Z

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

=Photon Energy deposition in a pure Ge crystal=

Here is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Here is the total photon attenuation for Germanium --- [[File:photoattenuationGe.png | 200 px]]

Here is the range of electrons in Ge as a function of energy --- [[File:eRangeGe.png | 200 px]]

Determine total energy scoring in the Germanium block.

1) 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.

2) 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.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==
1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==

Start with<math> E_{gamma}</math> = 100 eV

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

==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

=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.

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first hit position distribution is consistent with published attenuation coefficient.

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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]]

Initial verification - The first gamma interaction distribution is consistent with the published Ge mass attenuation coefficient at 500keV (0.08212) and a Ge density of 5.323.

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

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

Are the 2-D plots energy weights (normalized)? Set the Z-axis color range by hand in cold.

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction is consistent with published attenuation coefficient.

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

Initial verification - first interaction distribution is consistent with published attenuation coefficient.

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:32:15Z

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

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|2/16/16||2.5 hour
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:32:01Z

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

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}
|2/16/16||2.5 hour
|-

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:31:21Z

Milncody: /* 500 keV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

File:CM 3MeVPOSZvsEVENTSNORM 2-14-2016.png

2016-02-17T05:30:55Z

Milncody:

File:CM 1MeVPOSZvsEVENTSNORM 2-14-2016.png

2016-02-17T05:30:41Z

Milncody:

File:CM 500keVeLoss 2-7-2016.png

2016-02-17T05:30:16Z

Milncody: uploaded a new version of "File:CM 500keVeLoss 2-7-2016.png"

File:CM 500keVPOSZvsPOSXColz 2-14-2016.gif

2016-02-17T05:29:31Z

Milncody: uploaded a new version of "File:CM 500keVPOSZvsPOSXColz 2-14-2016.gif"

File:CM 100keVPOSZvsEVENTSNORM 2-14-2016.png

2016-02-17T05:28:49Z

Milncody:

CodyMilne G4Proj

2016-02-17T05:00:50Z

Milncody: /* 3 MeV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]
[[File:CM_3MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:00:33Z

Milncody: /* 1 MeV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

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

[[File:CM_1MeVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:00:18Z

Milncody: /* 100 keV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

[[File:CM_100keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T05:00:07Z

Milncody: /* 10 keV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_10keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T04:59:50Z

Milncody: /* 1 keV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_1keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-17T04:59:36Z

Milncody: /* 1 keV */

=Photon Energy deposition in a pure Ge crystal=

Here are is the mass attenuation cross sections for Germanium --- [[File:massattenuation.gif | 200 px]]

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

[[File:CM_1keVeLoss_2-7-2016.png | 200 px]]
[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

File:CM 500keVPOSZvsPOSXColz 2-14-2016.gif

2016-02-14T23:34:26Z

Milncody:

File:CM 500keVPOSZvsEVENTSNORM 2-14-2016.png

2016-02-14T23:34:03Z

Milncody:

CodyMilne G4Proj

2016-02-14T23:33:49Z

Milncody: /* 500 keV */

=Photon Energy deposition in a pure Ge crystal=

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

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

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

[[File:CM_500keVPOSZvsEVENTSNORM_2-14-2016.png | 200 px]]

[[File:CM_500keVPOSZvsPOSXColz_2-14-2016.gif | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-14T23:32:52Z

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

=Photon Energy deposition in a pure Ge crystal=

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|2/14/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

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

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-12T05:02:26Z

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

=Photon Energy deposition in a pure Ge crystal=

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|2/11/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

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

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

CodyMilne G4Proj

2016-02-12T05:02:02Z

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

=Photon Energy deposition in a pure Ge crystal=

Determine the deposition of photons ranging in energy from 10 eV to 3 MeV in a pure Germanium crystal assuming an infinite block.

Attenuation coefficient about 1/5.323 cm = 0.19 cm or about 2 mm to maybe 2 cm?

Create file containing Photon Energy, ProcessID, step #, Energy Lost, Position for each step a photon takes through the Germanium crystal.

{| border="1" |cellpadding="20" cellspacing="0
|-
| Date || Time
|-
|1/26/16||1 hour
|-
|1/27/16||1 hour
|-
|1/29/16||1 hour
|-
|2/4/16||1 hour
|-
|2/5/16||3 hour
|-
|2/8/16||2 hour
|-
|2/9/16||2 hour
|-
|2/10/16||1 hour
|-
|-2/11/16||2 hour
|-
|}

=Results of energy loss in Germanium=

Why do SN>1 photons loose energy at Z = -9000?

==10 eV==

1-D plot of events for each step as function of Z

1-D plot of the energy deposited for each step as function of Z

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

==100 eV==
Start with<math> E_{gamma}</math> = 100 eV

[[File:CM_100eVeLoss_2-7-2016.png | 200 px]]
==1 keV==
Start with<math> E_{gamma}</math> = 1 keV

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

==10 keV==
Start with<math> E_{gamma}</math> = 10 keV

[[File:CM_10keVeLoss_2-7-2016.png | 200 px]]
==100 keV==
Start with<math> E_{gamma}</math> = 100 keV

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

==500 keV==

Start with<math> E_{gamma}</math> = 500 keV

Change units on Y-axis to Energy (eV?)

"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_500keVPOSZvsEVENTScut_2-7-2016.png | 200 px]]

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

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

==1 MeV==
Start with<math> E_{gamma}</math> = 1 MeV

[[File:CM_1MeVeLoss_2-7-2016.png | 200 px]]
==3 MeV==
Start with<math> E_{gamma}</math> = 3 MeV

[[File:CM_13MeVeLoss_2-8-2016.png | 200 px]]

==Runing Types==

nohup ../exampleN02 run1.mac > /dev/null &

Divide canvas into 2 x 5

[[User_talk:Foretony#CodyMilne_G4Proj]]

File:CM 500keVPOSZvsPEL 2-7-2016.png

2016-02-12T05:00:40Z

Milncody: uploaded a new version of "File:CM 500keVPOSZvsPEL 2-7-2016.png"

File:CM 500keVPOSZvsEVENTScut 2-7-2016.png

2016-02-12T05:00:13Z

Milncody: