Changes in order to save the look-up table (3D) in the OCDB. Otherwise,
[u/mrichter/AliRoot.git] / STEER / AliESDCaloCells.cxx
CommitLineData
e649177a 1/**************************************************************************
2 * Copyright(c) 1998-2007, ALICE Experiment at CERN, All rights reserved. *
3 * *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
6 * *
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 **************************************************************************/
15
16/* $Id$ */
17
18//-------------------------------------------------------------------------
19// ESD class to store calorimeter cell data
20// Clone of AliAODCaloCells made by Markus Oldenburg, CERN
21// Author: Gustavo Conesa Balbastre INFN-LNF
22//-------------------------------------------------------------------------
23
24#include "AliESDCaloCells.h"
25
26ClassImp(AliESDCaloCells)
27
28//_______________________________________________________________________
29AliESDCaloCells::AliESDCaloCells() :
c8fe2783 30 AliVCaloCells(), fNCells(0), fCellNumber(0), fAmplitude(0), fTime(0), fIsSorted(kTRUE), fType(kUndef)
e649177a 31{
32 // default constructor
33}
34//_______________________________________________________________________
c8fe2783 35 AliESDCaloCells::AliESDCaloCells(const char* name, const char* title, VCells_t ttype) :
36 AliVCaloCells(name, title), fNCells(0), fCellNumber(0), fAmplitude(0), fTime(0), fIsSorted(kTRUE), fType(ttype)
37{
38 // AliVCaloCells constructor
e649177a 39 }
40
41//_______________________________________________________________________
42AliESDCaloCells::AliESDCaloCells(const AliESDCaloCells& c) :
c8fe2783 43 AliVCaloCells(c), fNCells(c.fNCells), fCellNumber(0), fAmplitude(0), fTime(0), fIsSorted(c.fIsSorted), fType(c.fType)
e649177a 44{
45 // copy constructor
46
c8fe2783 47 fCellNumber = new Short_t[fNCells];
48 fAmplitude = new Double32_t[fNCells];
49
e649177a 50 for(Int_t i = 0; i < fNCells; i++){
c8fe2783 51 fCellNumber[i] = c.fCellNumber[i];
52 fAmplitude[i] = c.fAmplitude[i];
53 fTime[i] = c.fTime[i];
54 fAmplitude[i] = c.fAmplitude[i];
55
e649177a 56 }
57}
58
59//_______________________________________________________________________
60AliESDCaloCells & AliESDCaloCells::operator =(const AliESDCaloCells& source)
61{
62 // assignment operator
63
64 if(&source == this) return *this;
e649177a 65
732a24fe 66 if(fNCells != source.fNCells){
67 DeleteContainer();
68 CreateContainer(source.fNCells);
69 }
70
e649177a 71 fNCells = source.fNCells;
72 fIsSorted = source.fIsSorted;
73 fType = source.fType;
732a24fe 74
75
76
e649177a 77 for(Int_t i = 0; i < fNCells; i++){
c8fe2783 78 fCellNumber[i] = source.fCellNumber[i];
79 fAmplitude[i] = source.fAmplitude[i];
80 fTime[i] = source.fTime[i];
e649177a 81 }
82
83 return *this;
84
85}
86
732a24fe 87
88void AliESDCaloCells::Copy(TObject &obj) const {
89
90 // this overwrites the virtual TOBject::Copy()
91 // to allow run time copying without casting
92 // in AliESDEvent
93
94 if(this==&obj)return;
95 AliESDCaloCells *robj = dynamic_cast<AliESDCaloCells*>(&obj);
96 if(!robj)return; // not an AliESDCaloCells
97 *robj = *this;
98
99}
100
e649177a 101//_______________________________________________________________________
102AliESDCaloCells::~AliESDCaloCells()
103{
104 // destructor
105
106 DeleteContainer();
107}
108
109//_______________________________________________________________________
110void AliESDCaloCells::CreateContainer(Short_t nCells)
111{
112 // function that creates container to store calorimeter cell data
113
114 DeleteContainer();
115
116 if (nCells <= 0) {
117 fNCells = 0;
118 return;
119 }
120
121 fNCells = nCells;
122
123 fCellNumber = new Short_t[fNCells];
c8fe2783 124 fAmplitude = new Double32_t[fNCells];
125 fTime = new Double32_t[fNCells];
126
127 // set to zero
128 for(int i = 0;i<fNCells;++i){
129 fAmplitude[i] = fCellNumber[i] = 0 ;
130 }
e649177a 131}
132
133//_______________________________________________________________________
134void AliESDCaloCells::DeleteContainer()
135{
136 // deletes allocated memory
137
138 if (fCellNumber)
139 {
140 delete[] fCellNumber;
141 fCellNumber = 0;
142 }
143
144 if (fAmplitude)
145 {
146 delete[] fAmplitude;
c8fe2783 147 fAmplitude = NULL;
e649177a 148 }
149
150 if (fTime)
151 {
152 delete[] fTime;
c8fe2783 153 fTime = NULL;
e649177a 154 }
c8fe2783 155
e649177a 156 fNCells = 0;
157 fIsSorted = kFALSE;
158}
159
160//_______________________________________________________________________
161void AliESDCaloCells::Sort()
162{
163 // sort the cell array by cell number
164
165 Int_t *idxArray = new Int_t[fNCells];
166 TMath::Sort(fNCells,fCellNumber,idxArray,kFALSE);
167
c8fe2783 168 Short_t *newIndex = new Short_t[fNCells];
e649177a 169 Double32_t *newAmplitude = new Double32_t[fNCells];
c8fe2783 170 Double32_t *newTime = new Double32_t[fNCells];
e649177a 171 for (Int_t i=0; i < fNCells; i++) {
c8fe2783 172 newIndex[i] = fCellNumber[idxArray[i]];
e649177a 173 newAmplitude[i] = fAmplitude[idxArray[i]];
c8fe2783 174 newTime[i] = fTime[idxArray[i]];
e649177a 175 }
176 delete [] fCellNumber;
177 delete [] fAmplitude;
178 delete [] fTime;
c8fe2783 179
e649177a 180 fCellNumber = newIndex;
c8fe2783 181 fAmplitude = newAmplitude;
182 fTime = newTime;
e649177a 183
184 delete [] idxArray;
185
186 fIsSorted = kTRUE;
187}
188
189//_______________________________________________________________________
190Bool_t AliESDCaloCells::SetCell(Short_t pos, Short_t cellNumber, Double32_t amplitude, Double32_t time)
191{
192 // Sets a cell at the given position
193
194 if (pos>=0 && pos < fNCells) {
195 fCellNumber[pos] = cellNumber;
c8fe2783 196 fAmplitude[pos] = amplitude;
197 fTime[pos] = time;
198
e649177a 199 fIsSorted = kFALSE;
200 return kTRUE;
201 } else {
202 return kFALSE;
203 }
204}