-// $Id: tpc_clusters.C 23497 2008-01-23 20:43:14Z mtadel $
+// $Id$
// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
/**************************************************************************
#include "AliRunLoader.h"
#include "AliCluster.h"
-#include "AliTRDcluster.h"
+#include "TRD/AliTRDcluster.h"
#endif
-TEvePointSet* trd_clusters(TEveElement *cont = 0){
- const Int_t kMaxClusters = 18 * 6 * 24 *10;
- AliEveEventManager::AssertGeometry();
+TEvePointSet* trd_clusters(TEveElement *cont = 0)
+{
+ const Int_t kMaxClusters = 18 * 6 * 24 *10;
- TEvePointSet *clusters = new TEvePointSet(kMaxClusters);
- clusters->SetOwnIds(kTRUE);
+ AliEveEventManager::AssertGeometry();
- AliRunLoader *rl = AliEveEventManager::AssertRunLoader();
- rl->LoadRecPoints("TRD");
+ AliRunLoader *rl = AliEveEventManager::AssertRunLoader();
+ rl->LoadRecPoints("TRD");
+
+ TTree *recPoints = rl->GetTreeR("TRD", kFALSE);
+ if (recPoints == 0)
+ return 0;
+
+ TObjArray *TRDcluster = 0x0;
+ recPoints->SetBranchAddress("TRDcluster", &TRDcluster);
+
+ TEvePointSet *clusters = new TEvePointSet(kMaxClusters);
+ clusters->SetOwnIds(kTRUE);
- TObjArray *TRDcluster = 0x0;
- TTree *recPoints = rl->GetTreeR("TRD", kFALSE);
- recPoints->SetBranchAddress("TRDcluster", &TRDcluster);
Int_t nentr=(Int_t)recPoints->GetEntries();
for (Int_t i=0; i<nentr; i++) {
AliTRDcluster *c = (AliTRDcluster*)TRDcluster->UncheckedAt(ncl);
Float_t g[3]; //global coordinates
c->GetGlobalXYZ(g);
- clusters->SetNextPoint(g[0], g[1], g[2]);
- AliCluster *atp = new AliCluster(*c);
- clusters->SetPointId(atp);
+ clusters->SetNextPoint(g[0], g[1], g[2]);
+ AliCluster *atp = new AliCluster(*c);
+ clusters->SetPointId(atp);
}
TRDcluster->Clear();
}
+ rl->UnloadRecPoints("TRD");
+
if(clusters->Size() == 0 && gEve->GetKeepEmptyCont() == kFALSE) {
- Warning("trd_clusters", "No TRD clusters");
+ Warning("trd_clusters.C", "No TRD clusters");
delete clusters;
return 0;
}
- clusters->SetMarkerStyle(2);
- clusters->SetMarkerSize(0.2);
- clusters->SetMarkerColor(4);
+ clusters->SetName("TRD Clusters");
+
+ clusters->SetTitle(Form("N=%d", clusters->Size()));
- char form[1000];
- sprintf(form,"TRD Clusters");
- clusters->SetName(form);
+ const TString viz_tag("REC Clusters TRD");
+
+ clusters->ApplyVizTag(viz_tag, "Clusters");
- char tip[1000];
- sprintf(tip,"N=%d", clusters->Size());
- clusters->SetTitle(tip);
gEve->AddElement(clusters, cont);
+
gEve->Redraw3D();
return clusters;