hadronization_kernel

Function hadronization_kernel 

Source
pub fn hadronization_kernel(
    energy_density: &[EnergyDensity],
    threshold: f64,
    dim: usize,
) -> Result<Vec<PhysicalVector>, PhysicsError>
Expand description

Simulates phenomenological hadronization: $E_{\text{density}} \rightarrow \vec{p}_{\text{jet}}$.

This kernel implements a simplified string fragmentation model where local energy density fluctuations exceeding a critical threshold $E_c$ condense into particle momentum vectors (jets).

Under the ‘Local Parton-Hadron Duality’ (LPHD) hypothesis, the energy density maps directly to the momentum magnitude of the resulting jet, assuming a massless limit ($E \approx |p|c$). The resulting momentum is assigned to the primary spatial component ($x$-axis/index 1) of the CausalMultiVector, representing a collimated jet in the local coordinate frame.

§Arguments

  • energy_density - A slice of scalar energy density fields (0-forms).
  • threshold - The critical deconfinement energy density $E_c$.
  • dim - The spatial dimension $d$ of the resulting manifold (requires $d > 0$).

§Returns

  • Result<Vec<PhysicalVector>, PhysicsError> - A collection of generated momentum vectors.

§Errors

  • DimensionMismatch - If dim == 0.
  • PhysicalInvariantBroken - If threshold < 0.
  • NumericalInstability - If vector creation fails.