For primary tracks, when mother-particle is defined, add label and
[u/mrichter/AliRoot.git] / STEER / AliESDFMD.cxx
CommitLineData
9da38871 1/**************************************************************************
2 * Copyright(c) 2004, 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
16/* $Id$ */
17
18//____________________________________________________________________
19//
20// ESD information from the FMD
21// Contains information on:
22// Charged particle multiplicty per strip (rough estimate)
23// Psuedo-rapdity per strip
24// Latest changes by Christian Holm Christensen
25//
26#include "AliESDFMD.h" // ALIFMDESD_H
27#include "AliLog.h" // ALILOG_H
28#include "Riostream.h" // ROOT_Riostream
29
30//____________________________________________________________________
31ClassImp(AliESDFMD)
32#if 0
33 ; // This is here to keep Emacs for indenting the next line
34#endif
35
36
37//____________________________________________________________________
38AliESDFMD::AliESDFMD()
39 : fMultiplicity(AliFMDFloatMap::kMaxDetectors,
40 AliFMDFloatMap::kMaxRings,
41 AliFMDFloatMap::kMaxSectors,
42 AliFMDFloatMap::kMaxStrips),
43 fEta(AliFMDFloatMap::kMaxDetectors,
44 AliFMDFloatMap::kMaxRings,
45 1,
c05d076f 46 AliFMDFloatMap::kMaxStrips),
0cd61c1d 47 fNoiseFactor(0),
c05d076f 48 fAngleCorrected(kFALSE)
9da38871 49{
50 // Default CTOR
51}
52
53//____________________________________________________________________
54AliESDFMD::AliESDFMD(const AliESDFMD& other)
55 : TObject(other),
56 fMultiplicity(other.fMultiplicity),
0cd61c1d 57 fEta(other.fEta),
58 fNoiseFactor(other.fNoiseFactor),
59 fAngleCorrected(other.fAngleCorrected)
9da38871 60{
61 // Default CTOR
62}
63
64//____________________________________________________________________
65AliESDFMD&
66AliESDFMD::operator=(const AliESDFMD& other)
67{
68 // Default CTOR
316c6cd9 69 if(this!=&other){
70 TObject::operator=(other);
71 fMultiplicity = other.fMultiplicity;
72 fEta = other.fEta;
73 }
9da38871 74 return *this;
75}
76
77//____________________________________________________________________
78void
c05d076f 79AliESDFMD::CheckNeedUShort(TFile* file)
80{
81 fMultiplicity.CheckNeedUShort(file);
82 fEta.CheckNeedUShort(file);
83}
84
85//____________________________________________________________________
86void
9da38871 87AliESDFMD::Clear(Option_t* )
88{
89 fMultiplicity.Reset(kInvalidMult);
90 fEta.Reset(kInvalidEta);
91}
92
93
94//____________________________________________________________________
95Float_t
96AliESDFMD::Multiplicity(UShort_t detector, Char_t ring, UShort_t sector,
97 UShort_t strip) const
98{
99 // Return rough estimate of charged particle multiplicity in the
100 // strip FMD<detector><ring>[<sector>,<strip>].
101 //
102 // Note, that this should at most be interpreted as the sum
103 // multiplicity of secondaries and primaries.
104 return fMultiplicity(detector, ring, sector, strip);
105}
106
107//____________________________________________________________________
108Float_t
109AliESDFMD::Eta(UShort_t detector, Char_t ring, UShort_t /* sector */,
110 UShort_t strip) const
111{
112 // Return pseudo-rapidity of the strip
113 // FMD<detector><ring>[<sector>,<strip>]. (actually, the sector
114 // argument is ignored, as it is assumed that the primary vertex is
115 // a (x,y) = (0,0), and that the modules are aligned with a
116 // precision better than 2 degrees in the azimuthal angle).
117 //
118 return fEta(detector, ring, 0, strip);
119}
120
121//____________________________________________________________________
122void
123AliESDFMD::SetMultiplicity(UShort_t detector, Char_t ring, UShort_t sector,
124 UShort_t strip, Float_t mult)
125{
126 // Return rough estimate of charged particle multiplicity in the
127 // strip FMD<detector><ring>[<sector>,<strip>].
128 //
129 // Note, that this should at most be interpreted as the sum
130 // multiplicity of secondaries and primaries.
131 fMultiplicity(detector, ring, sector, strip) = mult;
132}
133
134//____________________________________________________________________
135void
136AliESDFMD::SetEta(UShort_t detector, Char_t ring, UShort_t /* sector */,
137 UShort_t strip, Float_t eta)
138{
139 // Set pseudo-rapidity of the strip
140 // FMD<detector><ring>[<sector>,<strip>]. (actually, the sector
141 // argument is ignored, as it is assumed that the primary vertex is
142 // a (x,y) = (0,0), and that the modules are aligned with a
143 // precision better than 2 degrees in the azimuthal angle).
144 //
145 fEta(detector, ring, 0, strip) = eta;
146}
147
148//____________________________________________________________________
149void
150AliESDFMD::Print(Option_t* /* option*/) const
151{
6169f936 152 // Print all information to standard output.
9da38871 153 std::cout << "AliESDFMD:" << std::endl;
9eeb02aa 154 for (UShort_t det = 1; det <= fMultiplicity.MaxDetectors(); det++) {
155 for (UShort_t ir = 0; ir < fMultiplicity.MaxRings(); ir++) {
9da38871 156 Char_t ring = (ir == 0 ? 'I' : 'O');
157 std::cout << "FMD" << det << ring << ":" << std::endl;
9eeb02aa 158 for (UShort_t sec = 0; sec < fMultiplicity.MaxSectors(); sec++) {
9da38871 159 std::cout << " Sector # " << sec << ":" << std::flush;
9eeb02aa 160 for (UShort_t str = 0; str < fMultiplicity.MaxStrips(); str++) {
9da38871 161 if (str % 6 == 0) std::cout << "\n " << std::flush;
162 Float_t m = fMultiplicity(det, ring, sec, str);
163 Float_t e = fEta(det, ring, 0, str);
164 if (m == kInvalidMult && e == kInvalidEta) break;
165 if (m == kInvalidMult) std::cout << " ---- ";
166 else std::cout << Form("%6.3f", m);
167 std::cout << "/";
168 if (e == kInvalidEta) std::cout << " ---- ";
169 else std::cout << Form("%-6.3f", e);
170 std::cout << " " << std::flush;
171 }
172 std::cout << std::endl;
173 }
174 }
175 }
176}
177
178
179//____________________________________________________________________
180//
181// EOF
182//