3 /* Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
4 * See cxx source for full Copyright notice */
8 //-------------------------------------------------------------------------
9 // AOD Pid object for additional pid information
10 // Author: Annalisa Mastroserio, CERN
11 //-------------------------------------------------------------------------
17 class AliAODPid : public TObject {
22 AliAODPid(const AliAODPid& pid);
23 AliAODPid& operator=(const AliAODPid& pid);
28 void SetITSsignal(Double_t its) {fITSsignal=its;}
29 void SetITSdEdxSamples(const Double_t s[4]);
30 void SetTPCsignal(Double_t tpc) {fTPCsignal=tpc;}
31 void SetTPCsignalN(UShort_t tpcN) {fTPCsignalN=(UChar_t)((tpcN<160)?tpcN:160);}
32 void SetTPCmomentum(Double_t tpcMom) {fTPCmomentum=tpcMom;}
33 void SetTPCdEdxInfo(AliTPCdEdxInfo * dEdxInfo);
34 void SetTRDslices(Int_t nslices, const Double_t * const trdslices);
35 void SetTRDsignal(Double_t sig) {fTRDsignal = sig;}
36 void SetTRDmomentum(Int_t nplane, Float_t trdMom) {fTRDmomentum[nplane]=trdMom;}
37 void SetTRDncls(UChar_t ncls, Int_t layer = -1);
38 void SetTRDntrackletsPID(UChar_t ntls) {fTRDntls = ntls;}
39 void SetTRDChi2(Double_t chi2) {fTRDChi2 = chi2;}
40 void SetTOFsignal(Double_t tof) {fTOFesdsignal=tof;}
41 void SetTOFpidResolution(Double_t tofPIDres[5]);
42 void SetIntegratedTimes(Double_t timeint[5]);
44 Double_t GetITSsignal() const {return fITSsignal;}
45 void GetITSdEdxSamples(Double_t *s) const;
46 Double_t GetITSdEdxSample(Int_t i) const {
47 if(i>=0 && i<4) return fITSdEdxSamples[i];
50 Double_t GetTPCsignal() const {return fTPCsignal;}
51 UShort_t GetTPCsignalN() const {return (UShort_t)fTPCsignalN;}
52 AliTPCdEdxInfo * GetTPCdEdxInfo()const{return fTPCdEdxInfo;}
54 Double_t GetTPCmomentum() const {return fTPCmomentum;}
55 Int_t GetTRDnSlices() const {return fTRDnSlices/6;}
56 Double_t GetTRDsignal() const {return fTRDsignal;}
57 Double_t* GetTRDslices() const {return fTRDslices;}
58 Double_t GetTRDChi2() const {return fTRDChi2;}
59 const Double_t* GetTRDmomentum() const {return fTRDmomentum;}
60 UChar_t GetTRDncls(UChar_t layer) const { if(layer > 5) return 0; return fTRDncls[layer];}
61 UChar_t GetTRDncls() const;
62 UChar_t GetTRDntrackletsPID() const {return fTRDntls;}
63 Double_t GetTOFsignal() const {return fTOFesdsignal;}
65 void GetIntegratedTimes(Double_t timeint[5]) const;
66 void GetTOFpidResolution (Double_t tofRes[5]) const;
69 Double32_t fITSsignal; //[0.,0.,10] detector raw signal
70 Double32_t fITSdEdxSamples[4];//[0.,0.,10] ITS dE/dx samples
72 Double32_t fTPCsignal; //[0.,0.,10] detector raw signal
73 UChar_t fTPCsignalN; // number of points used for TPC dE/dx
74 Double32_t fTPCmomentum; //[0.,0.,20] momentum at the inner wall of TPC;
76 Int_t fTRDnSlices; // N slices used for PID in the TRD (as number of slices per tracklet * number of layers)
77 UChar_t fTRDntls; // number of tracklets used for PID calculation
78 UChar_t fTRDncls[6]; // number of clusters used for dE/dx calculation
79 Double32_t* fTRDslices; //[fTRDnSlices][0.,0.,10]
80 Double32_t fTRDsignal; //[0.,0.,10] TRD signal
81 Double32_t fTRDmomentum[6]; //[0.,0.,10] momentum at the TRD layers
82 Double32_t fTRDChi2; //[0.,0.,10] TRD chi2
84 Double32_t fTOFesdsignal; //[0.,0.,20] TOF signal - t0 (T0 interaction time)
85 Double32_t fTOFpidResolution[5]; //[0.,0.,20] TOF pid resolution for each mass hypotesys
86 Double32_t fIntTime[5]; //[0.,0.,20] track time hypothesis
88 AliTPCdEdxInfo * fTPCdEdxInfo; // object containing dE/dx information for different pad regions
90 ClassDef(AliAODPid, 14);
93 //_____________________________________________________________
94 inline void AliAODPid::SetTRDslices(Int_t nslices, const Double_t * const trdslices) {
96 // Set TRD dE/dx slices and the number of dE/dx slices per track
98 if(fTRDslices && fTRDnSlices != nslices) {
99 delete [] fTRDslices; fTRDslices = NULL;
101 if(!fTRDslices) fTRDslices = new Double32_t[nslices];
102 fTRDnSlices = nslices;
103 for(Int_t is = 0; is < fTRDnSlices; is++) fTRDslices[is] = trdslices[is];
106 //_____________________________________________________________
107 inline void AliAODPid::SetTRDncls(UChar_t ncls, Int_t layer) {
109 // Set the number of clusters / tracklet
110 // If no layer is specified the full number of clusters will be put in layer 0
112 if(layer > 5) return;
113 if(layer < 0) fTRDncls[0] = ncls;
114 else fTRDncls[layer] = ncls;
117 //_____________________________________________________________
118 inline UChar_t AliAODPid::GetTRDncls() const {
120 // Get number of clusters per track
121 // Calculated as sum of the number of clusters per tracklet
124 for(Int_t ily = 0; ily < 6; ily++) ncls += fTRDncls[ily];