Difference between revisions of "Occupancy for Sector 1"

From New IAC Wiki
Jump to navigation Jump to search
Line 74: Line 74:
  
 
=Occupancy=
 
=Occupancy=
 +
LH2_NOSol_0Tor_11GeV_IsotropicPhi_v2_6_ShieldOut
 +
 +
Run
 +
<pre>
 +
./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.....
 +
 +
</pre>
 +
 +
Create hipo file
 +
 +
 +
Move hipo file to clas12mon folder
 +
<pre>
 +
mv LH2_NOSol_0Tor_11GeV_IsotropicPhi_v2_6_ShieldOut.hipo ~/clas12mon
 +
</pre>
 +
 +
Run monitor program
 +
<pre>
 +
./README
 +
</pre>
 +
 +
Load hipo file
 +
<pre>
 +
"Press H for hipo"
 +
"Press play"
 +
"Switch to
 +
</pre>
 +
 
[[File:clas12monNoSolNoShield.png]]
 
[[File:clas12monNoSolNoShield.png]]

Revision as of 17:57, 12 April 2018

[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.pngMolThetaCM LUND DC limits.png

Applying the weight

MolThetaLab DClimits integral.pngCosMolThetaLab weightedDClimits.png


MolThetaCM DClimits weighted rebin integral.pngCosMolThetaCM weightedDClimits.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

CrossSectionMathematica1.png


"Integrating" with Cosine term

CrossSectionMathematica2.png

Finding the Cross Section

CrossSectionMathematicaProof.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
+------------------------------------------------------+
|                                                      |
|  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.....

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 

Clas12monNoSolNoShield.png