Theta-diagram technique for 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.