def AromaticAtoms(m):
    """
    This function calculates the number of aromatic atoms in a molecule.
    The input argument, m, is a rdkit molecular object
    """
    
    aromatic_atoms = [m.GetAtomWithIdx(i).GetIsAromatic() for i in range(m.GetNumAtoms())]
    aa_count = []
    for i in aromatic_atoms:
        if i==True:
            aa_count.append(1)
        
    sum_aa_count = sum(aa_count)
    
    return sum_aa_count