Difference between revisions of "Lead Shield Cone"
(→Target) |
|||
| Line 95: | Line 95: | ||
NOTE: This will need to be rotated 180 degrees so that the large end is nearest to the particle vertex. | NOTE: This will need to be rotated 180 degrees so that the large end is nearest to the particle vertex. | ||
| + | |||
| + | [[ExN02DetectorConstruction.cc]] | ||
| + | |||
| + | |||
| + | [[ExN02SteppingVerbose.cc]] | ||
=Particle Gun= | =Particle Gun= | ||
Revision as of 03:03, 23 January 2018
Target
We need to edit the EXn02DetectorConstruction file to allow for a target that is a lead cone to simulate the Moller Shield. The header file for an elliptical Cone gives the parameters need to call with.
// $Id: G4EllipticalCone.hh 67011 2013-01-29 16:17:41Z gcosmo $ // // // -------------------------------------------------------------------- // GEANT 4 class header file // // G4EllipticalCone // // Class description: // // G4EllipticalCone is a full cone with elliptical base which can be cut in Z. // // Member Data: // // xSemiAxis semi-axis, x, without dimentions // ySemiAxis semi-axis, y, without dimentions // zheight height, z // zTopCut upper cut plane level, z // // The height in Z corresponds to where the elliptical cone hits the // Z-axis if it had no Z cut. Also the cone is centered at zero having a // base at zTopCut and another at -zTopCut. The semi-major axes at the Z=0 // plane are given by xSemiAxis*zheight and ySemiAxis*zheight so that the // curved surface of our cone satisfies the equation: // // *************************************************************************** // * * // * (x/xSemiAxis)^2 + (y/ySemiAxis)^2 = (zheight - z)^2 * // * * // *************************************************************************** // // In case you want to construct G4EllipticalCone from : // 1. halflength in Z = zTopCut // 2. Dx and Dy = halflength of ellipse axis at z = -zTopCut // 3. dx and dy = halflength of ellipse axis at z = zTopCut // ! Attention : dx/dy=Dx/Dy // // You need to find xSemiAxis,ySemiAxis and zheight: // // xSemiAxis = (Dx-dx)/(2*zTopCut) // ySemiAxis = (Dy-dy)/(2*zTopCut) // zheight = (Dx+dx)/(2*xSemiAxis)
The geometry looks like
Solving the variables:
1. zTopCut
The cone shape starts 380 mm from the vertex point. It extends 1325.9 beyond the starting position. This gives the full length from where it would cross the z axis as 1705.9 mm. Since this cone is centered at the origin with zTopCut positive and negative placed an equal distance from the vertex, 1325.9/2=662.95 mm.
2. Dx and Dy
The length of the x and y components at the larger end of the cone are found from the radius at this position. Using geometry, for a right triangle with it's apex at the vertex and a height of 1705.9, with an interior angle of 5 degrees.
3. dx and dy
Using geometry, for a right triangle with it's apex at the vertex and a height of 380, with an interior angle of 5 degrees.
Solving for the variable to pass to the function:
The height in z is defined at the the point where the cone vertex crosses the z=0 plane in x&y. From the drawing, the vertex is originally placed at 0,0,0. Since the GEANT4 cone is centered at a new vertex with the circular cross sections of the cone at +/- 662.95, and the distance from the first smaller cross section is 380, this gives the intersection at 1042.95 cm with respect to the new vertex. Verifying this:
NOTE: This will need to be rotated 180 degrees so that the large end is nearest to the particle vertex.
Particle Gun
mkdir ParticleGunInputs
root -l LUND_Spread_LH2_IsotropicPhi.C
split -d -l 3000 -a 3 LUND_Spread_LH2_IsotropicPhi.LUND LUND_Spread_LH2_IsotropicPhi_
prename 's/(LUND_Spread_LH2_IsotropicPhi_\d{3})/$1.LUND/' LUND_Spread_LH2_IsotropicPhi_*
Using the file LUND2ParticleGun.C 106 macros for the GEANT4 ParticleGun are created that mirror the LUND files. The same random number is initially input for each macro as was used for the LUND file.
Running Simulations
==
Event # 902
Step# X Y Z Px Py Pz KineE dEStep StepLeng TrakLeng Volume Process Theta(pos) theta(mom) Phi(pos)
0 0 fm 0 fm -2 m 9.61695 keV 5.55235 keV 122.019 keV 14.48357 keV 0 eV 0 fm 0 fm World initStep NA NA NA
1 10.09023 cm 5.825597 cm -71.97609 cm 9.61695 keV 5.55235 keV 122.019 keV 14.48357 keV 4.258733e-16 eV 1.28553 m 1.28553 m World Transportation 5.2 5.2 30
2 10.09023 cm 5.825597 cm -71.97609 cm -31.23412 keV 91.41135 keV 72.1564 keV 14.03255 keV 451.0172 eV 49.21704 nm 1.28553 m Target msc 5.2 53.24 108.9
3 10.09023 cm 5.825598 cm -71.97609 cm -59.67862 keV 79.92509 keV 67.73916 keV 14.03255 keV 0 eV 2.216658 nm 1.28553 m Target Transportation 5.2 55.82 126.7
4 10.09023 cm 5.825598 cm -71.97609 cm -59.67862 keV 79.92509 keV 67.73916 keV 14.03255 keV 0 eV 0 fm 1.28553 m World Transportation 5.2 55.82 126.7
5 10.09022 cm 5.825604 cm -71.97608 cm -0 eV 0 eV 0 eV 0 eV 14.03255 keV 1.171968 um 1.285531 m Target eIoni 5.2 -nan -nan
Event # 903
Step# X Y Z Px Py Pz KineE dEStep StepLeng TrakLeng Volume Process Theta(pos) theta(mom) Phi(pos)
0 0 fm 0 fm -2 m 9.43847 keV -5.4493 keV 117.483 keV 13.44453 keV 0 eV 0 fm 0 fm World initStep NA NA NA
1 15.99911 cm -9.237085 cm -8.55066 mm 9.43847 keV -5.4493 keV 117.483 keV 13.44453 keV 7.06826e-16 eV 2 m 2 m World Transportation 5.3 5.3 -30
2 16.00511 cm -9.240549 cm -7.803867 mm 9.43847 keV -5.4493 keV 117.483 keV 13.44453 keV 2.650598e-19 eV 750 um 2.00075 m Tracker StepLimiter 5.3 5.3 -30
The Moller particles interact with the target up to 5.3 degrees in Theta.