Difference between revisions of "Occupancy for Sector 1"

From New IAC Wiki
Jump to navigation Jump to search
Line 75: Line 75:
  
 
[[File:IntegralDCLimitsdSigmaCM.png]][[File:IntegralDCLimitsdSigmaLab.png]]
 
[[File:IntegralDCLimitsdSigmaCM.png]][[File:IntegralDCLimitsdSigmaLab.png]]
 +
 +
=GEMC Cross Section=
  
 
=Using the Cross Section=
 
=Using the Cross Section=

Revision as of 16:16, 17 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


Performing a Riemann sum for [math]-30^{\circ} \lt \phi \lt 30^{\circ}[/math]


CrossSection60deg.png


MolThetaCMdsigmaIntegral.pngMolThetaLabdSigmaIntegral.png

AssociatedWeights2.pngDSigmaCMLab.png

Adjust for DC Sector 1 Limits

IntegralDCLimitsdSigmaCM.pngIntegralDCLimitsdSigmaLab.png

GEMC 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 

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 

Clas12monNoSolNoShield.png

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]