| 26 |
|
~Particle(){ |
| 27 |
|
}; |
| 28 |
|
|
| 29 |
< |
LorentzVector v4(){ |
| 29 |
> |
LorentzVector v4() const{ |
| 30 |
|
LorentzVector v4; |
| 31 |
|
v4.SetPt(m_pt); |
| 32 |
|
v4.SetEta(m_eta); |
| 35 |
|
return v4; |
| 36 |
|
}; |
| 37 |
|
|
| 38 |
< |
float charge(){return m_charge;} |
| 39 |
< |
float pt(){return m_pt;} |
| 40 |
< |
float eta(){return m_eta;} |
| 41 |
< |
float phi(){return m_phi;} |
| 42 |
< |
float energy(){return m_energy;} |
| 38 |
> |
float charge() const{return m_charge;} |
| 39 |
> |
float pt() const {return m_pt;} |
| 40 |
> |
float eta() const{return m_eta;} |
| 41 |
> |
float phi() const{return m_phi;} |
| 42 |
> |
float energy() const{return m_energy;} |
| 43 |
|
|
| 44 |
|
void set_charge(float charge){m_charge=charge;} |
| 45 |
|
void set_pt(float pt){m_pt=pt;} |
| 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 |
|
} |
