Add basic tests

Author: AntonReinhard

src/compton/scattering_process.jl

@@ -72,17 +72,17 @@ Return the implementation of [`AbstractPolarization`](@ref) or [`AbstractSpin`]
 function _spin_or_pol end
 
 function _spin_or_pol(process::Compton, ::Electron, ::Incoming)
-    return process.spin[1]
+    return process.spins[1]
 end
 
 function _spin_or_pol(process::Compton, ::Electron, ::Outgoing)
-    return process.spin[2]
+    return process.spins[2]
 end
 
 function _spin_or_pol(process::Compton, ::Photon, ::Incoming)
-    return process.polarization[1]
+    return process.polarizations[1]
 end
 
 function _spin_or_pol(process::Compton, ::Photon, ::Outgoing)
-    return process.polarization[2]
+    return process.polarizations[2]
 end

test/compton.jl

@@ -2,4 +2,54 @@ using Random
 using QEDbase
 using QEDprocesses
 
+@testset "constructors" begin
+    @testset "default" begin
+        proc = Compton()
+        @test QEDprocesses._spin_or_pol(proc, Photon(), Incoming()) == AllPol()
+        @test QEDprocesses._spin_or_pol(proc, Electron(), Incoming()) == AllSpin()
+        @test QEDprocesses._spin_or_pol(proc, Photon(), Outgoing()) == AllPol()
+        @test QEDprocesses._spin_or_pol(proc, Electron(), Outgoing()) == AllSpin()
+    end
+
+    @testset "non-default" begin
+        proc = Compton((PolX(), PolY()), (SpinUp(), SpinDown()))
+        @test QEDprocesses._spin_or_pol(proc, Photon(), Incoming()) == PolX()
+        @test QEDprocesses._spin_or_pol(proc, Electron(), Incoming()) == SpinUp()
+        @test QEDprocesses._spin_or_pol(proc, Photon(), Outgoing()) == PolY()
+        @test QEDprocesses._spin_or_pol(proc, Electron(), Outgoing()) == SpinDown()
+    end
+end
+
+@testset "scattering process interface" begin
+    proc = Compton()
+
+    @test incoming_particles(proc) == (Photon(), Electron())
+    @test outgoing_particles(proc) == (Photon(), Electron())
+
+    @test number_incoming_particles(proc) == 2
+    @test number_outgoing_particles(proc) == 2
+end
+
+@testset "invalid inputs" begin
+    model = PerturbativeQED()
+    proc = Compton()
+
+    momenta_2 = [zero(SFourMomentum) for _ in 1:2]
+    momenta_3 = [zero(SFourMomentum) for _ in 1:3]
+
+    momenta_2_2 = Matrix{SFourMomentum}(undef, 2, 2)
+    momenta_3_2 = Matrix{SFourMomentum}(undef, 3, 2)
+
+    # try compute single input with incorrect dimensions
+    @test_throws DimensionMismatch differential_cross_section(proc, model, momenta_3, momenta_2)
+    @test_throws "incoming" differential_cross_section(proc, model, momenta_3, momenta_2)
+    @test_throws DimensionMismatch differential_cross_section(proc, model, momenta_2, momenta_3)
+    @test_throws "outgoing" differential_cross_section(proc, model, momenta_2, momenta_3)
+
+    # try compute multiple inputs with incorrect dimensions
+    @test_throws DimensionMismatch differential_cross_section(proc, model, momenta_3_2, momenta_2_2)
+    @test_throws "incoming" differential_cross_section(proc, model, momenta_3_2, momenta_2_2)
+    @test_throws DimensionMismatch differential_cross_section(proc, model, momenta_2_2, momenta_3_2)
+    @test_throws "outgoing" differential_cross_section(proc, model, momenta_2_2, momenta_3_2)
+end
 

test/perturbative_qed.jl

@@ -0,0 +1,7 @@
+using Random
+using QEDbase
+using QEDprocesses
+
+@testset "interaction type" begin
+    @test fundamental_interaction_type(PerturbativeQED()) == :electromagnetic
+end

test/runtests.jl

@@ -15,6 +15,9 @@ begin
     end
 
     # modules
+    @time @safetestset "perturbative qed" begin
+        include("perturbative_qed.jl")
+    end
     @time @safetestset "compton" begin
         include("compton.jl")
     end