GitHub

RemboOnDiet

RemboOnDiet.jl implements a RAMBO-on-diet style phase-space generator for relativistic multi-particle final states. It uses FourVectors.jl and LorentzVectorBase.jl for four-vector storage and operations, and it exposes a single PhaseSpaceGenerator object that can be sampled with rand(generator).

The package is organized around three ideas:

  • a fixed-size map from 3n - 4 unit-hypercube variables to an n-body event,
  • exact four-momentum conservation by construction through sequential two-body decays,
  • a phase-space weight that is constant in the massless case and becomes the corrective Jacobian for massive final states.

The implementation always uses the same main sampler. The 2-body and 3-body cases are validation anchors because their expected geometry is especially transparent, not because they use a separate generator.

Quick Start

using RemboOnDiet

generator = PhaseSpaceGenerator([0.1, 0.2, 0.3, 0.4], 5.0)
point = rand(generator)

point.momenta
phase_space_weight(point)

Main API

  • PhaseSpaceGenerator(masses, sqrt_s)
  • PhaseSpaceGenerator(masses, total::FourVector)
  • rand(generator)
  • phase_space_weight(point)
  • total_momentum(point.momenta)

Guide