Fix to allow new initialization of AliReconstruction (C.H.Christinsen)
[u/mrichter/AliRoot.git] / FMD / AliFMDEdepMap.cxx
CommitLineData
e802be3e 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 **************************************************************************/
e802be3e 15/* $Id$ */
c2fc1258 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
02a27b50 20 @ingroup FMD_sim
c2fc1258 21*/
e802be3e 22//____________________________________________________________________
23//
02a27b50 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.
e802be3e 29//
30//
31#include "AliFMDEdepMap.h" // ALIFMDEDEPMAP_H
32
33//____________________________________________________________________
925e6570 34ClassImp(AliFMDEdepMap)
1a1fdef7 35#if 0
36 ; // This is here to keep Emacs for indenting the next line
37#endif
e802be3e 38
39//____________________________________________________________________
40AliFMDEdepMap::AliFMDEdepMap(const AliFMDEdepMap& other)
41 : AliFMDMap(other.fMaxDetectors, other.fMaxRings, other.fMaxSectors,
42 other.fMaxStrips),
021f1396 43 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
e802be3e 44 fData(0)
45{
088f8e79 46 // Copy constructor
021f1396 47 if (fTotal == 0) fTotal = 51200;
bfdc7f5d 48 fData = new AliFMDEdepHitPair[fTotal];
0ed9abeb 49 for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
e802be3e 50}
51
52
021f1396 53//____________________________________________________________________
54AliFMDEdepMap::AliFMDEdepMap()
55 : AliFMDMap(),
56 fTotal(0),
57 fData(0)
58{
59 // Construct a map
60 //
61 // Parameters:
62 // None
63}
e802be3e 64
65//____________________________________________________________________
6e79feeb 66AliFMDEdepMap::AliFMDEdepMap(UShort_t maxDet,
67 UShort_t maxRing,
68 UShort_t maxSec,
69 UShort_t maxStr)
e802be3e 70 : AliFMDMap(maxDet, maxRing, maxSec, maxStr),
021f1396 71 fTotal(fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips),
e802be3e 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
021f1396 81 if (fTotal == 0) fTotal = 51200;
bfdc7f5d 82 fData = new AliFMDEdepHitPair[fTotal];
e802be3e 83}
84
85//____________________________________________________________________
86AliFMDEdepMap&
87AliFMDEdepMap::operator=(const AliFMDEdepMap& other)
88{
088f8e79 89 // Assignment operator
d015ecfe 90 if (&other == this) return *this;
e802be3e 91 fMaxDetectors = other.fMaxDetectors;
92 fMaxRings = other.fMaxRings;
93 fMaxSectors = other.fMaxSectors;
94 fMaxStrips = other.fMaxStrips;
95 if (fData) delete [] fData;
bfdc7f5d 96 fTotal = fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips;
021f1396 97 if (fTotal == 0) fTotal = 51200;
bfdc7f5d 98 fData = new AliFMDEdepHitPair[fTotal];
0ed9abeb 99 for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
e802be3e 100 return *this;
101}
102
103//____________________________________________________________________
104void
088f8e79 105AliFMDEdepMap::Reset()
106{
107 // Reset to zero
0ed9abeb 108 for (Int_t i = 0; i < fTotal; i++) {
83ad576a 109 fData[i].fEdep = 0;
110 fData[i].fN = 0;
111 fData[i].fNPrim = 0;
8d00dfa3 112 fData[i].fLabels.Reset();
bfdc7f5d 113 };
088f8e79 114}
115
116//____________________________________________________________________
117void
69b696b9 118AliFMDEdepMap::Reset(const AliFMDEdepHitPair& val)
e802be3e 119{
088f8e79 120 // Reset to val
0ed9abeb 121 for (Int_t i = 0; i < fTotal; i++) {
8d00dfa3 122 fData[i].fEdep = val.fEdep;
123 fData[i].fN = val.fN;
124 fData[i].fNPrim = val.fNPrim;
125 fData[i].fLabels = val.fLabels;
bfdc7f5d 126 };
e802be3e 127}
128
129//____________________________________________________________________
130AliFMDEdepHitPair&
83ad576a 131AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
132 UShort_t sec, UShort_t str)
e802be3e 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&
83ad576a 149AliFMDEdepMap::operator()(UShort_t det, Char_t ring,
150 UShort_t sec, UShort_t str) const
e802be3e 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//