- dirac_trace() handles arbitrary expressions (even unexpanded), including
[ginac.git] / ginac / clifford.h
index 8ac696528873244d1ee6463e933425f256353973..af3cb7883568dea7fda15ade6d1abf191fd4ca99 100644 (file)
@@ -55,6 +55,10 @@ protected:
        unsigned return_type(void) const { return return_types::noncommutative; }
        unsigned return_type_tinfo(void) const { return TINFO_clifford + representation_label; }
 
        unsigned return_type(void) const { return return_types::noncommutative; }
        unsigned return_type_tinfo(void) const { return TINFO_clifford + representation_label; }
 
+       // non-virtual functions in this class
+public:
+       unsigned char get_representation_label(void) const {return representation_label;}
+
        // member variables
 private:
        unsigned char representation_label; /**< Representation label to distinguish independent spin lines */
        // member variables
 private:
        unsigned char representation_label; /**< Representation label to distinguish independent spin lines */
@@ -84,7 +88,8 @@ public:
 };
 
 
 };
 
 
-/** This class represents the Dirac gamma5 object. */
+/** This class represents the Dirac gamma5 object which anticommutes with
+ *  all other gammas. */
 class diracgamma5 : public tensor
 {
        GINAC_DECLARE_REGISTERED_CLASS(diracgamma5, tensor)
 class diracgamma5 : public tensor
 {
        GINAC_DECLARE_REGISTERED_CLASS(diracgamma5, tensor)
@@ -121,6 +126,21 @@ ex dirac_gamma(const ex & mu, unsigned char rl = 0);
  *  @return newly constructed object */
 ex dirac_gamma5(unsigned char rl = 0);
 
  *  @return newly constructed object */
 ex dirac_gamma5(unsigned char rl = 0);
 
+/** Create a term of the form e_mu * gamma~mu with a unique index mu.
+ *
+ *  @param dim Dimension of index
+ *  @param rl Representation label */
+ex dirac_slash(const ex & e, const ex & dim, unsigned char rl = 0);
+
+/** Calculate the trace of an expression containing gamma objects with
+ *  a specified representation label. The computed trace is a linear
+ *  functional that is equal to the usual trace only in D = 4 dimensions.
+ *  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);
+
 
 } // namespace GiNaC
 
 
 } // namespace GiNaC