1 /**************************************************************************
2 * Copyright(c) 2004, ALICE Experiment at CERN, All rights reserved. *
4 * Author: The ALICE Off-line Project. *
5 * Contributors are mentioned in the code where appropriate. *
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 **************************************************************************/
18 //____________________________________________________________________
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
26 #include "AliESDFMD.h" // ALIFMDESD_H
27 #include "AliLog.h" // ALILOG_H
28 #include "Riostream.h" // ROOT_Riostream
30 //____________________________________________________________________
33 ; // This is here to keep Emacs for indenting the next line
37 //____________________________________________________________________
38 AliESDFMD::AliESDFMD()
39 : fMultiplicity(AliFMDFloatMap::kMaxDetectors,
40 AliFMDFloatMap::kMaxRings,
41 AliFMDFloatMap::kMaxSectors,
42 AliFMDFloatMap::kMaxStrips),
43 fEta(AliFMDFloatMap::kMaxDetectors,
44 AliFMDFloatMap::kMaxRings,
46 AliFMDFloatMap::kMaxStrips),
48 fAngleCorrected(kFALSE)
53 //____________________________________________________________________
54 AliESDFMD::AliESDFMD(const AliESDFMD& other)
56 fMultiplicity(other.fMultiplicity),
58 fNoiseFactor(other.fNoiseFactor),
59 fAngleCorrected(other.fAngleCorrected)
64 //____________________________________________________________________
66 AliESDFMD::operator=(const AliESDFMD& other)
70 TObject::operator=(other);
71 fMultiplicity = other.fMultiplicity;
77 void AliESDFMD::Copy(TObject &obj) const{
78 // this overwrites the virtual TOBject::Copy()
79 // to allow run time copying without casting
83 AliESDFMD *robj = dynamic_cast<AliESDFMD*>(&obj);
84 if(!robj)return; // not an AliESDFMD
88 //____________________________________________________________________
90 AliESDFMD::CheckNeedUShort(TFile* file)
92 fMultiplicity.CheckNeedUShort(file);
93 fEta.CheckNeedUShort(file);
96 //____________________________________________________________________
98 AliESDFMD::Clear(Option_t* )
100 fMultiplicity.Reset(kInvalidMult);
101 fEta.Reset(kInvalidEta);
105 //____________________________________________________________________
107 AliESDFMD::Multiplicity(UShort_t detector, Char_t ring, UShort_t sector,
108 UShort_t strip) const
110 // Return rough estimate of charged particle multiplicity in the
111 // strip FMD<detector><ring>[<sector>,<strip>].
113 // Note, that this should at most be interpreted as the sum
114 // multiplicity of secondaries and primaries.
115 return fMultiplicity(detector, ring, sector, strip);
118 //____________________________________________________________________
120 AliESDFMD::Eta(UShort_t detector, Char_t ring, UShort_t /* sector */,
121 UShort_t strip) const
123 // Return pseudo-rapidity of the strip
124 // FMD<detector><ring>[<sector>,<strip>]. (actually, the sector
125 // argument is ignored, as it is assumed that the primary vertex is
126 // a (x,y) = (0,0), and that the modules are aligned with a
127 // precision better than 2 degrees in the azimuthal angle).
129 return fEta(detector, ring, 0, strip);
132 //____________________________________________________________________
134 AliESDFMD::SetMultiplicity(UShort_t detector, Char_t ring, UShort_t sector,
135 UShort_t strip, Float_t mult)
137 // Return rough estimate of charged particle multiplicity in the
138 // strip FMD<detector><ring>[<sector>,<strip>].
140 // Note, that this should at most be interpreted as the sum
141 // multiplicity of secondaries and primaries.
142 fMultiplicity(detector, ring, sector, strip) = mult;
145 //____________________________________________________________________
147 AliESDFMD::SetEta(UShort_t detector, Char_t ring, UShort_t /* sector */,
148 UShort_t strip, Float_t eta)
150 // Set pseudo-rapidity of the strip
151 // FMD<detector><ring>[<sector>,<strip>]. (actually, the sector
152 // argument is ignored, as it is assumed that the primary vertex is
153 // a (x,y) = (0,0), and that the modules are aligned with a
154 // precision better than 2 degrees in the azimuthal angle).
156 fEta(detector, ring, 0, strip) = eta;
159 //____________________________________________________________________
161 AliESDFMD::Print(Option_t* /* option*/) const
163 // Print all information to standard output.
164 std::cout << "AliESDFMD:" << std::endl;
165 for (UShort_t det = 1; det <= fMultiplicity.MaxDetectors(); det++) {
166 for (UShort_t ir = 0; ir < fMultiplicity.MaxRings(); ir++) {
167 Char_t ring = (ir == 0 ? 'I' : 'O');
168 std::cout << "FMD" << det << ring << ":" << std::endl;
169 for (UShort_t sec = 0; sec < fMultiplicity.MaxSectors(); sec++) {
170 std::cout << " Sector # " << sec << ":" << std::flush;
171 for (UShort_t str = 0; str < fMultiplicity.MaxStrips(); str++) {
172 if (str % 6 == 0) std::cout << "\n " << std::flush;
173 Float_t m = fMultiplicity(det, ring, sec, str);
174 Float_t e = fEta(det, ring, 0, str);
175 if (m == kInvalidMult && e == kInvalidEta) break;
176 if (m == kInvalidMult) std::cout << " ---- ";
177 else std::cout << Form("%6.3f", m);
179 if (e == kInvalidEta) std::cout << " ---- ";
180 else std::cout << Form("%-6.3f", e);
181 std::cout << " " << std::flush;
183 std::cout << std::endl;
190 //____________________________________________________________________