def _get_catalysts_in_reaction(reaction: Reaction) -> Set[BaseAbundance]:
return {<EOL>reactant<EOL>for reactant in reaction.reactants<EOL>if reactant in reaction.products<EOL>}<EOL>
Return nodes that are both in reactants and reactions in a reaction.