Parton jets in the hot and dense medium of a Quark Gluon Plasma (QGP) can undergo multiple processes of scatterings off medium particles as well as processes of coherent medium induced radiations. A Monte-Carlo algorithm and resulting program is presented that allows to obtain jets that were formed by these two types of processes from an initial highly energetic quark or gluon. The program accounts for the increase in the momentum components of jet-particles transverse to the jet-axis due to processes of scattering as well as medium induced radiations in addition to energy-loss due to the medium induced radiations. Program summaryProgram Title:TMDICECPC Library link to program files:https://doi.org/10.17632/jbbmgmtyrr.1Code Ocean capsule:https://codeocean.com/capsule/6186360Licensing provisions: GPLv3Programming language: C++, BashNature of problem : In order to describe the fragmentation of parton cascades/jets in the medium processes of coherent medium induced radiation, where a particle emission is formed simultaneously to multiple scatterings off medium particles, need to be considered (in addition to scatterings off medium particles without emissions) [1,2,3,4,5,6,7]. A description of jet-fragmentation in the medium needs to be found that uses the effective kernels for coherent medium induced radiation and scattering [8,9] and which provides distributions of jet-particles as a function of the time the jet needs for passing the medium.Solution method: A Monte-Carlo method is presented that allows to obtain a set of jet particles from a predefined initial particle. To this end, the variables relevant for the description of jet particles (such as the time of emission, momentum fraction and momentum component transverse to jet axis, and if the parton is a quark/antiquark or a gluon) are sampling from probability density functions that were obtained from the kernels for coherent medium induced radiation and scattering off medium particles in [8,9]. This is achieved in a two step process: First, before obtaining jet particles, the corresponding cumulative distribution functions are calculated and its inverse of the cumulative distribution functions obtained numerically. Then, samples are obtained by random selection from the inverse of the cumulative distribution functions.Additional comments including restrictions and unusual features: The kernels [8,9] for the coherent medium induced radiations were derived within the eikonal approximation that only momentum components transverse to the incoming particles are affected by medium transfers. Furthermore, these kernels do not depend on the time of emission (thus, neglecting effects of the finite size of the medium within the emission and scattering kernels). For simplicity so far only a medium with constant parameters for the jet-medium interactions have been assumed.