remoe duplicate QA initialisation and do ESD QA for same detectors as RecPoint QA
[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
49419555 503 TString fileName(gSystem->ExpandPathName(env.String("PadStoreFileName","$HOME/padstore.root")));
0145e89a 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;
49419555 526 fgInstance->fEnv->Set("PadStoreFileName",fileName.Data());
0145e89a 527 fgInstance->Save();
528
529 }
530 else
531 {
532 AliDebugClass(1,"Reading instance");
533 TFile f(fileName.Data());
534 fgInstance = static_cast<AliMUONPainterHelper*>(f.Get("AliMUONPainterHelper"));
535
536 TIter next(fgInstance->fPainterMatrices);
537 AliMUONPainterMatrix* matrix;
538 while ( ( matrix = static_cast<AliMUONPainterMatrix*>(next()) ) )
539 {
540 AliMUONPainterRegistry::Instance()->Register(matrix);
541 }
542 fgInstance->fPainterMatrices->SetOwner(kFALSE);
543 fgInstance->fEnv = new AliMUONPainterEnv;
544 }
545 return fgInstance;
546}
547
548//_____________________________________________________________________________
549void
550AliMUONPainterHelper::Global2Local(Int_t detElemId,
551 Double_t xg, Double_t yg, Double_t zg,
552 Double_t& xl, Double_t& yl, Double_t& zl) const
553{
554 /// Local to global transformation of coordinates
555
556 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fExplodedGlobalTransformations->GetValue(detElemId));
557 Double_t pg[3] = { xg, yg, zg };
558 Double_t pl[3] = { 0., 0., 0. };
559 matrix->MasterToLocal(pg, pl);
560 xl = pl[0];
561 yl = pl[1];
562 zl = pl[2];
563}
564
565//_____________________________________________________________________________
566void
567AliMUONPainterHelper::Global2LocalReal(Int_t detElemId,
568 Double_t xg, Double_t yg, Double_t zg,
569 Double_t& xl, Double_t& yl, Double_t& zl) const
570{
571 /// Local to global transformation of coordinates
572
573 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fRealGlobalTransformations->GetValue(detElemId));
574 Double_t pg[3] = { xg, yg, zg };
575 Double_t pl[3] = { 0., 0., 0. };
576 matrix->MasterToLocal(pg, pl);
577 xl = pl[0];
578 yl = pl[1];
579 zl = pl[2];
580}
581
582//_____________________________________________________________________________
583void
584AliMUONPainterHelper::Local2Global(Int_t detElemId,
585 Double_t xl, Double_t yl, Double_t zl,
586 Double_t& xg, Double_t& yg, Double_t& zg) const
587{
588 /// Local to (exploded) global transformation of coordinates
589
590 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fExplodedGlobalTransformations->GetValue(detElemId));
591 Double_t pl[3] = { xl, yl, zl };
592 Double_t pg[3] = { 0., 0., 0. };
593 matrix->LocalToMaster(pl, pg);
594 xg = pg[0];
595 yg = pg[1];
596 zg = pg[2];
597}
598
599//_____________________________________________________________________________
600void
601AliMUONPainterHelper::Local2GlobalReal(Int_t detElemId,
602 Double_t xl, Double_t yl, Double_t zl,
603 Double_t& xg, Double_t& yg, Double_t& zg) const
604{
605 /// Local to (real) global transformation of coordinates
606
607 TGeoHMatrix* matrix = static_cast<TGeoHMatrix*>(fRealGlobalTransformations->GetValue(detElemId));
608 Double_t pl[3] = { xl, yl, zl };
609 Double_t pg[3] = { 0., 0., 0. };
610 matrix->LocalToMaster(pl, pg);
611 xg = pg[0];
612 yg = pg[1];
613 zg = pg[2];
614}
615
616//_____________________________________________________________________________
617Int_t
618AliMUONPainterHelper::ColorFromValue(Double_t value, Double_t min, Double_t max) const
619{
620 /// Convert a value into a color, fitting within a given range
621
622 Int_t rv;
623
624 if (value > max) rv = 1;
625 else if (value <= min) rv = 0;
626 else
627 {
628 if ( max == min ) return gStyle->GetColorPalette(1);
629 Double_t range = max - min;
630 Double_t offset = value - min;
631 rv = gStyle->GetColorPalette( 1 + int( offset*(gStyle->GetNumberOfColors()-2)/range - 0.5 ) );
632 }
633 return rv;
634}
635
636//_____________________________________________________________________________
637AliMUONPainterContour*
638AliMUONPainterHelper::MergeContours(const TObjArray& contours,
639 const char* contourName)
640{
641 /// Merge a set of contours (delegating to the contour maker)
642 if (!fContourMaker)
643 {
644 fContourMaker = new AliMUONPainterContourMaker(fExplodedGlobalTransformations);
645 }
646
647 AliMUONPainterContour* contour = fContourMaker->MergeContours(contours,
648 contourName);
649
650 if (contour)
651 {
652 RegisterContour(contour);
653 }
654 return contour;
655}
656
657
658//_____________________________________________________________________________
659void
660AliMUONPainterHelper::Print(Option_t* opt) const
661{
662 /// Printout
663 TString sopt(opt);
664 sopt.ToUpper();
665
666 if ( sopt.Length() == 0 )
667 {
668 cout << Form("ExplodeFactor=%e,%e",fExplodeFactor[0],fExplodeFactor[1]) << endl;
669 cout << Form("PadStore=%x",fPadStore);
670 if ( fPadStore ) cout << Form(" with %d pads",fPadStore->GetSize());
671 cout << endl;
672 cout << Form("GlobalTransformations=%x",fExplodedGlobalTransformations);
673 if ( fExplodedGlobalTransformations ) cout << Form(" with %d transformations",fExplodedGlobalTransformations->GetSize());
674 cout << endl;
675 if ( fContourMaker )
676 {
677 cout << Form(" with %d contours",fContourMaker->Size());
678 }
679 else
680 {
681 cout << "No contour";
682 }
683 cout << endl;
684 cout << "Modified=";
685 if ( IsModified() )
686 {
687 cout << "YES";
688 }
689 else
690 {
691 cout << "NO";
692 }
693 cout << endl;
694 }
695
696 if ( sopt.Contains("CONTOUR") || sopt.Contains("FULL") )
697 {
698 fContourMaker->Print(opt);
699 }
700
701 if ( sopt.Contains("MATRI") || sopt.Contains("FULL") )
702 {
703 fPainterMatrices->Print(opt);
704 }
705}
706
707//_____________________________________________________________________________
708void
709AliMUONPainterHelper::RegisterContour(AliMUONPainterContour* contour)
710{
711 /// contour is adopted by contourMaker
712 AliCodeTimerAuto("")
713 AliDebug(1,contour->GetName());
714 if ( fContourMaker->HasContour(contour->GetName()) )
715 {
716 AliError(Form("Contour with name %s is already there",contour->GetName()));
717// Print("CONTOUR");
718 return;
719 }
720 fContourMaker->Add(contour);
721 Modified(kTRUE);
722}
723
724//_____________________________________________________________________________
725void
726AliMUONPainterHelper::Save()
727{
728 /// Save to disk
729
730 if (!IsModified()) return;
731
732 Modified(kFALSE);
733
734 AliInfo("");
735
736 fgInstance->Print();
737
49419555 738 fgInstance->Env()->Save();
739
740 TString fileName(gSystem->ExpandPathName(fgInstance->Env()->String("PadStoreFileName")));
0145e89a 741
742 AliInfo(Form("Saving to %s",fileName.Data()));
743
744 TFile f(fileName,"RECREATE");
745
746 fgInstance->Write("");
747
748 f.Close();
749}
750
751//_____________________________________________________________________________
752AliMpPad
753AliMUONPainterHelper::PadByExplodedPosition(Int_t detElemId, Int_t manuId,
754 Double_t x, Double_t y) const
755{
756 /// Find a pad by exploded position. FIXME: not really used nor tested !
757
758 Double_t xr, yr, zr;
759
760// Local2Global(detElemId,0.0,0.0,0.0,dummy,dummy,z); // to find z
761
762 AliDebug(1,Form("DE %04d ManuID %04d x %7.3f y %7.3f",detElemId,manuId,x,y));
763
764 Exploded2Real(detElemId,x,y,0,xr,yr,zr);
765
766 AliDebug(1,Form("xr %7.3f yr %7.3f zr %7.3f",xr,yr,zr));
767
768 Double_t xl,yl,zl;
769
770 Global2LocalReal(detElemId,xr,yr,zr,xl,yl,zl);
771
772 AliDebug(1,Form("xl %7.3f yl %7.3f zl %7.3f",xl,yl,zl));
773
774 const AliMpVSegmentation* seg = AliMpSegmentation::Instance()->GetMpSegmentationByElectronics(detElemId,manuId);
775
776 AliDebug(1,Form("dx,dy=%7.3f,%7.3f",seg->Dimensions().X(),seg->Dimensions().Y()));
777
778 return seg->PadByPosition(TVector2(xl,yl));
779}
780
781//_____________________________________________________________________________
782void
783AliMUONPainterHelper::Exploded2Real(Int_t detElemId,
784 Double_t xe, Double_t ye, Double_t ze,
785 Double_t& xr, Double_t& yr, Double_t& zr) const
786{
787 /// Convert exploded coordinates into real ones. FIXME: not really used nor tested !
788
789 // first go back to local
790
791 Double_t xl,yl,zl;
792
793 Global2Local(detElemId,xe,ye,ze,xl,yl,zl);
794
795 // and then back to global but not exploded
796
797 Local2GlobalReal(detElemId,xl,yl,zl,xr,yr,zr);
798}
799
800//_____________________________________________________________________________
801TString
802AliMUONPainterHelper::ChamberName(Int_t chamberId) const
803{
804 /// Build a name for one chamber
805 return Form("Chamber%1d",chamberId);
806}
807
808//_____________________________________________________________________________
809TString
810AliMUONPainterHelper::StationName(Int_t stationId) const
811{
812 /// Build a name for one station
813 return Form("Station%1d",stationId+1);
814}
815
816//_____________________________________________________________________________
817TString
818AliMUONPainterHelper::DEName(Int_t detElemId) const
819{
820 /// Build a name for one detection element
821 return Form("DE%04d",detElemId);
822}
823
824//_____________________________________________________________________________
825TString
826AliMUONPainterHelper::ManuName(Int_t manuId) const
827{
828 /// Build a name for one manu
829 return Form("MANU%04d",manuId);
830}
831
832//_____________________________________________________________________________
833TString
834AliMUONPainterHelper::BusPatchName(Int_t busPatchId) const
835{
836 /// Build a name for one buspatch
837 return Form("BUSPATCH%04d",busPatchId);
838}
839
840//_____________________________________________________________________________
841TString
842AliMUONPainterHelper::PCBName(Int_t pcbNumber) const
843{
844 /// Build a name for one pcb
845 return Form("PCB%1d",pcbNumber);
846}
847
848//_____________________________________________________________________________
849TString
850AliMUONPainterHelper::ChamberPathName(Int_t chamberId) const
851{
852 /// Build a name for one chamber
853 return Form("%s/%s",StationName(chamberId/2).Data(),ChamberName(chamberId).Data());
854}
855
856//_____________________________________________________________________________
857TString
858AliMUONPainterHelper::StationPathName(Int_t stationId) const
859{
860 /// Build a name for one station
861 return StationName(stationId);
862}
863
864//_____________________________________________________________________________
865TString
866AliMUONPainterHelper::DEPathName(Int_t detElemId) const
867{
868 /// Build a name for one detection element
869
870 Int_t chamberId = AliMpDEManager::GetChamberId(detElemId);
871
872 return Form("%s/%s/%s",
873 StationName(chamberId/2).Data(),
874 ChamberName(chamberId).Data(),
875 DEName(detElemId).Data());
876}
877
878//_____________________________________________________________________________
879TString
880AliMUONPainterHelper::ManuPathName(Int_t detElemId, Int_t manuId) const
881{
882 /// Build a name for one manu
883 return Form("%s/%s",DEPathName(detElemId).Data(),ManuName(manuId).Data());
884}
885
886//_____________________________________________________________________________
887TString
888AliMUONPainterHelper::BusPatchPathName(Int_t busPatchId) const
889{
890 /// Build a name for one buspatch
891 Int_t detElemId = AliMpDDLStore::Instance()->GetDEfromBus(busPatchId);
892
893 return Form("%s/%s",DEPathName(detElemId).Data(),BusPatchName(busPatchId).Data());
894}
895
896//_____________________________________________________________________________
897TString
898AliMUONPainterHelper::PCBPathName(Int_t detElemId, Int_t pcbNumber) const
899{
900 /// Build a name for one pcb
901 return Form("%s/%s",DEPathName(detElemId).Data(),PCBName(pcbNumber).Data());
902}
903
904//_____________________________________________________________________________
905TString
906AliMUONPainterHelper::FormatValue(const char* name, Double_t value) const
907{
908 /// Format a double value to be displayed
909 /// FIXME: should insure we have the right number of significant digits here...
910
911 return Form("%s = %e",name,value);
912}