G4Beamline PbBi

From New IAC Wiki
Jump to navigation Jump to search

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 [math]\pm[/math] 43 1984
2.5 2062
3 [math]\pm[/math] 13 1986
3.5 1938
4 [math]\pm[/math] 39 1858
5 1646
6 [math] \pm[/math] 37 1541
10 1216

Comparison of G4Beamline and MCNPX


Comparison.png


Energy Distribution

TF PosE 04-28-15.png Positrons2.png

Angular distribution of positrons

TF Theta 04-28-15.png


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


4-30-2015 PositronMomentum 2mm.png4-30-2015 ElectronMomentum 2mm.png

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.

4-30-2015 BeamPosDelta.png4-30-2015 ElectronMomentum.png


The positron and electron momentum distributions after the PbBi converter are shown below


4-30-2015 PositronMomentum 2mmDelta.png4-30-2015 ElectronMomentum 2mmDelta.png

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

Solenoid

Positrons#Simulations