2 #include "AliMultiplicity.h"
5 ClassImp(AliMultiplicity)
7 //______________________________________________________________________
8 AliMultiplicity::AliMultiplicity():
20 fFastOrFiredChips(1200)
22 // Default Constructor
25 for(Int_t ilayer = 0; ilayer < 6; ilayer++)fITSClusters[ilayer] = 0;
28 //______________________________________________________________________
29 AliMultiplicity::AliMultiplicity(Int_t ntr, Float_t *th, Float_t *ph, Float_t *dth, Float_t *dph, Int_t *labels, Int_t* labelsL2, Int_t ns, Float_t *ts, Float_t *ps, Short_t nfcL1, Short_t nfcL2, TBits fFastOr):
41 fFastOrFiredChips(1200)
43 // Standard constructor
45 fLabels = new Int_t[ntr];
46 fLabelsL2 = new Int_t[ntr];
47 fTh = new Double_t [ntr];
48 fPhi = new Double_t [ntr];
49 fDeltTh = new Double_t [ntr];
50 fDeltPhi = new Double_t [ntr];
51 for(Int_t i=0;i<fNtracks;i++){
56 fLabels[i] = labels[i];
57 fLabelsL2[i] = labelsL2[i];
61 fThsingle = new Double_t [ns];
62 fPhisingle = new Double_t [ns];
63 for(Int_t i=0;i<fNsingle;i++){
68 fFiredChips[0] = nfcL1;
69 fFiredChips[1] = nfcL2;
70 fFastOrFiredChips = fFastOr;
71 for(Int_t ilayer = 0; ilayer < 6; ilayer++)fITSClusters[ilayer] = 0;
74 //______________________________________________________________________
75 AliMultiplicity::AliMultiplicity(const AliMultiplicity& m):
87 fFastOrFiredChips(1200)
93 //______________________________________________________________________
94 AliMultiplicity &AliMultiplicity::operator=(const AliMultiplicity& m){
95 // assignment operator
96 if(this == &m)return *this;
97 ((TObject *)this)->operator=(m);
99 if(fTh)delete [] fTh;fTh = 0;
100 if(fPhi)delete [] fPhi;fPhi = 0;
101 if(fDeltTh)delete [] fDeltTh;fDeltTh= 0;
102 if(fDeltPhi)delete [] fDeltPhi;fDeltPhi = 0;
103 if(fLabels)delete [] fLabels;fLabels = 0;
104 if(fLabelsL2)delete [] fLabelsL2;fLabelsL2 = 0;
105 if(fThsingle)delete [] fThsingle;fThsingle = 0;
106 if(fPhisingle)delete [] fPhisingle;fPhisingle = 0;
112 void AliMultiplicity::Copy(TObject &obj) const {
114 // this overwrites the virtual TOBject::Copy()
115 // to allow run time copying without casting
118 if(this==&obj)return;
119 AliMultiplicity *robj = dynamic_cast<AliMultiplicity*>(&obj);
120 if(!robj)return; // not an AliMultiplicity
126 //______________________________________________________________________
127 void AliMultiplicity::Duplicate(const AliMultiplicity& m){
128 // used by copy constructor and assignment operator
129 fNtracks = m.fNtracks;
131 fTh = new Double_t[fNtracks];
132 fPhi = new Double_t[fNtracks];
133 fDeltTh = new Double_t[fNtracks];
134 fDeltPhi = new Double_t[fNtracks];
135 fLabels = new Int_t[fNtracks];
136 fLabelsL2 = new Int_t[fNtracks];
146 fNsingle = m.fNsingle;
148 fThsingle = new Double_t[fNsingle];
149 fPhisingle = new Double_t[fNsingle];
155 if(m.fTh)memcpy(fTh,m.fTh,fNtracks*sizeof(Double_t));
156 if(m.fPhi)memcpy(fPhi,m.fPhi,fNtracks*sizeof(Double_t));
157 if(m.fDeltTh)memcpy(fDeltTh,m.fDeltTh,fNtracks*sizeof(Double_t));
158 if(m.fDeltPhi)memcpy(fDeltPhi,m.fDeltPhi,fNtracks*sizeof(Double_t));
159 if(m.fLabels)memcpy(fLabels,m.fLabels,fNtracks*sizeof(Int_t));
160 if(m.fLabelsL2)memcpy(fLabelsL2,m.fLabelsL2,fNtracks*sizeof(Int_t));
161 if(m.fThsingle)memcpy(fThsingle,m.fThsingle,fNsingle*sizeof(Double_t));
162 if(m.fPhisingle)memcpy(fPhisingle,m.fPhisingle,fNsingle*sizeof(Double_t));
164 fFiredChips[0] = m.fFiredChips[0];
165 fFiredChips[1] = m.fFiredChips[1];
166 for(Int_t ilayer = 0; ilayer < 6; ilayer++){
167 fITSClusters[ilayer] = m.fITSClusters[ilayer];
172 fFastOrFiredChips = m.fFastOrFiredChips;
175 //______________________________________________________________________
176 AliMultiplicity::~AliMultiplicity(){
178 if(fTh)delete [] fTh;fTh = 0;
179 if(fPhi)delete [] fPhi;fPhi = 0;
180 if(fDeltTh)delete [] fDeltTh;fDeltTh = 0;
181 if(fDeltPhi)delete [] fDeltPhi;fDeltPhi = 0;
182 if(fLabels)delete [] fLabels;fLabels = 0;
183 if(fLabelsL2)delete [] fLabelsL2;fLabelsL2 = 0;
184 if(fThsingle)delete [] fThsingle;fThsingle = 0;
185 if(fPhisingle)delete [] fPhisingle;fPhisingle = 0;
189 //______________________________________________________________________
190 void AliMultiplicity::SetLabel(Int_t i, Int_t layer, Int_t label)
192 if(i>=0 && i<fNtracks) {
196 } else if (layer == 1) {
198 fLabelsL2[i] = label;
203 Error("SetLabel","Invalid track number %d or layer %d",i,layer);
206 //______________________________________________________________________
207 UInt_t AliMultiplicity::GetNumberOfITSClusters(Int_t layMin, Int_t layMax) const {
209 if(layMax < layMin) {
210 AliError("layer min > layer max");
214 AliError("layer min < 0");
218 AliError("layer max > 0");
223 for (Int_t i=layMin; i<=layMax; i++) sum+=fITSClusters[i];