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
|
// -*- C++ -*-
//
// DiagramBase.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 DiagramBase class.
//
#include "DiagramBase.h"
#include "ThePEG/Utilities/ClassDescription.h"
#include "ThePEG/Utilities/DescriptionList.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
using namespace ThePEG;
DiagramBase::~DiagramBase() {}
struct ParticleOrdering {
bool operator()(tcPDPtr p1, tcPDPtr p2) {
return abs(p1->id()) > abs(p2->id()) ||
( abs(p1->id()) == abs(p2->id()) && p1->id() > p2->id() ) ||
( p1->id() == p2->id() && p1->fullName() > p2->fullName() );
}
};
string DiagramBase::getTag() const {
string tag;
if ( !done() ) throw DiagramBaseSetupException(*this);
for ( int i = 0; i < nIncoming(); ++i ) {
if ( i ) tag += ",";
tag += partons()[i]->PDGName();
}
tag += "->";
multiset<tcPDPtr,ParticleOrdering> out;
for ( int i = nIncoming(), N = partons().size(); i < N; ++i )
out.insert(partons()[i]);
for ( multiset<tcPDPtr,ParticleOrdering>::iterator i = out.begin();
i != out.end(); ++i ) {
if ( i != out.begin() ) tag += ",";
tag += (**i).PDGName();
}
return tag;
}
void DiagramBase::persistentOutput(PersistentOStream & os) const {
os << theNIncoming << thePartons << theId;
}
void DiagramBase::persistentInput(PersistentIStream & is, int) {
is >> theNIncoming >> thePartons >> theId;
}
AbstractClassDescription<DiagramBase> DiagramBase::initDiagramBase;
// Definition of the static class description member.
void DiagramBase::Init() {}
DiagramBaseSetupException::DiagramBaseSetupException(const DiagramBase & db) {
const ClassDescriptionBase * cd = DescriptionList::find(typeid(db));
if ( !cd )
theMessage << "Tried to use an unknown sub class of DiagramBase.";
else
theMessage << "The '" << cd->name()
<< "' sub class did not setup the DiagramBase class correctly.";
severity(abortnow);
}
|
