import happi
import math
import h5py
import pint
import matplotlib
import pylab
import vtk
S = happi.Open("/home/whharrison/smilei-master-ed2612321e2266a41c3d247518a2877e2c5961d4/thermalplasma_toVTK") # Open a simulation
S.TrackParticles("proton", axes=["x","y","z","px","py","pz"]).getData()
tracked_particles = S.TrackParticles("proton", axes=["x","y","z","px","py","pz"],timesteps=[1,10])
tracked_particles.toVTK(rendering="cloud",data_format="xml");
print("End!")