26 |
|
~Particle(){ |
27 |
|
}; |
28 |
|
|
29 |
< |
LorentzVector v4(){ |
29 |
> |
LorentzVector v4() const{ |
30 |
|
LorentzVector v4; |
31 |
|
v4.SetPt(m_pt); |
32 |
|
v4.SetEta(m_eta); |
47 |
|
void set_phi(float phi){m_phi=phi;} |
48 |
|
void set_energy(float energy){m_energy=energy;} |
49 |
|
|
50 |
< |
double deltaPhi(Particle p2){ |
50 |
> |
void set_v4(LorentzVector v4){ |
51 |
> |
set_pt(v4.Pt()); |
52 |
> |
set_eta(v4.Eta()); |
53 |
> |
set_phi(v4.Phi()); |
54 |
> |
set_energy(v4.E()); |
55 |
> |
} |
56 |
> |
|
57 |
> |
double deltaPhi(const Particle p2) const{ |
58 |
|
double deltaphi = fabs(this->phi() - p2.phi()); |
59 |
|
if(deltaphi > PI) deltaphi = 2* PI - deltaphi; |
60 |
|
return deltaphi; |
61 |
|
} |
62 |
< |
double deltaR(Particle p2){ |
62 |
> |
double deltaR(const Particle p2) const{ |
63 |
|
double deltaeta = this->eta() - p2.eta(); |
64 |
|
return sqrt(deltaeta*deltaeta+deltaPhi(p2)*deltaPhi(p2)); |
65 |
|
} |