* @return newly constructed object */
ex dirac_gamma5(unsigned char rl = 0);
+/** This returns (dirac_ONE(rl) + dirac_gamma5(rl)). */
+ex dirac_gamma6(unsigned char rl = 0);
+
+/** This returns (dirac_ONE(rl) - dirac_gamma5(rl)). */
+ex dirac_gamma7(unsigned char rl = 0);
+
/** Create a term of the form e_mu * gamma~mu with a unique index mu.
*
* @param dim Dimension of index
* In particular, the functional is not always cyclic in D != 4 dimensions
* when gamma5 is involved.
*
- * @param rl Representation label */
-ex dirac_trace(const ex & e, unsigned char rl = 0);
+ * @param e Expression to take the trace of
+ * @param rl Representation label
+ * @param trONE Expression to be returned as the trace of the unit matrix */
+ex dirac_trace(const ex & e, unsigned char rl = 0, const ex & trONE = 4);
} // namespace GiNaC