Difference between revisions of "July2012PosSimulation"
| Line 136: | Line 136: | ||
The reason for the sharp drop in the positron distribution: there is very small probability for positrons to travel in the radial direction of the tungsten target and make it to the detector. | The reason for the sharp drop in the positron distribution: there is very small probability for positrons to travel in the radial direction of the tungsten target and make it to the detector. | ||
| − | [[File:hrrl_pos_sim_pos_sharp_drop_geometry.png | 400 px]] | + | [[File:hrrl_pos_sim_pos_sharp_drop_geometry.png | 400 px]] [[File:hrrl_pos_sim_s2_beam_DDNT1_xy_zoom_s1.png | 400 px]] |
| − | |||
[[File:hrrl_pos_sim_s2_beam_DDNT1_y_s1.png | 400 px]] [[File:hrrl_pos_sim_s2_beam_DDNT1_xy_s1.png | 400 px]] | [[File:hrrl_pos_sim_s2_beam_DDNT1_y_s1.png | 400 px]] [[File:hrrl_pos_sim_s2_beam_DDNT1_xy_s1.png | 400 px]] | ||
Revision as of 05:51, 17 April 2013
2 MeV Positrons
Measured Electron energy distribution at 10 MeV on 1.mm target
Simulation steps
1.) GEANT4 Simulated Beam energy distribution
2.) GEANT4 Simulated Positrons emitted from 1.x mm thick target
3.) Now use above Positron distribution as the particle source for G4beamline
Positrons hitting Tungsten Converter target
4.) Use GEANT4 to determine 511s from the positrons distribution impinging converter target
Oct 16th 2012 BenchMark
stages:
1. Program to input emittance, output beamsize, beam, divergence and beam energy.
2. e- on W, outcome e+. Incident electron distribution on T1. with general particle source with step 1 histogram. check graph x, y, z theta, En_dis out_put is the input. (directory for each).
3. insert T1 and generate positron distribution.
4. Acceleractor code to transport positron along the beamline. Out put positron theta, beamsize (x,y), energy distribution, out puts transported positron.
5. Geant4 takes step 4 output generates gamma (and other e+,e-) and look at those goes to detectors. beamline, plus shielding, and detectors.
Simulations
Partical Data Ground ID: PDGid=11 is electron. PDGid=-11 is positron. PDGid=22 is photon.
Parameters
Dipole vacuum chamber width is mm
The cavity exit diameter is about 7.3 mm.
Energy Spread Two Skewed Gaussian Fit
Energy spread fitted with two skewed Gaussian.
Beam Distributions Beyond RMS: Media:Beam_Distributions_Beyond_RMS.pdf
Amplitude = 2.13894, mean = 12.07181, sigma_L = 4.46986, sigma_R = 1.20046
Sigma = 2.83516, Es = -0.57658
rms = 2.56472, skew = -0.94853
Amplitude2 = 10.88318, mean2 = 12.32332, sigma_L2 = 0.69709, sigma_R2 = 0.45170
Sigma2 = 0.57440, Es2 = -0.21360
rms2 = 0.56719, skew2 = -0.34231
Electron Distribution Upstream and Positron Distribution Downstream of Target1
simulation: Electron Distribution Upstream and Positron Distribution Downstream of Target1
3067274 Electrons fired.
T1_s=943.5
FDT1_s=$T1_s-26.52
BDT1_s=$T1_s+26.52
Thickness of T1 is 1.016 mm, radius of T1 is 15.875 mm (0.625 inches).
Step 1: Generating positron from electron beam
Detector is after T1 to saple positron distribution. 1,379,974,500 electrons shot at the tungsten target to generate positrons.
step 1: ratio
DDNT1 = Detector DowN from T1
DQ4 = Detector at entrance Quad 4 located right after T1
DD1 = Detector at entrance of the first Dipole.
e-(DDT1)/e+(DDNT1) = 1248
e-(DDT1)/e+(DQ4) = 123698
e-(DDT1)/e+(DD1) = 936846
e+(DDNT1)/e+(DQ4) = 99
e+(DDNT1)/e+(DD1) = 750
e+(DQ4)/e+(DD1) = 7.6
Step 1: Generated Beam on DDNT1
Detector down T1:
Positrons on detector right before Q4.
X profile 
Y profile (why is the distribution clipped on the left side in the figure) 
X divergence profile 
Y divergence profile 
Momentum distribution 
Momentum distribution 
Sharp drop issue
The reason for the sharp drop in the positron distribution: there is very small probability for positrons to travel in the radial direction of the tungsten target and make it to the detector.
DDNT1 is at 26.52 mm (53.04 mm)downstream of T1.
T1 radius is 25.4 mm (d=50.8 mm)
DDNT1 radius is 31.75 mm (d=63.5 mm).
Step 1: Generated Beam on DQ4
Positrons on detector right before Q4.
X profile 
Y profile 
X divergence profile 
Y divergence profile 
Momentum distribution 
Momentum distribution 
Step 1: Generated Beam on DD1
Positrons on detector right before D1 - first dipole.
X profile 
Y profile 
X divergence profile 
Y divergence profile 
Momentum distribution 
Momentum distribution 
Step 2: generate positron
The positrons after T1 is detected on a virtual detector. The positrons beam size, divergence and momentum distributions are extracted and created new positron beam.
Generate positron beam from the detector after T1 and transport it to right before D1.
221032400*16=3,536,518,400 positrons generated.
step 2: ratio
e+(DDNT1)/e+(DQ4) = 3536518400/47788670 = 74
e+(DDNT1)/e+(DD1) = 3536518400/5464220 = 647
e+(DQ4)/e+(DD1) = 8.7
Generated Beam on DQ4
Positrons on detector right before Q4.
X profile 
Y profile 
X divergence profile 
Y divergence profile 
Momentum distribution 
Momentum distribution 
Generated Beam on DD1
Positrons on detector right before D1 - first dipole.
X profile 
Y profile 
X divergence profile 
Y divergence profile 
Momentum distribution 
Momentum distribution 
Beam Compare Step1 and Step2
Compare Ratio
| S1 | S2 | |
| e-(DDT1)/e+(DDNT1) | 1248 | |
| e-(DDT1)/e+(DQ4) | 123698 | |
| e-(DDT1)/e+(DD1) | 936846 | |
| e+(DDNT1)/e+(DQ4) | 99 | 74 |
| e+(DDNT1)/e+(DD1) | 750 | 647 |
| e+(DQ4)/e+(DD1) | 7.6 | 8.7 |
Positrons at DDNT1
| Step 1 | Step 2 | |
| X profile | |
|
| Y profile | |
|
| X divergence profile | |
|
| Y divergence profile | |
|
| Momentum distribution | |
|
| Momentum distribution | |
Compare at DQ4
| Step 1 | Step 2 | |
| X profile | |
|
| Y profile | |
|
| X divergence profile | |
|
| Y divergence profile | |
|
| Momentum distribution | |
|
| Momentum distribution | |
|
Compare at DD1
| Step 1 | Step 2 | |
| X profile | |
|
| Y profile | |
|
| X divergence profile | |
|
| Y divergence profile | |
|
| Momentum distribution | |
|
| Momentum distribution | |
|