]>
Commit | Line | Data |
---|---|---|
df9db588 | 1 | #ifndef AliAODJet_H |
2 | #define AliAODJet_H | |
3 | /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. * | |
4 | * See cxx source for full Copyright notice */ | |
5 | ||
6 | /* $Id$ */ | |
7 | ||
8 | //------------------------------------------------------------------------- | |
87102d3c | 9 | // AOD jet class |
10 | // Author: Andreas Morsch, CERN | |
df9db588 | 11 | //------------------------------------------------------------------------- |
12 | ||
87102d3c | 13 | #include <TLorentzVector.h> |
ff7c57dd | 14 | #include "AliVParticle.h" |
2e7293a2 | 15 | #include <TArrayI.h> |
df9db588 | 16 | #include "AliAODVertex.h" |
17 | ||
87102d3c | 18 | |
ff7c57dd | 19 | class AliAODJet : public AliVParticle { |
df9db588 | 20 | |
21 | public: | |
87102d3c | 22 | AliAODJet(); |
23 | AliAODJet(Double_t px, Double_t py, Double_t pz, Double_t e); | |
a1708071 | 24 | AliAODJet(TLorentzVector & p); |
87102d3c | 25 | virtual ~AliAODJet(); |
26 | AliAODJet(const AliAODJet& jet); | |
27 | AliAODJet& operator=(const AliAODJet& jet); | |
ff7c57dd | 28 | // AliVParticle methods |
87102d3c | 29 | virtual Double_t Px() const { return fMomentum->Px(); } |
30 | virtual Double_t Py() const { return fMomentum->Py(); } | |
31 | virtual Double_t Pz() const { return fMomentum->Pz(); } | |
32 | virtual Double_t Pt() const { return fMomentum->Pt(); } | |
33 | virtual Double_t P() const { return fMomentum->P(); } | |
34 | virtual Double_t OneOverPt() const { return 1. / fMomentum->Pt(); } | |
c683ddc2 | 35 | virtual Bool_t PxPyPz(Double_t p[3]) const { p[0] = Px(); p[1] = Py(); p[2] = Pz(); return kTRUE; } |
4124c6c5 | 36 | virtual Double_t Phi() const; |
87102d3c | 37 | virtual Double_t Theta() const { return fMomentum->Theta(); } |
38 | virtual Double_t E() const { return fMomentum->E(); } | |
39 | virtual Double_t M() const { return fMomentum->M(); } | |
40 | virtual Double_t Eta() const { return fMomentum->Eta(); } | |
41 | virtual Double_t Y() const { return fMomentum->Rapidity();} | |
c683ddc2 | 42 | virtual Double_t Xv() const {return -999.;} // put reasonable values here |
43 | virtual Double_t Yv() const {return -999.;} // | |
44 | virtual Double_t Zv() const {return -999.;} // | |
45 | virtual Bool_t XvYvZv(Double_t x[3]) const { x[0] = Xv(); x[1] = Yv(); x[2] = Zv(); return kTRUE; } | |
a5817b6b | 46 | virtual Bool_t IsTriggeredEMCAL(){return (fTrigger&kEMCALTriggered)==kEMCALTriggered;} |
47 | virtual Bool_t IsTriggeredTRD(){return (fTrigger&kTRDTriggered)==kTRDTriggered;} | |
48 | virtual UChar_t Trigger(){return fTrigger;} | |
4124c6c5 | 49 | |
086d6981 | 50 | virtual void AddTrack(TObject *tr); |
51 | ||
2e7293a2 | 52 | TObject* GetTrack(Int_t i) {return fRefTracks->At(i);} |
87102d3c | 53 | virtual void SetBgEnergy(Double_t bgEnCh, Double_t bgEnNe) |
54 | {fBackgEnergy[0] = bgEnCh; fBackgEnergy[1] = bgEnNe;} | |
a5817b6b | 55 | virtual void SetEffArea(Double_t effACh, Double_t effANe, Double_t effAErrCh = 0, Double_t effAErrNe = 0) |
56 | { | |
57 | fEffectiveArea[0] = effACh; fEffectiveArea[1] = effANe; | |
58 | fEffectiveAreaError[0] = effAErrCh; | |
59 | fEffectiveAreaError[1] = effAErrNe; | |
60 | } | |
944efc80 | 61 | virtual void SetPxPyPzE(Double_t px, Double_t py, Double_t pz, Double_t e); |
a5817b6b | 62 | virtual void SetTrigger(UChar_t f){fTrigger |= f;} |
63 | virtual void ResetTrigger(UChar_t f){fTrigger &= ~f;} | |
944efc80 | 64 | |
2e7293a2 | 65 | virtual TRefArray* GetRefTracks() const { return fRefTracks;} |
66 | virtual Double_t ChargedBgEnergy() const { return fBackgEnergy[0];} | |
67 | virtual Double_t NeutralBgEnergy() const { return fBackgEnergy[1];} | |
68 | virtual Double_t TotalBgEnergy() const { return (fBackgEnergy[0] + fBackgEnergy[1]);} | |
87102d3c | 69 | |
2e7293a2 | 70 | virtual Double_t EffectiveAreaCharged() const { return fEffectiveArea[0];} |
71 | virtual Double_t EffectiveAreaNeutral() const { return fEffectiveArea[1];} | |
a5817b6b | 72 | virtual Double_t ErrorEffectiveAreaCharged() const { return fEffectiveAreaError[0];} |
73 | virtual Double_t ErrorEffectiveAreaNeutral() const { return fEffectiveAreaError[1];} | |
a19f293d | 74 | virtual Double_t DeltaR(const AliVParticle* part); |
87102d3c | 75 | |
a19f293d | 76 | |
66ac9c97 | 77 | TLorentzVector* MomentumVector() const {return fMomentum;} |
a19f293d | 78 | virtual void Print(Option_t* /*option*/) const; |
87102d3c | 79 | |
80 | // Dummy | |
81 | virtual Short_t Charge() const { return 0;} | |
82 | virtual const Double_t* PID() const { return NULL;} | |
3c43fb2b | 83 | virtual Int_t GetLabel() const { return -1;} |
87102d3c | 84 | // |
3a4dcedf | 85 | |
a5817b6b | 86 | |
87 | ||
88 | ||
3a4dcedf | 89 | /** Compare this class with an other instance of this class |
90 | * used in a TClonesArray::Sort() | |
91 | * @param obj ptr to other instance | |
92 | * @return Returns 0 when equal, 1 when this is smaller | |
93 | * and -1 when bigger -- sorts descending | |
94 | */ | |
95 | Int_t Compare( const TObject* obj) const; | |
87102d3c | 96 | |
97 | ||
3a4dcedf | 98 | /** Defines this class as being sortable in a TClonesArray |
99 | * @return always kTRUE; | |
100 | */ | |
101 | Bool_t IsSortable() const { return kTRUE; } | |
102 | ||
a5817b6b | 103 | // first only one bit for EMCAL and TRD, leave space for more |
104 | // trigger types and/or other detectors | |
105 | enum {kEMCALTriggered = 1<<0, | |
106 | kTRDTriggered = 4<<0}; | |
107 | ||
9333290e | 108 | |
a5817b6b | 109 | private: |
086d6981 | 110 | Double32_t fBackgEnergy[2]; // Subtracted background energy |
111 | Double32_t fEffectiveArea[2]; // Effective jet area used for background subtraction | |
a5817b6b | 112 | Double32_t fEffectiveAreaError[2]; //[0,1,10] relative error of jet areas, 10 bit precision |
086d6981 | 113 | Double32_t fNeutralFraction; //[0,1,12] Neutral fraction between 0 and 1 12 bit precision; |
114 | UChar_t fTrigger; // Bit mask to flag jets triggered by a certain detector | |
115 | TLorentzVector* fMomentum; // Jet 4-momentum vector | |
116 | TRefArray* fRefTracks; // array of references to the tracks belonging to the jet | |
9333290e | 117 | |
a5817b6b | 118 | ClassDef(AliAODJet,5); |
4124c6c5 | 119 | |
df9db588 | 120 | }; |
121 | ||
4124c6c5 | 122 | inline Double_t AliAODJet::Phi() const |
123 | { | |
124 | // Return phi | |
125 | Double_t phi = fMomentum->Phi(); | |
126 | if (phi < 0.) phi += 2. * TMath::Pi(); | |
127 | return phi; | |
128 | } | |
129 | ||
df9db588 | 130 | #endif |