1 /**************************************************************************
2 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
7 * Permission to use, copy, modify and distribute this software and its *
8 * documentation strictly for non-commercial purposes is hereby granted *
9 * without fee, provided that the above copyright notice appears in all *
10 * copies and that both the copyright notice and this permission notice *
11 * appear in the supporting documentation. The authors make no claims *
12 * about the suitability of this software for any purpose. It is *
13 * provided "as is" without express or implied warranty. *
14 **************************************************************************/
18 Revision 1.3 2000/06/30 12:07:49 kowal2
19 Updated from the TPC-PreRelease branch
21 Revision 1.2.4.2 2000/06/26 07:39:42 kowal2
22 Changes to obey the coding rules
24 Revision 1.2.4.1 2000/06/25 08:38:41 kowal2
25 Splitted from AliTPCtracking
27 Revision 1.2 2000/04/17 09:37:33 kowal2
28 removed obsolete AliTPCDigitsDisplay.C
30 Revision 1.1.4.2 2000/04/10 11:32:37 kowal2
32 "ROOT"-based class with some extra functionality
37 //----------------------------------------------------------------------------
38 // Author: Marian Ivanov
40 // Implementation of class AliH2F
42 //-----------------------------------------------------------------------------
45 #include "TClonesArray.h"
48 #include "AliTPCClusterFinder.h"
52 //***********************************************************************
53 //***********************************************************************
54 //***********************************************************************
55 //***********************************************************************
56 AliH2F::AliH2F():TH2F()
61 AliH2F::AliH2F(const Text_t *name,const Text_t *title,
62 Int_t nbinsx,Axis_t xlow,Axis_t xup
63 ,Int_t nbinsy,Axis_t ylow,Axis_t yup):
64 TH2F(name,title,nbinsx,xlow,xup
76 AliH2F::AliH2F(const AliH2F &his)
82 AliH2F & AliH2F::operator = (const AliH2F & his)
88 TClonesArray * AliH2F::FindPeaks(Float_t threshold, Float_t noise)
90 //find peaks and write it in form of AliTPCcluster to array
92 //firstly we need to create object for cluster finding
93 //and fill it with contents of histogram
94 AliTPCClusterFinder cfinder;
95 cfinder.SetThreshold(threshold);
96 cfinder.SetNoise(noise);
97 cfinder.GetHisto(this);
98 return cfinder.FindPeaks3();
101 void AliH2F::ClearSpectrum()
104 Int_t dimx = fXaxis.GetNbins();
105 Int_t dimy = fYaxis.GetNbins();
106 for (Int_t i = 0 ;i<dimx;i++)
107 for (Int_t j = 0 ;j<dimy;j++)
109 SetCellContent(i,j,0);
115 void AliH2F::AddNoise(Float_t sn)
117 // add gauss noise with sigma sn
118 Int_t dimx = fXaxis.GetNbins();
119 Int_t dimy = fYaxis.GetNbins();
120 for (Int_t i = 0 ;i<dimx;i++)
121 for (Int_t j = 0 ;j<dimy;j++)
123 Float_t noise = gRandom->Gaus(0,sn);
124 Float_t oldv =GetCellContent(i,j);
125 Float_t olds =GetCellError(i,j);
128 SetCellContent(i,j,noise+oldv);
129 SetCellError(i,j,TMath::Sqrt((noise*noise+olds*olds)));
133 void AliH2F::AddGauss(Float_t x, Float_t y,
134 Float_t sx, Float_t sy, Float_t max)
136 //transform to histogram coordinata
137 Int_t dimx = fXaxis.GetNbins();
138 Int_t dimy = fYaxis.GetNbins();
139 Float_t dx =(GetXaxis()->GetXmax()-GetXaxis()->GetXmin())/Float_t(dimx);
140 Float_t dy =(GetYaxis()->GetXmax()-GetYaxis()->GetXmin())/Float_t(dimy);
141 // x=(x-GetXaxis()->GetXmin())/dx;
142 //y=(y-GetYaxis()->GetXmin())/dy;
147 for (Int_t i = 0 ;i<dimx;i++)
148 for (Int_t j = 0 ;j<dimy;j++)
150 Float_t x2 =GetXaxis()->GetBinCenter(i+1);
151 Float_t y2 =GetYaxis()->GetBinCenter(j+1);
152 Float_t dx2 = (x2-x)*(x2-x);
153 Float_t dy2 = (y2-y)*(y2-y);
154 Float_t amp =max*exp(-(dx2/(2*sx*sx)+dy2/(2*sy*sy)));
155 //Float_t oldv =GetCellContent(i+1,j+1);
156 // SetCellContent(i+1,j+1,amp+oldv);
161 void AliH2F::ClearUnderTh(Int_t threshold)
163 //clear histogram for bin under threshold
164 Int_t dimx = fXaxis.GetNbins();
165 Int_t dimy = fYaxis.GetNbins();
166 for (Int_t i = 0 ;i<=dimx;i++)
167 for (Int_t j = 0 ;j<=dimy;j++)
169 Float_t oldv =GetCellContent(i,j);
171 SetCellContent(i,j,0);
177 //round float to integer
178 Int_t dimx = fXaxis.GetNbins();
179 Int_t dimy = fYaxis.GetNbins();
180 for (Int_t i = 0 ;i<=dimx;i++)
181 for (Int_t j = 0 ;j<=dimy;j++)
183 Float_t oldv =GetCellContent(i,j);
185 SetCellContent(i,j,oldv);
191 AliH2F *AliH2F::GetSubrange2d(Float_t xmin, Float_t xmax,
192 Float_t ymin, Float_t ymax)
194 //this function return pointer to the new created
195 //histogram which is subhistogram of the
197 //subhistogram range must be inside histogram
200 xmin=fXaxis.GetXmin();
201 xmax=fXaxis.GetXmax();
204 ymin=fYaxis.GetXmin();
205 ymax=fYaxis.GetXmax();
208 Int_t nx = Int_t((xmax-xmin)/(fXaxis.GetXmax()-fXaxis.GetXmin()) *
209 Float_t(fXaxis.GetNbins()));
210 Int_t ny = Int_t((ymax-ymin)/(fYaxis.GetXmax()-fYaxis.GetXmin()) *
211 Float_t(fYaxis.GetNbins()));
213 TString t2 = fTitle ;
216 const Text_t *ktt1 = t1;
217 const Text_t *ktt2 = t2;
219 AliH2F * sub = new AliH2F(ktt1,ktt2,nx,xmin,xmax,ny,ymin,ymax);
221 Int_t i1 = Int_t( Float_t(fXaxis.GetNbins())*(xmin-fXaxis.GetXmin())/
222 (fXaxis.GetXmax()-fXaxis.GetXmin()) ) ;
223 Int_t i2 = Int_t( Float_t(fYaxis.GetNbins())*(ymin-fYaxis.GetXmin())/
224 (fYaxis.GetXmax()-fYaxis.GetXmin()) ) ;
225 for (Int_t i=0;i<nx;i++)
226 for (Int_t j=0;j<ny;j++)
228 Int_t index1 = GetBin(i1+i,i2+j);
229 // Int_t index2 = sub->GetBin(i,j);
230 Float_t val = GetBinContent(index1);
231 // sub->SetBinContent(index2,val);
232 // Float_t err = GetBinError(index1);
233 //sub->SetBinError(index2,GetBinError(index1));
234 sub->SetCellContent(i,j,val);
239 TH1F *AliH2F::GetAmplitudes(Float_t zmin, Float_t zmax, Float_t th, Float_t xmin, Float_t xmax,
240 Float_t ymin, Float_t ymax)
242 //this function return pointer to the new created
243 //histogram which is subhistogram of the
245 //subhistogram range must be inside histogram
248 xmin=fXaxis.GetXmin();
249 xmax=fXaxis.GetXmax();
252 ymin=fYaxis.GetXmin();
253 ymax=fYaxis.GetXmax();
255 Int_t nx = Int_t((xmax-xmin)/(fXaxis.GetXmax()-fXaxis.GetXmin()) *
256 Float_t(fXaxis.GetNbins()));
257 Int_t ny = Int_t((ymax-ymin)/(fYaxis.GetXmax()-fYaxis.GetXmin()) *
258 Float_t(fYaxis.GetNbins()));
260 TString t2 = fTitle ;
263 const Text_t *ktt1 = t1;
264 const Text_t *ktt2 = t2;
266 TH1F * h = new TH1F(ktt1,ktt2,100,zmin,zmax);
268 Int_t i1 = Int_t( Float_t(fXaxis.GetNbins())*(xmin-fXaxis.GetXmin())/
269 (fXaxis.GetXmax()-fXaxis.GetXmin()) ) ;
270 Int_t i2 = Int_t( Float_t(fYaxis.GetNbins())*(ymin-fYaxis.GetXmin())/
271 (fYaxis.GetXmax()-fYaxis.GetXmin()) ) ;
272 for (Int_t i=0;i<nx;i++)
273 for (Int_t j=0;j<ny;j++)
275 Int_t index1 = GetBin(i1+i,i2+j);
276 Float_t val = GetBinContent(index1);
277 if (val>th) h->Fill(val);
282 Float_t AliH2F::GetOccupancy(Float_t th , Float_t xmin, Float_t xmax,
283 Float_t ymin, Float_t ymax)
285 //this function return pointer to the new created
286 //histogram which is subhistogram of the
288 //subhistogram range must be inside histogram
291 xmin=fXaxis.GetXmin();
292 xmax=fXaxis.GetXmax();
295 ymin=fYaxis.GetXmin();
296 ymax=fYaxis.GetXmax();
298 Int_t nx = Int_t((xmax-xmin)/(fXaxis.GetXmax()-fXaxis.GetXmin()) *
299 Float_t(fXaxis.GetNbins()));
300 Int_t ny = Int_t((ymax-ymin)/(fYaxis.GetXmax()-fYaxis.GetXmin()) *
301 Float_t(fYaxis.GetNbins()));
304 Int_t i1 = Int_t( Float_t(fXaxis.GetNbins())*(xmin-fXaxis.GetXmin())/
305 (fXaxis.GetXmax()-fXaxis.GetXmin()) ) ;
306 Int_t i2 = Int_t( Float_t(fYaxis.GetNbins())*(ymin-fYaxis.GetXmin())/
307 (fYaxis.GetXmax()-fYaxis.GetXmin()) ) ;
308 for (Int_t i=0;i<nx;i++)
309 for (Int_t j=0;j<ny;j++)
311 Int_t index1 = GetBin(i1+i,i2+j);
312 Float_t val = GetBinContent(index1);
316 if (all>0) return Float_t(over)/Float_t(all);