]> git.uio.no Git - u/mrichter/AliRoot.git/blame - MUON/AliMUONPainterHelper.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
effc++ warnings
[u/mrichter/AliRoot.git] / MUON / AliMUONPainterHelper.cxx
CommitLineData
0145e89a 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
16// $Id$
17
18#include "AliMUONPainterHelper.h"
19
20#include "AliMUONGeometryDetElement.h"
21#include "AliMUONGeometryTransformer.h"
22#include "AliMUONPainterContour.h"
23#include "AliMUONPainterContourMaker.h"
24#include "AliMUONPainterEnv.h"
25#include "AliMUONPainterMatrix.h"
26#include "AliMUONPainterPadStore.h"
27#include "AliMUONPainterRegistry.h"
28#include "AliMUONVCalibParam.h"
29#include "AliMUONVDigit.h"
30#include "AliMUONVTrackerData.h"
31#include "AliMpCDB.h"
32#include "AliMpConstants.h"
33#include "AliMpDDLStore.h"
34#include "AliMpDEIterator.h"
35#include "AliMpDEManager.h"
36#include "AliMpExMap.h"
37#include "AliMpMotifMap.h"
38#include "AliMpMotifPosition.h"
39#include "AliMpPCB.h"
40#include "AliMpPad.h"
41#include "AliMpSector.h"
42#include "AliMpSectorSegmentation.h"
43#include "AliMpSegmentation.h"
44#include "AliMpSlat.h"
45#include "AliMpSlatSegmentation.h"
46#include "AliMpStationType.h"
47#include "AliMpVPadIterator.h"
48#include "AliCodeTimer.h"
49#include "AliLog.h"
50#include <Riostream.h>
51#include <TArrayI.h>
52#include <TCanvas.h>
53#include <TClass.h>
54#include <TCollection.h>
55#include <TFile.h>
56#include <TGLabel.h>
57#include <TGeoMatrix.h>
58#include <TGMsgBox.h>
59#include <TLine.h>
60#include <TList.h>
61#include <TMap.h>
62#include <TObjArray.h>
63#include <TObjString.h>
64#include <TStyle.h>
65#include <TSystem.h>
66#include <TVector2.h>
67#include <TVector3.h>
68#include <TVirtualPad.h>
69#include <TVirtualX.h>
70
71///\class AliMUONPainterHelper
72///
73/// Helper class for painters
74///
75///\author Laurent Aphecetche, Subatech
76
77///\cond CLASSIMP
78ClassImp(AliMUONPainterHelper)
79///\endcond
80
81AliMUONPainterHelper* AliMUONPainterHelper::fgInstance(0x0);
82
83//_____________________________________________________________________________
84AliMUONPainterHelper::AliMUONPainterHelper() :
85 TObject(),
86 fPadStore(0x0),
87 fExplodedGlobalTransformations(0x0),
88 fRealGlobalTransformations(0x0),
89 fIsModified(kFALSE),
90 fContourMaker(0x0),
91 fPainterMatrices(0x0),
92 fEnv(0x0)
93{
94 /// ctor
95 fExplodeFactor[0] = 1.00;
96 fExplodeFactor[1] = 1.50;
97
98 if ( ! AliMpCDB::LoadMpSegmentation() )
99 {
100 AliFatal("Could not access mapping from OCDB !");
101 }
102
103 // Load DDL store
104 if ( ! AliMpCDB::LoadDDLStore() )
105 {
106 AliFatal("Could not access DDL Store from OCDB !");
107 }
108}
109
110//_____________________________________________________________________________
111AliMUONPainterHelper::~AliMUONPainterHelper()
112{
113 /// dtor
114 if ( fIsModified ) Save();
115 delete fExplodedGlobalTransformations;
116 delete fRealGlobalTransformations;
117 delete fPadStore;
118 delete fContourMaker;
119 delete fPainterMatrices;
120 fgInstance = 0;
121}
122
123//_____________________________________________________________________________
124AliMUONPainterContour*
125AliMUONPainterHelper::GetContour(const char* contourName) const
126{
127 /// Get a contour by name
128
129 AliCodeTimerAuto("")
130
131 if ( fContourMaker )
132 {
133 return fContourMaker->GetContour(contourName);
134 }
135 return 0x0;
136}
137
138//_____________________________________________________________________________
139Int_t
140AliMUONPainterHelper::FindPadID(const TArrayI& pads, Double_t x, Double_t y) const
141{
142 /// Find a pad by position
143
144 return fPadStore->FindPadID(pads,x,y);
145}
146
147//_____________________________________________________________________________
148void
149AliMUONPainterHelper::GenerateDefaultMatrices()
150{
151 /// Kind of bootstrap method to trigger the generation of all contours
152
153 fPainterMatrices = new TObjArray;
154 fPainterMatrices->SetOwner(kFALSE);
155
156 TObjArray attributes;
157
158 AliMUONAttPainter att;
159
160 att.SetViewPoint(kTRUE,kFALSE);
161 att.SetPlane(kFALSE,kFALSE);
162 att.SetCathode(kTRUE,kFALSE);
163
164 AliWarningClass("Should generate back views as well here");
165
166 attributes.Add(new AliMUONAttPainter(att));
167 att.SetCathode(kFALSE,kTRUE);
168 attributes.Add(new AliMUONAttPainter(att));
169 att.SetCathode(kFALSE,kFALSE);
170 att.SetPlane(kTRUE,kFALSE);
171 attributes.Add(new AliMUONAttPainter(att));
172 att.SetPlane(kFALSE,kTRUE);
173 attributes.Add(new AliMUONAttPainter(att));
174
175 TIter next(&attributes);
176 AliMUONAttPainter* a;
177
178 while ( ( a = static_cast<AliMUONAttPainter*>(next()) ) )
179 {
180 AliMUONPainterMatrix* matrix = new AliMUONPainterMatrix("Tracker",5,2);
181
182 for ( Int_t i = 0; i < 10; ++i )
183 {
184 AliMUONVPainter* painter = AliMUONVPainter::CreatePainter("AliMUONChamberPainter",*a,i,-1);
185
186 painter->SetResponder("Chamber");
187
188 painter->SetOutlined("*",kFALSE);
189
190 painter->SetOutlined("MANU",kTRUE);
191
192 for ( Int_t j = 0; j < 3; ++j )
193 {
194 painter->SetLine(j,1,4-j);
195 }
196
197 matrix->Adopt(painter);
198 }
199 AliMUONPainterRegistry::Instance()->Register(matrix);
200 fPainterMatrices->Add(matrix);
201 }
202}
203
204//_____________________________________________________________________________
205void
206AliMUONPainterHelper::GenerateGeometry()
207{
208 /// Generate the geometry (FIXME: using transform.dat for the moment)
209 /// The geometry is not the "normal" one as we "explode" it to avoid
210 /// having overlapping detection elements as in the reality, which
211 /// would be inconvenient for a display ;-)
212
213 AliDebug(1,Form(" with explodeFactor=%e,%e",fExplodeFactor[0],fExplodeFactor[1]));
214
215 AliMUONGeometryTransformer transformer;
216 transformer.LoadGeometryData("transform.dat");
217// transformer.LoadGeometryData("geometry.root"); //FIXME: add a protection if geometry.root file does not exist
218 fExplodedGlobalTransformations = new AliMpExMap(true);
219 fRealGlobalTransformations = new AliMpExMap(true);
220 AliMpDEIterator deIt;
221 deIt.First();
222 while ( !deIt.IsDone() )
223 {
224 Int_t detElemId = deIt.CurrentDEId();
225 const AliMUONGeometryDetElement* de = transformer.GetDetElement(detElemId);
226
227 fRealGlobalTransformations->Add(detElemId,de->GetGlobalTransformation()->Clone());
228
229 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(de->GetGlobalTransformation()->Clone());
230 Double_t* translation = matrix->GetTranslation();
231
232 AliDebug(1,Form("Initial translation for DE %04d is %7.3f, %7.3f",
233 detElemId,translation[0],translation[1]));
234
235 if ( AliMpDEManager::GetStationType(detElemId) == AliMp::kStation345 )
236 {
237 translation[0] *= fExplodeFactor[0];
238 translation[1] *= fExplodeFactor[1];
239 }
240 else
241 {
242 Double_t shift = 5; // cm
243 Double_t xshift[] = { shift, -shift, -shift, shift };
244 Double_t yshift[] = { shift, shift, -shift, -shift };
245 Int_t ishift = detElemId % 100;
246
247 translation[0] += xshift[ishift];
248 translation[1] += yshift[ishift];
249 }
250 matrix->SetTranslation(translation);
251 fExplodedGlobalTransformations->Add(detElemId,matrix);
252 deIt.Next();
253 }
254}
255
256//_____________________________________________________________________________
257AliMUONPainterContour*
258AliMUONPainterHelper::GenerateManuContour(Int_t detElemId,
259 Int_t manuId,
260 AliMUONAttPainter viewType,
261 const char* contourName)
262{
263 /// Generate the contour of the list of pads
264
265 if (!fContourMaker) fContourMaker = new AliMUONPainterContourMaker(fExplodedGlobalTransformations);
266
267 AliMUONPainterContour* contour =
268 fContourMaker->GenerateManuContour(contourName,detElemId,manuId,viewType);
269
270 if (contour)
271 {
272 RegisterContour(contour);
273 }
274
275 return contour;
276}
277
278//_____________________________________________________________________________
279void
280AliMUONPainterHelper::GeneratePadStore()
281{
282 /// Generate the pad store
283
284 AliCodeTimerAuto("")
285 AliDebugClass(1,"Generating pad store");
286 fPadStore = new AliMUONPainterPadStore();
287
288 AliMpDEIterator deIt;
289
290 deIt.First();
291 while ( !deIt.IsDone() )
292 {
293 Int_t detElemId = deIt.CurrentDEId();
294 if ( AliMpDEManager::GetStationType(detElemId) != AliMp::kStationTrigger )
295 {
296 GeneratePadStore(detElemId);
297 }
298 deIt.Next();
299 }
300}
301
302//_____________________________________________________________________________
303void
304AliMUONPainterHelper::GeneratePadStore(Int_t detElemId)
305{
306 /// Generate part of the padstore for one detection element
307
308 AliMp::CathodType cathode[] = { AliMp::kCath0, AliMp::kCath1 };
309
310 for ( Int_t i = 0; i < 2; ++i )
311 {
312 const AliMpVSegmentation* seg =
313 AliMpSegmentation::Instance()->GetMpSegmentation(detElemId,cathode[i]);
314 AliMpVPadIterator* it = seg->CreateIterator();
315 it->First();
316
317 while ( !it->IsDone() )
318 {
319 AliMpPad pad = it->CurrentItem();
320
321 TVector2 localPosition(pad.Position());
322 Double_t x,y,z;
323 Local2Global(detElemId,localPosition.X(),localPosition.Y(),0,
324 x,y,z);
325 Int_t manuId = pad.GetLocation().GetFirst();
326 Int_t manuChannel = pad.GetLocation().GetSecond();
327 AliMUONVCalibParam* param = fPadStore->Get(detElemId,manuId);
328 param->SetValueAsDouble(manuChannel,0,x);
329 param->SetValueAsDouble(manuChannel,1,y);
330 param->SetValueAsDouble(manuChannel,2,pad.Dimensions().X());
331 param->SetValueAsDouble(manuChannel,3,pad.Dimensions().Y());
332 it->Next();
333 }
334 delete it;
335 }
336}
337
338//_____________________________________________________________________________
339void
340AliMUONPainterHelper::GetBoundaries(const TArrayI& pads, Double_t& xmin, Double_t& ymin,
341 Double_t& xmax, Double_t& ymax) const
342{
343 /// Get the area covered by an array of pads
344
345 return fPadStore->GetBoundaries(pads,xmin,ymin,xmax,ymax);
346}
347
348//_____________________________________________________________________________
349AliMp::CathodType
350AliMUONPainterHelper::GetCathodeType(Int_t detElemId, Int_t manuId) const
351{
352 /// Get the cathode type of a given manu
353
354 AliMp::PlaneType planeType(AliMp::kBendingPlane);
355 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) )
356 {
357 planeType = AliMp::kNonBendingPlane;
358 }
359 return AliMpDEManager::GetCathod(detElemId,planeType);
360}
361
362//_____________________________________________________________________________
363AliMUONPainterContour*
364AliMUONPainterHelper::GetLocalManuContour(Int_t detElemId, Int_t manuId) const
365{
366 /// Retrieve a manu contour (in local coordinates)
367 return fContourMaker->FindLocalManuContour(detElemId,manuId);
368}
369
370//_____________________________________________________________________________
371AliMpMotifPosition*
372AliMUONPainterHelper::GetMotifPosition(Int_t detElemId, Int_t manuId) const
373{
374 /// Get a given motif position
375 AliMp::StationType stationType = AliMpDEManager::GetStationType(detElemId);
376 if ( stationType == AliMp::kStation345 )
377 {
378 AliMp::PlaneType planeType(AliMp::kBendingPlane);
379 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) )
380 {
381 planeType = AliMp::kNonBendingPlane;
382 }
383 const AliMpSlat* slat = GetSlat(detElemId,planeType);
384 return slat->FindMotifPosition(manuId);
385 }
386 else if ( stationType != AliMp::kStationTrigger )
387 {
388 AliMp::PlaneType planeType(AliMp::kBendingPlane);
389 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) )
390 {
391 planeType = AliMp::kNonBendingPlane;
392 }
393 const AliMpSector* sector = GetSector(detElemId,planeType);
394 return sector->GetMotifMap()->FindMotifPosition(manuId);
395 }
396 AliFatalClass("Not supposed to work with trigger");
397 return 0x0;
398}
399
400
401//_____________________________________________________________________________
402AliMpPCB*
403AliMUONPainterHelper::GetPCB(Int_t detElemId, AliMp::CathodType cathodeType,
404 Int_t pcbNumber) const
405{
406 /// Get a given PCB
407 const AliMpSlat* slat = GetSlat(detElemId,cathodeType);
408 if ( slat ) return slat->GetPCB(pcbNumber);
409 return 0x0;
410}
411
412//_____________________________________________________________________________
413AliMpPCB*
414AliMUONPainterHelper::GetPCB(Int_t detElemId, AliMp::PlaneType planeType,
415 Int_t pcbNumber) const
416{
417 /// Get a given PCB
418 AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId,
419 planeType);
420 return GetPCB(detElemId,cathodeType,pcbNumber);
421}
422
423//_____________________________________________________________________________
424AliMp::PlaneType
425AliMUONPainterHelper::GetPlaneType(Int_t manuId) const
426{
427 /// Get the planeType of a given manu
428
429 if ( manuId & AliMpConstants::ManuMask(AliMp::kNonBendingPlane) )
430 {
431 return AliMp::kNonBendingPlane;
432 }
433 return AliMp::kBendingPlane;
434}
435
436//_____________________________________________________________________________
437const AliMpSlat*
438AliMUONPainterHelper::GetSlat(Int_t detElemId, AliMp::PlaneType planeType) const
439{
440 /// Get a given slat
441
442 AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId,
443 planeType);
444
445 return GetSlat(detElemId,cathodeType);
446}
447
448//_____________________________________________________________________________
449const AliMpSector*
450AliMUONPainterHelper::GetSector(Int_t detElemId, AliMp::PlaneType planeType) const
451{
452 /// Get a given sector
453 AliMp::StationType stationType = AliMpDEManager::GetStationType(detElemId);
454 if ( stationType != AliMp::kStation1 && stationType != AliMp::kStation2 ) return 0x0;
455
456 AliMp::CathodType cathodeType = AliMpDEManager::GetCathod(detElemId,planeType);
457
458 const AliMpVSegmentation* seg =
459 AliMpSegmentation::Instance()->GetMpSegmentation(detElemId,cathodeType);
460
461 const AliMpSectorSegmentation* sectorSeg = static_cast<const AliMpSectorSegmentation*>(seg);
462 return sectorSeg->GetSector();
463
464}
465
466//_____________________________________________________________________________
467const AliMpSlat*
468AliMUONPainterHelper::GetSlat(Int_t detElemId, AliMp::CathodType cathodeType) const
469{
470 /// Get a given slat
471 AliMp::StationType stationType = AliMpDEManager::GetStationType(detElemId);
472 if ( stationType != AliMp::kStation345 ) return 0x0;
473
474 const AliMpVSegmentation* seg =
475 AliMpSegmentation::Instance()->GetMpSegmentation(detElemId,cathodeType);
476
477 const AliMpSlatSegmentation* slatSeg = static_cast<const AliMpSlatSegmentation*>(seg);
478 return slatSeg->Slat();
479}
480
481//_____________________________________________________________________________
482const AliMpSlat*
483AliMUONPainterHelper::GetSlat(Int_t detElemId, Int_t manuId) const
484{
485 /// Get a given slat
486 const AliMpVSegmentation* seg =
487 AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
488
489 const AliMpSlatSegmentation* slatSeg = static_cast<const AliMpSlatSegmentation*>(seg);
490 return slatSeg->Slat();
491}
492
493//_____________________________________________________________________________
494AliMUONPainterHelper*
495AliMUONPainterHelper::Instance()
496{
497 /// Return the global and unique instance of this class
498
499 if (fgInstance) return fgInstance;
500
501 AliMUONPainterEnv env;
502
503 TString fileName(gSystem->ExpandPathName(env.String("PadStoreFileName","padstore.root")));
504
505 if ( gSystem->AccessPathName(fileName.Data(),kFileExists) ) // mind the strange return value of that method...
506 {
507 // file does NOT exist yet. Create it
508 AliDebugClass(1,"Generating instance");
509
510 Int_t ret;
511
512 new TGMsgBox(gClient->GetRoot(),gClient->GetRoot(),"",
513 Form("File %s not found.\nI will generate it, and this will take a while.\n"
514 "Click OK (and grab a cup of coffee ;-) ) to proceed,\n or Cancel to quit.",fileName.Data()),
515 kMBIconQuestion,
516 kMBOk | kMBCancel,
517 &ret);
518 if ( ret == kMBCancel ) exit(1);
519
520 fgInstance = new AliMUONPainterHelper;
521 fgInstance->GenerateGeometry();
522 fgInstance->GeneratePadStore();
523 fgInstance->GenerateDefaultMatrices();
524 fgInstance->Modified(kTRUE);
525 fgInstance->fEnv = new AliMUONPainterEnv;
526 fgInstance->Save();
527
528 }
529 else
530 {
531 AliDebugClass(1,"Reading instance");
532 TFile f(fileName.Data());
533 fgInstance = static_cast<AliMUONPainterHelper*>(f.Get("AliMUONPainterHelper"));
534
535 TIter next(fgInstance->fPainterMatrices);
536 AliMUONPainterMatrix* matrix;
537 while ( ( matrix = static_cast<AliMUONPainterMatrix*>(next()) ) )
538 {
539 AliMUONPainterRegistry::Instance()->Register(matrix);
540 }
541 fgInstance->fPainterMatrices->SetOwner(kFALSE);
542 fgInstance->fEnv = new AliMUONPainterEnv;
543 }
544 return fgInstance;
545}
546
547//_____________________________________________________________________________
548void
549AliMUONPainterHelper::Global2Local(Int_t detElemId,
550 Double_t xg, Double_t yg, Double_t zg,
551 Double_t& xl, Double_t& yl, Double_t& zl) const
552{
553 /// Local to global transformation of coordinates
554
555 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fExplodedGlobalTransformations->GetValue(detElemId));
556 Double_t pg[3] = { xg, yg, zg };
557 Double_t pl[3] = { 0., 0., 0. };
558 matrix->MasterToLocal(pg, pl);
559 xl = pl[0];
560 yl = pl[1];
561 zl = pl[2];
562}
563
564//_____________________________________________________________________________
565void
566AliMUONPainterHelper::Global2LocalReal(Int_t detElemId,
567 Double_t xg, Double_t yg, Double_t zg,
568 Double_t& xl, Double_t& yl, Double_t& zl) const
569{
570 /// Local to global transformation of coordinates
571
572 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fRealGlobalTransformations->GetValue(detElemId));
573 Double_t pg[3] = { xg, yg, zg };
574 Double_t pl[3] = { 0., 0., 0. };
575 matrix->MasterToLocal(pg, pl);
576 xl = pl[0];
577 yl = pl[1];
578 zl = pl[2];
579}
580
581//_____________________________________________________________________________
582void
583AliMUONPainterHelper::Local2Global(Int_t detElemId,
584 Double_t xl, Double_t yl, Double_t zl,
585 Double_t& xg, Double_t& yg, Double_t& zg) const
586{
587 /// Local to (exploded) global transformation of coordinates
588
589 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fExplodedGlobalTransformations->GetValue(detElemId));
590 Double_t pl[3] = { xl, yl, zl };
591 Double_t pg[3] = { 0., 0., 0. };
592 matrix->LocalToMaster(pl, pg);
593 xg = pg[0];
594 yg = pg[1];
595 zg = pg[2];
596}
597
598//_____________________________________________________________________________
599void
600AliMUONPainterHelper::Local2GlobalReal(Int_t detElemId,
601 Double_t xl, Double_t yl, Double_t zl,
602 Double_t& xg, Double_t& yg, Double_t& zg) const
603{
604 /// Local to (real) global transformation of coordinates
605
606 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fRealGlobalTransformations->GetValue(detElemId));
607 Double_t pl[3] = { xl, yl, zl };
608 Double_t pg[3] = { 0., 0., 0. };
609 matrix->LocalToMaster(pl, pg);
610 xg = pg[0];
611 yg = pg[1];
612 zg = pg[2];
613}
614
615//_____________________________________________________________________________
616Int_t
617AliMUONPainterHelper::ColorFromValue(Double_t value, Double_t min, Double_t max) const
618{
619 /// Convert a value into a color, fitting within a given range
620
621 Int_t rv;
622
623 if (value > max) rv = 1;
624 else if (value <= min) rv = 0;
625 else
626 {
627 if ( max == min ) return gStyle->GetColorPalette(1);
628 Double_t range = max - min;
629 Double_t offset = value - min;
630 rv = gStyle->GetColorPalette( 1 + int( offset*(gStyle->GetNumberOfColors()-2)/range - 0.5 ) );
631 }
632 return rv;
633}
634
635//_____________________________________________________________________________
636AliMUONPainterContour*
637AliMUONPainterHelper::MergeContours(const TObjArray& contours,
638 const char* contourName)
639{
640 /// Merge a set of contours (delegating to the contour maker)
641 if (!fContourMaker)
642 {
643 fContourMaker = new AliMUONPainterContourMaker(fExplodedGlobalTransformations);
644 }
645
646 AliMUONPainterContour* contour = fContourMaker->MergeContours(contours,
647 contourName);
648
649 if (contour)
650 {
651 RegisterContour(contour);
652 }
653 return contour;
654}
655
656
657//_____________________________________________________________________________
658void
659AliMUONPainterHelper::Print(Option_t* opt) const
660{
661 /// Printout
662 TString sopt(opt);
663 sopt.ToUpper();
664
665 if ( sopt.Length() == 0 )
666 {
667 cout << Form("ExplodeFactor=%e,%e",fExplodeFactor[0],fExplodeFactor[1]) << endl;
668 cout << Form("PadStore=%x",fPadStore);
669 if ( fPadStore ) cout << Form(" with %d pads",fPadStore->GetSize());
670 cout << endl;
671 cout << Form("GlobalTransformations=%x",fExplodedGlobalTransformations);
672 if ( fExplodedGlobalTransformations ) cout << Form(" with %d transformations",fExplodedGlobalTransformations->GetSize());
673 cout << endl;
674 if ( fContourMaker )
675 {
676 cout << Form(" with %d contours",fContourMaker->Size());
677 }
678 else
679 {
680 cout << "No contour";
681 }
682 cout << endl;
683 cout << "Modified=";
684 if ( IsModified() )
685 {
686 cout << "YES";
687 }
688 else
689 {
690 cout << "NO";
691 }
692 cout << endl;
693 }
694
695 if ( sopt.Contains("CONTOUR") || sopt.Contains("FULL") )
696 {
697 fContourMaker->Print(opt);
698 }
699
700 if ( sopt.Contains("MATRI") || sopt.Contains("FULL") )
701 {
702 fPainterMatrices->Print(opt);
703 }
704}
705
706//_____________________________________________________________________________
707void
708AliMUONPainterHelper::RegisterContour(AliMUONPainterContour* contour)
709{
710 /// contour is adopted by contourMaker
711 AliCodeTimerAuto("")
712 AliDebug(1,contour->GetName());
713 if ( fContourMaker->HasContour(contour->GetName()) )
714 {
715 AliError(Form("Contour with name %s is already there",contour->GetName()));
716// Print("CONTOUR");
717 return;
718 }
719 fContourMaker->Add(contour);
720 Modified(kTRUE);
721}
722
723//_____________________________________________________________________________
724void
725AliMUONPainterHelper::Save()
726{
727 /// Save to disk
728
729 if (!IsModified()) return;
730
731 Modified(kFALSE);
732
733 AliInfo("");
734
735 fgInstance->Print();
736
737 TString fileName(gSystem->ExpandPathName(fgInstance->Env()->String("PadStoreFileName","padstore.root")));
738
739 AliInfo(Form("Saving to %s",fileName.Data()));
740
741 TFile f(fileName,"RECREATE");
742
743 fgInstance->Write("");
744
745 f.Close();
746}
747
748//_____________________________________________________________________________
749AliMpPad
750AliMUONPainterHelper::PadByExplodedPosition(Int_t detElemId, Int_t manuId,
751 Double_t x, Double_t y) const
752{
753 /// Find a pad by exploded position. FIXME: not really used nor tested !
754
755 Double_t xr, yr, zr;
756
757// Local2Global(detElemId,0.0,0.0,0.0,dummy,dummy,z); // to find z
758
759 AliDebug(1,Form("DE %04d ManuID %04d x %7.3f y %7.3f",detElemId,manuId,x,y));
760
761 Exploded2Real(detElemId,x,y,0,xr,yr,zr);
762
763 AliDebug(1,Form("xr %7.3f yr %7.3f zr %7.3f",xr,yr,zr));
764
765 Double_t xl,yl,zl;
766
767 Global2LocalReal(detElemId,xr,yr,zr,xl,yl,zl);
768
769 AliDebug(1,Form("xl %7.3f yl %7.3f zl %7.3f",xl,yl,zl));
770
771 const AliMpVSegmentation* seg = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
772
773 AliDebug(1,Form("dx,dy=%7.3f,%7.3f",seg->Dimensions().X(),seg->Dimensions().Y()));
774
775 return seg->PadByPosition(TVector2(xl,yl));
776}
777
778//_____________________________________________________________________________
779void
780AliMUONPainterHelper::Exploded2Real(Int_t detElemId,
781 Double_t xe, Double_t ye, Double_t ze,
782 Double_t& xr, Double_t& yr, Double_t& zr) const
783{
784 /// Convert exploded coordinates into real ones. FIXME: not really used nor tested !
785
786 // first go back to local
787
788 Double_t xl,yl,zl;
789
790 Global2Local(detElemId,xe,ye,ze,xl,yl,zl);
791
792 // and then back to global but not exploded
793
794 Local2GlobalReal(detElemId,xl,yl,zl,xr,yr,zr);
795}
796
797//_____________________________________________________________________________
798TString
799AliMUONPainterHelper::ChamberName(Int_t chamberId) const
800{
801 /// Build a name for one chamber
802 return Form("Chamber%1d",chamberId);
803}
804
805//_____________________________________________________________________________
806TString
807AliMUONPainterHelper::StationName(Int_t stationId) const
808{
809 /// Build a name for one station
810 return Form("Station%1d",stationId+1);
811}
812
813//_____________________________________________________________________________
814TString
815AliMUONPainterHelper::DEName(Int_t detElemId) const
816{
817 /// Build a name for one detection element
818 return Form("DE%04d",detElemId);
819}
820
821//_____________________________________________________________________________
822TString
823AliMUONPainterHelper::ManuName(Int_t manuId) const
824{
825 /// Build a name for one manu
826 return Form("MANU%04d",manuId);
827}
828
829//_____________________________________________________________________________
830TString
831AliMUONPainterHelper::BusPatchName(Int_t busPatchId) const
832{
833 /// Build a name for one buspatch
834 return Form("BUSPATCH%04d",busPatchId);
835}
836
837//_____________________________________________________________________________
838TString
839AliMUONPainterHelper::PCBName(Int_t pcbNumber) const
840{
841 /// Build a name for one pcb
842 return Form("PCB%1d",pcbNumber);
843}
844
845//_____________________________________________________________________________
846TString
847AliMUONPainterHelper::ChamberPathName(Int_t chamberId) const
848{
849 /// Build a name for one chamber
850 return Form("%s/%s",StationName(chamberId/2).Data(),ChamberName(chamberId).Data());
851}
852
853//_____________________________________________________________________________
854TString
855AliMUONPainterHelper::StationPathName(Int_t stationId) const
856{
857 /// Build a name for one station
858 return StationName(stationId);
859}
860
861//_____________________________________________________________________________
862TString
863AliMUONPainterHelper::DEPathName(Int_t detElemId) const
864{
865 /// Build a name for one detection element
866
867 Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
868
869 return Form("%s/%s/%s",
870 StationName(chamberId/2).Data(),
871 ChamberName(chamberId).Data(),
872 DEName(detElemId).Data());
873}
874
875//_____________________________________________________________________________
876TString
877AliMUONPainterHelper::ManuPathName(Int_t detElemId, Int_t manuId) const
878{
879 /// Build a name for one manu
880 return Form("%s/%s",DEPathName(detElemId).Data(),ManuName(manuId).Data());
881}
882
883//_____________________________________________________________________________
884TString
885AliMUONPainterHelper::BusPatchPathName(Int_t busPatchId) const
886{
887 /// Build a name for one buspatch
888 Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(busPatchId);
889
890 return Form("%s/%s",DEPathName(detElemId).Data(),BusPatchName(busPatchId).Data());
891}
892
893//_____________________________________________________________________________
894TString
895AliMUONPainterHelper::PCBPathName(Int_t detElemId, Int_t pcbNumber) const
896{
897 /// Build a name for one pcb
898 return Form("%s/%s",DEPathName(detElemId).Data(),PCBName(pcbNumber).Data());
899}
900
901//_____________________________________________________________________________
902TString
903AliMUONPainterHelper::FormatValue(const char* name, Double_t value) const
904{
905 /// Format a double value to be displayed
906 /// FIXME: should insure we have the right number of significant digits here...
907
908 return Form("%s = %e",name,value);
909}