Occupancy for Sector 1

[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.

Applying the weight

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

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
+------------------------------------------------------+
|                                                      |
|  GEMC Simulations with root  V1.0                    |
|  Daniel VanWasshenova                                |
|                                                      |
+------------------------------------------------------+
********************************************************
Target material (LH2, NH3, NO) ?
LH2
Solenoid Field (-1 to 1, NO) ?
NO
Torus Field (-1 to 1) ?
-1
Incident Electron Energy (GeV) ?
11
Phi Distribution (Isotropic or set value) ?
Isotropic
GEMC Version (4 or 6) ?
6
Shield Status (In/Out) ?
Out
Creating directory LH2_NOSol_-1Tor_11GeV_IsotropicPhi_v2_6_ShieldOut
*********************************************************
Creating LUND files......
root [0] 
Processing LUND_Spread_IsotropicPhi.C...
******************************************************************************
*Tree    :Moller    : Moller data from ascii file                            *
*Entries :        0 : Total =             289 bytes  File  Size =          0 *
*        :          : Tree compression factor =   1.00                       *
******************************************************************************
*********************************************************
Building specific gcards.....
root [0] 
Processing gcard_builder.C("LH2","NO","-1","11","Isotropic","6","Out","x","x","x","x")...
Creating submit files........
root [0] 
Processing submit_builder.C("LH2_NOSol_-1Tor_11GeV_IsotropicPhi_v2_6_ShieldOut")...
root [0] 
Processing Block_builder.C...
Changed directory to:/home/lds/src/CLAS/GEMC/experiments/eg12/MolrBckGrd/DV/Isotropic_study/LH2_NOSol_-1Tor_11GeV_IsotropicPhi_v2_6_ShieldOut/submit

*********************************************************
CREATING SCREEN: Name
screen: zeroth
screen: first
screen: second
screen: third
screen: fourth
screen: fifth
screen: sixth
screen: seventh
screen: eighth
screen: ninth
screen: tenth
Creating evio files!
[0 %] finished....[0 %] finished.....[0 %] finished.....[0 %] finished.....[0 %] finished.....[0 %] finished.....[0 %] finished.....[0 %] finished.....[10.0 %] finished.....[10.0 %] finished.....[10.0 %] finished.....[10.0 %] finished.....[20.0 %] finished.....[20.0 %] finished.....[20.0 %] finished.....[20.0 %] finished.....[30.0 %] finished.....[30.0 %] finished.....[30.0 %] finished.....[30.0 %] finished.....[30.0 %] finished.....[30.0 %] finished.....[40.0 %] finished.....[40.0 %] finished.....[40.0 %] finished.....[40.0 %] finished.....[50.0 %] finished.....[50.0 %] finished.....[50.0 %] finished.....[50.0 %] finished.....[50.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[60.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[70.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....[80.0 %] finished.....

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]t_{simulated}=\frac{N_{events}}{\sigma_{events} \Phi \rho \ell}=\frac{1000000\ barn \cdot s}{7.87\times 10^{-2} \cdot 1.33\times 10^{11}\ barn}=9.54\times 10^{-5} s[/math]

[math]N_0=\Delta t \cdot R_{events}=\Delta t \cdot \frac{N_{events}}{t_{simulated}}=500\times 10^{-9}\ s \cdot \frac{1\times 10^6}{0.0954\times 10^{-3}\ s}=5240[/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]