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