$$\theta $$-diagram technique for $${\mathcal {N}}=1$$, $$d=4$$ superfields

We describe a diagrammatic procedure to carry out the Grassmann integration in super-Feynman diagrams of 4d theories expressed in terms of N=1 superfields. This method is alternative to the well known D-algebra approach. We develop it in detail for theories containing vector, chiral and anti-chiral superfields, with the type of interactions which occur in N=2 SYM theories with massless matter, but it would be possible to extend it to other cases. The main advantage is that this method is algorithmic; we implemented it as a Mathematica program that, given the description of a super Feynman diagram in momentum space, returns directly the polynomial in the momenta produced by the Grassmann integration.