[u/mrichter/AliRoot.git] / TPC / AliTPCclusterKr.cxx

/**************************************************************************\r
 * Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
 *                                                                        *\r
 * Author: The ALICE Off-line Project.                                    *\r
 * Contributors are mentioned in the code where appropriate.              *\r
 *                                                                        *\r
 * Permission to use, copy, modify and distribute this software and its   *\r
 * documentation strictly for non-commercial purposes is hereby granted   *\r
 * without fee, provided that the above copyright notice appears in all   *\r
 * copies and that both the copyright notice and this permission notice   *\r
 * appear in the supporting documentation. The authors make no claims     *\r
 * about the suitability of this software for any purpose. It is          *\r
 * provided "as is" without express or implied warranty.                  *\r
 **************************************************************************/\r
\r
/* $Id: AliTPCclusterKr.cxx,v 1.7 2008/01/22 17:24:53 matyja Exp $ */\r
\r
//-----------------------------------------------------------------\r
//           Implementation of the TPC Kr cluster class\r
//\r
// Origin: Adam Matyja, INP PAN, adam.matyja@ifj.edu.pl\r
//-----------------------------------------------------------------\r
\r
#include "AliTPCclusterKr.h"\r
#include "AliCluster.h"\r
#include "AliTPCvtpr.h"\r
#include "TObjArray.h"\r
//#include "TH1F.h"\r
#include "TMath.h"\r
#include "TArrayI.h"\r
\r
ClassImp(AliTPCclusterKr)\r
\r
\r
AliTPCclusterKr::AliTPCclusterKr()\r
:AliCluster(),\r
 fMax(),\r
 fADCcluster(0),\r
 fSec(0),\r
 fNPads(0),\r
 fNRows(0),\r
 fTimebins1D(0),\r
 fPads1D(0),\r
 fPadRMS(0),\r
 fRowRMS(0),\r
 fTimebinRMS(0),\r
 fSize(0),\r
 fCenterX(0),\r
 fCenterY(0),\r
 fCenterT(0),\r
 fCluster(0)\r
{\r
//\r
// default constructor\r
//\r
  fCluster=new TObjArray();\r
}\r
\r
AliTPCclusterKr::AliTPCclusterKr(const AliTPCclusterKr &param)\r
:AliCluster(param),\r
 fMax(),\r
 fADCcluster(0),\r
 fSec(0),\r
 fNPads(0),\r
 fNRows(0),\r
 fTimebins1D(0),\r
 fPads1D(0),\r
 fPadRMS(0),\r
 fRowRMS(0),\r
 fTimebinRMS(0),\r
 fSize(0),\r
 fCenterX(0),\r
 fCenterY(0),\r
 fCenterT(0),\r
 fCluster(0)\r
{\r
//\r
// copy constructor\r
//\r
  fADCcluster = param.fADCcluster;\r
  fSec  = param.fSec ;\r
  fNPads = param.fNPads;\r
  fNRows = param.fNRows;\r
  fMax = param.fMax;\r
  //  fCluster = param.fCluster;\r
  fCenterX = param.fCenterX;\r
  fCenterY = param.fCenterY;\r
  fCenterT = param.fCenterT;\r
  fCluster=new TObjArray(*(param.fCluster));\r
  fSize = param.fSize;\r
  fTimebins1D = param.fTimebins1D;\r
  fPads1D = param.fPads1D;\r
  fPadRMS = param.fPadRMS;\r
  fRowRMS = param.fRowRMS;\r
  fTimebinRMS = param.fTimebinRMS;\r
} \r
\r
AliTPCclusterKr &AliTPCclusterKr::operator = (const AliTPCclusterKr & param)\r
{\r
  //\r
  // assignment operator\r
  // \r
  (AliCluster&)(*this) = (AliCluster&)param;\r
  fADCcluster = param.fADCcluster;\r
  fSec  = param.fSec ;\r
  fNPads = param.fNPads;\r
  fNRows = param.fNRows;\r
  fMax = param.fMax;\r
  //  fCluster=param.fCluster;\r
  fCenterX = param.fCenterX;\r
  fCenterY = param.fCenterY;\r
  fCenterT = param.fCenterT;\r
  delete fCluster;\r
  fCluster=new TObjArray(*(param.fCluster));\r
  fSize=param.fSize;\r
  fTimebins1D = param.fTimebins1D;\r
  fPads1D = param.fPads1D;\r
  fPadRMS = param.fPadRMS;\r
  fRowRMS = param.fRowRMS;\r
  fTimebinRMS = param.fTimebinRMS;\r
  return (*this);\r
}\r
\r
AliTPCclusterKr::~AliTPCclusterKr()\r
{\r
  //\r
  // destructor\r
  //\r
  if(fCluster) {\r
    fCluster->SetOwner(kTRUE);\r
    fCluster->Delete();\r
    delete fCluster;\r
  }\r
  fCluster=0;\r
}\r
\r
////____________________________________________________________________________\r
void AliTPCclusterKr::SetCenter(){\r
  //\r
  // calculate geometrical center of the cluster\r
  //\r
  Double_t rX=0;\r
  Double_t rY=0;\r
  Double_t rT=0;\r
\r
  Short_t adc;\r
  fADCcluster=0;\r
  for(Int_t iter = 0; iter < fCluster->GetEntriesFast(); ++iter) {\r
    AliTPCvtpr *iclus=(AliTPCvtpr *)fCluster->At(iter);\r
\r
    //for( std::vector<AliTPCvtpr*>::iterator iclus  = fCluster.begin();\r
    //iclus != fCluster.end(); ++iclus ) {\r
    adc = (iclus)->GetAdc();\r
    fADCcluster+=adc;\r
    rX += ((iclus)->GetX() * adc);\r
    rY += ((iclus)->GetY() * adc);\r
    rT += ((iclus)->GetT() * adc);\r
  }\r
  fCenterX=rX/fADCcluster;\r
  fCenterY=rY/fADCcluster;\r
  fCenterT=rT/fADCcluster;\r
\r
  return;\r
}\r
\r
void AliTPCclusterKr::SetPadRMS(){\r
  //\r
  // calculate RMS in pad direction\r
  //\r
  //  TH1F *histo= new TH1F("","",200,0,200);\r
  TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
  for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
    {\r
      array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetPad(),i);\r
      //histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetPad() );\r
    }\r
  //  fPadRMS=histo->GetRMS();\r
  fPadRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
  //  delete histo;\r
  delete array;\r
  return;\r
}\r
\r
void AliTPCclusterKr::SetRowRMS(){\r
  //\r
  // calculate RMS in row direction\r
  //\r
  TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
  //  TH1F *histo= new TH1F("","",120,0,120);\r
  for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
    {\r
      array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetRow(),i);\r
      //      histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetRow() );\r
    }\r
  //  fRowRMS=histo->GetRMS();\r
  fRowRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
  //  delete histo;\r
  delete array;\r
  return;\r
}\r
\r
void AliTPCclusterKr::SetTimebinRMS(){\r
  //\r
  // calculate RMS in timebin direction\r
  //\r
  TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
  //  TH1F *histo= new TH1F("","",1000,0,1000);\r
  for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
    {\r
      array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetTime(),i);\r
      //      histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetTime() );\r
    }\r
  fTimebinRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
  //histo->GetRMS();\r
  //  delete histo;\r
  delete array;\r
  return;\r
}\r
\r
void AliTPCclusterKr::SetRMS(){\r
  //\r
  // calculate RMS in pad,row,timebin direction\r
  //\r
  TArrayI *arrayPad = new TArrayI(fCluster->GetEntriesFast());\r
  TArrayI *arrayRow = new TArrayI(fCluster->GetEntriesFast());\r
  TArrayI *arrayTime= new TArrayI(fCluster->GetEntriesFast());\r
  //  TH1F *histoPad= new TH1F("p","p",200,0,200);\r
  //  TH1F *histoRow= new TH1F("r","r",120,0,120);\r
  //  TH1F *histoTime= new TH1F("t","t",1000,0,1000);\r
  for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
    {\r
      arrayPad->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetPad(),i);\r
      arrayRow->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetRow(),i);\r
      arrayTime->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetTime(),i);\r
\r
      //histoPad->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetPad() );\r
      //histoRow->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetRow() );\r
      //histoTime->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetTime() );\r
    }\r
  //  fPadRMS=histoPad->GetRMS();\r
  fPadRMS=TMath::RMS(arrayPad->GetSize(),arrayPad->GetArray());\r
  fRowRMS=TMath::RMS(arrayRow->GetSize(),arrayRow->GetArray());\r
    //histoRow->GetRMS();\r
  fTimebinRMS=TMath::RMS(arrayTime->GetSize(),arrayTime->GetArray());\r
    //histoTime->GetRMS();\r
\r
  delete arrayPad;\r
  delete arrayRow;\r
  delete arrayTime;\r
  //  delete histoPad;\r
  //  delete histoRow;\r
  //  delete histoTime;\r
\r
  return;\r
}\r
\r
\r
void AliTPCclusterKr::Set1D(){\r
  //\r
  //\r
  //\r
  Short_t maxTime=0;\r
  Short_t minTime=1000;\r
  Short_t maxPad=0;\r
  Short_t minPad=1000;\r
 \r
  for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
    {\r
      if(((AliTPCvtpr *)(fCluster->At(i)))->GetPad()>maxPad)maxPad   =((AliTPCvtpr *)(fCluster->At(i)))->GetPad();\r
      if(((AliTPCvtpr *)(fCluster->At(i)))->GetPad()<minPad)minPad   =((AliTPCvtpr *)(fCluster->At(i)))->GetPad();\r
      if(((AliTPCvtpr *)(fCluster->At(i)))->GetTime()>maxTime)maxTime=((AliTPCvtpr *)(fCluster->At(i)))->GetTime();\r
      if(((AliTPCvtpr *)(fCluster->At(i)))->GetTime()<minTime)minTime=((AliTPCvtpr *)(fCluster->At(i)))->GetTime();\r
    }\r
  fPads1D=maxPad-minPad+1;\r
  fTimebins1D=maxTime-minTime+1;\r
  return;\r
}\r
Commit	Line	Data
e2a4d72c	1	/**************************************************************************\r
	2	* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *\r
	3	* *\r
	4	* Author: The ALICE Off-line Project. *\r
	5	* Contributors are mentioned in the code where appropriate. *\r
	6	* *\r
	7	* Permission to use, copy, modify and distribute this software and its *\r
	8	* documentation strictly for non-commercial purposes is hereby granted *\r
	9	* without fee, provided that the above copyright notice appears in all *\r
	10	* copies and that both the copyright notice and this permission notice *\r
	11	* appear in the supporting documentation. The authors make no claims *\r
	12	* about the suitability of this software for any purpose. It is *\r
	13	* provided "as is" without express or implied warranty. *\r
	14	**************************************************************************/\r
	15	\r
	16	/* $Id: AliTPCclusterKr.cxx,v 1.7 2008/01/22 17:24:53 matyja Exp $ */\r
	17	\r
	18	//-----------------------------------------------------------------\r
	19	// Implementation of the TPC Kr cluster class\r
	20	//\r
	21	// Origin: Adam Matyja, INP PAN, adam.matyja@ifj.edu.pl\r
	22	//-----------------------------------------------------------------\r
	23	\r
	24	#include "AliTPCclusterKr.h"\r
	25	#include "AliCluster.h"\r
	26	#include "AliTPCvtpr.h"\r
	27	#include "TObjArray.h"\r
	28	//#include "TH1F.h"\r
	29	#include "TMath.h"\r
	30	#include "TArrayI.h"\r
	31	\r
	32	ClassImp(AliTPCclusterKr)\r
	33	\r
	34	\r
	35	AliTPCclusterKr::AliTPCclusterKr()\r
	36	:AliCluster(),\r
	37	fMax(),\r
	38	fADCcluster(0),\r
	39	fSec(0),\r
	40	fNPads(0),\r
	41	fNRows(0),\r
	42	fTimebins1D(0),\r
	43	fPads1D(0),\r
	44	fPadRMS(0),\r
	45	fRowRMS(0),\r
	46	fTimebinRMS(0),\r
	47	fSize(0),\r
	48	fCenterX(0),\r
	49	fCenterY(0),\r
	50	fCenterT(0),\r
	51	fCluster(0)\r
	52	{\r
	53	//\r
	54	// default constructor\r
	55	//\r
	56	fCluster=new TObjArray();\r
	57	}\r
	58	\r
	59	AliTPCclusterKr::AliTPCclusterKr(const AliTPCclusterKr &param)\r
	60	:AliCluster(param),\r
	61	fMax(),\r
	62	fADCcluster(0),\r
	63	fSec(0),\r
	64	fNPads(0),\r
65	fNRows(0),\r
66	fTimebins1D(0),\r
67	fPads1D(0),\r
68	fPadRMS(0),\r
69	fRowRMS(0),\r
70	fTimebinRMS(0),\r
71	fSize(0),\r
72	fCenterX(0),\r
73	fCenterY(0),\r
74	fCenterT(0),\r
75	fCluster(0)\r
76	{\r
77	//\r
78	// copy constructor\r
79	//\r
80	fADCcluster = param.fADCcluster;\r
81	fSec = param.fSec ;\r
82	fNPads = param.fNPads;\r
83	fNRows = param.fNRows;\r
84	fMax = param.fMax;\r
85	// fCluster = param.fCluster;\r
86	fCenterX = param.fCenterX;\r
87	fCenterY = param.fCenterY;\r
88	fCenterT = param.fCenterT;\r
89	fCluster=new TObjArray(*(param.fCluster));\r
90	fSize = param.fSize;\r
91	fTimebins1D = param.fTimebins1D;\r
92	fPads1D = param.fPads1D;\r
93	fPadRMS = param.fPadRMS;\r
94	fRowRMS = param.fRowRMS;\r
95	fTimebinRMS = param.fTimebinRMS;\r
96	} \r
97	\r
98	AliTPCclusterKr &AliTPCclusterKr::operator = (const AliTPCclusterKr & param)\r
99	{\r
100	//\r
101	// assignment operator\r
102	// \r
103	(AliCluster&)(*this) = (AliCluster&)param;\r
104	fADCcluster = param.fADCcluster;\r
105	fSec = param.fSec ;\r
106	fNPads = param.fNPads;\r
107	fNRows = param.fNRows;\r
108	fMax = param.fMax;\r
109	// fCluster=param.fCluster;\r
110	fCenterX = param.fCenterX;\r
111	fCenterY = param.fCenterY;\r
112	fCenterT = param.fCenterT;\r
113	delete fCluster;\r
114	fCluster=new TObjArray(*(param.fCluster));\r
115	fSize=param.fSize;\r
116	fTimebins1D = param.fTimebins1D;\r
117	fPads1D = param.fPads1D;\r
118	fPadRMS = param.fPadRMS;\r
119	fRowRMS = param.fRowRMS;\r
120	fTimebinRMS = param.fTimebinRMS;\r
121	return (*this);\r
122	}\r
123	\r
124	AliTPCclusterKr::~AliTPCclusterKr()\r
125	{\r
126	//\r
127	// destructor\r
128	//\r
129	if(fCluster) {\r
130	fCluster->SetOwner(kTRUE);\r
131	fCluster->Delete();\r
132	delete fCluster;\r
133	}\r
134	fCluster=0;\r
135	}\r
136	\r
137	////____________________________________________________________________________\r
138	void AliTPCclusterKr::SetCenter(){\r
139	//\r
140	// calculate geometrical center of the cluster\r
141	//\r
142	Double_t rX=0;\r
143	Double_t rY=0;\r
144	Double_t rT=0;\r
145	\r
146	Short_t adc;\r
147	fADCcluster=0;\r
148	for(Int_t iter = 0; iter < fCluster->GetEntriesFast(); ++iter) {\r
149	AliTPCvtpr iclus=(AliTPCvtpr )fCluster->At(iter);\r
150	\r
151	//for( std::vector<AliTPCvtpr*>::iterator iclus = fCluster.begin();\r
152	//iclus != fCluster.end(); ++iclus ) {\r
153	adc = (iclus)->GetAdc();\r
154	fADCcluster+=adc;\r
155	rX += ((iclus)->GetX() * adc);\r
156	rY += ((iclus)->GetY() * adc);\r
157	rT += ((iclus)->GetT() * adc);\r
158	}\r
159	fCenterX=rX/fADCcluster;\r
160	fCenterY=rY/fADCcluster;\r
161	fCenterT=rT/fADCcluster;\r
162	\r
163	return;\r
164	}\r
165	\r
166	void AliTPCclusterKr::SetPadRMS(){\r
167	//\r
168	// calculate RMS in pad direction\r
169	//\r
170	// TH1F *histo= new TH1F("","",200,0,200);\r
171	TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
172	for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
173	{\r
174	array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetPad(),i);\r
175	//histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetPad() );\r
176	}\r
177	// fPadRMS=histo->GetRMS();\r
178	fPadRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
179	// delete histo;\r
180	delete array;\r
181	return;\r
182	}\r
183	\r
184	void AliTPCclusterKr::SetRowRMS(){\r
185	//\r
186	// calculate RMS in row direction\r
187	//\r
188	TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
189	// TH1F *histo= new TH1F("","",120,0,120);\r
190	for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
191	{\r
192	array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetRow(),i);\r
193	// histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetRow() );\r
194	}\r
195	// fRowRMS=histo->GetRMS();\r
196	fRowRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
197	// delete histo;\r
198	delete array;\r
199	return;\r
200	}\r
201	\r
202	void AliTPCclusterKr::SetTimebinRMS(){\r
203	//\r
204	// calculate RMS in timebin direction\r
205	//\r
206	TArrayI *array= new TArrayI(fCluster->GetEntriesFast());\r
207	// TH1F *histo= new TH1F("","",1000,0,1000);\r
208	for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
209	{\r
210	array->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetTime(),i);\r
211	// histo->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetTime() );\r
212	}\r
213	fTimebinRMS=TMath::RMS(array->GetSize(),array->GetArray());\r
214	//histo->GetRMS();\r
215	// delete histo;\r
216	delete array;\r
217	return;\r
218	}\r
219	\r
220	void AliTPCclusterKr::SetRMS(){\r
221	//\r
222	// calculate RMS in pad,row,timebin direction\r
223	//\r
224	TArrayI *arrayPad = new TArrayI(fCluster->GetEntriesFast());\r
225	TArrayI *arrayRow = new TArrayI(fCluster->GetEntriesFast());\r
226	TArrayI *arrayTime= new TArrayI(fCluster->GetEntriesFast());\r
227	// TH1F *histoPad= new TH1F("p","p",200,0,200);\r
228	// TH1F *histoRow= new TH1F("r","r",120,0,120);\r
229	// TH1F *histoTime= new TH1F("t","t",1000,0,1000);\r
230	for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
231	{\r
232	arrayPad->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetPad(),i);\r
233	arrayRow->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetRow(),i);\r
234	arrayTime->SetAt(((AliTPCvtpr *)(fCluster->At(i)))->GetTime(),i);\r
235	\r
236	//histoPad->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetPad() );\r
237	//histoRow->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetRow() );\r
238	//histoTime->Fill( ((AliTPCvtpr *)(fCluster->At(i)))->GetTime() );\r
239	}\r
240	// fPadRMS=histoPad->GetRMS();\r
241	fPadRMS=TMath::RMS(arrayPad->GetSize(),arrayPad->GetArray());\r
242	fRowRMS=TMath::RMS(arrayRow->GetSize(),arrayRow->GetArray());\r
243	//histoRow->GetRMS();\r
244	fTimebinRMS=TMath::RMS(arrayTime->GetSize(),arrayTime->GetArray());\r
245	//histoTime->GetRMS();\r
246	\r
247	delete arrayPad;\r
248	delete arrayRow;\r
249	delete arrayTime;\r
250	// delete histoPad;\r
251	// delete histoRow;\r
252	// delete histoTime;\r
253	\r
254	return;\r
255	}\r
256	\r
257	\r
258	void AliTPCclusterKr::Set1D(){\r
259	//\r
260	//\r
261	//\r
262	Short_t maxTime=0;\r
263	Short_t minTime=1000;\r
264	Short_t maxPad=0;\r
265	Short_t minPad=1000;\r
266	\r
267	for(Int_t i=0;i<fCluster->GetEntriesFast();i++)\r
268	{\r
269	if(((AliTPCvtpr )(fCluster->At(i)))->GetPad()>maxPad)maxPad =((AliTPCvtpr )(fCluster->At(i)))->GetPad();\r
270	if(((AliTPCvtpr )(fCluster->At(i)))->GetPad()<minPad)minPad =((AliTPCvtpr )(fCluster->At(i)))->GetPad();\r
271	if(((AliTPCvtpr )(fCluster->At(i)))->GetTime()>maxTime)maxTime=((AliTPCvtpr )(fCluster->At(i)))->GetTime();\r
272	if(((AliTPCvtpr )(fCluster->At(i)))->GetTime()<minTime)minTime=((AliTPCvtpr )(fCluster->At(i)))->GetTime();\r
273	}\r
274	fPads1D=maxPad-minPad+1;\r
275	fTimebins1D=maxTime-minTime+1;\r
276	return;\r
277	}\r