We have developed a Monte Carlo method to estimate the transport of different groups of particles for plasmas in toroidal geometries. This method can determine the important transport mechanisms driving the anomalous transport by comparing the numerical results with the experimental data. The important groups of particles whose transport can be estimated by this method include runaway electrons, thermal electrons, both passing and trapped diagnostic beam ions etc. The three basic mechanisms driving the anomalous transport are: spatial variation of magnetic field strength, spatial variation of electrostatic potential within the flux surfaces, and the loss of flux surfaces. The equation of motion are obtained from the drift hamiltonian. The equations of motion are developed in the canonical and in the non-canonical, practical co-ordinates as well. The effects of collisions are represented by appropriate stochastic changes in the constants of motion at each time-step. Here we present the results of application of this method to three cases: superathermal alphas in the rippled field of tokamaks, motion in the magnetic turbulence of takapole II, and transport in the stochastic fields of ZT40. This work is supported by DOE OFE and ORAU HBCU program.