add calculation and histograms for MC cross section
[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   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 //____________________________________________________________________
104 void
105 AliFMDEdepMap::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 //____________________________________________________________________
117 void
118 AliFMDEdepMap::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 //____________________________________________________________________
130 AliFMDEdepHitPair& 
131 AliFMDEdepMap::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 //____________________________________________________________________
148 const AliFMDEdepHitPair& 
149 AliFMDEdepMap::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 //