From b784720859c969196fbbef5dff4b7cd03b305c55 Mon Sep 17 00:00:00 2001 From: Anton Reinhard Date: Thu, 4 Jul 2024 20:40:22 +0200 Subject: [PATCH] Minor fixes and commints --- examples/congruent_in_ph.jl | 10 ++-------- 1 file changed, 2 insertions(+), 8 deletions(-) diff --git a/examples/congruent_in_ph.jl b/examples/congruent_in_ph.jl index 1f8367d..24b430e 100644 --- a/examples/congruent_in_ph.jl +++ b/examples/congruent_in_ph.jl @@ -26,15 +26,12 @@ end function congruent_input_momenta(processDescription::GenericQEDProcess, omega::Number) # generate an input sample for given e + nk -> e' + k' process, where the nk are equal - massSum = 0 inputMasses = Vector{Float64}() for particle in incoming_particles(processDescription) - massSum += mass(particle) push!(inputMasses, mass(particle)) end outputMasses = Vector{Float64}() for particle in outgoing_particles(processDescription) - massSum += mass(particle) push!(outputMasses, mass(particle)) end @@ -49,7 +46,7 @@ function congruent_input_momenta(processDescription::GenericQEDProcess, omega::N return (tuple(initial_momenta...), tuple(final_momenta...)) end -# cos_theta ∈ [-1, 1] and phi ∈ [0, 2π] +# theta ∈ [0, 2π] and phi ∈ [0, 2π] function congruent_input_momenta_scenario_2( processDescription::GenericQEDProcess, omega::Number, @@ -60,15 +57,12 @@ function congruent_input_momenta_scenario_2( # same as above # generate an input sample for given e + nk -> e' + k' process, where the nk are equal - massSum = 0 inputMasses = Vector{Float64}() for particle in incoming_particles(processDescription) - massSum += mass(particle) push!(inputMasses, mass(particle)) end outputMasses = Vector{Float64}() for particle in outgoing_particles(processDescription) - massSum += mass(particle) push!(outputMasses, mass(particle)) end @@ -163,7 +157,7 @@ for photons in [6] end end - +# scenario n = 1000000 # n is the number of incoming photons