2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
4 /**************************************************************************
5 * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
6 * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for *
7 * full copyright notice. *
8 **************************************************************************/
13 #include "TLinearFitter.h"
15 #include "TGeoMatrix.h"
16 #include "TGeoManager.h"
18 #include "TEveTrans.h"
19 #include "TEveManager.h"
21 #include "EveBase/AliEveEventManager.h"
23 #include "AliEveTRDData.h"
24 #include "AliEveTRDModuleImp.h"
25 #include "AliEveTRDLoader.h"
26 #include "AliEveTRDLoaderImp.h"
28 #include "AliGeomManager.h"
29 #include "AliESDtrack.h"
32 #include "AliTrackPointArray.h"
33 #include "AliRieman.h"
35 #include "AliTRDhit.h"
36 #include "AliTRDcluster.h"
37 #include "AliTRDseedV1.h"
38 #include "AliTRDtrackletMCM.h"
39 #include "AliTRDtrackletWord.h"
40 #include "AliTRDmcmSim.h"
41 #include "AliTRDtrackV1.h"
42 #include "AliTRDtrackerV1.h"
43 #include "AliTRDpadPlane.h"
44 #include "AliTRDdigitsManager.h"
45 #include "AliTRDmcmSim.h"
46 #include "AliTRDarrayADC.h"
47 #include "AliTRDSignalIndex.h"
48 #include "AliTRDgeometry.h"
49 #include "AliTRDtransform.h"
50 #include "AliTRDReconstructor.h"
51 #include "AliTRDrecoParam.h"
53 ClassImp(AliEveTRDHits)
54 ClassImp(AliEveTRDDigits)
55 ClassImp(AliEveTRDClusters)
56 ClassImp(AliEveTRDTracklet)
57 ClassImp(AliEveTRDTrack)
58 ClassImp(AliEveTRDTrackletOnline)
59 ClassImp(AliEveTRDmcm)
61 ///////////////////////////////////////////////////////////
62 ///////////// AliEveTRDDigits /////////////////////
63 ///////////////////////////////////////////////////////////
65 //______________________________________________________________________________
66 AliEveTRDDigits::AliEveTRDDigits(AliEveTRDChamber *p)
67 :TEveQuadSet("digits", "")
72 gStyle->SetPalette(1, 0);
73 SetPalette(new TEveRGBAPalette(0, 512));
74 Reset(TEveQuadSet::kQT_RectangleYZ, kFALSE, 32);
77 //______________________________________________________________________________
78 AliEveTRDDigits::~AliEveTRDDigits()
80 // AliInfo(GetTitle());
83 //______________________________________________________________________________
84 void AliEveTRDDigits::SetData(AliTRDdigitsManager *digits)
88 Int_t det(fParent->GetID());
89 AliTRDarrayADC *data = digits->GetDigits(det);
90 if(!data->GetDim()) return;
93 AliTRDSignalIndex *indexes = digits->GetIndexes(det);
94 if(!indexes->IsAllocated()) digits->BuildIndexes(det);
96 Double_t scale, dy, dz;
97 Int_t ly = AliTRDgeometry::GetLayer(det),
98 stk = AliTRDgeometry::GetStack(det),
99 sec = AliTRDgeometry::GetSector(det),
100 vid = AliGeomManager::LayerToVolUID(AliGeomManager::kTRD1 + ly, stk + AliTRDgeometry::Nstack() * sec);
101 SetNameTitle(Form("digits%03d", det), Form("D-%03d [%02d_%d_%d]", det, sec, stk, ly));
102 Short_t sig[7]={0,0,0,10,0,0,0};
104 AliTRDtransform transform(det);
105 AliTRDpadPlane *pp(fParent->fGeo->GetPadPlane(ly, stk));
109 indexes->ResetCounters();
110 while (indexes->NextRCIndex(row, col)){
111 dz = pp->GetRowSize(row);
112 dy = pp->GetColSize(col);
113 Short_t *const adc = data->GetDataAddress(row,col);
114 for (Int_t time(0); time<30; time++){
115 if(data->IsPadCorrupted(row, col, time)) break;
116 if(adc[time] <= 1) continue;
117 new (&c) AliTRDcluster(det, col, row, time, sig, vid);
118 transform.Transform(&c);
120 scale = adc[time] < 512 ? adc[time]/512. : 1.;
121 AddQuad(c.GetY()-0.5*dy, c.GetZ()-0.5*dz*scale, c.GetX(), dy*0.95, dz*scale);
122 QuadValue(Int_t(adc[time]));
123 QuadId(new TNamed(Form("ADC %d", adc[time]), Form("det[%3d(%02d_%d_%d)] col[%3d] row[%2d] tb[%2d]", det, sec, stk, ly, col, row, time)));
127 // rotate to global coordinates
129 TEveTrans& t = RefMainTrans();
130 t.SetRotByAngles((sec+.5)*AliTRDgeometry::GetAlpha(), 0.,0.);
134 // //______________________________________________________________________________
135 // void AliEveTRDDigits::Paint(Option_t *option)
137 // // Paint the object.
139 // if(fParent->GetDigitsBox()) fBoxes.Paint(option);
140 // else TEveQuadSet::Paint(option);
143 // //______________________________________________________________________________
144 // void AliEveTRDDigits::Reset()
146 // // Reset raw and visual data.
148 // TEveQuadSet::Reset(TEveQuadSet::kQT_RectangleYZ, kTRUE, 64);
149 // // MT fBoxes.fBoxes.clear();
153 ///////////////////////////////////////////////////////////
154 ///////////// AliEveTRDHits /////////////////////
155 ///////////////////////////////////////////////////////////
157 //______________________________________________________________________________
158 AliEveTRDHits::AliEveTRDHits() : TEvePointSet("hits", 20)
162 SetMarkerColor(kGreen);
166 //______________________________________________________________________________
167 AliEveTRDHits::~AliEveTRDHits()
169 //AliInfo(GetTitle());
172 //______________________________________________________________________________
173 void AliEveTRDHits::PointSelected(Int_t n)
175 // Handle an individual point selection from GL.
178 if(!(h = dynamic_cast<AliTRDhit*>(GetPointId(n)))) return;
179 printf("Id[%3d] Det[%3d] Reg[%c] TR[%c] Q[%3d] MC[%d] t[%f]\n",
181 h->FromAmplification() ? 'A' : 'D',
182 h->FromTRphoton() ? 'y' : 'n',
183 h->GetCharge(), h->GetTrack(), h->GetTime());
187 ///////////////////////////////////////////////////////////
188 ///////////// AliEveTRDClusters /////////////////////
189 ///////////////////////////////////////////////////////////
191 //______________________________________________________________________________
192 AliEveTRDClusters::AliEveTRDClusters():AliEveTRDHits()
199 SetMarkerColor(kGray);
203 //______________________________________________________________________________
204 void AliEveTRDClusters::PointSelected(Int_t n)
206 // Handle an individual point selection from GL.
208 AliTRDcluster *c = dynamic_cast<AliTRDcluster*>(GetPointId(n));
211 Emit("PointSelected(Int_t)", n);
212 // Bool_t AliCluster::GetGlobalCov(Float_t* cov) const
213 // Bool_t AliCluster::GetGlobalXYZ(Float_t* xyz) const
214 // Float_t AliCluster::GetSigmaY2() const
215 // Float_t AliCluster::GetSigmaYZ() const
216 // Float_t AliCluster::GetSigmaZ2() const
219 //______________________________________________________________________________
220 void AliEveTRDClusters::Print(Option_t *o) const
222 AliTRDcluster *c = NULL;
224 for(Int_t n = GetN(); n--;){
225 if(!(c = dynamic_cast<AliTRDcluster*>(GetPointId(n)))) continue;
230 //______________________________________________________________________________
231 void AliEveTRDClusters::Load(const Char_t *w) const
234 if(strcmp(w, "hit")==0) typ = 0;
235 else if(strcmp(w, "dig")==0) typ = 1;
236 else if(strcmp(w, "cls")==0) typ = 2;
237 else if(strcmp(w, "tlt")==0) typ = 3;
238 else if(strcmp(w, "all")==0) typ = 4;
240 AliInfo("The following arguments are accepted:");
241 AliInfo(" \"hit\" : loading of MC hits");
242 AliInfo(" \"dig\" : loading of digits");
243 AliInfo(" \"cls\" : loading of reconstructed clusters");
244 AliInfo(" \"tlt\" : loading of online tracklets");
245 AliInfo(" \"all\" : loading everything available");
249 AliTRDcluster *c = NULL;
251 while((n = GetN() && !(c = dynamic_cast<AliTRDcluster*>(GetPointId(n))))) n++;
254 Int_t det = c->GetDetector();
255 AliEveTRDLoader *loader(NULL);
258 loader = new AliEveTRDLoader("Hits");
259 if(!loader->Open("TRD.Hits.root")){
265 loader = new AliEveTRDLoader("Digits");
266 if(!loader->Open("TRD.Digits.root")){
272 loader = new AliEveTRDLoader("Clusters");
273 if(!loader->Open("TRD.RecPoints.root")){
279 loader = new AliEveTRDLoader("Tracklets");
280 if(!loader->Open("TRD.Tracklets.root")){
286 loader = new AliEveTRDLoaderSim("MC"); // Assume MC data. Try to load via run loader
287 if(!loader->Open("galice.root")){
288 delete loader; loader=NULL;
289 // Assume measured data. Try to load data types individualy
291 loader = new AliEveTRDLoader("Data");
292 if(!loader->Open("TRD.Digits.root")) failed++;
293 if(!loader->Open("TRD.RecPoints.root")) failed++;
294 if(!loader->Open("TRD.Tracklets.root")) failed++;
295 if(failed==3){delete loader;loader=NULL;}
296 }else loader->SetDataType(AliEveTRDLoader::kTRDHits | AliEveTRDLoader::kTRDDigits | AliEveTRDLoader::kTRDClusters);
302 loader->AddChambers(AliTRDgeometry::GetSector(det),AliTRDgeometry::GetStack(det), AliTRDgeometry::GetLayer(det));
304 if(!loader->GoToEvent(AliEveEventManager::GetCurrent()->GetEventId())) AliWarning(Form("No data loaded for event %d", AliEveEventManager::GetCurrent()->GetEventId()));
305 // register loader with alieve
306 gEve->AddElement(loader->GetChamber(det), *(BeginParents()));
307 //loader->SpawnEditor();
311 ///////////////////////////////////////////////////////////
312 ///////////// AliEveTRDTracklet /////////////////////
313 ///////////////////////////////////////////////////////////
315 //______________________________________________________________________________
316 AliEveTRDTracklet::AliEveTRDTracklet(AliTRDseedV1 *trklt):TEveLine()
322 if(!gGeoManager) AliEveEventManager::AssertGeometry();
324 // init tracklet line
325 Int_t sec = AliTRDgeometry::GetSector(trklt->GetDetector());
326 Float_t alpha((0.5+sec)*AliTRDgeometry::GetAlpha()),
327 cphi(TMath::Cos(alpha)),
328 sphi(TMath::Sin(alpha));
330 Float_t x0 = trklt->GetX0();
331 Float_t y0 = trklt->GetYref(0);
332 Float_t z0 = trklt->GetZref(0);
333 Float_t dydx = trklt->GetYref(1);
334 Float_t dzdx = trklt->GetZref(1);
335 Float_t tilt = trklt->GetTilt();
337 AliTRDcluster *c = NULL;
338 for(Int_t ic=0; ic<AliTRDseedV1::kNclusters; ic++){
339 if(!(c = trklt->GetClusters(ic))) continue;
340 if(!fClusters) AddElement(fClusters = new AliEveTRDClusters());
342 //Float_t yt = y0 - dx*dydx;
343 Float_t zt = z0 - dx*dzdx;
344 // backup yc - for testing purposes
345 Float_t yc = c->GetY();
346 c->SetY(yc-tilt*(c->GetZ()-zt));
348 // Ben: Temporary(?) fix -> Issue with wrongly set flag "fIsMisaligned" of AliCluster
349 // Changes can be undone ofter the issue with the flag has been solved
350 //c->GetGlobalXYZ(g);
352 // The following code for the trasformation local-to-global is just adapted from
353 // AliTRDgeometry::RotateBack(...)
355 // tracking to global coordinates transformation
356 Float_t loc[3] = { c->GetX(), c->GetY(), c->GetZ() };
357 g[0] = loc[0] * cphi - loc[1] * sphi;
358 g[1] = loc[0] * sphi + loc[1] * cphi;
361 Int_t id = fClusters->SetNextPoint(g[0], g[1], g[2]);
363 fClusters->SetPointId(id, new AliTRDcluster(*c));
366 fClusters->SetName("TC clusters");
367 fClusters->SetTitle(Form("N[%d]", trklt->GetN()));
368 fClusters->SetMarkerColor(kMagenta);
371 SetTitle(Form("Det[%d] RC[%c] Layer[%d] P[%7.3f]", trklt->GetDetector(), trklt->IsRowCross()?'y':'n', trklt->GetPlane(), trklt->GetMomentum()));
376 y0 = trklt->GetYfit(0);
377 dydx = trklt->GetYfit(1);
378 Double_t xg(x0 * cphi - y0 * sphi),
379 yg(x0 * sphi + y0 * cphi);
380 SetPoint(0, xg, yg, z0);
381 //SetPoint(0, x0, y0, z0);
384 dx = .5*AliTRDgeometry::CamHght()+AliTRDgeometry::CdrHght();
388 xg = x0 * cphi - y0 * sphi;
389 yg = x0 * sphi + y0 * cphi;
390 SetPoint(1, xg, yg, z0);
391 //SetPoint(1, x0, y0, z0);
395 AliEveTRDTracklet::~AliEveTRDTracklet()
400 //______________________________________________________________________________
401 void AliEveTRDTracklet::Print(Option_t *o) const
403 AliTRDseedV1 *tracklet = (AliTRDseedV1*)GetUserData();
404 if(!tracklet) return;
408 ///////////////////////////////////////////////////////////
409 ///////////// AliEveTRDTrack /////////////////////
410 ///////////////////////////////////////////////////////////
412 //______________________________________________________________________________
413 AliEveTRDTrack::AliEveTRDTrack(AliTRDtrackV1 *trk)
425 // AliTRDtrackerV1::SetNTimeBins(24);
427 fRim = new AliRieman(trk->GetNumberOfClusters());
428 AliTRDseedV1 *tracklet = NULL;
429 for(Int_t il=0; il<AliTRDgeometry::kNlayer; il++){
430 if(!(tracklet = trk->GetTracklet(il))) continue;
431 if(!tracklet->IsOK()) continue;
432 AddElement(new AliEveTRDTracklet(tracklet));
434 // tracklet->ResetClusterIter(kFALSE);
435 // while((c = tracklet->PrevCluster())){
436 // AliTRDcluster *c(NULL);
437 /* for(Int_t ic=AliTRDseedV1::kNtb; ic--;){
438 if(!(c=tracklet->GetClusters(ic))) continue;
439 Float_t xc = c->GetX();
440 Float_t yc = c->GetY();
441 Float_t zc = c->GetZ();
442 Float_t zt = tracklet->GetZref(0) - (tracklet->GetX0()-xc)*tracklet->GetZref(1);
443 yc -= tracklet->GetTilt()*(zc-zt);
444 fRim->AddPoint(xc, yc, zc, .05, 2.3);
447 if(trk->GetNumberOfTracklets()>1) fRim->Update();
448 SetStatus(fTrackState);
451 //______________________________________________________________________________
452 AliEveTRDTrack::~AliEveTRDTrack()
454 if(fPoints) delete [] fPoints; fPoints = NULL;
455 //delete dynamic_cast<AliTRDtrackV1*>(GetUserData());
458 //______________________________________________________________________________
459 void AliEveTRDTrack::Print(Option_t *o) const
461 AliTRDtrackV1 *track = (AliTRDtrackV1*)GetUserData();
466 //______________________________________________________________________________
467 void AliEveTRDTrack::SetStatus(UChar_t s)
469 // nothing to be done
470 if(fPoints && fTrackState == s) return;
472 const Int_t nc = AliTRDtrackV1::kMAXCLUSTERSPERTRACK;
473 AliTRDtrackV1 *trk(NULL);
474 if(!(trk=static_cast<AliTRDtrackV1*>(GetUserData()))) {
475 AliError("Failed casting data to TRD track.");
479 Bool_t BUILD = kFALSE;
481 fPoints = new AliTrackPoint[nc];
483 // define the radial span of the track in the TRD
484 Double_t xmin = -1., xmax = -1.;
486 AliTRDseedV1 *trklt = NULL;
487 for(Int_t ily=0; ily<AliTRDgeometry::kNlayer; ily++){
488 if(!(trklt = trk->GetTracklet(ily))) continue;
489 if(xmin<0.) xmin = trklt->GetX0() - AliTRDgeometry::CamHght() - AliTRDgeometry::CdrHght();
490 if(trklt->GetX0()>xmax) xmax = trklt->GetX0();
491 det = trklt->GetDetector();
493 Int_t sec = det/AliTRDgeometry::kNdets;
494 fAlpha = AliTRDgeometry::GetAlpha() * (sec<9 ? sec + .5 : sec - 17.5); //trk->GetAlpha()
496 Double_t dx =(xmax - xmin)/nc;
497 for(Int_t ip=0; ip<nc; ip++){
498 fPoints[ip].SetXYZ(xmin, 0., 0.);
504 // select track model
505 if(BUILD || ((s&12) != (fTrackState&12))){
506 if(TESTBIT(s, kTrackCosmics)){
507 //printf("Straight track\n");
508 AliTRDtrackerV1::FitLine(trk, NULL, kFALSE, nc, fPoints);
510 if(TESTBIT(s, kTrackModel)){
511 //printf("Kalman track\n");
512 if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitKalman(trk, NULL, kFALSE, nc, fPoints);
514 //printf("Rieman track rim[%p] nc[%d]\n", (void*)fRim, nc);
515 //if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitRiemanTilt(trk, NULL, kTRUE, nc, fPoints);
517 for(Int_t ip = nc; ip--;){
518 x = fPoints[ip].GetX();
519 //printf("%2d x[%f] y[%f] z[%f]\n", ip, x, fRim->GetYat(x), fRim->GetZat(x));
520 fPoints[ip].SetXYZ(x, fRim->GetYat(x), fRim->GetZat(x));
526 for(Int_t ip=0; ip<nc; ip++){
527 fPoints[ip].Rotate(-fAlpha).GetXYZ(global);
528 SetPoint(ip, global[0], global[1], global[2]);
529 //printf("*** %2d x[%f] y[%f] z[%f]\n", ip, global[0], global[1], global[2]);
536 if(BUILD || ((s&3) != (fTrackState&3))){
537 if(TESTBIT(s, kSource)){
538 //printf("Source color\n");
539 if(fESDStatus&AliESDtrack::kTRDin){
540 SetMarkerColor(kGreen);
541 SetLineColor(kGreen);
543 SetMarkerColor(kMagenta);
544 SetLineColor(kMagenta);
547 if(TESTBIT(s, kPID) == AliTRDpidUtil::kLQ){
548 //printf("PID color kLQPID\n");
549 //trk->GetReconstructor()->SetOption("!nn");
551 //printf("PID color kNNPID\n");
552 //trk->GetReconstructor()->SetOption("nn");
556 Int_t species = 0; Float_t pid = 0.;
557 for(Int_t is=0; is<AliPID::kSPECIES; ++is)
558 if(trk->GetPID(is) > pid){
559 pid = trk->GetPID(is);
560 species = (AliPID::EParticleType) is;
563 case AliPID::kElectron:
564 SetMarkerColor(kRed);
568 SetMarkerColor(kBlue);
576 const Char_t *model = "line";
577 if(!TESTBIT(s, kTrackCosmics)){
578 if(TESTBIT(s, kTrackModel)) model = "kalman";
579 else model = "rieman";
581 Int_t species = 0; Float_t pid = 0.;
582 for(Int_t is=0; is<AliPID::kSPECIES; is++)
583 if(trk->GetPID(is) > pid){
584 pid = trk->GetPID(is);
589 "Tracklets[%d] Clusters[%d]\n"
590 "Reconstruction Source[%s]\n"
591 "PID[%4.1f %4.1f %4.1f %4.1f %4.1f]\n"
592 "MC[%d]", trk->GetNumberOfTracklets(), trk->GetNumberOfClusters(), fESDStatus&AliESDtrack::kTRDin ? "barrel" : "sa",
593 1.E2*trk->GetPID(0), 1.E2*trk->GetPID(1),
594 1.E2*trk->GetPID(2), 1.E2*trk->GetPID(3), 1.E2*trk->GetPID(4), trk->GetLabel()));
597 char id[6]; snprintf(id, 6, "%s", GetName());
598 SetName(Form("%s %s", id, AliPID::ParticleName(species)));
605 //______________________________________________________________________________
606 void AliEveTRDTrack::Load(const Char_t *what) const
608 // Spread downwards to tracklets the command "what"
610 const AliEveTRDTracklet* trklt(NULL);
611 TEveElement::List_ci itrklt=BeginChildren();
612 while(itrklt!=EndChildren()){
613 if((trklt = dynamic_cast<const AliEveTRDTracklet*>(*itrklt))) trklt->Load(what);
618 //______________________________________________________________________________
619 AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletMCM *tracklet) :
625 AliTRDtrackletMCM *trkl = new AliTRDtrackletMCM(*tracklet);
628 fDetector = trkl->GetDetector();
629 fROB = trkl->GetROB();
630 fMCM = trkl->GetMCM();
632 SetName("sim. tracklet");
633 SetTitle(Form("Det: %i, ROB: %i, MCM: %i, Label: %i\n0x%08x",
634 trkl->GetDetector(), trkl->GetROB(), trkl->GetMCM(), trkl->GetLabel(),
635 trkl->GetTrackletWord()));
636 SetLineColor(kGreen);
640 TGeoHMatrix *matrix = geo.GetClusterMatrix(trkl->GetDetector());
642 fDetector = trkl->GetDetector();
643 fROB = trkl->GetROB();
644 fMCM = trkl->GetMCM();
650 x[0] = AliTRDgeometry::AnodePos();
652 x[2] = trkl->GetLocalZ();
654 matrix->LocalToMaster(x, p);
655 geo.RotateBack(trkl->GetDetector(), p, p2);
656 SetPoint(0, p2[0], p2[1], p2[2]);
659 x[1] -= length * trkl->GetdYdX();
660 matrix->LocalToMaster(x, p);
661 p[2] *= p[0] / (p[0] + length);
662 geo.RotateBack(trkl->GetDetector(), p, p2);
663 SetPoint(1, p2[0], p2[1], p2[2]);
666 AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletWord *tracklet) :
672 AliTRDtrackletWord *trkl = new AliTRDtrackletWord(*tracklet);
675 fDetector = trkl->GetDetector();
676 fROB = trkl->GetROB();
677 fMCM = trkl->GetMCM();
679 SetName("raw tracklet");
680 SetTitle(Form("Det: %i, ROB: %i, MCM: %i, Label: %i\n0x%08x",
681 trkl->GetDetector(), fROB, fMCM, -1,
682 trkl->GetTrackletWord()));
687 TGeoHMatrix *matrix = geo.GetClusterMatrix(trkl->GetDetector());
693 x[0] = AliTRDgeometry::AnodePos();
695 x[2] = trkl->GetLocalZ();
697 matrix->LocalToMaster(x, p);
698 geo.RotateBack(trkl->GetDetector(), p, p2);
699 SetPoint(0, p2[0], p2[1], p2[2]);
702 x[1] -= length * trkl->GetdYdX();
703 matrix->LocalToMaster(x, p);
704 p[2] *= p[0] / (p[0] + length);
705 geo.RotateBack(trkl->GetDetector(), p, p2);
706 SetPoint(1, p2[0], p2[1], p2[2]);
709 AliEveTRDTrackletOnline::~AliEveTRDTrackletOnline()
711 delete ((AliTRDtrackletBase*) GetUserData());
715 void AliEveTRDTrackletOnline::ShowMCM(Option_t *opt) const
717 if (fDetector < 0 || fROB < 0 || fMCM < 0)
720 AliEveTRDmcm *evemcm = new AliEveTRDmcm();
721 evemcm->Init(fDetector, fROB, fMCM);
722 evemcm->LoadDigits();
725 TEveElementList *mcmlist = NULL;
726 if (gEve->GetCurrentEvent())
727 mcmlist = (TEveElementList*) gEve->GetCurrentEvent()->FindChild("TRD MCMs");
729 mcmlist = new TEveElementList("TRD MCMs");
730 gEve->AddElement(mcmlist);
732 gEve->AddElement(evemcm, mcmlist);
736 AliEveTRDmcm::AliEveTRDmcm() :
739 fMCM(new AliTRDmcmSim())
742 SetTitle("Unknown MCM");
745 AliEveTRDmcm::~AliEveTRDmcm()
750 Bool_t AliEveTRDmcm::Init(Int_t det, Int_t rob, Int_t mcm)
752 SetName(Form("MCM: Det. %i, ROB %i, MCM %i", det, rob, mcm));
753 SetTitle(Form("MCM: Det. %i, ROB %i, MCM %i", det, rob, mcm));
754 fMCM->Init(det, rob, mcm);
759 Bool_t AliEveTRDmcm::LoadDigits()
761 AliRunLoader *rl = AliEveEventManager::AssertRunLoader();
762 return fMCM->LoadMCM(rl, fMCM->GetDetector(),
763 fMCM->GetRobPos(), fMCM->GetMcmPos());
766 Bool_t AliEveTRDmcm::Filter()
772 Bool_t AliEveTRDmcm::Tracklet()
778 void AliEveTRDmcm::Draw(Option_t* option)
780 const char *mcmname = Form("mcm_%i_%i_%i", fMCM->GetDetector(),
781 fMCM->GetRobPos(), fMCM->GetMcmPos());
783 TCanvas *c = dynamic_cast<TCanvas*> (gROOT->FindObject(mcmname));
785 c = gEve->AddCanvasTab("TRD MCM");
786 c->SetTitle(Form("MCM %i on ROB %i of det. %i",
787 fMCM->GetMcmPos(), fMCM->GetRobPos(), fMCM->GetDetector()));
793 Bool_t AliEveTRDmcm::AssignPointer(const char* ptrname)
795 gROOT->ProcessLine(Form("AliTRDmcmSim* %s = (AliTRDmcmSim *)%p", ptrname, (void*)fMCM));