species in an initial H, O, C gas mixture, a 0D collisionalradiative model (CRM) was built and experimentally validated [1]. This initial CRM involves 1050 reactions (electrons impacting atoms in their fundamental, excited and ionized states, atom de-excitation, three body recombination, ionization stepwise, etc) and 79 species (electrons, atoms, ions and excited states) [1]. On the other hand, the electron beam can be simulated based on the LSP particle-in-cell (PIC) code [7] to provide a detailed understanding of the beam dynamics [2]. In order to investigate the plasma effects on the electron beam, the plasma kinetic must be implemented in the LSP-PIC