]> git.uio.no Git - u/mrichter/AliRoot.git/blobdiff - EVE/EveDet/AliEveTRDData.cxx
add identifiers to digits for readable histos
[u/mrichter/AliRoot.git] / EVE / EveDet / AliEveTRDData.cxx
index f10a5ec91d1dfc4a44a14d94243b3f3671c41de9..ed88bec10573fb4337c509df35db37710a2d6ab7 100644 (file)
@@ -8,9 +8,12 @@
  **************************************************************************/
 
 #include "TROOT.h"
+#include "TStyle.h"
 #include "TVector.h"
 #include "TLinearFitter.h"
 #include "TCanvas.h"
+#include "TGeoMatrix.h"
+#include "TGeoManager.h"
 
 #include "TEveTrans.h"
 #include "TEveManager.h"
@@ -22,6 +25,7 @@
 #include "AliEveTRDLoader.h"
 #include "AliEveTRDLoaderImp.h"
 
+#include "AliGeomManager.h"
 #include "AliESDtrack.h"
 #include "AliLog.h"
 #include "AliPID.h"
@@ -59,10 +63,15 @@ ClassImp(AliEveTRDmcm)
 ///////////////////////////////////////////////////////////
 
 //______________________________________________________________________________
-AliEveTRDDigits::AliEveTRDDigits(AliEveTRDChamber *p) :
-  TEveQuadSet("digits", ""), fParent(p), fBoxes(), fData()
+AliEveTRDDigits::AliEveTRDDigits(AliEveTRDChamber *p) 
+  :TEveQuadSet("digits", "")
+  ,fParent(p)
 {
   // Constructor.
+  SetOwnIds(kTRUE);
+  gStyle->SetPalette(1, 0);
+  SetPalette(new TEveRGBAPalette(0, 512));
+  Reset(TEveQuadSet::kQT_RectangleYZ, kFALSE, 32);
 }
 
 //______________________________________________________________________________
@@ -71,72 +80,12 @@ AliEveTRDDigits::~AliEveTRDDigits()
 //  AliInfo(GetTitle());
 }
 
-//______________________________________________________________________________
-void AliEveTRDDigits::ComputeRepresentation()
-{
-  // Calculate digits representation according to user settings. The
-  // user can set the following parameters:
-  // - digits scale (log/lin)
-  // - digits threshold
-  // - digits apparence (quads/boxes)
-
-  if(!fData.HasData()){
-    return;
-  }
-
-  TEveQuadSet::Reset(TEveQuadSet::kQT_RectangleYZ, kTRUE, 64);
-
-  Double_t scale, dy, dz;
-  Int_t q, color;
-  Int_t nrows = fData.GetNrow(),
-        ncols = fData.GetNcol(),
-        ntbs  = fData.GetNtime(),
-        det   = fParent->GetID();
-  Float_t threshold = fParent->GetDigitsThreshold();
-
-  AliTRDtransform transform(det);
-  AliTRDgeometry *geo = fParent->fGeo;
-  AliTRDpadPlane *pp = geo->GetPadPlane(geo->GetLayer(det), geo->GetStack(det));
-
-  // express position in tracking coordinates
-  fData.Expand();
-  for (Int_t ir = 0; ir < nrows; ir++) {
-    dz = pp->GetRowSize(ir);
-    for (Int_t ic = 0; ic < ncols; ic++) {
-      dy = pp->GetColSize(ic);
-      for (Int_t it = 0; it < ntbs; it++) {
-        q = fData.GetData(ir, ic, it);
-        if (q < threshold) continue;
-
-        Double_t x[6] = {0., 0., Double_t(q), 0., 0., 0.}; 
-        //Int_t  roc[3] = {ir, ic, 0}; 
-        //Bool_t    out = kTRUE;
-        //transform.Transform(&x[0], &roc[0], UInt_t(it), out, 0);
-
-        scale = q < 512 ? q/512. : 1.;
-        color  = 50+int(scale*50.);
-    
-        AddQuad(x[1]-.45*dy, x[2]-.5*dz*scale, x[0], .9*dy, dz*scale);
-        QuadValue(q);
-        QuadColor(color);
-        QuadId(new TNamed(Form("Charge%d", q), "dummy title"));
-      }  // end time loop
-    }  // end col loop
-  }  // end row loop
-  fData.Compress();
-  
-  // rotate to global coordinates
-  //RefitPlex();
-  TEveTrans& t = RefMainTrans();
-  t.SetRotByAngles((geo->GetSector(det)+.5)*AliTRDgeometry::GetAlpha(), 0.,0.);
-}
-
 //______________________________________________________________________________
 void AliEveTRDDigits::SetData(AliTRDdigitsManager *digits)
 {
   // Set data source.
 
-  Int_t det = fParent->GetID();
+  Int_t det(fParent->GetID());
   AliTRDarrayADC *data = digits->GetDigits(det);
   if(!data->GetDim()) return;
   data->Expand();
@@ -144,47 +93,62 @@ void AliEveTRDDigits::SetData(AliTRDdigitsManager *digits)
   AliTRDSignalIndex *indexes = digits->GetIndexes(det);
   if(!indexes->IsAllocated()) digits->BuildIndexes(det);
 
-  if(!fData.HasData()) fData.Allocate(data->GetNrow(), data->GetNcol(), data->GetNtime());
-  fData.Expand();
+  Double_t scale, dy, dz;
+  Int_t ly    = AliTRDgeometry::GetLayer(det),
+        stk   = AliTRDgeometry::GetStack(det),
+        sec   = AliTRDgeometry::GetSector(det),
+        vid   = AliGeomManager::LayerToVolUID(AliGeomManager::kTRD1 + ly, stk + AliTRDgeometry::Nstack() * sec);
+  SetNameTitle(Form("digits%03d", det), Form("det %03d [%02d_%d_%d]", sec, stk, ly));
+  Short_t sig[7]={0,0,0,10,0,0,0};
 
-  Int_t row, col, time, adc;
+  AliTRDtransform transform(det);
+  AliTRDpadPlane *pp(fParent->fGeo->GetPadPlane(ly, stk));
+
+  Int_t row, col;
+  AliTRDcluster c;
   indexes->ResetCounters();
   while (indexes->NextRCIndex(row, col)){
-    indexes->ResetTbinCounter();
-    while (indexes->NextTbinIndex(time)){
-      if(data->IsPadCorrupted(row, col, time)){
-        // we should mark this position
-        break;
-      }
-      adc = data->GetData(row, col, time);
-      if(adc <= 1) continue;
-      fData.SetData(row, col, time, adc);
-      //fIndex->AddIndexTBin(row,col,time);
-      //printf("\tr[%d] c[%d] t[%d] ADC[%d]\n", row, col, time, adc);
+    dz = pp->GetRowSize(row);
+    dy = pp->GetColSize(col);
+    Short_t *const adc = data->GetDataAddress(row,col);
+    for (Int_t time(0); time<30; time++){     
+      if(data->IsPadCorrupted(row, col, time)) break;
+      if(adc[time] <= 1) continue;
+      new (&c) AliTRDcluster(det, col, row, time, sig, vid);
+      transform.Transform(&c);
+
+      scale = adc[time] < 512 ? adc[time]/512. : 1.;
+      AddQuad(c.GetY()-0.5*dy, c.GetZ()-0.5*dz*scale, c.GetX(), dy*0.95, dz*scale);
+      QuadValue(Float_t(adc[time]));
+      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)));
     } 
   }
-  fData.Compress();
-}
-
-
-//______________________________________________________________________________
-void AliEveTRDDigits::Paint(Option_t *option)
-{
-  // Paint the object.
 
-  if(fParent->GetDigitsBox()) fBoxes.Paint(option);
-  else TEveQuadSet::Paint(option);
+  // rotate to global coordinates
+  RefitPlex();
+  TEveTrans& t = RefMainTrans();
+  t.SetRotByAngles((sec+.5)*AliTRDgeometry::GetAlpha(), 0.,0.);
 }
 
-//______________________________________________________________________________
-void AliEveTRDDigits::Reset()
-{
-  // Reset raw and visual data.
 
-  TEveQuadSet::Reset(TEveQuadSet::kQT_RectangleYZ, kTRUE, 64);
-  // MT fBoxes.fBoxes.clear();
-  fData.Reset();
-}
+// //______________________________________________________________________________
+// void AliEveTRDDigits::Paint(Option_t *option)
+// {
+//   // Paint the object.
+// 
+//   if(fParent->GetDigitsBox()) fBoxes.Paint(option);
+//   else TEveQuadSet::Paint(option);
+// }
+
+// //______________________________________________________________________________
+// void AliEveTRDDigits::Reset()
+// {
+//   // Reset raw and visual data.
+// 
+//   TEveQuadSet::Reset(TEveQuadSet::kQT_RectangleYZ, kTRUE, 64);
+//   // MT fBoxes.fBoxes.clear();
+//   fData.Reset();
+// }
 
 ///////////////////////////////////////////////////////////
 /////////////   AliEveTRDHits         /////////////////////
@@ -210,7 +174,7 @@ void AliEveTRDHits::PointSelected(Int_t n)
 {
   // Handle an individual point selection from GL.
 
-  AliTRDhit *h = 0x0;
+  AliTRDhit *h = NULL;
   if(!(h = dynamic_cast<AliTRDhit*>(GetPointId(n)))) return;
   printf("Id[%3d] Det[%3d] Reg[%c] TR[%c] Q[%3d] MC[%d] t[%f]\n", 
     n, h->GetDetector(), 
@@ -242,16 +206,9 @@ void AliEveTRDClusters::PointSelected(Int_t n)
   // Handle an individual point selection from GL.
 
   AliTRDcluster *c = dynamic_cast<AliTRDcluster*>(GetPointId(n));
-  printf("\nDetector             : %d\n", c->GetDetector());
-  printf("Charge               : %f\n", c->GetQ());
-  printf("Sum S                : %4.0f\n", c->GetSumS());
-  printf("Time bin             : %d\n", c->GetLocalTimeBin());
-  printf("Signals              : ");
-  Short_t *cSignals = c->GetSignals();
-  for(Int_t ipad=0; ipad<7; ipad++) printf("%d ", cSignals[ipad]); printf("\n");
-  printf("Central pad          : %d\n", c->GetPadCol());
-  printf("MC track labels      : ");
-  for(Int_t itrk=0; itrk<3; itrk++) printf("%d ", c->GetLabel(itrk)); printf("\n");
+  if(!c) return;
+  c->Print();
+  Emit("PointSelected(Int_t)", n);
   // Bool_t    AliCluster::GetGlobalCov(Float_t* cov) const
   // Bool_t    AliCluster::GetGlobalXYZ(Float_t* xyz) const
   // Float_t   AliCluster::GetSigmaY2() const
@@ -262,7 +219,7 @@ void AliEveTRDClusters::PointSelected(Int_t n)
 //______________________________________________________________________________
 void AliEveTRDClusters::Print(Option_t *o) const
 {
-  AliTRDcluster *c = 0x0;
+  AliTRDcluster *c = NULL;
 
   for(Int_t n = GetN(); n--;){
     if(!(c = dynamic_cast<AliTRDcluster*>(GetPointId(n)))) continue;
@@ -271,7 +228,7 @@ void AliEveTRDClusters::Print(Option_t *o) const
 }
 
 //______________________________________________________________________________
-void AliEveTRDClusters::Load(Char_t *w, Bool_t stk) const
+void AliEveTRDClusters::Load(Char_t *w) const
 {
   Int_t typ = -1;
   if(strcmp(w, "hit")==0) typ = 0;
@@ -287,13 +244,13 @@ void AliEveTRDClusters::Load(Char_t *w, Bool_t stk) const
     return;
   }
 
-  AliTRDcluster *c = 0x0;
+  AliTRDcluster *c = NULL;
   Int_t n = 0;
   while((n = GetN() && !(c = dynamic_cast<AliTRDcluster*>(GetPointId(n))))) n++;
   if(!c) return;
 
   Int_t det = c->GetDetector();
-  AliEveTRDLoader *loader = 0x0;
+  AliEveTRDLoader *loader = NULL;
   switch(typ){
   case 0:  
     loader = new AliEveTRDLoader("Hits");
@@ -330,12 +287,12 @@ void AliEveTRDClusters::Load(Char_t *w, Bool_t stk) const
   default: return;
   }
 
-  loader->AddChambers(AliTRDgeometry::GetSector(det),AliTRDgeometry::GetStack(det), stk ? -1 : AliTRDgeometry::GetLayer(det));
+  loader->AddChambers(AliTRDgeometry::GetSector(det),AliTRDgeometry::GetStack(det), AliTRDgeometry::GetLayer(det));
   // load first event
   loader->GoToEvent(AliEveEventManager::GetCurrent()->GetEventId());
   
   // register loader with alieve
-  gEve->AddElement(loader);
+  gEve->AddElement(loader->GetChamber(det), *(BeginParents()));
   //loader->SpawnEditor();
   gEve->Redraw3D();
 }
@@ -346,11 +303,15 @@ void AliEveTRDClusters::Load(Char_t *w, Bool_t stk) const
 
 //______________________________________________________________________________
 AliEveTRDTracklet::AliEveTRDTracklet(AliTRDseedV1 *trklt):TEveLine()
-  ,fClusters(0x0)
+  ,fClusters(NULL)
 {
   // Constructor.
   SetName("tracklet");
   
+  if(!gGeoManager){ 
+    AliEveEventManager::AssertGeometry();
+    AliInfo(Form("gGeo[%p] Closed[%c]", (void*)gGeoManager, gGeoManager->IsClosed()?'y':'n'));
+  }
   SetUserData(trklt);
   Float_t dx;
   Float_t x0   = trklt->GetX0();
@@ -358,9 +319,9 @@ AliEveTRDTracklet::AliEveTRDTracklet(AliTRDseedV1 *trklt):TEveLine()
   Float_t z0   = trklt->GetZref(0);
   Float_t dydx = trklt->GetYref(1);
   Float_t dzdx = trklt->GetZref(1);
-  Float_t tilt = trklt->GetTilt();
+  Float_t  tilt = trklt->GetTilt();
   Float_t g[3];
-  AliTRDcluster *c = 0x0;
+  AliTRDcluster *c = NULL;
   for(Int_t ic=0; ic<AliTRDseedV1::kNclusters; ic++){
     if(!(c = trklt->GetClusters(ic))) continue;
     if(!fClusters) AddElement(fClusters = new AliEveTRDClusters());
@@ -371,17 +332,17 @@ AliEveTRDTracklet::AliEveTRDTracklet(AliTRDseedV1 *trklt):TEveLine()
     Float_t yc = c->GetY(); 
     c->SetY(yc-tilt*(c->GetZ()-zt));
     c->GetGlobalXYZ(g); 
-    Int_t id = fClusters->SetNextPoint(g[0], g[1], g[2]);    
-    //Int_t id = fClusters->SetNextPoint(c->GetX(), c->GetY(), c->GetZ());    
+    Int_t id = fClusters->SetNextPoint(g[0], g[1], g[2]);
     c->SetY(yc);
     fClusters->SetPointId(id, new AliTRDcluster(*c));
   } 
   if(fClusters){
-    fClusters->SetTitle(Form("N[%d]", trklt->GetN2()));
+    fClusters->SetName("TC clusters");
+    fClusters->SetTitle(Form("N[%d]", trklt->GetN()));
     fClusters->SetMarkerColor(kMagenta);
   }
 
-  SetTitle(Form("Det[%d] Plane[%d] P[%7.3f]", trklt->GetDetector(), trklt->GetPlane(), trklt->GetMomentum()));
+  SetTitle(Form("Det[%d] RC[%c] Layer[%d] P[%7.3f]", trklt->GetDetector(), trklt->IsRowCross()?'y':'n', trklt->GetPlane(), trklt->GetMomentum()));
   SetLineColor(kRed);
   //SetOwnIds(kTRUE);
   
@@ -432,26 +393,26 @@ AliEveTRDTrack::AliEveTRDTrack(AliTRDtrackV1 *trk)
   ,fTrackState(0)
   ,fESDStatus(0)
   ,fAlpha(0.)
-  ,fPoints(0x0)
-  ,fRim(0x0)
+  ,fPoints(NULL)
+  ,fRim(NULL)
 {
   // Constructor.
   SetUserData(trk);
   SetName("");
 
-  AliTRDtrackerV1::SetNTimeBins(24);
+//  AliTRDtrackerV1::SetNTimeBins(24);
 
   fRim = new AliRieman(trk->GetNumberOfClusters());
-  AliTRDseedV1 *tracklet = 0x0;
+  AliTRDseedV1 *tracklet = NULL;
   for(Int_t il=0; il<AliTRDgeometry::kNlayer; il++){
     if(!(tracklet = trk->GetTracklet(il))) continue;
     if(!tracklet->IsOK()) continue;
     AddElement(new AliEveTRDTracklet(tracklet));
 
-    AliTRDcluster *c = 0x0;
 //     tracklet->ResetClusterIter(kFALSE);
 //     while((c = tracklet->PrevCluster())){
-    for(Int_t ic=AliTRDseedV1::kNtb; ic--;){
+//    AliTRDcluster *c(NULL);
+/*    for(Int_t ic=AliTRDseedV1::kNtb; ic--;){
       if(!(c=tracklet->GetClusters(ic))) continue;
       Float_t xc = c->GetX();
       Float_t yc = c->GetY();
@@ -459,7 +420,7 @@ AliEveTRDTrack::AliEveTRDTrack(AliTRDtrackV1 *trk)
       Float_t zt = tracklet->GetZref(0) - (tracklet->GetX0()-xc)*tracklet->GetZref(1); 
       yc -= tracklet->GetTilt()*(zc-zt);
       fRim->AddPoint(xc, yc, zc, .05, 2.3);
-    }
+    }*/
   }
   if(trk->GetNumberOfTracklets()>1) fRim->Update();
   SetStatus(fTrackState);
@@ -468,7 +429,7 @@ AliEveTRDTrack::AliEveTRDTrack(AliTRDtrackV1 *trk)
 //______________________________________________________________________________
 AliEveTRDTrack::~AliEveTRDTrack()
 {
-  if(fPoints) delete [] fPoints; fPoints = 0x0;
+  if(fPoints) delete [] fPoints; fPoints = NULL;
   //delete dynamic_cast<AliTRDtrackV1*>(GetUserData());
 }
 
@@ -485,9 +446,13 @@ void AliEveTRDTrack::SetStatus(UChar_t s)
 {
   // nothing to be done
   if(fPoints && fTrackState == s) return;
-
+  //return;
   const Int_t nc = AliTRDtrackV1::kMAXCLUSTERSPERTRACK;
-  AliTRDtrackV1 *trk = (AliTRDtrackV1*)GetUserData();
+  AliTRDtrackV1 *trk(NULL);
+  if(!(trk=static_cast<AliTRDtrackV1*>(GetUserData()))) {
+    AliError("Failed casting data to TRD track.");
+    return;
+  }
 
   Bool_t BUILD = kFALSE;
   if(!fPoints){ 
@@ -496,11 +461,11 @@ void AliEveTRDTrack::SetStatus(UChar_t s)
     // define the radial span of the track in the TRD
     Double_t xmin = -1., xmax = -1.;
     Int_t det = 0;
-    AliTRDseedV1 *trklt = 0x0;
-    for(Int_t ily=AliTRDgeometry::kNlayer; ily--;){
+    AliTRDseedV1 *trklt = NULL;
+    for(Int_t ily=0; ily<AliTRDgeometry::kNlayer; ily++){
       if(!(trklt = trk->GetTracklet(ily))) continue;
       if(xmin<0.) xmin = trklt->GetX0() - AliTRDgeometry::CamHght() - AliTRDgeometry::CdrHght();
-      if(xmax<0.) xmax = trklt->GetX0();
+      if(trklt->GetX0()>xmax) xmax = trklt->GetX0();
       det = trklt->GetDetector();
     }
     Int_t sec = det/AliTRDgeometry::kNdets;
@@ -518,27 +483,29 @@ void AliEveTRDTrack::SetStatus(UChar_t s)
   if(BUILD || ((s&12) != (fTrackState&12))){
     if(TESTBIT(s, kTrackCosmics)){
       //printf("Straight track\n");
-      AliTRDtrackerV1::FitLine(trk, 0x0, kFALSE, nc, fPoints);
+      AliTRDtrackerV1::FitLine(trk, NULL, kFALSE, nc, fPoints);
     } else {
       if(TESTBIT(s, kTrackModel)){
         //printf("Kalman track\n");
-        if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitKalman(trk, 0x0, kFALSE, nc, fPoints);
+        if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitKalman(trk, NULL, kFALSE, nc, fPoints);
       } else { 
-        //printf("Rieman track\n");
-        // if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitRiemanTilt(trk, 0x0, kTRUE, nc, fPoints);
-
-/*        Float_t x = 0.;
+        //printf("Rieman track rim[%p] nc[%d]\n", (void*)fRim, nc);
+        //if(trk->GetNumberOfTracklets() >=4) AliTRDtrackerV1::FitRiemanTilt(trk, NULL, kTRUE, nc, fPoints);
+        Float_t x = 0.;
         for(Int_t ip = nc; ip--;){
           x = fPoints[ip].GetX();
+          //printf("%2d x[%f] y[%f] z[%f]\n", ip, x, fRim->GetYat(x), fRim->GetZat(x));
           fPoints[ip].SetXYZ(x, fRim->GetYat(x), fRim->GetZat(x));
-        }*/
+        }
       }
     }
   
     Float_t global[3];
-    for(Int_t ip=0; ip<0/*nc*/; ip++){
+    for(Int_t ip=0; ip<nc; ip++){
       fPoints[ip].Rotate(-fAlpha).GetXYZ(global);
       SetPoint(ip, global[0], global[1], global[2]);
+      //printf("*** %2d x[%f] y[%f] z[%f]\n", ip, global[0], global[1], global[2]);
+    
     }
     SetSmooth(kTRUE);
   }
@@ -562,7 +529,7 @@ void AliEveTRDTrack::SetStatus(UChar_t s)
         //printf("PID color kNNPID\n");
         //trk->GetReconstructor()->SetOption("nn");
       }
-      trk->CookPID();
+      //trk->CookPID();
   
       Int_t species = 0; Float_t pid = 0.;
       for(Int_t is=0; is<AliPID::kSPECIES; is++) 
@@ -613,6 +580,16 @@ void AliEveTRDTrack::SetStatus(UChar_t s)
   fTrackState = s;
 }
 
+//______________________________________________________________________________
+void AliEveTRDTrack::Load(Char_t *what) const
+{
+  TEveElement::List_ci itrklt=BeginChildren();
+  while(itrklt!=EndChildren()){
+    dynamic_cast<const AliEveTRDTracklet*>(*itrklt)->Load(what);
+    itrklt++;
+  }
+}
+
 
 AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletMCM *tracklet) :
   TEveLine(),
@@ -627,41 +604,38 @@ AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletMCM *tracklet) :
   fROB = trkl->GetROB();
   fMCM = trkl->GetMCM();
 
-  SetName("TRD tracklet");
+  SetName("sim. tracklet");
   SetTitle(Form("Det: %i, ROB: %i, MCM: %i, Label: %i\n0x%08x", 
                 trkl->GetDetector(), trkl->GetROB(), trkl->GetMCM(), trkl->GetLabel(),
                 trkl->GetTrackletWord()));
   SetLineColor(kGreen);
+  SetLineWidth(3);
 
-  AliTRDgeometry *geo = new AliTRDgeometry();
-//  TGeoHMatrix *matrix = geo->GetClusterMatrix(trkl->GetDetector());
+  AliTRDgeometry geo;
+  TGeoHMatrix *matrix = geo.GetClusterMatrix(trkl->GetDetector());
 
+  fDetector = trkl->GetDetector();
+  fROB = trkl->GetROB();
+  fMCM = trkl->GetMCM();
+  
   Float_t length = 3.;
   Double_t x[3];
   Double_t p[3];
-  x[0] = trkl->GetX();
+  Double_t p2[3];
+  x[0] = AliTRDgeometry::AnodePos(); 
   x[1] = trkl->GetY();
-  x[2] = trkl->GetZ();
+  x[2] = trkl->GetLocalZ();
 
-  fDetector = trkl->GetDetector();
-  AliTRDpadPlane *pp = geo->GetPadPlane(geo->GetLayer(fDetector), geo->GetStack(fDetector));
-  fROB = 2 * (trkl->GetZbin() / 4) + (trkl->GetY() > 0 ? 1 : 0);
-  fMCM = (((Int_t) ((trkl->GetY()) / pp->GetWidthIPad()) + 72) / 18) % 4 
-    + 4 * (trkl->GetZbin() % 4) ;
-  AliInfo(Form("From position/tracklet: ROB: %i/%i, MCM: %i/%i", 
-               fROB, trkl->GetROB(), fMCM, trkl->GetMCM()));  
-  
-  geo->RotateBack(trkl->GetDetector(), x, p);
-//  matrix->LocalToMaster(x, p);
-  SetPoint(0, p[0], p[1], p[2]);
+  matrix->LocalToMaster(x, p);
+  geo.RotateBack(trkl->GetDetector(), p, p2);
+  SetPoint(0, p2[0], p2[1], p2[2]);
 
   x[0] -= length;
   x[1] -= length * trkl->GetdYdX();
-  x[2] *= x[0] / (x[0] + length);
-  geo->RotateBack(trkl->GetDetector(), x, p);
-//  matrix->LocalToMaster(x, p);
-  SetPoint(1, p[0], p[1], p[2]);
-  delete geo;
+  matrix->LocalToMaster(x, p);
+  p[2] *= p[0] / (p[0] + length);
+  geo.RotateBack(trkl->GetDetector(), p, p2);
+  SetPoint(1, p2[0], p2[1], p2[2]);
 }
 
 AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletWord *tracklet) :
@@ -673,47 +647,44 @@ AliEveTRDTrackletOnline::AliEveTRDTrackletOnline(AliTRDtrackletWord *tracklet) :
   AliTRDtrackletWord *trkl = new AliTRDtrackletWord(*tracklet);
   SetUserData(trkl);
 
-  AliTRDgeometry *geo = new AliTRDgeometry();
   fDetector = trkl->GetDetector();
-  AliTRDpadPlane *pp = geo->GetPadPlane(geo->GetLayer(fDetector), geo->GetStack(fDetector));
-  fROB = 2 * (trkl->GetZbin() / 4) + (trkl->GetY() > 0 ? 1 : 0);
-  fMCM = (((Int_t) ((trkl->GetY()) / pp->GetWidthIPad()) + 72) / 18) % 4 
-    + 4 * (trkl->GetZbin() % 4) ;
+  fROB = trkl->GetROB(); 
+  fMCM = trkl->GetMCM(); 
 
-  SetName("TRD tracklet");
+  SetName("raw tracklet");
   SetTitle(Form("Det: %i, ROB: %i, MCM: %i, Label: %i\n0x%08x", 
                 trkl->GetDetector(), fROB, fMCM, -1,
                 trkl->GetTrackletWord()));
   SetLineColor(kRed);
+  SetLineWidth(3);
 
-//  AliTRDgeometry *geo = new AliTRDgeometry();
+  AliTRDgeometry geo;
+  TGeoHMatrix *matrix = geo.GetClusterMatrix(trkl->GetDetector());
 
   Float_t length = 3.;
   Double_t x[3];
   Double_t p[3];
-  x[0] = trkl->GetX();
+  Double_t p2[3];
+  x[0] = AliTRDgeometry::AnodePos();
   x[1] = trkl->GetY();
-  x[2] = trkl->GetZ();
+  x[2] = trkl->GetLocalZ();
   
-  geo->RotateBack(trkl->GetDetector(), x, p);
-  SetPoint(0, p[0], p[1], p[2]);
+  matrix->LocalToMaster(x, p);
+  geo.RotateBack(trkl->GetDetector(), p, p2);
+  SetPoint(0, p2[0], p2[1], p2[2]);
 
   x[0] -= length;
   x[1] -= length * trkl->GetdYdX();
-  x[2] *= x[0] / (x[0] + length);
-  geo->RotateBack(trkl->GetDetector(), x, p);
-  SetPoint(1, p[0], p[1], p[2]);
-  delete geo;
+  matrix->LocalToMaster(x, p);
+  p[2] *= p[0] / (p[0] + length);
+  geo.RotateBack(trkl->GetDetector(), p, p2);
+  SetPoint(1, p2[0], p2[1], p2[2]);
 }
 
 AliEveTRDTrackletOnline::~AliEveTRDTrackletOnline() 
 {
-  AliTRDtrackletMCM *trkl = dynamic_cast<AliTRDtrackletMCM*> ((AliTRDtrackletBase*) GetUserData());
-  printf("trkl: %p\n", (void*)trkl);
-//  delete trkl;
-  AliTRDtrackletWord *trklWord = dynamic_cast<AliTRDtrackletWord*> ((AliTRDtrackletBase*) GetUserData());
-  printf("trklWord: %p\n", (void*)trklWord);
-//  delete trklWord;
+  delete ((AliTRDtrackletBase*) GetUserData());
+  SetUserData(0x0);
 }
 
 void AliEveTRDTrackletOnline::ShowMCM(Option_t *opt) const
@@ -726,7 +697,7 @@ void AliEveTRDTrackletOnline::ShowMCM(Option_t *opt) const
   evemcm->LoadDigits();
   evemcm->Draw(opt);
 
-  TEveElementList *mcmlist = 0x0;
+  TEveElementList *mcmlist = NULL;
   if (gEve->GetCurrentEvent()) 
     mcmlist = (TEveElementList*) gEve->GetCurrentEvent()->FindChild("TRD MCMs");
   if (!mcmlist) {
@@ -796,6 +767,6 @@ void AliEveTRDmcm::Draw(Option_t* option)
 
 Bool_t AliEveTRDmcm::AssignPointer(const char* ptrname)
 {
-  gROOT->ProcessLine(Form("AliTRDmcmSim* %s = (AliTRDmcmSim *) 0x%x", ptrname, fMCM));
+  gROOT->ProcessLine(Form("AliTRDmcmSim* %s = (AliTRDmcmSim *)%p", ptrname, (void*)fMCM));
   return kTRUE;
 }