/data/Documents/eclipse/workspace/ColorFull/Col_str.h

Go to the documentation of this file.

// -*- C++ -*-
/*
 *  Col_str.h
 *      Contains declaration of the class Col_str and associated types and operators.
 *  Created on: Jul 7, 2010
 *  Author: Malin Sjodahl
 */

#ifndef COLORFULL_Col_str_h
#define COLORFULL_Col_str_h

#include "Quark_line.h"

namespace ColorFull {

typedef std::vector < Quark_line > col_str;


class Col_str {
public:

        Col_str(){};

        Col_str( const std::string str );

        Col_str( Quark_line Ql ) {cs.push_back(Ql);}

        col_str cs;

        Polynomial Poly;

        const Quark_line & at( int i ) const {return cs.at(i);}

        Quark_line & at( int i ) {return cs.at(i);}

        int at( int i, int j ) const;

        uint size() const{ return cs.size(); }

        bool empty() const { return cs.empty(); }

        void clear() { cs.clear(); }

        void erase( int i );

        void erase( int i, int j );

        void erase(std::pair<int, int> place);

        void append( Quark_line Ql ) { cs.push_back( Ql ); }

        void insert( int i, int j, int part_num );

        void append( col_str cs_in );

        void read_in_Col_str( std::string filename );

        void write_out_Col_str( std::string filename ) const;

        std::pair<int, int> find_parton( int part_num ) const;

        bool neighbor( int p1, int p2 ) const;

        bool right_neighbor( int p1, int p2 ) const;

        bool left_neighbor( int p1, int p2 ) const;

        void replace( int old_ind, int new_ind );

        std::string find_kind( int p ) const;

        bool gluons_only() const;

        int n_gluon() const;

        int n_quark() const;

        void normal_order();

        int smallest( const Col_str & Cs1, const Col_str & Cs2 ) const;

        int longest_quark_line() const;

        void remove_1_rings();

        void remove_0_rings();

        void simplify();

        void conjugate();

        void contract_next_neighboring_gluons( );

        void contract_2_rings( );

        void contract_quarks( const Col_str Cs1, const Col_str Cs2 );


private:

        int compare_quark_lines( int i1, int i2 ) const;

        void Col_str_of_str( const std::string str  );

        void col_str_of_str( std::string );

}; //end class Col_str


bool operator==(const col_str & cs1, const col_str & cs2);

bool operator!=(const col_str & cs1, const col_str & cs2);

std::ostream& operator<<( std::ostream& out, const col_str & cs );

Col_str operator*( const Col_str & Cs, const int i );

Col_str operator*( const int i, const Col_str & Cs );

Col_str operator*( const Col_str & Cs, const double d );

Col_str operator*( const double d, const Col_str & Cs );

Col_str operator*( const Col_str & Cs, const cnum c);

Col_str operator*( const cnum c, const Col_str & Cs );

Col_str operator*( const Col_str & Cs, const Monomial & Mon );

Col_str operator*( const Monomial & Mon, const Col_str & Cs);

Col_str operator*( const Col_str & Cs, const Polynomial & Poly );

Col_str operator*( const Polynomial & Poly, const Col_str & Cs );

Col_str operator*( const Col_str & Cs, const Quark_line & Ql);

Col_str operator*( const Quark_line & Ql, const Col_str & Cs );

Col_str operator*( const Col_str & Cs1, const Col_str & Cs2 );

Col_str operator*( const Quark_line & Ql1, const Quark_line & Ql2 );

std::ostream& operator<<(std::ostream& out, const Col_str & Cs);

bool operator==( const Col_str & Cs1, const Col_str & Cs2 );

bool operator!=( const Col_str & Cs1, const Col_str & Cs2 );

}

#endif /* COLORFULL_Col_str_h */