polarization.h

Go to the documentation of this file.

// This file is part of MARTY.
// 
// MARTY is free software: you can redistribute it and/or modify
// it under the terms of the GNU General Public License as published by
// the Free Software Foundation, either version 3 of the License, or
// (at your option) any later version.
// 
// MARTY is distributed in the hope that it will be useful,
// but WITHOUT ANY WARRANTY; without even the implied warranty of
// MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
// GNU General Public License for more details.
// 
// You should have received a copy of the GNU General Public License
// along with MARTY. If not, see <https://www.gnu.org/licenses/>.

#ifndef POLARIZATION_H_INCLUDED
#define POLARIZATION_H_INCLUDED

#include <csl.h>
#include "quantumField.h"

namespace mty {

class PolarizationField: public mty::QuantumField {

    public:

        void test(int);

    PolarizationField();

    PolarizationField(csl::Tensor const& impulsion,
                      csl::Parent const&  t_parent);

    PolarizationField(csl::Tensor const&            impulsion,
                      csl::Parent const&             t_parent,
                      const std::vector<csl::Index>& indices);

    PolarizationField(csl::Tensor const&        impulsion,
                      csl::Parent const&         t_parent,
                      const csl::IndexStructure& indices);

    PolarizationField(PolarizationField const&) = default;

    bool getCommutable() const override;

    std::optional<csl::Expr> getComplexConjugate() const override;

    bool hasContractionProperty(csl::Expr_info other) const override;

    csl::Expr contraction(csl::Expr_info other) const override;

    csl::unique_Expr copy_unique() const override;

    void printCode(
            int mode = 0,
            std::ostream &out = std::cout
            ) const override;

    void printProp(std::ostream& out = std::cout) const override;

    inline 
    Particle getParticle() const {
        return std::static_pointer_cast<QuantumFieldParent>(parent);
    }

    bool isConjugationLocked() const {
        return lockConjugation;
    }

    void setConjugationLock(bool t_lock) {
        lockConjugation = t_lock;
    }

    void setParticle(bool) override;
    void setIncoming(bool) override;

    friend
    std::ostream& operator<<(std::ostream& fout, const PolarizationField& pol);

    private:

    PolarizationField complexConjugatedField(
            bool keepFermionOrder = false
            ) const;

    bool hasFieldChargeConjugation(
            mty::PolarizationField const *other
            ) const;

    bool spaceIndexContraction(csl::Expr_info other) const;
    bool spinIndexContraction(csl::Expr_info other) const;

    bool isContractibleWith(csl::Expr_info otherPolar) const;

    csl::Expr sumPolarization(csl::Expr_info other) const;

    csl::Expr matrixChargeConjugation(csl::Expr_info other) const override;
    csl::Expr fieldChargeConjugation(mty::PolarizationField const *other) const;

    void updateComplexConjugation();

    private:

    bool lockConjugation;
};

} // End of namespace mty

#endif