Aims. In this work, we present an extended study of the dust environment of a sample of short period comets and their dynamical history. With this aim, we characterized the dust tails when the comets are active, and we made a statistical study to determine their dynamical evolution. The targets selected were 22P/Kopff, 30P/Reinmuth 1, 78P/Gehrels 2, 115P/Maury, 118P/Shoemaker-Levy 4, 123P/West-Hartley, 157P/Tritton, 185/Petriew, and P/2011 W2 (Rinner). Methods. We use two different observational data: a set of images taken at the Observatorio de Sierra Nevada and the Afrho curves provided by the amateur astronomical association Cometas-Obs. To model these observations, we use our Monte Carlo dust tail code. From this analysis, we derive the dust parameters, which best describe the dust environment: dust loss rates, ejection velocities, and size distribution of particles. On the other hand, we use a numerical integrator to study the dynamical history of the comets, which allows us to determine with a 90% of confidence level the time spent by these objects in the region of Jupiter Family Comets. Results. From the Monte Carlo dust tail code, we derived three categories attending to the amount of dust emitted: Weakly active (115P, 157P, and Rinner), moderately active (30P, 123P, and 185P), and highly active (22P, 78P, and 118P). The dynamical studies showed that the comets of this sample are young in the Jupiter Family region, where the youngest ones are 22P (100 yr), 78P (500 yr), and 118P (600 yr). The study points to a certain correlation between comet activity and time spent in the Jupiter Family region, although this trend is not always fulfilled. The largest particle sizes are not tightly constrained, so that the total dust mass derived should be regarded as lower limits.}

We present optical observations and Monte Carlo models of the dust coma, tail, and trail structures of comet 22P/Kopff during the 2002 and 2009 apparitions. Dust loss rates, ejection velocities, and power-law size distribution functions are derived as functions of the heliocentric distance using pre- and post-perihelion imaging observations during both apparitions. The 2009 post-perihelion images can be accurately fitted by an isotropic ejection model. On the other hand, strong dust ejection anisotropies are required to fit the near-coma regions at large heliocentric distances (both inbound at $r_h$=2.5 AU and outbound at $r_h$=2.6 AU) for the 2002 apparition. These asymmetries are compatible with a scenario where dust ejection is mostly seasonally-driven, coming mainly from regions near subsolar latitudes at far heliocentric distances inbound and outbound. At intermediate to near-perihelion heliocentric distances, the outgassing would affect much more extended latitude regions, the emission becoming almost isotropic near perihelion. We derived a maximum dust production rate of 260 kg s$^{-1}$ at perihelion, and an averaged production rate over one orbit of 40 kg s$^{-1}$. An enhanced emission rate, accompanied also by a large ejection velocity, is predicted at $r_h>$2.5 pre-perihelion. The model has also been extended to the thermal infrared in order to be applied to available trail observations with IRAS and ISO spacecrafts of this comet. The modeled trail intensities are in good agreement with those observations, which is remarkable taking into account that those data are sensitive to dust ejection patterns corresponding to several orbits before the 2002 and 2009 apparitions.