Lead Shield Cone

From New IAC Wiki
Jump to navigation Jump to search

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

HTCC section with ic v3.png

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.


[math]1705.9\ mm\ Tan\ 5^\circ=149.247\ mm[/math]



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.


[math]380\ mm\ Tan\ 5^\circ=33.246\ mm[/math]



Solving for the variable to pass to the function:


[math] xSemiAxis = \frac{Dx-dx}{2*zTopCut}= \frac{149.25-33.25}{2*662.95}= \frac{116}{1325.9}=.087=ySemiAxis [/math]


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:


[math]zheight = \frac{Dx+dx}{2*xSemiAxis}=\frac{149.247+33.25}{2*0.087}=\frac{182.5}{.175}=1042.95[/math]


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

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

GEANT4submit_builder.C

Moller Shield effects

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.

Effects of DC/Tracker on Moller particles

The DC chamber is filled with a mixture of 90%/10% of Ar/CO2. The density of this gas can be calculated can be calculated using the ideal gas law, expressed as a function of temperature and pressure:


[math]\rho = \frac{p}{R_{\rm specific} T} [/math]

where:

[math]\rho = density\ (g/cm3)[/math]


[math]p = absolute\ pressure\ (Pa)[/math]


[math]T = absolute\ temperature\ (K)[/math]


[math]R_{\rm specific} = specific\ gas\ constant\ (J/(g·K))[/math]


The specific gas constant of a gas or a mixture of gases ([math]R_{specific}[/math]) is given by the molar gas constant divided by the molar mass (M) of the gas or mixture, or in turn by Boltzmann constant over the molecular mass of the gas.


[math]R_{\rm specific}=\frac {R}{M}=\frac {k_{\rm {B}}}{m}[/math]


The average molar mass of mixtures [math] \bar {M} [/math] can be calculated from the mole fractions [math] x_{i} [/math] of the components and their molar masses [math] M_{i}[/math]:

[math] \bar {M}=\sum _{i}x_{i}M_{i}[/math]

M(Ar) = 39.948 × 1 g/mol = 39.948 g/mol M(CO2) = 12.01+16.00*2× 1 g/mol = 44.01 g/mol


[math] \bar {M}=0.9*39.948 g/mol+0.1*44.01 g/mol=35.9532+4.401=40.354 g/mol[/math]


[math]R_{\rm specific}=\frac {R}{M}[/math]

where:

R=8.314 J / mol. K.


[math]R_{\rm specific}=\frac {8.314 J / mol. K.}{40.354 g/mol}=0.2060 \frac{J}{g \cdot K}[/math]


[math]\rho = \frac{p}{R_{\rm specific} T} [/math]

Using the normal temperature and pressure of 20 °C (293.15 K, 68 °F) and an absolute pressure of 1 atm (14.696 psi, 101.325 kPa) 1Pa=1 kg⋅m−1⋅s−2 1 J=kg⋅m2⋅s−2


[math]\rho = \frac{101325 kg⋅m^{−1}⋅s^{−2}}{0.2060 \frac{J}{g \cdot K} 293.15 K} = \frac{101325kg⋅m^{−1}⋅s^{−2}}{60.40 \frac{kg⋅m^2⋅s^{−2}}{g}} [/math]


[math]\rho =\frac{1678g}{m^3}\frac{1m^3}{1000000cm^3}\frac{1000mg}{1g}=1.678mg/cm^3[/math]