1 #ifndef ALIMULTIPLICITY_H
2 #define ALIMULTIPLICITY_H
8 ////////////////////////////////////////////////////////
9 //// Class containing multiplicity information //
10 //// to stored in the ESD //
11 ////////////////////////////////////////////////////////
13 class AliMultiplicity : public TObject {
17 AliMultiplicity(); // default constructor
18 // standard constructor
19 AliMultiplicity(Int_t ntr,Float_t *th, Float_t *ph, Float_t *dth, Float_t *dph, Int_t *labels,
20 Int_t* labelsL2, Int_t ns, Float_t *ts, Float_t *ps, Short_t nfcL1, Short_t nfcL2, TBits fFastOrFiredChips);
21 AliMultiplicity(const AliMultiplicity& m);
22 AliMultiplicity& operator=(const AliMultiplicity& m);
23 virtual void Copy(TObject &obj) const;
24 virtual ~AliMultiplicity();
25 // methods to access tracklet information
26 Int_t GetNumberOfTracklets() const {return fNtracks;}
27 Double_t GetTheta(Int_t i) const { if(i>=0 && i<fNtracks) {return fTh[i];}
28 else {Error("GetTheta","Invalid track number %d",i); return -9999.;}}
29 Double_t GetEta(Int_t i) const { if(i>=0 && i<fNtracks) {return -TMath::Log(TMath::Tan(fTh[i]/2.));}
30 else {Error("GetEta","Invalid track number %d",i); return -9999.;}}
31 Double_t GetPhi(Int_t i) const { if(i>=0 && i<fNtracks) {return fPhi[i];}
32 else {Error("GetPhi","Invalid track number %d",i); return -9999.;}}
33 Double_t GetDeltaTheta(Int_t i) const {if(fDeltTh && i>=0 && i<fNtracks) {return fDeltTh[i];}
34 else {Error("GetDeltaTheta","DeltaTheta not available in data or Invalid track number %d(max %d)",i, fNtracks); return -9999.;}}
35 Double_t GetDeltaPhi(Int_t i) const {if(i>=0 && i<fNtracks) {return fDeltPhi[i];}
36 else {Error("GetDeltaPhi","Invalid track number %d",i); return -9999.;}}
38 Int_t GetLabel(Int_t i, Int_t layer) const;
39 void SetLabel(Int_t i, Int_t layer, Int_t label);
41 // methods to access single cluster information
42 Int_t GetNumberOfSingleClusters() const {return fNsingle;}
43 Double_t GetThetaSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fThsingle[i];}
44 else {Error("GetThetaSingle","Invalid cluster number %d",i); return -9999.;}}
45 Double_t GetPhiSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fPhisingle[i];}
46 else {Error("GetPhisingle","Invalid cluster number %d",i); return -9999.;}}
48 Short_t GetNumberOfFiredChips(Int_t layer) const { return fFiredChips[layer]; }
49 void SetFiredChips(Int_t layer, Short_t firedChips) { fFiredChips[layer] = firedChips; }
51 void SetFastOrFiredChips(UInt_t chipKey){fFastOrFiredChips.SetBitNumber(chipKey);}
52 TBits GetFastOrFiredChips() const {return fFastOrFiredChips;}
53 Bool_t TestFastOrFiredChips(UInt_t chipKey) const {return fFastOrFiredChips.TestBitNumber(chipKey);}
56 void Duplicate(const AliMultiplicity &m); // used by copy ctr.
58 Int_t fNtracks; // Number of tracklets
59 Int_t fNsingle; // Number of clusters on SPD layer 1, not associated
60 // with a tracklet on SPD layer 2
61 Int_t *fLabels; //[fNtracks] array with labels of cluster in L1 used for tracklet
62 Int_t *fLabelsL2; //[fNtracks] array with labels of cluster in L2 used for tracklet
63 Double32_t *fTh; //[fNtracks] array with theta values
64 Double32_t *fPhi; //[fNtracks] array with phi values
65 Double32_t *fDeltTh; //[fNtracks] array with delta theta values
66 Double32_t *fDeltPhi; //[fNtracks] array with delta phi values
67 Double32_t *fThsingle; //[fNsingle] array with theta values of L1 clusters
68 Double32_t *fPhisingle; //[fNsingle] array with phi values of L2 clusters
69 Short_t fFiredChips[2]; // Number of fired chips in the two SPD layers
71 TBits fFastOrFiredChips; // Map of FastOr fired chips
73 ClassDef(AliMultiplicity,9);
76 inline Int_t AliMultiplicity::GetLabel(Int_t i, Int_t layer) const
78 if(i>=0 && i<fNtracks) {
81 } else if (layer == 1) {
85 Warning("GetLabel", "No information for layer 2 available !");
89 Error("GetLabel","Invalid layer number %d",layer); return -9999;
92 Error("GetLabel","Invalid track number %d",i); return -9999;