Changes in STEER:
[u/mrichter/AliRoot.git] / FMD / AliFMDEdepMap.cxx
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 //____________________________________________________________________
34 ClassImp(AliFMDEdepMap)
35 #if 0
36   ; // This is here to keep Emacs for indenting the next line
37 #endif
38
39 //____________________________________________________________________
40 AliFMDEdepMap::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 //____________________________________________________________________
54 AliFMDEdepMap::AliFMDEdepMap()
55   : AliFMDMap(), 
56     fTotal(0),
57     fData(0)
58 {
59   // Construct a map
60   //
61   // Parameters:
62   //   None
63 }
64
65 //____________________________________________________________________
66 AliFMDEdepMap::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 //____________________________________________________________________
86 AliFMDEdepMap&
87 AliFMDEdepMap::operator=(const AliFMDEdepMap& other) 
88 {
89   // Assignment operator
90   fMaxDetectors = other.fMaxDetectors;
91   fMaxRings     = other.fMaxRings;
92   fMaxSectors   = other.fMaxSectors;
93   fMaxStrips    = other.fMaxStrips;
94   if (fData) delete [] fData;
95   fTotal = fMaxDetectors * fMaxRings * fMaxSectors * fMaxStrips;
96   if (fTotal == 0) fTotal = 51200;
97   fData  = new AliFMDEdepHitPair[fTotal];
98   for (Int_t i = 0; i < fTotal; i++) fData[i] = other.fData[i];
99   return *this;
100 }
101
102 //____________________________________________________________________
103 void
104 AliFMDEdepMap::Reset() 
105 {
106   // Reset to zero
107   for (Int_t i = 0; i < fTotal; i++) { 
108     fData[i].fEdep  = 0; 
109     fData[i].fN     = 0; 
110     fData[i].fNPrim = 0;
111     fData[i].fLabels.Reset();
112   };
113 }
114
115 //____________________________________________________________________
116 void
117 AliFMDEdepMap::Reset(const AliFMDEdepHitPair& val) 
118 {
119   // Reset to val
120   for (Int_t i = 0; i < fTotal; i++) { 
121     fData[i].fEdep   = val.fEdep; 
122     fData[i].fN      = val.fN; 
123     fData[i].fNPrim  = val.fNPrim;
124     fData[i].fLabels = val.fLabels;
125   };
126 }
127
128 //____________________________________________________________________
129 AliFMDEdepHitPair& 
130 AliFMDEdepMap::operator()(UShort_t det, Char_t ring, 
131                           UShort_t sec, UShort_t str) 
132 {
133   // Get data 
134   // 
135   // Parameters: 
136   //     det       Detector #
137   //     ring      Ring ID
138   //     sec       Sector # 
139   //     str       Strip # 
140   //
141   // Returns appropriate data
142   //
143   return fData[CalcIndex(det, ring, sec, str)];
144 }
145
146 //____________________________________________________________________
147 const AliFMDEdepHitPair& 
148 AliFMDEdepMap::operator()(UShort_t det, Char_t ring, 
149                           UShort_t sec, UShort_t str) const
150 {
151   // Get data 
152   // 
153   // Parameters: 
154   //     det       Detector #
155   //     ring      Ring ID
156   //     sec       Sector # 
157   //     str       Strip # 
158   //
159   // Returns appropriate data
160   //
161   return fData[CalcIndex(det, ring, sec, str)];
162 }
163
164
165 //___________________________________________________________________
166 //
167 // EOF
168 //