1 #include "isajet/pilot.h"
3 C CALCULATE TOTALLY ANTISYMMETRIC TENSOR EPSILON CONTRACTED
5 #include "isajet/itapes.inc"
6 DIMENSION A(4),B(4),C(4),D(4)
7 #if defined(CERNLIB_DOUBLE)
9 DOUBLE PRECISION A,B,C,D,CD,BCD
11 CD(I,J)=C(I)*D(J)-C(J)*D(I)
12 BCD(I,J,K)=B(I)*CD(J,K)-B(J)*CD(I,K)+B(K)*CD(I,J)
13 EPF=A(1)*BCD(2,3,4)-A(2)*BCD(1,3,4)+A(3)*BCD(1,2,4)