1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
120
121
122
123
124
125
126
127
128
129
130
131
132
133
134
135
136
137
138
139
140
141
142
143
144
145
146
147
148
149
150
151
152
153
154
155
156
157
158
159
|
// -*- C++ -*-
//
// ColourLine.cc is a part of ThePEG - Toolkit for HEP Event Generation
// Copyright (C) 1999-2011 Leif Lonnblad
//
// ThePEG is licenced under version 2 of the GPL, see COPYING for details.
// Please respect the MCnet academic guidelines, see GUIDELINES for details.
//
//
// This is the implementation of the non-inlined, non-templated member
// functions of the ColourLine class.
//
#include "ColourLine.h"
#include "ThePEG/EventRecord/Particle.h"
#include "ThePEG/EventRecord/Event.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Config/algorithm.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
using namespace ThePEG;
tColinePtr ColourLine::create(tPPtr col, tPPtr anti) {
if ( col->colourLine() || anti->antiColourLine() ) return tColinePtr();
ColinePtr l = new_ptr(ColourLine());
l->addColoured(col);
l->addAntiColoured(anti);
return l;
}
tColinePtr ColourLine::create(tPPtr col, bool anti) {
if ( col->colourLine(anti) ) return tColinePtr();
ColinePtr l = new_ptr(ColourLine());
l->addColoured(col, anti);
return l;
}
tColinePtr ColourLine::create(tColinePtr son1, tColinePtr son2,
tColinePtr sin1, tColinePtr sin2) {
if ( !son1 || !son2 || !sin1 || !sin2 ) return tColinePtr();
ColinePtr l = new_ptr(ColourLine());
l->theSourceNeighbours = make_pair(son1, son2);
son1->theSourceNeighbours = make_pair(son2, l);
son2->theSourceNeighbours = make_pair(l, son1);
l->theSinkNeighbours = make_pair(sin1, sin2);
sin1->theSinkNeighbours = make_pair(sin2, l);
sin2->theSinkNeighbours = make_pair(l, sin1);
son1->orphanedConnectors.push_back(l);
return l;
}
ColourLine::~ColourLine() {}
tPPtr ColourLine::startParticle() const {
if ( sourceNeighbours().first ) return tPPtr();
for ( tPVector::const_reverse_iterator it = antiColoured().rbegin();
it != antiColoured().rend(); ++it )
if ( !(**it).outgoingAntiColour() ) return *it;
return tPPtr();
}
tPPtr ColourLine::endParticle() const {
if ( sinkNeighbours().first ) return tPPtr();
for ( tPVector::const_reverse_iterator it = coloured().rbegin();
it != coloured().rend(); ++it )
if ( !(**it).outgoingColour() ) return *it;
return tPPtr();
}
void ColourLine::addAntiColouredIndexed(tPPtr p, int index) {
theAntiColoured.push_back(p);
Ptr<MultiColour>::pointer colour =
dynamic_ptr_cast<Ptr<MultiColour>::pointer>
(p->colourInfo());
colour->antiColourLine(this, index);
}
void ColourLine::addColouredIndexed(tPPtr p, int index, bool anti) {
if ( anti ) addAntiColouredIndexed(p, index);
else {
theColoured.push_back(p);
Ptr<MultiColour>::pointer colour =
dynamic_ptr_cast<Ptr<MultiColour>::pointer>
(p->colourInfo());
colour->colourLine(this, index);
}
}
void ColourLine::addAntiColoured(tPPtr p) {
theAntiColoured.push_back(p);
p->colourInfo()->antiColourLine(this);
}
void ColourLine::addColoured(tPPtr p, bool anti) {
if ( anti ) addAntiColoured(p);
else {
theColoured.push_back(p);
p->colourInfo()->colourLine(this);
}
}
void ColourLine::removeAntiColoured(tPPtr p) {
tPVector::iterator cp=find(range(theAntiColoured), p);
if(cp!=theAntiColoured.end()) theAntiColoured.erase(cp);
p->colourInfo()->removeAntiColourLine(this);
}
void ColourLine::removeColoured(tPPtr p, bool anti) {
if ( anti ) removeAntiColoured(p);
else {
tPVector::iterator cp=find(range(theColoured), p);
if(cp!=theColoured.end()) theColoured.erase(cp);
p->colourInfo()->removeColourLine(this);
}
}
bool ColourLine::join(ColinePtr line) {
if ( !startParticle() || startParticle() != line->endParticle() )
return false;
while ( line->coloured().size() ) {
tPPtr p = line->coloured()[0];
line->removeColoured(p);
theColoured.insert(theColoured.begin(), p);
p->colourInfo()->colourLine(this);
}
while ( line->antiColoured().size() ) {
tPPtr p = line->antiColoured()[0];
line->removeAntiColoured(p);
theAntiColoured.push_back(p);
p->colourInfo()->antiColourLine(this);
}
return true;
}
void ColourLine::write(ostream & os, tcEventPtr event, bool anti) const {
os << ( anti? '-': '+' );
int index = event->colourLineIndex(this);
if ( sourceNeighbours().first && sourceNeighbours().second )
os << '(' << event->colourLineIndex(sourceNeighbours().first)
<< '*' << event->colourLineIndex(sourceNeighbours().second) << ')';
os << index;
if ( sinkNeighbours().first && sinkNeighbours().second )
os << '(' << event->colourLineIndex(sinkNeighbours().first)
<< '.' << event->colourLineIndex(sinkNeighbours().second) << ')';
}
void ColourLine::persistentOutput(PersistentOStream & os) const {
os << theColoured << theAntiColoured << theSourceNeighbours
<< theSinkNeighbours << orphanedConnectors;
}
void ColourLine::persistentInput(PersistentIStream & is, int) {
is >> theColoured >> theAntiColoured >> theSourceNeighbours
>> theSinkNeighbours >> orphanedConnectors;
}
ClassDescription<ColourLine> ColourLine::initColourLine;
// Definition of the static class description member.
|
