]>
Commit | Line | Data |
---|---|---|
7df864a3 | 1 | // $Id$ |
7efbea04 | 2 | |
7df864a3 | 3 | #include "AliESDJet.h" |
4 | ||
7efbea04 | 5 | ClassImp(AliESDJet) |
6 | ||
7 | //__________________________________________________________________________________________________ | |
8 | AliESDJet::AliESDJet(Double_t px, Double_t py, Double_t pz) | |
9 | : AliVParticle(), | |
10 | fPt(TMath::Sqrt(px*px+py*py)), | |
11 | fEta(TMath::ASinH(pz/fPt)), | |
12 | fPhi(0), fM(0), fNEF(0), fArea(0), fNch(0), fNn(0), | |
13 | fMaxCPt(0), fMaxNPt(0) | |
14 | { | |
15 | // Constructor. | |
16 | ||
17 | if (fPt != 0) { | |
18 | fPhi = TMath::ATan2(py, px); | |
19 | if (fPhi<0.) | |
20 | fPhi += 2. * TMath::Pi(); | |
21 | } | |
22 | } | |
23 | ||
24 | //_________________________________________________________________________________________________ | |
25 | AliESDJet::AliESDJet(const AliESDJet &jet) : | |
26 | AliVParticle(jet), | |
27 | fPt(jet.fPt), fEta(jet.fEta), fPhi(jet.fPhi), | |
28 | fM(jet.fM), fNEF(jet.fNEF), fArea(jet.fArea), | |
29 | fNch(jet.fNch), fNn(jet.fNn), | |
30 | fMaxCPt(jet.fMaxCPt), fMaxNPt(jet.fMaxNPt) | |
7df864a3 | 31 | { |
7efbea04 | 32 | // Constructor. |
33 | } | |
34 | ||
35 | //_________________________________________________________________________________________________ | |
36 | AliESDJet &AliESDJet::operator=(const AliESDJet &jet) | |
37 | { | |
38 | // Assignment operator. | |
39 | ||
40 | if (this!=&jet) { | |
41 | AliVParticle::operator=(jet); | |
42 | fPt = jet.fPt; | |
43 | fEta = jet.fEta; | |
44 | fPhi = jet.fPhi; | |
45 | fM = jet.fM; | |
46 | fNEF = jet.fNEF; | |
47 | fArea = jet.fArea; | |
48 | fNch = jet.fNch; | |
49 | fNn = jet.fNn; | |
50 | fMaxCPt = jet.fMaxCPt; | |
51 | fMaxNPt = jet.fMaxNPt; | |
52 | } | |
53 | ||
54 | return *this; | |
7df864a3 | 55 | } |
7efbea04 | 56 | |
57 | //__________________________________________________________________________________________________ | |
58 | void AliESDJet::GetMom(TLorentzVector &vec) const | |
59 | { | |
60 | // Return momentum as four vector. | |
61 | ||
62 | Double_t p = fPt *TMath::CosH(fEta); | |
63 | vec.SetPtEtaPhiE(fPt,fEta,fPhi,TMath::Sqrt(p*p+fM*fM)); | |
64 | } | |
65 | ||
66 | //__________________________________________________________________________________________________ | |
67 | void AliESDJet::Print(Option_t* /*option*/) const | |
68 | { | |
69 | // Print jet information. | |
70 | ||
71 | printf("Jet pt=%.2f, eta=%.2f, phi=%.2f, area=%.2f, NEF=%.2f\n", fPt, fEta, fPhi, fArea, fNEF); | |
72 | } | |
73 |