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 **************************************************************************/
16 /** @file AliFMDEdepMap.cxx
17 @author Christian Holm Christensen <cholm@nbi.dk>
18 @date Mon Mar 27 12:39:50 2006
19 @brief Per strip map of energy deposited and number of hits
22 //____________________________________________________________________
24 // Contains a pair of energy deposited fEdep and number of hits
25 // fN, fEdep is the summed energy deposition, and fN is the
26 // number of hits. The map contains one such object or each strip.
27 // It is used to cache the data in the digitization classes
28 // AliFMDBaseDigitizer and so on.
31 #include "AliFMDEdepMap.h" // ALIFMDEDEPMAP_H
33 //____________________________________________________________________
34 ClassImp(AliFMDEdepMap)
36 ; // This is here to keep Emacs for indenting the next line
39 //____________________________________________________________________
40 AliFMDEdepMap::AliFMDEdepMap(const AliFMDEdepMap& other)
41 : AliFMDMap(other.fMaxDetectors, other.fMaxRings, other.fMaxSectors,
43 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
47 if (fTotal == 0) fTotal = 51200;
48 fData = new AliFMDEdepHitPair[fTotal];
49 for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
53 //____________________________________________________________________
54 AliFMDEdepMap::AliFMDEdepMap()
65 //____________________________________________________________________
66 AliFMDEdepMap::AliFMDEdepMap(UShort_t maxDet,
70 : AliFMDMap(maxDet, maxRing, maxSec, maxStr),
71 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
77 // maxDet Maximum # of detectors
78 // maxRinf Maximum # of rings
79 // maxSec Maximum # of sectors
80 // maxStr Maximum # of strips
81 if (fTotal == 0) fTotal = 51200;
82 fData = new AliFMDEdepHitPair[fTotal];
85 //____________________________________________________________________
87 AliFMDEdepMap::operator=(const AliFMDEdepMap& other)
89 // Assignment operator
90 if (&other == this) return *this;
91 fMaxDetectors = other.fMaxDetectors;
92 fMaxRings = other.fMaxRings;
93 fMaxSectors = other.fMaxSectors;
94 fMaxStrips = other.fMaxStrips;
95 if (fData) delete [] fData;
96 fTotal = fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips;
97 if (fTotal == 0) fTotal = 51200;
98 fData = new AliFMDEdepHitPair[fTotal];
99 for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
103 //____________________________________________________________________
105 AliFMDEdepMap::Reset()
108 for (Int_t i = 0; i < fTotal; i++) {
112 fData[i].fLabels.Reset();
116 //____________________________________________________________________
118 AliFMDEdepMap::Reset(const AliFMDEdepHitPair& val)
121 for (Int_t i = 0; i < fTotal; i++) {
122 fData[i].fEdep = val.fEdep;
123 fData[i].fN = val.fN;
124 fData[i].fNPrim = val.fNPrim;
125 fData[i].fLabels = val.fLabels;
129 //____________________________________________________________________
131 AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
132 UShort_t sec, UShort_t str)
142 // Returns appropriate data
144 return fData[CalcIndex(det, ring, sec, str)];
147 //____________________________________________________________________
148 const AliFMDEdepHitPair&
149 AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
150 UShort_t sec, UShort_t str) const
160 // Returns appropriate data
162 return fData[CalcIndex(det, ring, sec, str)];
166 //___________________________________________________________________