Difference between revisions of "G4Beamline PbBi"
Line 72: | Line 72: | ||
=== 2mm thick PbBi, 10 MeV, 1 cm cylindrical incident electron distribution=== | === 2mm thick PbBi, 10 MeV, 1 cm cylindrical incident electron distribution=== | ||
+ | |||
+ | [[File:Setup.png| 400 px]] | ||
[[File:ElectronTracks.png| 200 px]][[File:PhotonTracks.png| 200 px]] | [[File:ElectronTracks.png| 200 px]][[File:PhotonTracks.png| 200 px]] |
Revision as of 14:58, 1 May 2015
Development of a Positron source using a PbBi converter and a Solenoid
Converter target properties
Definition of Lead Bismuth
1cm diameter target
2 mm thick PbBi
0.5 Tesla solenoid
G4BeamLine and MCNPX
Target thickness optimization
First simple test is to send 1 million, 10 MeV electrons towards a PbBi target and count how many positrons leave the downstream side
The Random number seed is set by Time in G4beamline to use a different set of pseudo random numbers each time it is run
The G4Beamlin incident electron beam has the following properties
beam gaussian particle=e- nEvents=1000000 beamZ=0.0 sigmaX=1.0 sigmaY=1.0 sigmaXp=0.100 sigmaYp=0.100 meanMomentum=10.0 sigmaP=4.0 meanT=0.0 sigmaT=0.0
PbBi Thickness (mm) | #positrons/million electrons (G4Beamline) | #positrons/million electrons (MCNPX) |
1 | 1091 | |
1.5 | 1728 | |
2 | 1873,1945 | 431984 |
2.5 | 2062 | |
3 | 13 | 1986 |
3.5 | 1938 | |
4 | 39 | 1858 |
5 | 1646 | |
6 | 37 | 1541 |
10 | 1216 |
Comparison of G4Beamline and MCNPX
Energy Distribution
Angular distribution of positrons
2mm thick PbBi, 10 MeV, 1 cm cylindrical incident electron distribution
I was unable to do anything other than a gaussian beam right now, I will try to do one later
For now I have a gaussian with an 8mm RMS and 10 MeV incident electrons as shown below.
The positron and electron momentum distributions after the PbBi converter are shown below
A comma delimited text file with the above events in the format of
x,y,z,Px,Py,Pz
in units of cm for distance and MeV for momentum is located at
for positrons
http://www2.cose.isu.edu/~foretony/Positrons_2mm10MeV.dat
and
http://www2.cose.isu.edu/~foretony/Electrons_2mm10MeV.dat
for electrons
2mm thick PbBi, 10 MeV, point source
I was unable to do anything other than a gaussian beam right now, I will try to do one later
For now I have a gaussian with an 1mm RMS and 10 MeV incident electrons as shown below.
The positron and electron momentum distributions after the PbBi converter are shown below
A comma delimited text file with the above events in the format of
x,y,z,Px,Py,Pz
in units of cm for distance and MeV for momentum is located at
for positrons
http://www2.cose.isu.edu/~foretony/Positrons_2mm10MeVDelta.dat
and
http://www2.cose.isu.edu/~foretony/Electrons_2mm10MeVDelta.dat
for electrons