Structure and dynamics of size-selected charged pyrrole clusters have been studied by means of molecular beam scattering experiments and ab initio calculations. Small neutral Pyn clusters were produced in Py/He mixture expansions, and the scattering experiment with a secondary beam of He-atoms was exploited to select the neutral clusters of different sizes. The complete size-selected fragmentation patterns for the neutral dimer to the tetramer after an electron impact ionization at 70 eV were obtained from the measurements of the angular and velocity distributions at different fragment masses. All the investigated cluster sizes decay mainly to the monomer ions Py+1 (from 60 to 80% of the corresponding neutral size) and to the dimer ion Py+2 (20-30%). The trimer ions Py+3 are generated to less than 10% from the neutral trimer and tetramer. To explain the observed results, we have calculated the structures and energetics of pyrrole clusters up to the trimer for the neutral and the ionic state using DFT and PMP2 methods. The ab initio calculations show that ionized pyrrole clusters are formed with a dimeric core that is solvated by neutral pyrrole molecules. In addition, the ground and ionic state of Py-Ar complexes were calculated at CCSD(T) level with extended basis in relevance to the mixed clusters produced in supersonic expansions of Py seeded in Ar. The calculated dissociation energies of the Py-Ar and (Py-Ar)+ complexes indicate that Ar atoms are able to rapidly evaporate after ionization. The combined analysis of the fragmentation probabilities, and calculations allowed us to estimate the distribution of energy deposited in the clusters after the electron impact, which peaks above 1 eV and has a tail up to 5 eV.