Difference between revisions of "LB January Flux Simulations"
| Line 33: | Line 33: | ||
[[File:LB JanFluxSim RotatedView.png|200px]] | [[File:LB JanFluxSim RotatedView.png|200px]] | ||
| − | + | [[File:LB NickelFlux G4Sim FullView Jan2018.png|200px]] | |
| + | |||
| + | [[File:LB SeFlux G4Sim FullView Jan2018.png|200px]] | ||
| + | |||
| + | [[File:LB NickelFlux G4Sim Jan2018.png|200px]] | ||
| + | |||
| + | [[File:LB SeFlux G4Sim Jan2018.png|200px]] | ||
| + | |||
| − | |||
There were some strange issues while performing this simulation. The first issue that stuck out was that when trying to create a G4Orb for the selenium pellet, the program would compile, but the orb would not be present in the visualizer and furthermore would not have particles written to the outfile. A similar issue happened when trying to curve the nickel foil and the aluminum cylinder. If the nickel foil and the aluminum cylinder were curved, perhaps the number of particles passing through them would be lower which would raise the ratio of the number of particles. | There were some strange issues while performing this simulation. The first issue that stuck out was that when trying to create a G4Orb for the selenium pellet, the program would compile, but the orb would not be present in the visualizer and furthermore would not have particles written to the outfile. A similar issue happened when trying to curve the nickel foil and the aluminum cylinder. If the nickel foil and the aluminum cylinder were curved, perhaps the number of particles passing through them would be lower which would raise the ratio of the number of particles. | ||
Revision as of 23:59, 24 January 2018
Simulations were run using geant4 to try and see if the ratio of the flux between the nickel foil and the inner selenium pellet would explain why the ratios of the rates are not correct. A general particle source was used to fire electrons in a beam at a titanium beam window, which then proceeds to an aluminum/tungsten radiator. After that the particles will encounter a nickel foil followed by an aluminum sheet and finally a selenium pellet that is -0.75cm off of the beam axis (to try and replicate the location of the foil relative to the selenium pellet). The run.mac file is shown below
/gps/particle e-
/gps/position 0 0 -5 m
/gps/pos/type Plane
/gps/pos/shape Circle
/gps/pos/radius 0.5 cm
/tracking/verbose 2
/gps/direction 0 0 1
/gps/verbose 2
/gps/pos/sigma_r 0.5 cm
/gps/energy 32 MeV
/run/beamOn 50000
The dimensions of the pieces of interest are:
Se Pellet: 0.25cm radius Nickel Foil: 6cm in x, 2cm in y
A picture of the setup is shown below:
There were some strange issues while performing this simulation. The first issue that stuck out was that when trying to create a G4Orb for the selenium pellet, the program would compile, but the orb would not be present in the visualizer and furthermore would not have particles written to the outfile. A similar issue happened when trying to curve the nickel foil and the aluminum cylinder. If the nickel foil and the aluminum cylinder were curved, perhaps the number of particles passing through them would be lower which would raise the ratio of the number of particles.