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

Go to the documentation of this file.

// -*- C++ -*-
/*
 * Col_functions.h
 * Contains declarations of the class Col_str and associated types and operators
 * Author: Malin Sjodahl
 */

#ifndef COLORFULL_Col_functions_h
#define COLORFULL_Col_functions_h


#include "Col_amp.h"
#include "Poly_matr.h"
#include <list>
#include <map>
#include <memory>

namespace ColorFull {

using std::shared_ptr;

class Col_functions {

private:

double Nc;

double TR;

double CF;

bool full_CF;


public:

Col_functions()
  : Nc(3.0), TR(0.5), CF(4.0/3.0), full_CF( false )  {}

void set_Nc( double n) {
  Nc = n;
  CF = TR*(Nc*Nc-1.)/Nc;
}

void set_TR( double tr) {
  CF *= tr/TR;
  TR = tr;
}

void set_CF( double cf) {
  CF = cf;
}

void set_full_CF( bool is_full ) { full_CF = is_full; }

double get_Nc() const { return Nc; }

double get_TR() const { return TR; }

double get_CF() const { return CF; }

bool get_full_CF() const { return full_CF; }


/****************** Functions for leading terms *************************/
// The functions called leading(...) depend on the variable full_CF.
// As it would be messy to let each Polynomial carry around
// its own full_CF, these functions are kept here.

int  leading_Nc_pow( const Polynomial & Poly ) const;

int  leading_Nc_pow( const Poly_vec & Pv ) const;

/*
int  leading_Nc_pow( const std::vector< shared_ptr<Polynomial> > & Pvp) const;
*/

Polynomial leading( const Polynomial & Poly ) const;

// Used only by Poly_matr version of leading
Poly_vec leading( const Poly_vec & Pv ) const;

// Used only once in Col_basis
Poly_matr leading( const Poly_matr & Pm ) const;

/*
// Currently never used
Poly_vec leading( const std::vector<shared_ptr<Polynomial> > & Pvp) const;

// used only by scalar_product_matrix_mem in Col_functions
dmatr leading( const std::vector< std::vector< shared_ptr<Polynomial> > > & Pm ) const;
*/

/********************* Functions for numerical evaluation *************************/
// These functions has to be kept in Col_functions class as they need numerical
// values for evaluation. Letting each Polynomial carry around its own Nc etc.
// would be messy.


cnum cnum_num( const Monomial & Mon ) const;

cnum cnum_num( const Polynomial & Poly ) const;

cvec cnum_num( const Poly_vec & Pv ) const;

cmatr cnum_num( const Poly_matr & Pm ) const;

double double_num( const Monomial & Mon ) const;

double double_num( const Polynomial & Poly ) const;

dvec double_num( const Poly_vec & Pv ) const;

dmatr double_num( const Poly_matr & Pm ) const;

/*
dvec double_num( const std::vector<std::shared_ptr<Polynomial> > & Pv ) const;

// (not used 14 06 08)
dmatr double_num( const std::vector<std::vector<std::shared_ptr<Polynomial> > > & Pm ) const;

// (not used 14 06 08)
std::map< std::string, double > double_num( std::map< std::string, std::shared_ptr<Polynomial> > mem_map ) const;

// (not used 14 06 08)
std::map< std::string, double > double_num( std::map< std::string, Polynomial > mem_map ) const;
*/

Polynomial Polynomial_cnum_num( const Polynomial & Poly ) const;

Poly_vec Poly_vec_cnum_num( const Poly_vec & Pv ) const;

Poly_matr Poly_matr_cnum_num( const Poly_matr & Pm ) const;


/****************** Functions for scalar products *************************/

Polynomial scalar_product( const Col_amp & Ca1, const Col_amp & Ca2 ) const;

Polynomial scalar_product( const Col_str & Cs1, const Col_str & Cs2 ) const;


/****************** Functions for gluon emission exchange, and splitting *************/

Col_amp emit_gluon( const Col_str & in_Col_str, int emitter, int g_new ) const;

Col_amp emit_gluon( const Col_amp & Ca_in, int emitter, int g_new ) const;

Col_amp split_gluon( const Col_str & in_Col_str, int g_old, int q_new, int qbar_new ) const;

Col_amp split_gluon( const Col_amp & in_Col_amp, int g_old, int q_new, int qbar_new ) const;


Col_amp exchange_gluon( const Col_str & Cs, int p1, int p2 ) const;

Col_amp exchange_gluon( const Col_amp & Ca, int p1, int p2 ) const;

Polynomial color_correlator( const Col_amp Ca, int i, int j ) const;


/********************* Other functions *************************/

// As dvec and dmatr are not classes some read and write functions
// are contained here.

dvec read_in_dvec( std::string filename ) const;

dmatr read_in_dmatr( std::string filename ) const;

void write_out_dvec( const dvec & dv, std::string filename ) const;

void write_out_dmatr( const dmatr & matr, std::string filename ) const;

int factorial( int i ) const;

std::map<int, int> default_parton_numbers( const Col_str &, int g_old, int q_new, int qbar_new ) const;


Col_str rename_partons( const Col_str &, const std::map <int, int> replacements ) const;


Col_amp rename_partons( const Col_amp &, const std::map <int, int> replacements ) const;


}; // end class Col_functions



std::list<int>::iterator operator+( std::list<int>::iterator x, int n );

std::list<int>::iterator operator-( std::list<int>::iterator x, int n );

std::list< Quark_line >::iterator operator+( std::list < Quark_line >::iterator x, int n );

col_str::iterator operator-( col_str::iterator x, int n );

std::ostream& operator<<( std::ostream& out, const std::vector<int> & vec );

std::ostream& operator<<( std::ostream& out, const cvec & cv );

std::ostream& operator<<( std::ostream& out, const dvec & dv );

std::ostream& operator<<( std::ostream& out, const cmatr & cm );

std::ostream& operator<<( std::ostream& out, const dmatr & matr );

std::ostream& operator<<( std::ostream& out, std::pair<int, int> pair );
} // end namespace ColorFull

#endif /* COLORFULL_Col_functions_h */