Defining Occupancy

General Occupancy

The occupancy measures the number of particles interactions per a detector cell per an event. For the CLAS12 drift chamber, there are 112 wires on each layer, with 12 layers within a region, giving 1344 cells. This can simply be defined as the "Unweighted Occupancy" for the CLAS12 DC and follows the equation:

Unweighted CLAS12 DC Occupancy[math]\equiv \frac{N_{hits}}{N_{evts}N_{cells}}[/math]

where

[math]N_{hits}\equiv [/math]The number of DC wires intersected by primary and secondary events throughout the drift chamber in Region 1

[math]N_{evt} \equiv \phi \times Prob(interacting)[/math]

[math]\phi \equiv [/math] Number of incident particles on the face of drift chamber per cm[math]^2[/math]

[math]N_{cells} \equiv 112 \frac{wires}{layer} \times 12 \frac{layers}{Region}[/math]

CLAS12 DC Occupancy

The registering of a "hit" takes a finite time in which the detector and its associated electronics are not able to register an additional signal if it occurs. This time window is known as the "dead time" during which only limited events are registered. For Region 1:

[math]\Delta t \equiv [/math] 250 ns: The time needed for events to be read by the electronics within Region 1

Since the events are simulated outside the dead time constraints of the DC, we can factor in the number of event windows that occur by dividing the dead time window per region by the time that would have been required to produce the number of incident electrons given a known current.

[math]t_{sim} \equiv [/math] Time of simulation = [math]\frac{N_{incident}}{I(A)}\frac{1A}{1C}\frac{1s}{}\frac{1.602E-19\ C}{1\ e^{-}}[/math]

When applying the Moller differential cross-section as a weight, this gives the CLAS12 occupancy as:

CLAS12 DC occupancy [math]\equiv \frac{N_{hits}}{N_{evt}}\frac{\Delta t}{t_{sim}}\frac{1}{112}\frac{1}{12}[/math]

Weighted Hits Occupancy

Lab Cross-Section

Using the definition of the cross-section:

[math]\sigma \equiv \frac{N_{scattered}}{\mathcal L t}=\frac{N_{scattered}}{\Phi \rho l t}[/math]

where the flux is defined as:

[math]\Phi \equiv \frac{Number\ of\ e^{-}}{s}[/math]

Making some assumptions that the flux can be taken over an same time range as the time found in the cross-section, which allows

[math]\Phi \approx N_{incident}[/math]

For a LH2 target of length 5cm.

[math]\rho_{target}\times l_{target}=\frac{70.85 kg}{1 m^3}\times \frac{1 mole}{2.02 g} \times \frac{1000g}{1 kg} \times \frac{6.022\times10^{23} molecules\ LH_2}{1 mole} \frac{2 atoms}{1\ molecule\ LH_2} \times \frac{1m^3}{(100 cm)^3} \times \frac{5 cm}{ }=\frac{2.11\times 10^{23}}{cm^2} \times \frac{1^{-24} cm^{2}}{barn}=0.211 barns^{-1}[/math]

Additionally, for Moller Scattering, we can assume that almost 100% of the scattered electrons occur as events.

[math]N_{scattered}=N_{evt}[/math]

This allows us to rewrite the cross-section expression as,

[math]\sigma \frac{N_{scattered}}{\Phi \rho l}\Rightarrow N_{evt}=\sigma N_{incident} \rho l[/math]

we can define a differential scattering cross-section, [math]\frac{d\sigma(\theta')}{d{\mit\Omega}'}[/math], in the laboratory frame, where [math]d{\mit\Omega}'=\sin\theta'\,d\theta'd\phi[/math] is an element of solid angle in this frame. Therefore, [math]\frac{d\sigma(\theta)}{d{\mit\Omega}'}\,d{\mit\Omega}'[/math] is the effective cross-sectional area in the laboratory frame for scattering into the range of scattering angles [math]\theta[/math] to [math]\theta+d\theta[/math] and transverse angles [math]\phi[/math] to [math]\phi+d\phi[/math]. Likewise, [math]\frac{d\sigma(\theta)}{d{\mit\Omega}}\,d{\mit\Omega}[/math]is the effective cross-sectional area in the CM frame for scattering into the range of scattering angles [math]\theta [/math] to [math]\theta+d\theta[/math] and transverse angles [math]\phi[/math] to [math]\phi+d\phi[/math]. By relativity, only the scattering angles are Lorentz contracted in the direction of motion, leaving the transverse angles invariant between frames of reference. Additionally, an effective cross-section with corresponding cross-sectional area is not changed when we transform between different inertial frames. The number of scattered particles and incident particles are invariant between frames, while the cross-sectional areas are also measuring the same relative space.

[math]\frac{d\sigma(\theta')}{d{\mit\Omega}'} \,d{\mit\Omega}' = \frac{d\sigma(\theta)}{d{\mit\Omega}}\,d{\mit\Omega}[/math]

[math]d\sigma(\theta') =\frac{d\sigma(\theta)}{d{\mit\Omega}}\, \frac{d{\mit\Omega}}{d{\mit\Omega}'}[/math]

[math]d\sigma(\theta') =\frac{d\sigma(\theta)}{d{\mit\Omega}}\, \sin\theta d\theta d\phi[/math]

Since the expression for the differential cross-section for Moller Scattering is well known in the CM, we can solve for the minimum angle detected by the DC (.55 degrees in Theta) in the lab frame.

[math]\sigma(\theta'=0.55^{\circ})=\frac{d\sigma(\theta)}{d{\mit\Omega}}\, \sin\theta\, d\theta\, d\phi[/math]

[math]\int_{\theta}^{\theta+d\theta} \int_{0}^{\d\phi} \frac{\left(3+\left(\cos{\theta}\right)^2\right)^2}{\left(\sin{\theta}\right)^4}\sin{\theta} d\theta[/math]

Here we will use a indefinite integral to find the cross-section in the lab. For a constant [math]\Delta \theta_{lab}=0.1^{\circ}[/math] this corresponds to a [math]\Delta \theta_{CM}[/math] that will become smaller as [math]\theta_{lab}[/math] increases. [math]\Delta \phi=0.2^{/circ}[/math] in both frames due to it being Lorentz invariant. Using the Lorentz boost vectors in the LUND file generator we find:

Labθ=5.0	Labdσ=0.000000020919615	ΔLabθ=0.1        CMθ=167.426100001950005	CMdσ=0.000000020919615	ΔCMθ=0.00429

Number of Incident Electrons

Solving for the number of incident electrons that would be required to produce 1 scattered electron (event),

[math]0.000000020919615=\frac{1\,e^{-}}{N_{in}\rho l}\rightarrow N_{in}=226,549,890.755\,e^{-}[/math] for a 5cm long LH2 target.

Near the maximum of the DC, this corresponds to more events.

Labθ=39.8	Labdσ=0.000000341148416	ΔLabθ=0.1        CMθ=178.674400001996105	CMdσ=0.000000341148416	ΔCMθ=0.00008

[math]0.000000341148416=\frac{N_{out}}{N_{in}\rho l}\rightarrow N_{out}=16.3\,e^{-}[/math]