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 **************************************************************************/
17 //======================================================================
19 // Fill Unit Array class
20 // Base class used by AliJetESDReader to fill a UnitArray from the information extracted
21 // from the particle tracks or emcal cells
22 // Author: magali.estienne@subatech.in2p3.fr
23 //======================================================================
26 #include "AliJetFillUnitArray.h"
28 // --- ROOT system ---
35 // --- AliRoot header files ---
38 class AliJetESDReader;
39 class AliJetESDReaderHeader;
40 class AliJetUnitArray;
42 ClassImp(AliJetFillUnitArray)
44 //_____________________________________________________________________________
45 AliJetFillUnitArray::AliJetFillUnitArray()
46 : TTask("AliJetFillUnitArray","Fill Unit Array with tpc/its and emcal information"),
83 AliJetFillUnitArray::AliJetFillUnitArray(const AliJetFillUnitArray& cpfrom)
84 : TTask("AliJetFillUnitArray","Fill Unit Array with tpc/its and emcal information"),
114 fEtaBinInEMCalAcc(0),
115 fPhiBinInEMCalAcc(0),
121 fNTracks = cpfrom.fNTracks;
122 fNTracksCut = cpfrom.fNTracksCut;
125 fDebug = cpfrom.fDebug;
126 fReaderHeader = cpfrom.fReaderHeader;
127 fMomentumArray = cpfrom.fMomentumArray;
128 fUnitArray = cpfrom.fUnitArray;
129 fRefArray = cpfrom.fRefArray;
131 fSignalFlag = cpfrom.fSignalFlag;
132 fCutFlag = cpfrom.fCutFlag;
133 fProcId = cpfrom.fProcId;
134 fTPCGrid = cpfrom.fTPCGrid;
135 fEMCalGrid = cpfrom.fEMCalGrid;
136 fGeom = cpfrom.fGeom;
137 fNphi = cpfrom.fNphi;
138 fNeta = cpfrom.fNeta;
139 fGrid = cpfrom.fGrid;
140 fPhi2 = cpfrom.fPhi2;
141 fEta2 = cpfrom.fEta2;
142 fIndex = cpfrom.fIndex;
143 fParams = cpfrom.fParams;
144 fPhiMin = cpfrom.fPhiMin;
145 fPhiMax = cpfrom.fPhiMax;
146 fEtaMin = cpfrom.fEtaMin;
147 fEtaMax = cpfrom.fEtaMax;
148 fEtaBinInTPCAcc = cpfrom.fEtaBinInTPCAcc;
149 fPhiBinInTPCAcc = cpfrom.fPhiBinInTPCAcc;
150 fEtaBinInEMCalAcc = cpfrom.fEtaBinInEMCalAcc;
151 fPhiBinInEMCalAcc = cpfrom.fPhiBinInEMCalAcc;
152 fNbinPhi = cpfrom.fNbinPhi;
155 AliJetFillUnitArray& AliJetFillUnitArray::operator=(const AliJetFillUnitArray& rhs)
158 // Assignment operator
161 fNTracks = rhs.fNTracks;
162 fNTracksCut = rhs.fNTracksCut;
166 fReaderHeader = rhs.fReaderHeader;
167 fMomentumArray = rhs.fMomentumArray;
168 fUnitArray = rhs.fUnitArray;
169 fRefArray = rhs.fRefArray;
171 fSignalFlag = rhs.fSignalFlag;
172 fCutFlag = rhs.fCutFlag;
173 fProcId = rhs.fProcId;
174 fTPCGrid = rhs.fTPCGrid;
175 fEMCalGrid = rhs.fEMCalGrid;
183 fParams = rhs.fParams;
184 fPhiMin = rhs.fPhiMin;
185 fPhiMax = rhs.fPhiMax;
186 fEtaMin = rhs.fEtaMin;
187 fEtaMax = rhs.fEtaMax;
188 fEtaBinInTPCAcc = rhs.fEtaBinInTPCAcc;
189 fPhiBinInTPCAcc = rhs.fPhiBinInTPCAcc;
190 fEtaBinInEMCalAcc = rhs.fEtaBinInEMCalAcc;
191 fPhiBinInEMCalAcc = rhs.fPhiBinInEMCalAcc;
192 fNbinPhi = rhs.fNbinPhi;
197 //_____________________________________________________________________________
198 AliJetFillUnitArray::~AliJetFillUnitArray()
203 //_____________________________________________________________________________
204 void AliJetFillUnitArray::GetEtaPhiFromIndex(Int_t index, Float_t &eta, Float_t &phi)
206 // Get the eta,phi position from the index
208 for(Int_t j=0; j<fNphi+1; j++) {
209 for(Int_t i=0; i<fNeta+1; i++) {
212 //-------------------------------------
214 if(j*(fNeta+1)+i == index) {
221 //-------------------------------------
224 if(i>0 && i<(fEtaBinInTPCAcc-fEtaBinInEMCalAcc)/2) ii = i;
225 if(i>=(fEtaBinInTPCAcc+fEtaBinInEMCalAcc)/2 && i<fNeta+1) ii = i-fEtaBinInEMCalAcc;
228 if(j<(fNbinPhi+1) && j*(fNeta+1)+i == index) {
233 if((j>=(fNbinPhi+1) && j<(fNbinPhi+1+fPhiBinInEMCalAcc)) &&
234 ((fNbinPhi+1)*(fNeta+1) + (j-fNbinPhi-1)*(fEtaBinInTPCAcc-fEtaBinInEMCalAcc) + ii)== index ) {
235 if(ii==0) {Int_t ind = 0; eta = fEta2->At(ind);}
236 else eta = fEta2->At(i);
240 if(j>=(fNbinPhi+1+fPhiBinInEMCalAcc) && ((fNbinPhi+1)*(fNeta+1)+fPhiBinInEMCalAcc*((fEtaBinInTPCAcc-fEtaBinInEMCalAcc))+(j-(fNbinPhi+1+fPhiBinInEMCalAcc))*(fNeta+1)+i == index)) {
249 //_____________________________________________________________________________
250 Float_t AliJetFillUnitArray::EtaToTheta(Float_t arg)
252 return 2.*atan(exp(-arg));