Fix to allow new initialization of AliReconstruction (C.H.Christinsen)
[u/mrichter/AliRoot.git] / FMD / AliFMDEdepMap.cxx
... / ...
CommitLineData
1/**************************************************************************
2 * Copyright(c) 1998-1999, 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/* $Id$ */
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
20 @ingroup FMD_sim
21*/
22//____________________________________________________________________
23//
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.
29//
30//
31#include "AliFMDEdepMap.h" // ALIFMDEDEPMAP_H
32
33//____________________________________________________________________
34ClassImp(AliFMDEdepMap)
35#if 0
36 ; // This is here to keep Emacs for indenting the next line
37#endif
38
39//____________________________________________________________________
40AliFMDEdepMap::AliFMDEdepMap(const AliFMDEdepMap& other)
41 : AliFMDMap(other.fMaxDetectors, other.fMaxRings, other.fMaxSectors,
42 other.fMaxStrips),
43 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
44 fData(0)
45{
46 // Copy constructor
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];
50}
51
52
53//____________________________________________________________________
54AliFMDEdepMap::AliFMDEdepMap()
55 : AliFMDMap(),
56 fTotal(0),
57 fData(0)
58{
59 // Construct a map
60 //
61 // Parameters:
62 // None
63}
64
65//____________________________________________________________________
66AliFMDEdepMap::AliFMDEdepMap(UShort_t maxDet,
67 UShort_t maxRing,
68 UShort_t maxSec,
69 UShort_t maxStr)
70 : AliFMDMap(maxDet, maxRing, maxSec, maxStr),
71 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
72 fData(0)
73{
74 // Construct a map
75 //
76 // Parameters:
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];
83}
84
85//____________________________________________________________________
86AliFMDEdepMap&
87AliFMDEdepMap::operator=(const AliFMDEdepMap& other)
88{
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];
100 return *this;
101}
102
103//____________________________________________________________________
104void
105AliFMDEdepMap::Reset()
106{
107 // Reset to zero
108 for (Int_t i = 0; i < fTotal; i++) {
109 fData[i].fEdep = 0;
110 fData[i].fN = 0;
111 fData[i].fNPrim = 0;
112 fData[i].fLabels.Reset();
113 };
114}
115
116//____________________________________________________________________
117void
118AliFMDEdepMap::Reset(const AliFMDEdepHitPair& val)
119{
120 // Reset to 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;
126 };
127}
128
129//____________________________________________________________________
130AliFMDEdepHitPair&
131AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
132 UShort_t sec, UShort_t str)
133{
134 // Get data
135 //
136 // Parameters:
137 // det Detector #
138 // ring Ring ID
139 // sec Sector #
140 // str Strip #
141 //
142 // Returns appropriate data
143 //
144 return fData[CalcIndex(det, ring, sec, str)];
145}
146
147//____________________________________________________________________
148const AliFMDEdepHitPair&
149AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
150 UShort_t sec, UShort_t str) const
151{
152 // Get data
153 //
154 // Parameters:
155 // det Detector #
156 // ring Ring ID
157 // sec Sector #
158 // str Strip #
159 //
160 // Returns appropriate data
161 //
162 return fData[CalcIndex(det, ring, sec, str)];
163}
164
165
166//___________________________________________________________________
167//
168// EOF
169//