EvtSpinAmp Class Reference

#include <EvtSpinAmp.hh>

Public Member Functions
	EvtSpinAmp ()
	EvtSpinAmp (const vector< EvtSpinType::spintype > &twospin)
	EvtSpinAmp (const vector< EvtSpinType::spintype > &twospin, const EvtComplex &val)
	EvtSpinAmp (const vector< EvtSpinType::spintype > &twospin, const vector< EvtComplex > &elem)
	EvtSpinAmp (const EvtSpinAmp &)
	~EvtSpinAmp ()
EvtComplex &	operator() (const vector< int > &)
const EvtComplex &	operator() (const vector< int > &) const
EvtComplex &	operator() (int,...)
const EvtComplex &	operator() (int,...) const
EvtSpinAmp &	operator= (const EvtSpinAmp &)
EvtSpinAmp	operator+ (const EvtSpinAmp &) const
EvtSpinAmp &	operator+= (const EvtSpinAmp &)
EvtSpinAmp	operator- (const EvtSpinAmp &) const
EvtSpinAmp &	operator-= (const EvtSpinAmp &)
EvtSpinAmp	operator* (const EvtSpinAmp &) const
EvtSpinAmp &	operator*= (const EvtSpinAmp &)
EvtSpinAmp &	operator*= (const EvtComplex &)
EvtSpinAmp &	operator/= (const EvtComplex &)
void	intcont (int, int)
void	extcont (const EvtSpinAmp &, int, int)
void	assign (const EvtComplex &val)
int	rank () const
const vector< int > &	dims () const
void	addspin (int twospin)
void	setelem (const vector< EvtComplex > &elem)
bool	iterate (vector< int > &index) const
vector< int >	iterinit () const
bool	allowed (const vector< int > &index) const
bool	iterateallowed (vector< int > &index) const
vector< int >	iterallowedinit () const

Friends
EvtSpinAmp	operator* (const EvtComplex &real, const EvtSpinAmp &cont)
EvtSpinAmp	operator* (const EvtSpinAmp &cont, const EvtComplex &real)
EvtSpinAmp	operator/ (const EvtSpinAmp &cont, const EvtComplex &real)
std::ostream &	operator<< (std::ostream &os, const EvtSpinAmp &amp)

Detailed Description

Definition at line 17 of file EvtSpinAmp.hh.

Constructor & Destructor Documentation

EvtSpinAmp() [1/5]

EvtSpinAmp::EvtSpinAmp ( )

inline

Definition at line 27 of file EvtSpinAmp.hh.

{};

EvtSpinAmp() [2/5]

EvtSpinAmp::EvtSpinAmp

(

const vector< EvtSpinType::spintype > &

twospin

)

Definition at line 63 of file EvtSpinAmp.cc.

{
    int num = 1;
    _type = type;
    _twospin=calctwospin( type );
    
    for( int i=0; i<_twospin.size(); ++i )
        num*=_twospin[i]+1;
 
    _elem=vector<EvtComplex>( num );
}

EvtSpinAmp() [3/5]

EvtSpinAmp::EvtSpinAmp	(	const vector< EvtSpinType::spintype > &	twospin,
		const EvtComplex &	val )

Definition at line 75 of file EvtSpinAmp.cc.

{
    int num = 1;
    _type = type;
    _twospin=calctwospin( type );
 
    for( int i=0; i<_twospin.size(); ++i )
        num*=_twospin[i]+1;
 
    _elem=vector<EvtComplex>( num, val );
}

EvtSpinAmp() [4/5]

EvtSpinAmp::EvtSpinAmp	(	const vector< EvtSpinType::spintype > &	twospin,
		const vector< EvtComplex > &	elem )

Definition at line 87 of file EvtSpinAmp.cc.

{  
    int num = 1;
 
    _type = type;
    _twospin=calctwospin( type );
    _elem=elem;
    
    for( int i=0; i<_twospin.size(); ++i )
        num*=_twospin[i]+1;
 
    if(_elem.size() != num ) {
        report(ERROR,"EvtGen")<<"Wrong number of elements input:"
            <<_elem.size()<<" vs. "<<num<<endl;
       ::abort(); 
    }
 
}

EvtSpinAmp() [5/5]

EvtSpinAmp::EvtSpinAmp

(

const EvtSpinAmp &

copy

)

Definition at line 107 of file EvtSpinAmp.cc.

{
    _twospin = copy._twospin;
    _elem = copy._elem;
    _type = copy._type;
}

~EvtSpinAmp()

EvtSpinAmp::~EvtSpinAmp ( )

inline

Definition at line 33 of file EvtSpinAmp.hh.

{};

Member Function Documentation

addspin()

void EvtSpinAmp::addspin ( int twospin )

inline

Definition at line 71 of file EvtSpinAmp.hh.

{ _twospin.push_back( twospin ); }

allowed()

bool EvtSpinAmp::allowed

(

const vector< int > &

index

)

const

Definition at line 386 of file EvtSpinAmp.cc.

{
    if( index.size() != _type.size() ) {
        report(ERROR,"EvtGen")
            <<"Wrong dimensino index input to allowed."<<endl;
        ::abort();
    }
    
    for(int i=0; i<index.size(); ++i) {
        switch(_type[i]) {
            case EvtSpinType::PHOTON: 
                if(abs(index[i])!=2) return false;
                break;
            case EvtSpinType::NEUTRINO:
                report(ERROR,"EvtGen")
                    <<"EvMultibody currently cannot handle neutrinos."<<endl;
                ::abort();
        default:
          break;
        }
    }
    
    return true;
}

Referenced by EvtMHelAmp::EvtMHelAmp(), iterallowedinit(), and iterateallowed().

assign()

void EvtSpinAmp::assign ( const EvtComplex & val )

inline

Definition at line 61 of file EvtSpinAmp.hh.

{ _elem.assign( _elem.size(), val ); }

dims()

const vector< int > & EvtSpinAmp::dims ( ) const

inline

Definition at line 67 of file EvtSpinAmp.hh.

{ return _twospin; } 

Referenced by EvtMultibody::decay().

extcont()

void EvtSpinAmp::extcont	(	const EvtSpinAmp &	cont,
		int	a,
		int	b )

Definition at line 489 of file EvtSpinAmp.cc.

{
    EvtSpinAmp ret = (*this)*cont;
    ret.intcont( a, rank()+b );
 
    *this=ret;
}

Referenced by EvtMHelAmp::amplitude(), and EvtMultibody::decay().

intcont()

void EvtSpinAmp::intcont	(	int	a,
		int	b )

Definition at line 431 of file EvtSpinAmp.cc.

{
    int newrank=rank()-2;
    if(newrank<=0) {
        report(ERROR,"EvtGen")<<"EvtSpinAmp can't handle no indices" << endl;
        ::abort();
    }
 
    if(_twospin[a]!=_twospin[b]) {
        report(ERROR,"EvtGen")
            <<"Contaction called on indices of different dimension" 
            <<endl;
        report(ERROR,"EvtGen")
            <<"Called on "<<_twospin[a]<<" and "<<_twospin[b]
            <<endl;
        ::abort();
    }
 
    vector<int> newtwospin( newrank );
    vector<EvtSpinType::spintype> newtype( newrank );
 
    for( int i=0, j=0; i<_twospin.size(); ++i ){
        if(i==a || i==b) continue;
 
        newtwospin[j] = _twospin[i];
        newtype[j] = _type[i];
        ++j;
    }
 
    EvtSpinAmp newamp( newtype );
    vector<int> index( rank() ), newindex = newamp.iterinit();
 
    for( int i=0; i<newrank; ++i )
        newindex[i] = -newtwospin[i];
 
    while(true) {
        
        for( int i=0, j=0; i<rank(); ++i ) {
            if(i==a || i==b) continue;
            index[i]=newindex[j];
            ++j;
        }
        
        index[b]=index[a]=-_twospin[a];
        newamp(newindex) = (*this)(index);
        for( int i=-_twospin[a]+2; i<=_twospin[a]; i+=2 ) {
            index[b]=index[a]=i;
            newamp(newindex) += (*this)(index);
        }
       
        if(!newamp.iterate(newindex)) break;
    }
    
    *this=newamp;
}

Referenced by extcont().

iterallowedinit()

vector< int > EvtSpinAmp::iterallowedinit

(

)

const

Definition at line 421 of file EvtSpinAmp.cc.

{
    vector<int> init = iterinit();
    while(!allowed(init)) {
        iterate(init);
    }
 
    return init;
}

Referenced by EvtMHelAmp::amplitude(), and EvtMultibody::decay().

iterate()

bool EvtSpinAmp::iterate

(

vector< int > &

index

)

const

Definition at line 369 of file EvtSpinAmp.cc.

{
    int last = _twospin.size() - 1;
 
    index[0]+=2;
    for( int j=0; j<last; ++j ) {
        if( index[j] > _twospin[j] ) {
            index[j] = -_twospin[j];
            index[j+1]+=2;
        }
    }
 
    return abs(index[last])<=_twospin[last];
}

Referenced by EvtMHelAmp::EvtMHelAmp(), intcont(), iterallowedinit(), iterateallowed(), and operator*().

iterateallowed()

bool EvtSpinAmp::iterateallowed

(

vector< int > &

index

)

const

Definition at line 411 of file EvtSpinAmp.cc.

{
    while(true) {
        if(!iterate( index ))
            return false;
        if(allowed( index )) 
            return true;
    }
}

Referenced by EvtMHelAmp::amplitude(), and EvtMultibody::decay().

iterinit()

vector< int > EvtSpinAmp::iterinit

(

)

const

Definition at line 359 of file EvtSpinAmp.cc.

{
    vector<int> init( _twospin.size() );
 
    for( int i=0; i<_twospin.size(); ++i )
        init[i]=-_twospin[i];
 
    return init;
}

Referenced by EvtMHelAmp::EvtMHelAmp(), intcont(), and iterallowedinit().

operator()() [1/4]

EvtComplex & EvtSpinAmp::operator()

(

const vector< int > &

index

)

Definition at line 168 of file EvtSpinAmp.cc.

{
    checkindexargs( index );
    
    int trueindex = findtrueindex(index);
    if(trueindex >= _elem.size()) {
        report(ERROR,"EvtGen")<<"indexing error "<<trueindex<<" "<<_elem.size()<<endl;
        for(int i=0; i<_twospin.size(); ++i) {
            report(ERROR,"")<<_twospin[i]<<" ";
        }
        report(ERROR,"")<<endl;
 
        for(int i=0; i<index.size(); ++i) {
            report(ERROR,"")<<index[i]<<" ";
        }
        report(ERROR,"")<<endl;
 
        ::abort();
    }
 
    return _elem[trueindex];
}

operator()() [2/4]

const EvtComplex & EvtSpinAmp::operator()

(

const vector< int > &

index

)

const

Definition at line 191 of file EvtSpinAmp.cc.

{
    checkindexargs( index );
    
    int trueindex = findtrueindex(index);
    if(trueindex >= _elem.size()) {
        report(ERROR,"EvtGen")<<"indexing error "<<trueindex<<" "<<_elem.size()<<endl;
        for(int i=0; i<_twospin.size(); ++i) {
            report(ERROR,"")<<_twospin[i]<<" ";
        }
        report(ERROR,"")<<endl;
 
        for(int i=0; i<index.size(); ++i) {
            report(ERROR,"")<<index[i]<<" ";
        }
        report(ERROR,"")<<endl;
 
        ::abort();
    }
    
    return _elem[trueindex];
}

operator()() [3/4]

EvtComplex & EvtSpinAmp::operator()	(	int	i,
			... )

Definition at line 214 of file EvtSpinAmp.cc.

{
    va_list ap;
    vector<int> index( _twospin.size() );
    
    va_start(ap, i);
    
    index[0]=i;
    for(int n=1; n<_twospin.size(); ++n )
        index[n]=va_arg( ap, int );
 
    va_end(ap);
 
    return (*this)( index );
}

operator()() [4/4]

const EvtComplex & EvtSpinAmp::operator()	(	int	i,
			... ) const

Definition at line 230 of file EvtSpinAmp.cc.

{
    vector<int> index( _twospin.size() );
    va_list ap;
 
    va_start(ap, i);
 
    index[0]=i;
    for(int n=1; n<_twospin.size(); ++n )
        index[n]=va_arg( ap, int );
 
    va_end(ap);
 
    return (*this)( index );
}

operator*()

EvtSpinAmp EvtSpinAmp::operator*

(

const EvtSpinAmp &

amp2

)

const

Definition at line 300 of file EvtSpinAmp.cc.

{
    vector<int> index(rank()+amp2.rank());
    vector<int> index1(rank()), index2(amp2.rank()); 
    EvtSpinAmp amp;
    
    amp._twospin=_twospin;
    amp._type=_type;
 
    for( int i=0; i<amp2._twospin.size(); ++i ) {
        amp._twospin.push_back( amp2._twospin[i] );
        amp._type.push_back( amp2._type[i] );
    }
    
    amp._elem = vector<EvtComplex>( _elem.size() * amp2._elem.size() );
 
    for( int i=0; i<index1.size(); ++i )
        index[i]=index1[i]=-_twospin[i];
    
    for( int i=0; i<index2.size(); ++i )
        index[i+rank()]=index2[i]=-amp2._twospin[i];
 
    while(true) {
        amp( index ) = (*this)( index1 )*amp2( index2 );
        if(!amp.iterate( index )) break;
        
        for( int i=0; i<index1.size(); ++i )
            index1[i]=index[i];
 
        for( int i=0; i<index2.size(); ++i )
            index2[i]=index[i+rank()];
    }
    
    return amp;
}

operator*=() [1/2]

EvtSpinAmp & EvtSpinAmp::operator*=

(

const EvtComplex &

real

)

Definition at line 343 of file EvtSpinAmp.cc.

{
    for( int i=0; i<_elem.size(); ++i )
        _elem[i] *= real;
 
    return *this;
}

operator*=() [2/2]

EvtSpinAmp & EvtSpinAmp::operator*=

(

const EvtSpinAmp &

cont

)

Definition at line 336 of file EvtSpinAmp.cc.

{
    EvtSpinAmp ret = (*this)*cont;
    *this=ret;
    return *this;
}

operator+()

EvtSpinAmp EvtSpinAmp::operator+

(

const EvtSpinAmp &

cont

)

const

Definition at line 255 of file EvtSpinAmp.cc.

{
    checktwospin( cont._twospin );
 
    EvtSpinAmp ret( cont );
    for( int i=0; i<ret._elem.size(); ++i ) {
        ret._elem[i]+=_elem[i];
    }
 
    return ret;
}

operator+=()

EvtSpinAmp & EvtSpinAmp::operator+=

(

const EvtSpinAmp &

cont

)

Definition at line 267 of file EvtSpinAmp.cc.

{
    checktwospin( cont._twospin );
 
    for( int i=0; i<_elem.size(); ++i )
        _elem[i]+=cont._elem[i];
 
    return *this;
}

operator-()

EvtSpinAmp EvtSpinAmp::operator-

(

const EvtSpinAmp &

cont

)

const

Definition at line 278 of file EvtSpinAmp.cc.

{
    checktwospin( cont._twospin );
 
    EvtSpinAmp ret( *this );
    for( int i=0; i<ret._elem.size(); ++i )
        ret._elem[i]-=cont._elem[i];
 
    return ret;
}

operator-=()

EvtSpinAmp & EvtSpinAmp::operator-=

(

const EvtSpinAmp &

cont

)

Definition at line 289 of file EvtSpinAmp.cc.

{
    checktwospin( cont._twospin );
 
    for( int i=0; i<_elem.size(); ++i )
        _elem[i]-=cont._elem[i];
 
    return *this;
}

operator/=()

EvtSpinAmp & EvtSpinAmp::operator/=

(

const EvtComplex &

real

)

Definition at line 351 of file EvtSpinAmp.cc.

{
    for( int i=0; i<_elem.size(); ++i )
        _elem[i] /= real;
    
    return *this;
}

operator=()

EvtSpinAmp & EvtSpinAmp::operator=

(

const EvtSpinAmp &

cont

)

Definition at line 246 of file EvtSpinAmp.cc.

{
    _twospin=cont._twospin;
    _elem=cont._elem;
    _type=cont._type;
 
    return *this;
}

rank()

int EvtSpinAmp::rank ( ) const

inline

Definition at line 64 of file EvtSpinAmp.hh.

{ return _twospin.size(); }

Referenced by EvtMultibody::decay(), extcont(), intcont(), and operator*().

setelem()

void EvtSpinAmp::setelem ( const vector< EvtComplex > & elem )

inline

Definition at line 72 of file EvtSpinAmp.hh.

{ _elem = elem; }

Friends And Related Symbol Documentation

operator* [1/2]

EvtSpinAmp operator* ( const EvtComplex & real, const EvtSpinAmp & cont )

friend

Definition at line 25 of file EvtSpinAmp.cc.

{
    EvtSpinAmp ret( cont );
 
    for( int i=0; i<ret._elem.size(); ++i ) {
        ret._elem[i] *= real;
    }
 
    return ret;
}

operator* [2/2]

EvtSpinAmp operator* ( const EvtSpinAmp & cont, const EvtComplex & real )

friend

Definition at line 36 of file EvtSpinAmp.cc.

{
    return real*cont;
}

operator/

EvtSpinAmp operator/ ( const EvtSpinAmp & cont, const EvtComplex & real )

friend

Definition at line 41 of file EvtSpinAmp.cc.

{
    EvtSpinAmp ret( cont );
 
    for( int i=0; i<ret._elem.size(); ++i ) {
        ret._elem[i] /= real;
    }
 
    return ret;
}

operator<<

std::ostream & operator<< ( std::ostream & os, const EvtSpinAmp & amp )

friend

Definition at line 8 of file EvtSpinAmp.cc.

{
    vector<int> index = amp.iterinit();
    
    os << ":";
    do {
        os<<"<";
        for(int i=0; i<index.size()-1; ++i) {
            os<<index[i];
        }
        os<<index[index.size()-1]<<">"<<amp(index)<<":";
    } while( amp.iterate( index ) );
 
    return os;
}

---

The documentation for this class was generated from the following files:

• EvtSpinAmp.hh
• EvtSpinAmp.cc