I am trying to calculate the sasa for each atom but I cannot get the arguments correct for rosetta.core.scoring.calc_per_atom_sasa.
I guess for "utility::vector1< Real > & rsd_sasa" I can use rosetta.utility.vector1_Real().
However, I have trouble with "id::AtomID_Map< Real > & atom_sasa". I tried rosetta.core.id.AtomID_Map_AtomID() but it doesn't match the C++ signature.
I would like help to understand how to properly set in pyrosetta "id::AtomID_Map< Real > atom_sasa" as well as "id::AtomID_Map< bool > atom_subset" and initialise it to true for heavy only atoms.
Many thanks for the help.