diff --git a/data/see/constant.dat b/data/see/constant.dat
new file mode 100644
index 0000000..060c8ad
--- /dev/null
+++ b/data/see/constant.dat
@@ -0,0 +1,4 @@
+#Relative energy (eV) yield ()
+0.000 1.000E-01
+1.000 1.000E-01
+
diff --git a/src/modules/mesh/moduleMeshBoundary.f90 b/src/modules/mesh/moduleMeshBoundary.f90
index 34dc8a7..ffd1af9 100644
--- a/src/modules/mesh/moduleMeshBoundary.f90
+++ b/src/modules/mesh/moduleMeshBoundary.f90
@@ -193,7 +193,7 @@ MODULE moduleMeshBoundary
 
       END SELECT
 
-      !Removes ionizing electron regardless the number of pair created
+      !Removes ionizing electron regardless the number of pairs created
       part%n_in = .FALSE.
 
     END SUBROUTINE ionization
@@ -214,6 +214,7 @@ MODULE moduleMeshBoundary
 
     SUBROUTINE secondaryEmission(edge, part)
       USE moduleSpecies
+      USE moduleRandom
       IMPLICIT NONE
 
       CLASS(meshEdge), INTENT(inout):: edge
@@ -270,6 +271,7 @@ MODULE moduleMeshBoundary
           secondaryElectron%n_in = .TRUE.
 
           !Assign velocity
+          secondaryElectron%v = 2.D0 * edge%normal + 1.D-1 * (/ randomMaxwellian(), randomMaxwellian(), randomMaxwellian() /)
 
           !Add particle to list
           CALL partSurfaces%setLock()
@@ -314,6 +316,7 @@ MODULE moduleMeshBoundary
 
       TYPE IS(boundarySEE)
         edge%fBoundary(s)%apply => secondaryEmission
+
       CLASS DEFAULT
         CALL criticalError("Boundary type not defined in this geometry", 'pointBoundaryFunction')