This is an on-going problem in the field. The main issue is how to divide/distribute the charge surrounding each atoms. This problem is easy for single atom systems but the difficulty scales rapidly as more atoms are introduced and electron exchange (bonding) happens. In this sense, you won't really know to which atom does a particular charge density belongs to.

There are some codes that have proposed some sort of charge calculation (using charge density). I would suggest exploring Bader Charge, Lowdin charge (the one in DOS calculation) and Mulliken charge. These implement different "slicing" methods for the charge density.

TLDR: Check Bader charge analysis, and Lowdin charge in DOS output

you can easily find those charge distributions in the Mulliken charge analysis using gaussian 09 and AIMAll quantum computation software.

Khatua, you need to first of all know which is the dependent and independent variable. Plot the dependent variable on the vertical and independent on the horizontal axis