1 #ifndef ALIAODHMPIDRINGS_H
2 #define ALIAODHMPIDRINGS_H
6 // Class to handle the AOD tracks with good HMPID rings
7 // Author: Levente Molnar
8 // levente.molnar@cern.ch , March 2012
19 class AliAODHMPIDrings : public TObject {
41 AliAODHMPIDrings(const AliAODHMPIDrings& hmpidAOD);//
42 AliAODHMPIDrings &operator=(const AliAODHMPIDrings& hmpidAOD);//
43 virtual ~AliAODHMPIDrings() {};
46 Int_t GetHmpTrkID() const { return fHmpidAODtrkId; }
48 Double32_t GetHmpMipCharge() const { return fHmpidAODqn%1000000; }
49 Double32_t GetHmpNumOfPhotonClusters() const { return fHmpidAODqn/1000000;}
51 Int_t GetHmpChamber() const { return fHmpidAODcluIdx/1000000; }
53 Int_t GetHmpCluIdx() const { return fHmpidAODcluIdx; }
55 Double_t GetHmpTrackTheta() const { return fHmpidAODtrkTheta;}
56 Double_t GetHmpTrackPhi() const { return fHmpidAODtrkPhi;}
58 Double_t GetHmpSignal() const { return fHmpidAODsignal;}
59 Double_t GetHmpOccupancy() const { return fHmpidAODocc;}
61 Double_t GetHmpChi2() const { return fHmpidAODchi2;}
63 Double_t GetHmpTrackX() const { return fHmpidAODtrkX;}
64 Double_t GetHmpTrackY() const { return fHmpidAODtrkY;}
66 Double_t GetHmpMipX() const { return fHmpidAODmipX;}
67 Double_t GetHmpMipY() const { return fHmpidAODmipY;}
69 Double_t GetHmpDX() const { return fHmpidAODmipX - fHmpidAODtrkX;}
70 Double_t GetHmpDY() const { return fHmpidAODmipY - fHmpidAODtrkY;}
71 Double_t GetHmpDist() const { return TMath::Sqrt((fHmpidAODmipX - fHmpidAODtrkX)*(fHmpidAODmipX - fHmpidAODtrkX) + (fHmpidAODmipY - fHmpidAODtrkY)*(fHmpidAODmipY - fHmpidAODtrkY));}
74 void GetHmpPidProbs(Double32_t *pid) const; //defined in cxx
75 void GetHmpMom(Double32_t *mom) const; //defined in cxx
79 void SetHmpMipCharge(Int_t q) { fHmpidAODqn = q; }
80 void SetHmpCluIdx(Int_t ch,Int_t idx) { fHmpidAODcluIdx=ch*1000000+idx;}
82 void SetHmpNumOfPhotonClusters(Int_t nph) { fHmpidAODqn = 1000000 * nph;}
84 void SetHmpTrackTheta(Double_t trkTheta) { fHmpidAODtrkTheta = trkTheta;}
85 void SetHmpTrackPhi(Double_t trkPhi) { fHmpidAODtrkPhi = trkPhi;}
87 void SetHmpSignal(Double_t thetaC) { fHmpidAODsignal = thetaC;}
88 void SetHmpOccupancy(Double_t occ) { fHmpidAODocc = occ;}
90 void SetHmpChi2(Double_t chi2) { fHmpidAODchi2 = chi2;}
92 void SetHmpTrackX(Double_t trkX) { fHmpidAODtrkX = trkX;}
93 void SetHmpTrackY(Double_t trkY) { fHmpidAODtrkY = trkY;}
95 void SetHmpMipX(Double_t mipX) { fHmpidAODmipX = mipX;}
96 void SetHmpMipY(Double_t mipY) { fHmpidAODmipY = mipY;}
98 void SetHmpPidProbs(Double_t *pid);
99 void SetHmpMom(Double_t *mom);
105 Int_t fHmpidAODtrkId; // Unique track id as in ESD
106 Int_t fHmpidAODqn; // 1000000*number of photon clusters + QDC
107 Int_t fHmpidAODcluIdx; // 1000000*chamber id + cluster idx of the assigned MIP cluster
109 Double32_t fHmpidAODtrkTheta; // [-2*pi,2*pi,16] theta of the track extrapolated to the HMPID, LORS
110 Double32_t fHmpidAODtrkPhi; // [-2*pi,2*pi,16] theta of the track extrapolated to the HMPID, LORS
111 Double32_t fHmpidAODsignal; // HMPID signal (Theta ckov, rad)
112 Double32_t fHmpidAODocc; // chamber occupancy where the track passed through: number of pads
113 Double32_t fHmpidAODchi2; // [0.,0.,8] chi2 in the HMPID
114 Double32_t fHmpidAODtrkX; // x of the track impact, LORS
115 Double32_t fHmpidAODtrkY; // y of the track impact, LORS
116 Double32_t fHmpidAODmipX; // x of the MIP in LORS
117 Double32_t fHmpidAODmipY; // y of the MIP in LORS
118 Double32_t fHmpidAODpid[AliPID::kSPECIES]; // [0.,0.,8] "detector response probabilities" (for the PID)
119 Double32_t fHMPIDmom[3]; // track momentum at the HMPID ring reconstruction
121 ClassDef(AliAODHMPIDrings,2)