The position of the trd clusters in EVE is transformed with different
authorabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 May 2012 13:01:51 +0000 (13:01 +0000)
committerabercuci <abercuci@f7af4fe6-9843-0410-8265-dc069ae4e863>
Thu, 10 May 2012 13:01:51 +0000 (13:01 +0000)
matrices depending on the flag fIsMisaligned. As figured out in the
past, the flag is not set correctly or this information is lost at some
point. Just a reminder:
https://indico.cern.ch/getFile.py/access?contribId=1&resId=0&materialId=slides&confId=190673

Adapted the trd_cluster.C macro to just skip the flag and take the
correct transformation. Now, in May 2012, the flag seems to be still
wrong (see slides 1-4 in):

https://indico.cern.ch/getFile.py/access?contribId=0&resId=0&materialId=slides&confId=190673

Ben

EVE/alice-macros/trd_clusters.C

index c8eb0d1..51a7c6f 100644 (file)
@@ -7,20 +7,18 @@
  * full copyright notice.                                                 *
  **************************************************************************/
 
-#if !defined(__CINT__) || defined(__MAKECINT__)
-#include <TObjArray.h>
-#include <TTree.h>
+#ifdef __CINT__
+class TEvePointSet;
+class TEveElement;
+#else
 #include <TEveManager.h>
-#include <TEveElement.h>
 #include <TEvePointSet.h>
+#include <TTree.h>
+#include <EveBase/AliEveEventManager.h>
 
-#include <AliCluster.h>
-#include <AliRunLoader.h>
-#include <AliTRDcluster.h>
-#include <AliEveEventManager.h>
-#else
-class TEvePointSet;
-class TEveElement;
+#include "AliRunLoader.h"
+#include "AliCluster.h"
+#include "TRD/AliTRDcluster.h"
 #endif
 
 TEvePointSet* trd_clusters(TEveElement *cont = 0)
@@ -42,7 +40,6 @@ TEvePointSet* trd_clusters(TEveElement *cont = 0)
   TEvePointSet *clusters = new TEvePointSet(kMaxClusters);
   clusters->SetOwnIds(kTRUE);
 
-
   Int_t nentr=(Int_t)recPoints->GetEntries();
   for (Int_t i=0; i<nentr; i++) {
     if (!recPoints->GetEvent(i)) continue;
@@ -51,9 +48,21 @@ TEvePointSet* trd_clusters(TEveElement *cont = 0)
 
     while (ncl--) {
       AliTRDcluster *c = (AliTRDcluster*)TRDcluster->UncheckedAt(ncl);
-      Float_t g[3]; //global coordinates
-      c->GetGlobalXYZ(g);
-      clusters->SetNextPoint(g[0], g[1], g[2]);
+      //Float_t g[3]; //global coordinates
+      //c->GetGlobalXYZ(g);
+      
+      Int_t fVolumeId = c->GetVolumeId();
+      const TGeoHMatrix *mt =AliGeomManager::GetTracking2LocalMatrix(fVolumeId);;
+      Double_t txyz[3] = {c->GetX(), c->GetY(), c->GetZ()};
+      Double_t lxyz[3] = {0, 0, 0};
+      mt->LocalToMaster(txyz,lxyz);
+   
+      TGeoHMatrix *mlIdeal = AliGeomManager::GetOrigGlobalMatrix(fVolumeId);
+      Double_t gxyzIdeal[3] = {0, 0, 0};
+      mlIdeal->LocalToMaster(lxyz,gxyzIdeal);
+      
+      clusters->SetNextPoint(gxyzIdeal[0], gxyzIdeal[1], gxyzIdeal[2]);
+      
       AliCluster *atp = new AliCluster(*c);
       clusters->SetPointId(atp);
     }