Col_str.h

// -*- 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 */