]>
Commit | Line | Data |
---|---|---|
32e449be | 1 | #ifndef ALIMULTIPLICITY_H |
2 | #define ALIMULTIPLICITY_H | |
3 | ||
4 | #include<TObject.h> | |
44ca59b0 | 5 | #include<TMath.h> |
32e449be | 6 | |
7 | //////////////////////////////////////////////////////// | |
8 | //// Class containing multiplicity information // | |
9 | //// to stored in the ESD // | |
10 | //////////////////////////////////////////////////////// | |
11 | ||
12 | class AliMultiplicity : public TObject { | |
13 | ||
14 | public: | |
15 | ||
16 | AliMultiplicity(); // default constructor | |
17 | // standard constructor | |
59ea252e | 18 | AliMultiplicity(Int_t ntr,Float_t *t, Float_t *ph, Float_t *df, Int_t *labels, |
0939e22a | 19 | Int_t* labelsL2, Int_t ns, Float_t *ts, Float_t *ps); |
32e449be | 20 | AliMultiplicity(const AliMultiplicity& m); |
21 | AliMultiplicity& operator=(const AliMultiplicity& m); | |
22 | virtual ~AliMultiplicity(); | |
968e8539 | 23 | // methods to access tracklet information |
32e449be | 24 | Int_t GetNumberOfTracklets() const {return fNtracks;} |
cd888a89 | 25 | Double_t GetTheta(Int_t i) const { if(i>=0 && i<fNtracks) {return fTh[i];} |
32e449be | 26 | else {Error("GetTheta","Invalid track number %d",i); return -9999.;}} |
44ca59b0 | 27 | Double_t GetEta(Int_t i) const { if(i>=0 && i<fNtracks) {return -TMath::Log(TMath::Tan(fTh[i]/2.));} |
28 | else {Error("GetTheta","Invalid track number %d",i); return -9999.;}} | |
cd888a89 | 29 | Double_t GetPhi(Int_t i) const { if(i>=0 && i<fNtracks) {return fPhi[i];} |
968e8539 | 30 | else {Error("GetPhi","Invalid track number %d",i); return -9999.;}} |
cd888a89 | 31 | Double_t GetDeltaPhi(Int_t i) const {if(i>=0 && i<fNtracks) {return fDeltPhi[i];} |
32e449be | 32 | else {Error("GetDeltaPhi","Invalid track number %d",i); return -9999.;}} |
0939e22a | 33 | Int_t GetLabel(Int_t i, Int_t layer) const {if(i>=0 && i<fNtracks) {return (layer == 0) ? fLabels[i] : fLabelsL2[i];} |
59ea252e | 34 | else {Error("GetLabel","Invalid track number %d",i); return -9999;}} |
968e8539 | 35 | // methods to access single cluster information |
36 | Int_t GetNumberOfSingleClusters() const {return fNsingle;} | |
cd888a89 | 37 | Double_t GetThetaSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fThsingle[i];} |
968e8539 | 38 | else {Error("GetThetaSingle","Invalid cluster number %d",i); return -9999.;}} |
cd888a89 | 39 | Double_t GetPhiSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fPhisingle[i];} |
968e8539 | 40 | else {Error("GetPhisingle","Invalid cluster number %d",i); return -9999.;}} |
32e449be | 41 | |
466ffa28 | 42 | Short_t GetFiredChips(Int_t layer) { return fFiredChips[layer]; } |
43 | void SetFiredChips(Int_t layer, Short_t firedChips) { fFiredChips[layer] = firedChips; } | |
44 | ||
32e449be | 45 | protected: |
46 | void Duplicate(const AliMultiplicity &m); // used by copy ctr. | |
32e449be | 47 | |
cd888a89 | 48 | Int_t fNtracks; // Number of tracklets |
49 | Int_t fNsingle; // Number of clusters on SPD layer 1, not associated | |
cd888a89 | 50 | // with a tracklet on SPD layer 2 |
0939e22a | 51 | Int_t *fLabels; //[fNtracks] array with labels of cluster in L1 used for tracklet |
52 | Int_t *fLabelsL2; //[fNtracks] array with labels of cluster in L2 used for tracklet | |
cd888a89 | 53 | Double32_t *fTh; //[fNtracks] array with theta values |
54 | Double32_t *fPhi; //[fNtracks] array with phi values | |
55 | Double32_t *fDeltPhi; //[fNtracks] array with delta phi values | |
56 | Double32_t *fThsingle; //[fNsingle] array with theta values of L1 clusters | |
57 | Double32_t *fPhisingle; //[fNsingle] array with phi values of L2 clusters | |
466ffa28 | 58 | Short_t fFiredChips[2]; // number of fired chips in the two SPD layers |
59 | ||
0939e22a | 60 | ClassDef(AliMultiplicity,7); |
32e449be | 61 | }; |
62 | ||
63 | #endif |