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
160
161
162
163
164
165
166
167
168
169
170
171
172
173
174
175
176
177
178
179
180
181
|
// -*- C++ -*-
//
// Strategy.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 Strategy class.
//
#include "Strategy.h"
#include "ThePEG/Persistency/PersistentOStream.h"
#include "ThePEG/Persistency/PersistentIStream.h"
#include "ThePEG/PDT/ParticleData.h"
#include "ThePEG/Interface/RefVector.h"
#include "ThePEG/Interface/Parameter.h"
#include "ThePEG/Interface/ParVector.h"
#include "ThePEG/Interface/ClassDocumentation.h"
#include "ThePEG/Repository/Repository.h"
#include "ThePEG/Utilities/Throw.h"
using namespace ThePEG;
string Strategy::localParticlesDir() const {
return theLocalParticlesDir;
}
bool Strategy::checkDir(string x) {
if ( x == "" ) return false;
if ( x[0] != '/' )
Throw<InterfaceException>()
<< "Directory name must start with a '/'."
<< Exception::setuperror;
if ( x[x.length() - 1] != '/' ) x += '/';
try {
Repository::CheckObjectDirectory(x);
}
catch ( Exception & e ) {
e.handle();
Throw<InterfaceException>()
<< "No directory named " << x << " in the repository."
<< Exception::setuperror;
}
return true;
}
void Strategy::setDefaultParticlesDirs(string dir, int i) {
if ( i < 0 || unsigned(i) >= theDefaultParticlesDirs.size() )
Throw<InterfaceException>()
<< "Index out of range in DefaultParticlesDirs"
<< Exception::setuperror;
if ( !checkDir(dir) )
Throw<InterfaceException>()
<< "Empty directory name not allowd in DefaultParticlesDirs"
<< Exception::setuperror;
theDefaultParticlesDirs[i] = dir;
}
void Strategy::insDefaultParticlesDirs(string dir, int i) {
if ( i < 0 || unsigned(i) > theDefaultParticlesDirs.size() )
Throw<InterfaceException>()
<< "Index out of range in DefaultParticlesDirs"
<< Exception::setuperror;
if ( !checkDir(dir) )
Throw<InterfaceException>()
<< "Empty directory name not allowd in DefaultParticlesDirs"
<< Exception::setuperror;
theDefaultParticlesDirs.insert(theDefaultParticlesDirs.begin() + i, dir);
}
void Strategy::setLocalParticlesDir(string x) {
if ( checkDir(x) ) theLocalParticlesDir = x;
}
IBPtr Strategy::clone() const {
return new_ptr(*this);
}
IBPtr Strategy::fullclone() const {
return new_ptr(*this);
}
void Strategy::persistentOutput(PersistentOStream & os) const {
os << theParticles << theLocalParticlesDir << theDefaultObjects
<< theDefaultParticlesDirs;
}
void Strategy::persistentInput(PersistentIStream & is, int) {
is >> theParticles >> theLocalParticlesDir >> theDefaultObjects
>> theDefaultParticlesDirs;
}
ClassDescription<Strategy> Strategy::initStrategy;
void Strategy::setLocalParticles(PDPtr pd, int) {
particles()[pd->id()] = pd;
}
void Strategy::insLocalParticles(PDPtr pd, int) {
particles()[pd->id()] = pd;
}
void Strategy::delLocalParticles(int place) {
ParticleMap::iterator it = particles().begin();
while ( place-- && it != particles().end() ) ++it;
if ( it != particles().end() ) particles().erase(it);
}
vector<PDPtr> Strategy::getLocalParticles() const {
vector<PDPtr> ret;
for ( ParticleMap::const_iterator it = particles().begin();
it != particles().end(); ++it ) ret.push_back(it->second);
return ret;
}
void Strategy::Init() {
static ClassDocumentation<Strategy> documentation
("Represents a general strategy to be assigned to an EventGenerator. "
"It contains a set of default ParticleData objects which takes "
"presedence over the ones in the Repository (although not over "
"the ones in the EventGenerator). It also contains a set of other "
"default objects which are automatically assigned to all Reference "
"and RefVector interfaces which have the "
"InterfaceBase::defaultIfNull() flag set.");
static RefVector<Strategy,ParticleData> interfaceLocalParticles
("LocalParticles",
"Special versions of ThePEG::ParticleData objects to be used. Note "
"that to delete an object, its number in the list "
"should be given, rather than its id number.",
0, 0, false, false, true, false,
&Strategy::setLocalParticles, &Strategy::insLocalParticles,
&Strategy::delLocalParticles, &Strategy::getLocalParticles);
static RefVector<Strategy,Interfaced> interfaceDefaultObjects
("DefaultObjects",
"A vector of pointers to default objects. In a ThePEG::Reference or "
"ThePEG::RefVector interface with the defaultIfNull() flag set, if a "
"null pointer is encountered this vector is gone through until an "
"acceptable object is found in which case the null pointer is replaced "
"by a pointer to this object. Note that the default objects given in the "
"ThePEG::EventGenerator are gone through first and are given precedence.",
&Strategy::theDefaultObjects, 0, true, false, true, false, false);
static Parameter<Strategy,string> interfaceLocalParticlesDir
("LocalParticlesDir",
"A directory in the repository which will be scanned for particles "
"which will be included as default particles in a run. These particles "
"will be overridden by particles specified in "
"<interface>LocalParticles</interface> and default particles "
"specified directly in the EventGenerator.",
&Strategy::theLocalParticlesDir, "",
true, false,
&Strategy::setLocalParticlesDir, (string(Strategy::*)()const)(0),
(string(Strategy::*)()const)(0));
static ParVector<Strategy,string> interfaceDefaultParticlesDirs
("DefaultParticlesDirs",
"By default all particles in the Repository are included in a run, "
"although only one particle object per PDG id number. If directories are "
"listed in DefaultParticlesDirs, only particles in these will be "
"considered for default inclusion in a run. Only particles which have "
"a PDG id which is not given by particles in "
"<interface>LocalParticlesDir</interface>, "
"<interface>LocalParticles</interface>, or in "
"<interface>EventGenerator::LocalParticles</interface> will be "
"considered.",
&Strategy::theDefaultParticlesDirs, -1, "", "", "",
true, false, Interface::nolimits,
&Strategy::setDefaultParticlesDirs, &Strategy::insDefaultParticlesDirs);
interfaceLocalParticles.rank(10);
interfaceLocalParticlesDir.rank(11);
interfaceDefaultObjects.rank(9);
interfaceDefaultParticlesDirs.rank(8);
}
|
