for pidCG in range(num_particles_CG):

    # we put CG molecule in first atom, later CG molecules will be positioned in the center
    # note that we put the CG molecule at the first atom's position. Later the CG molecule will be positioned in the center of mass of all it's atoms
    tmptuple = [pidCG+num_particles_AT+1]