New macro to keep track of timing performances of the segmentation methods (Laurent)
[u/mrichter/AliRoot.git] / MUON / mapping / AliMpSt345Reader.cxx
CommitLineData
dee1d5f1 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 purpeateose. It is *
13* provided "as is" without express or implied warranty. *
14**************************************************************************/
15
16// $Id$
13985652 17// $MpId: AliMpSt345Reader.cxx,v 1.11 2006/05/24 13:58:50 ivana Exp $
dee1d5f1 18
19#include "AliMpSt345Reader.h"
20
21#include "AliLog.h"
d5c9a22c 22#include "AliMpSlatMotifMap.h"
dee1d5f1 23#include "AliMpMotifReader.h"
24#include "AliMpFiles.h"
228fd720 25#include "AliMpDataStreams.h"
dee1d5f1 26#include "AliMpMotifType.h"
27#include "AliMpPCB.h"
28#include "AliMpSlat.h"
dee1d5f1 29#include "AliMpMotifPosition.h"
30#include "AliMpMotif.h"
31#include "AliMpHelper.h"
f0eac1a3 32#include "AliMpConstants.h"
dee1d5f1 33
34#include "Riostream.h"
dee1d5f1 35#include "TClass.h"
36#include "TObjString.h"
37#include "TString.h"
38
39#include <sstream>
40
3d1463c8 41
42//-----------------------------------------------------------------------------
85fec35d 43/// \class AliMpSt345Reader
44//
45/// Read slat and pcb ASCII files.
46///
d5c9a22c 47/// Basically this class provides two methods :
85fec35d 48/// - AliMpSlat* ReadSlat()
49/// - AliMpPCB ReadPCB()
50///
51/// \author Laurent Aphecetche
3d1463c8 52//-----------------------------------------------------------------------------
dee1d5f1 53
13985652 54/// \cond CLASSIMP
f0eac1a3 55ClassImp(AliMpSt345Reader)
13985652 56/// \endcond
f0eac1a3 57
dee1d5f1 58//_____________________________________________________________________________
183279c1 59AliMpSt345Reader::AliMpSt345Reader(const AliMpDataStreams& dataStreams, AliMpSlatMotifMap* motifMap)
d5c9a22c 60:
61TObject(),
7d5d0cc5 62fkDataStreams(dataStreams),
183279c1 63fMotifMap(motifMap)
dee1d5f1 64{
71a2d3aa 65 ///
66 /// Default ctor.
67 ///
dee1d5f1 68}
69
70//_____________________________________________________________________________
71AliMpSt345Reader::~AliMpSt345Reader()
72{
71a2d3aa 73 ///
74 /// Dtor.
75 ///
dee1d5f1 76}
77
78//_____________________________________________________________________________
79AliMpPCB*
80AliMpSt345Reader::ReadPCB(const char* pcbType)
81{
d5c9a22c 82 ///
83 /// Create a new AliMpPCB object, by reading it from file.
84 /// The returned object must be deleted by the client
dee1d5f1 85
228fd720 86 istream& in
7d5d0cc5 87 = fkDataStreams.
ab167304 88 CreateDataStream(AliMpFiles::SlatPCBFilePath(
228fd720 89 AliMp::kStation345, pcbType));
dee1d5f1 90
7d5d0cc5 91 AliMpMotifReader reader(fkDataStreams,
4e51cfd2 92 AliMp::kStation345, AliMq::kNotSt12, AliMp::kNonBendingPlane);
dee1d5f1 93 // note that the nonbending
94 // parameter is of no use for station345, as far as reading motif is
95 // concerned, as all motifs are supposed to be in the same directory
96 // (as they are shared by bending/non-bending planes).
97
98 char line[80];
99
84aac932 100 const TString kSizeKeyword("SIZES");
101 const TString kMotifKeyword("MOTIF");
dee1d5f1 102
103 AliMpPCB* pcb = 0;
104
105 while ( in.getline(line,80) )
106 {
107 if ( line[0] == '#' ) continue;
108
109 TString sline(line);
110
84aac932 111 if ( sline(0,kSizeKeyword.Length()) == kSizeKeyword )
dee1d5f1 112 {
84aac932 113 std::istringstream sin(sline(kSizeKeyword.Length(),
114 sline.Length()-kSizeKeyword.Length()).Data());
d5c9a22c 115 double padSizeX = 0.0;
116 double padSizeY = 0.0;
117 double pcbSizeX = 0.0;
118 double pcbSizeY = 0.0;
dee1d5f1 119 sin >> padSizeX >> padSizeY >> pcbSizeX >> pcbSizeY;
c918393d 120 if (pcb)
121 {
122 AliError("pcb not null as expected");
123 }
657b96e3 124 pcb = new AliMpPCB(fMotifMap,pcbType,padSizeX,padSizeY,pcbSizeX,pcbSizeY);
dee1d5f1 125 }
126
84aac932 127 if ( sline(0,kMotifKeyword.Length()) == kMotifKeyword )
dee1d5f1 128 {
84aac932 129 std::istringstream sin(sline(kMotifKeyword.Length(),
130 sline.Length()-kMotifKeyword.Length()).Data());
dee1d5f1 131 TString sMotifType;
132 int ix;
133 int iy;
134 sin >> sMotifType >> ix >> iy;
135
657b96e3 136 AliMpMotifType* motifType = fMotifMap->FindMotifType(sMotifType);
d5c9a22c 137 if (!motifType)
138 {
139 AliDebug(1,Form("Reading motifType %s from file",sMotifType.Data()));
140 motifType = reader.BuildMotifType(sMotifType.Data());
657b96e3 141 fMotifMap->AddMotifType(motifType);
d5c9a22c 142 }
143 else
144 {
145 AliDebug(1,Form("Got motifType %s from motifMap",sMotifType.Data()));
146 }
dee1d5f1 147
dee1d5f1 148 pcb->Add(motifType,ix,iy);
149 }
150 }
151
228fd720 152 delete &in;
dee1d5f1 153
dee1d5f1 154 return pcb;
155}
156
157//_____________________________________________________________________________
158AliMpSlat*
cddd101e 159AliMpSt345Reader::ReadSlat(const char* slatType, AliMp::PlaneType planeType)
dee1d5f1 160{
d5c9a22c 161 ///
162 /// Create a new AliMpSlat object, by reading it from file.
163 /// The returned object must be deleted by the client.
dee1d5f1 164
228fd720 165 istream& in
7d5d0cc5 166 = fkDataStreams.
ab167304 167 CreateDataStream(AliMpFiles::SlatFilePath(
228fd720 168 AliMp::kStation345, slatType, planeType));
169
dee1d5f1 170 char line[80];
171
84aac932 172 const TString kpcbKeyword("PCB");
dee1d5f1 173
aad465b5 174 AliMpSlat* slat = new AliMpSlat(slatType, planeType);
dee1d5f1 175
176 while ( in.getline(line,80) )
177 {
178 if ( line[0] == '#' ) continue;
179
180 TString sline(AliMpHelper::Normalize(line));
181
84aac932 182 if ( sline(0,kpcbKeyword.Length()) == kpcbKeyword )
dee1d5f1 183 {
84aac932 184 TString tmp(sline(kpcbKeyword.Length()+1,sline.Length()-kpcbKeyword.Length()));
dee1d5f1 185 Ssiz_t blankPos = tmp.First(' ');
186 if ( blankPos < 0 )
187 {
188 AliErrorClass("Syntax error in PCB file, should get a list of "
189 "manu ids after the pcbname");
190 delete slat;
191 return 0;
192 }
193
194 TString pcbName(tmp(0,blankPos));
195 TString manus(tmp(blankPos+1,tmp.Length()-blankPos));
196
d5c9a22c 197 AliMpPCB* pcbType = ReadPCB(pcbName.Data());
dee1d5f1 198 if (!pcbType)
199 {
200 AliErrorClass(Form("Cannot read pcbType=%s",pcbName.Data()));
201 delete slat;
202 return 0;
203 }
204
205 TArrayI manuList;
206 AliMpHelper::DecodeName(manus,';',manuList);
c6e05ab2 207 if ( manuList.GetSize() != Int_t(pcbType->GetSize()) )
dee1d5f1 208 {
209 AliErrorClass(Form("Wrong number of manu ids for this PCB ("
210 "%s) : %d out of %d",pcbName.Data(),
211 manuList.GetSize(),pcbType->GetSize()));
d5c9a22c 212 delete pcbType;
dee1d5f1 213 delete slat;
214 return 0;
215 }
f0eac1a3 216
217 for ( Int_t i = 0; i < manuList.GetSize(); ++i )
218 {
219 manuList[i] |= AliMpConstants::ManuMask(planeType);
220 }
d5c9a22c 221 slat->Add(*pcbType,manuList);
222 delete pcbType;
dee1d5f1 223 }
224 }
225
228fd720 226 delete &in;
dee1d5f1 227
228 return slat;
228fd720 229}
230
dee1d5f1 231