diff --git a/tests/test_field_quadrupole.py b/tests/test_field_quadrupole.py
index 65e93bd..b24c9aa 100644
--- a/tests/test_field_quadrupole.py
+++ b/tests/test_field_quadrupole.py
@@ -26,7 +26,7 @@ def test_field_vs_tm_quadrupole():
         ene=gamma_avg, ene_sp=ene_spread, s_t=s_t, xi_c=xi_avg,
         q_tot=q_bunch, n_part=n_part, name='elec_bunch')
 
-    bunch_2 = copy.deepcopy(bunch_1)
+    bunch_2 = bunch_1.copy()
 
     foc_strength = 100  # T/m
     quadrupole_length = 0.05  # m
diff --git a/tests/test_parabolic_profile.py b/tests/test_parabolic_profile.py
index ebd0859..ba0a095 100644
--- a/tests/test_parabolic_profile.py
+++ b/tests/test_parabolic_profile.py
@@ -55,8 +55,8 @@ def parabolic_coefficient(z):
     # Create identical bunches for each case.
     bunch = get_matched_bunch(
         1e-6, 1e-6, 200, 1, 3, laser.xi_c - 30e-6, 1e-6, 1e4, n_p=n_p)
-    bunch_1 = deepcopy(bunch)
-    bunch_2 = deepcopy(bunch)
+    bunch_1 = bunch.copy()
+    bunch_2 = bunch.copy()
 
     # Create single plasma stage (containing all sections).
     plasma_single = PlasmaStage(
diff --git a/tutorials/01_single_plasma_simulation.py b/tutorials/01_single_plasma_simulation.py
index 2588f31..df67f71 100644
--- a/tutorials/01_single_plasma_simulation.py
+++ b/tutorials/01_single_plasma_simulation.py
@@ -68,7 +68,7 @@
 
 plasma_target = PlasmaStage(
     length=1e-2, density=1e23, wakefield_model='quasistatic_2d',
-    xi_max=30e-6, xi_min=-30e-6, r_max=30e-6, n_xi=60, n_r=30, ppc=4)
+    xi_max=30e-6, xi_min=-30e-6, r_max=30e-6, n_xi=120, n_r=60, ppc=4)
 
 # %%
 # Once the target is defined, we can track the beam through it simply by doing:
diff --git a/wake_t/__init__.py b/wake_t/__init__.py
index a8c5eb4..ef1abb7 100644
--- a/wake_t/__init__.py
+++ b/wake_t/__init__.py
@@ -1,4 +1,4 @@
-__version__ = '0.7.1'
+__version__ = '0.8.0'
 
 
 from .beamline_elements import (PlasmaStage, PlasmaRamp, ActivePlasmaLens,