For band decomposed charge density you have to first identify the band and k-point corresponding to the VBM and CBM. Make an input file shown below and perform pp.x calculation.
I totally agree with your answer. But still I have some major doubts-
1. From where can I get to know the band number at which VBM or CBM occurs (if I am not wrong, Is it found from the xmgr file that is generated for plotting band structure).
2. Mine is a hexagonal system and the electronic band structure is plotted at M-K-G-M, then how to select the k-point number??? I mean suppose my VBM is somewhere between K-G and CBM is at M point. Then how should I select the kpoint value for the pp.x input file????
For identifying k-point, use "verbosity='high' " input tag and see the .band file form in which the eigenvalues at each k-point are written. From that, you can get a k-point and see the xmgrace file to identify the k-band corresponding to VBM and CBM.