[math]\textbf{Navigation}[/math]
[math]\vartriangleleft [/math]
[math]\triangle [/math]
[math]\vartriangleright [/math]
A bash script to run the GEMC simulations is created. tcsh scripts to run root2evio on lds2 is called using sshpass. The lds2 scripts use sshfs
The main script on lds3:
BUILD_GEMC_SIMULATION.sh
The 3 scripts on lds2:
first_commands.tcsh
second_commands.tcsh
last_commands.tcsh
LUND File Output
0.1 degree spacing in the Lab frame. CM Frame is not evenly spaced.
![MolThetaLab LUND DC limits.png](/./images/f/f1/MolThetaLab_LUND_DC_limits.png)
Applying the weight
![MolThetaLab DClimits integral.png](/./images/f/f3/MolThetaLab_DClimits_integral.png)
![MolThetaCM DClimits weighted rebin integral.png](/./images/a/a5/MolThetaCM_DClimits_weighted_rebin_integral.png)
Looking at the angles and the associated weight, we can find the sums
Once_Angles_and_weight=3399.930890560805437
Total_Angles_and_weight=1023379.198058736044914
Checking with Mathematica
"Integrating" with Cosine term
Finding the Cross Section
Performing a Riemann sum for [math]-30^{\circ} \lt \phi \lt 30^{\circ}[/math]
![MolThetaCMdsigmaIntegral.png](/./images/3/33/MolThetaCMdsigmaIntegral.png)
![AssociatedWeights2.png](/./images/b/b3/AssociatedWeights2.png)
Using the Cross Section
[math]\sigma \equiv \frac{N_{scattered}}{\mathcal{L}t}[/math]
[math]\sigma \mathcal{L} =\frac{N_{scattered}}{t}[/math]
For a Luminosity of [math]\mathcal{L}=\frac{1.32\times 10^{11}}{barn\cdot s}[/math]
[math].087\ barn \frac{1.32\times 10^{11}}{barn\cdot s} =\frac{N_{scattered}}{t}=\frac{1.15\times 10^{10}}{t}[/math]
Occupancy
LH2_NOSol_0Tor_11GeV_IsotropicPhi_v2_6_ShieldOut
Run
./BUILD_GEMC_SIMULATION.sh
DVMacro
Clas12Mon
Create hipo file
Move hipo file to clas12mon folder
mv LH2_NOSol_0Tor_11GeV_IsotropicPhi_v2_6_ShieldOut.hipo ~/clas12mon
Run monitor program
./README
Load hipo file
"Press H for hipo"
"Press play"
"Switch to
Calculating
[math]\sigma_{Moller\ (\theta_{lab} = 5^{\circ}-40^{\circ})}=0.86\ barn[/math]
[math]t_{simulated}=\frac{N_{events}}{\sigma_{events} \Phi \rho \ell}=\frac{96105\ barn \cdot s}{7.87\times 10^{-2} \cdot 1.33\times 10^{11}\ barn}=9.3\times 10^{-6} s[/math]
[math]N_0=\Delta t \cdot R_{events}=\Delta t \cdot \frac{N_{events}}{t_{simulated}}=250\times 10^{-9}\ s \cdot \frac{98181}{9.3\times 10^{-6}\ s}=2639[/math]
[math]Occupancy=\frac{N_{hits}}{N_0}=\frac{N_{hits}}{\Delta t \cdot R_{events}}=\frac{t_{simulated}\cdot N_{hits}}{N_{events}\cdot \Delta t}=[/math]