]> git.uio.no Git - u/mrichter/AliRoot.git/blame_incremental - STEER/AliMultiplicity.h
Memory leaks in TOF visualization code
[u/mrichter/AliRoot.git] / STEER / AliMultiplicity.h
... / ...
CommitLineData
1#ifndef ALIMULTIPLICITY_H
2#define ALIMULTIPLICITY_H
3
4#include<TObject.h>
5#include <TBits.h>
6#include<TMath.h>
7
8////////////////////////////////////////////////////////
9//// Class containing multiplicity information //
10//// to stored in the ESD //
11////////////////////////////////////////////////////////
12
13class AliMultiplicity : public TObject {
14
15 public:
16
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(i>=0 && i<fNtracks) {return fDeltTh[i];}
34 else {Error("GetDeltaTheta","Invalid track number %d",i); 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.;}}
37
38 Int_t GetLabel(Int_t i, Int_t layer) const;
39
40// methods to access single cluster information
41 Int_t GetNumberOfSingleClusters() const {return fNsingle;}
42 Double_t GetThetaSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fThsingle[i];}
43 else {Error("GetThetaSingle","Invalid cluster number %d",i); return -9999.;}}
44 Double_t GetPhiSingle(Int_t i) const { if(i>=0 && i<fNsingle) {return fPhisingle[i];}
45 else {Error("GetPhisingle","Invalid cluster number %d",i); return -9999.;}}
46
47 Short_t GetNumberOfFiredChips(Int_t layer) const { return fFiredChips[layer]; }
48 void SetFiredChips(Int_t layer, Short_t firedChips) { fFiredChips[layer] = firedChips; }
49
50 void SetFastOrFiredChips(UInt_t chipKey){fFastOrFiredChips.SetBitNumber(chipKey);}
51 TBits GetFastOrFiredChips() const {return fFastOrFiredChips;}
52 Bool_t TestFastOrFiredChips(UInt_t chipKey) const {return fFastOrFiredChips.TestBitNumber(chipKey);}
53
54 protected:
55 void Duplicate(const AliMultiplicity &m); // used by copy ctr.
56
57 Int_t fNtracks; // Number of tracklets
58 Int_t fNsingle; // Number of clusters on SPD layer 1, not associated
59 // with a tracklet on SPD layer 2
60 Int_t *fLabels; //[fNtracks] array with labels of cluster in L1 used for tracklet
61 Int_t *fLabelsL2; //[fNtracks] array with labels of cluster in L2 used for tracklet
62 Double32_t *fTh; //[fNtracks] array with theta values
63 Double32_t *fPhi; //[fNtracks] array with phi values
64 Double32_t *fDeltTh; //[fNtracks] array with delta theta values
65 Double32_t *fDeltPhi; //[fNtracks] array with delta phi values
66 Double32_t *fThsingle; //[fNsingle] array with theta values of L1 clusters
67 Double32_t *fPhisingle; //[fNsingle] array with phi values of L2 clusters
68 Short_t fFiredChips[2]; // Number of fired chips in the two SPD layers
69
70 TBits fFastOrFiredChips; // Map of FastOr fired chips
71
72 ClassDef(AliMultiplicity,9);
73};
74
75inline Int_t AliMultiplicity::GetLabel(Int_t i, Int_t layer) const
76{
77 if(i>=0 && i<fNtracks) {
78 if (layer == 0) {
79 return fLabels[i];
80 } else if (layer == 1) {
81 if (fLabelsL2) {
82 return fLabelsL2[i];
83 } else {
84 Warning("GetLabel", "No information for layer 2 available !");
85 return -9999;
86 }
87 } else {
88 Error("GetLabel","Invalid layer number %d",layer); return -9999;
89 }
90 } else {
91 Error("GetLabel","Invalid track number %d",i); return -9999;
92 }
93}
94#endif