]>
Commit | Line | Data |
---|---|---|
5d8eaff6 | 1 | #ifndef AliAODPID_H |
2 | #define AliAODPID_H | |
1c28d174 | 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 | //------------------------------------------------------------------------- | |
9 | // AOD Pid object for additional pid information | |
d7bdc804 | 10 | // Author: Annalisa Mastroserio, CERN |
1c28d174 | 11 | //------------------------------------------------------------------------- |
12 | ||
13 | #include <TObject.h> | |
14 | ||
af48efa8 | 15 | class AliTPCdEdxInfo; |
16 | ||
1c28d174 | 17 | class AliAODPid : public TObject { |
18 | ||
19 | public: | |
1c28d174 | 20 | AliAODPid(); |
21 | virtual ~AliAODPid(); | |
22 | AliAODPid(const AliAODPid& pid); | |
23 | AliAODPid& operator=(const AliAODPid& pid); | |
d7bdc804 | 24 | |
00a38d07 | 25 | enum{kTRDnPlanes=6}; |
d7bdc804 | 26 | |
373fc041 | 27 | //setters |
28 | void SetITSsignal(Double_t its) {fITSsignal=its;} | |
6d354ed1 | 29 | void SetITSdEdxSamples(const Double_t s[4]); |
373fc041 | 30 | void SetTPCsignal(Double_t tpc) {fTPCsignal=tpc;} |
c997f0f9 | 31 | void SetTPCsignalN(UShort_t tpcN) {fTPCsignalN=(UChar_t)((tpcN<160)?tpcN:160);} |
febefa51 | 32 | void SetTPCmomentum(Double_t tpcMom) {fTPCmomentum=tpcMom;} |
c997f0f9 | 33 | void SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo); |
6c52f2f2 | 34 | void SetTRDsignal(Int_t nslices, const Double_t * const trdslices); |
febefa51 | 35 | void SetTRDmomentum(Int_t nplane, Float_t trdMom) {fTRDmomentum[nplane]=trdMom;} |
6c52f2f2 | 36 | void SetTRDncls(UChar_t ncls, Int_t layer = -1); |
59a8e853 | 37 | void SetTRDntrackletsPID(UChar_t ntls) {fTRDntls = ntls;} |
820214a7 | 38 | void SetTRDChi2(Double_t chi2) {fTRDChi2 = chi2;} |
373fc041 | 39 | void SetTOFsignal(Double_t tof) {fTOFesdsignal=tof;} |
b13bcc35 | 40 | void SetTOFpidResolution(Double_t tofPIDres[5]); |
373fc041 | 41 | void SetIntegratedTimes(Double_t timeint[5]); |
d7bdc804 | 42 | |
febefa51 | 43 | Double_t GetITSsignal() const {return fITSsignal;} |
6d354ed1 | 44 | void GetITSdEdxSamples(Double_t *s) const; |
45 | Double_t GetITSdEdxSample(Int_t i) const { | |
46 | if(i>=0 && i<4) return fITSdEdxSamples[i]; | |
47 | else return 0.; | |
48 | } | |
febefa51 | 49 | Double_t GetTPCsignal() const {return fTPCsignal;} |
c997f0f9 | 50 | UShort_t GetTPCsignalN() const {return (UShort_t)fTPCsignalN;} |
af48efa8 | 51 | AliTPCdEdxInfo * GetTPCdEdxInfo()const{return fTPCdEdxInfo;} |
52 | ||
febefa51 | 53 | Double_t GetTPCmomentum() const {return fTPCmomentum;} |
de8b2440 | 54 | Int_t GetTRDnSlices() const {return fTRDnSlices/6;} |
febefa51 | 55 | Double_t* GetTRDsignal() const {return fTRDslices;} |
820214a7 | 56 | Double_t GetTRDChi2() const {return fTRDChi2;} |
c997f0f9 | 57 | const Double_t* GetTRDmomentum() const {return fTRDmomentum;} |
59a8e853 | 58 | UChar_t GetTRDncls(UChar_t layer) const { if(layer > 5) return 0; return fTRDncls[layer];} |
6c52f2f2 | 59 | UChar_t GetTRDncls() const; |
59a8e853 | 60 | UChar_t GetTRDntrackletsPID() const {return fTRDntls;} |
febefa51 | 61 | Double_t GetTOFsignal() const {return fTOFesdsignal;} |
febefa51 | 62 | |
63 | void GetIntegratedTimes(Double_t timeint[5]) const; | |
b13bcc35 | 64 | void GetTOFpidResolution (Double_t tofRes[5]) const; |
1c28d174 | 65 | |
66 | private : | |
febefa51 | 67 | Double32_t fITSsignal; //[0.,0.,10] detector raw signal |
6d354ed1 | 68 | Double32_t fITSdEdxSamples[4];//[0.,0.,10] ITS dE/dx samples |
c997f0f9 | 69 | |
febefa51 | 70 | Double32_t fTPCsignal; //[0.,0.,10] detector raw signal |
c997f0f9 | 71 | UChar_t fTPCsignalN; // number of points used for TPC dE/dx |
72 | Double32_t fTPCmomentum; //[0.,0.,20] momentum at the inner wall of TPC; | |
73 | ||
de8b2440 | 74 | Int_t fTRDnSlices; // N slices used for PID in the TRD (as number of slices per tracklet * number of layers) |
59a8e853 | 75 | UChar_t fTRDntls; // number of tracklets used for PID calculation |
76 | UChar_t fTRDncls[6]; // number of clusters used for dE/dx calculation | |
c997f0f9 | 77 | Double32_t* fTRDslices; //[fTRDnSlices][0.,0.,10] |
78 | Double32_t fTRDmomentum[6]; //[0.,0.,10] momentum at the TRD layers | |
7b7338b4 | 79 | Double32_t fTRDChi2; //[0.,0.,10] TRD chi2 |
c997f0f9 | 80 | |
81 | Double32_t fTOFesdsignal; //[0.,0.,20] TOF signal - t0 (T0 interaction time) | |
82 | Double32_t fTOFpidResolution[5]; //[0.,0.,20] TOF pid resolution for each mass hypotesys | |
83 | Double32_t fIntTime[5]; //[0.,0.,20] track time hypothesis | |
996918c5 | 84 | |
771ec43e | 85 | AliTPCdEdxInfo * fTPCdEdxInfo; // object containing dE/dx information for different pad regions |
febefa51 | 86 | |
7b7338b4 | 87 | ClassDef(AliAODPid, 13); |
1c28d174 | 88 | }; |
89 | ||
59a8e853 | 90 | //_____________________________________________________________ |
6c52f2f2 | 91 | inline void AliAODPid::SetTRDsignal(Int_t nslices, const Double_t * const trdslices) { |
59a8e853 | 92 | // |
93 | // Set TRD dE/dx slices and the number of dE/dx slices per track | |
94 | // | |
95 | if(fTRDslices && fTRDnSlices != nslices) { | |
96 | delete [] fTRDslices; fTRDslices = NULL; | |
97 | }; | |
98 | if(!fTRDslices) fTRDslices = new Double32_t[nslices]; | |
99 | fTRDnSlices = nslices; | |
100 | for(Int_t is = 0; is < fTRDnSlices; is++) fTRDslices[is] = trdslices[is]; | |
101 | } | |
102 | ||
103 | //_____________________________________________________________ | |
6c52f2f2 | 104 | inline void AliAODPid::SetTRDncls(UChar_t ncls, Int_t layer) { |
59a8e853 | 105 | // |
106 | // Set the number of clusters / tracklet | |
107 | // If no layer is specified the full number of clusters will be put in layer 0 | |
108 | // | |
109 | if(layer > 5) return; | |
110 | if(layer < 0) fTRDncls[0] = ncls; | |
111 | else fTRDncls[layer] = ncls; | |
112 | } | |
113 | ||
114 | //_____________________________________________________________ | |
6c52f2f2 | 115 | inline UChar_t AliAODPid::GetTRDncls() const { |
59a8e853 | 116 | // |
117 | // Get number of clusters per track | |
118 | // Calculated as sum of the number of clusters per tracklet | |
119 | // | |
120 | UChar_t ncls = 0; | |
121 | for(Int_t ily = 0; ily < 6; ily++) ncls += fTRDncls[ily]; | |
122 | return ncls; | |
123 | } | |
1c28d174 | 124 | #endif |