Newer
Older
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
#include "Event/Event.hh"
#include <vector>
Event::Event()
{}
Event::~Event()
{
std::vector<EvtPartData*>::iterator iter;
for (iter = particles.begin(); iter != particles.end(); ++iter)
delete *iter;
}
Vector4<float>* Event::p4(int i)
{
if (particles.size() > i)
return particles[i]->vector4;
else
ErrMsg(fatal) << "accessing " << i << "th of "
<< particles.size() << " particles" << endmsg;
return 0;
}
float* Event::pid(int i)
{
if (particles.size() > i)
return particles[i]->pidVector;
else
ErrMsg(fatal) << "accessing pid of " << i << "th of "
<< particles.size() << " particles" << endmsg;
}
void Event::addParticle(double e, double px, double py, double pz)
{
EvtPartData* evtData = new EvtPartData(e,px,py,pz);
particles.push_back(evtData);
return;
}
int Event::size()
{
return particles.size();
}