// $Id$
+// Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
-Reve::Track* esd_make_track(Reve::TrackRnrStyle* rnrStyle,
- Int_t index,
- AliESDtrack* at,
- AliExternalTrackParam* tp=0)
+/**************************************************************************
+ * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
+ * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for *
+ * full copyright notice. *
+ **************************************************************************/
+
+#if !defined(__CINT__) || defined(__MAKECINT__)
+
+#include "TString.h"
+#include "TMath.h"
+#include "TGListTree.h"
+#include "TEveVSDStructs.h"
+#include "TEveManager.h"
+#include "TEveTrackPropagator.h"
+
+#include "AliESDEvent.h"
+#include "AliESDtrackCuts.h"
+#include "AliESDtrack.h"
+#include "AliExternalTrackParam.h"
+
+#include "EVE/EveBase/AliEveTrack.h"
+#include "EVE/EveBase/AliEveMagField.h"
+#include "EVE/EveBase/AliEveEventManager.h"
+
+#endif
+
+
+// Use inner-tpc track params when its refit failed.
+Bool_t g_esd_tracks_use_ip_on_failed_its_refit = kFALSE;
+
+// Use magnetic-field as retrieved from GRP.
+Bool_t g_esd_tracks_true_field = kTRUE;
+
+// Use Runge-Kutta track stepper.
+Bool_t g_esd_tracks_rk_stepper = kFALSE;
+
+
+//==============================================================================
+
+void esd_track_propagator_setup(TEveTrackPropagator* trkProp,
+ Float_t magF, Float_t maxR)
+{
+ if (g_esd_tracks_true_field)
+ {
+ trkProp->SetMagFieldObj(new AliEveMagField);
+ }
+ else
+ {
+ trkProp->SetMagField(magF);
+ }
+ if (g_esd_tracks_rk_stepper)
+ {
+ trkProp->SetStepper(TEveTrackPropagator::kRungeKutta);
+ }
+ trkProp->SetMaxR(maxR);
+}
+
+//==============================================================================
+
+TString esd_track_title(AliESDtrack* t)
+{
+ TString s;
+
+ Int_t label = t->GetLabel(), index = t->GetID();
+ TString idx(index == kMinInt ? "<undef>" : Form("%d", index));
+ TString lbl(label == kMinInt ? "<undef>" : Form("%d", label));
+
+ Double_t p[3], v[3];
+ t->GetXYZ(v);
+ t->GetPxPyPz(p);
+ Double_t pt = t->Pt();
+ Double_t ptsig = TMath::Sqrt(t->GetSigma1Pt2());
+ Double_t ptsq = pt*pt;
+ Double_t ptm = pt / (1.0 + pt*ptsig);
+ Double_t ptM = pt / (1.0 - pt*ptsig);
+
+ s = Form("Index=%s, Label=%s\nChg=%d, Pdg=%d\n"
+ "pT = %.3f + %.3f - %.3f [%.3f]\n"
+ "P = (%.3f, %.3f, %.3f)\n"
+ "V = (%.3f, %.3f, %.3f)\n",
+ idx.Data(), lbl.Data(), t->Charge(), 0,
+ pt, ptM - pt, pt - ptm, ptsig*ptsq,
+ p[0], p[1], p[2],
+ v[0], v[1], v[2]);
+
+ Int_t o;
+ s += "Det (in,out,refit,pid):\n";
+ o = AliESDtrack::kITSin;
+ s += Form("ITS (%d,%d,%d,%d) ", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
+ o = AliESDtrack::kTPCin;
+ s += Form("TPC(%d,%d,%d,%d)\n", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
+ o = AliESDtrack::kTRDin;
+ s += Form("TRD(%d,%d,%d,%d) ", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
+ o = AliESDtrack::kTOFin;
+ s += Form("TOF(%d,%d,%d,%d)\n", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
+ o = AliESDtrack::kHMPIDout;
+ s += Form("HMPID(out=%d,pid=%d)\n", t->IsOn(o), t->IsOn(o<<1));
+ s += Form("ESD pid=%d", t->IsOn(AliESDtrack::kESDpid));
+
+ if (t->IsOn(AliESDtrack::kESDpid))
+ {
+ Double_t pid[5];
+ t->GetESDpid(pid);
+ s += Form("\n[%.2f %.2f %.2f %.2f %.2f]", pid[0], pid[1], pid[2], pid[3], pid[4]);
+ }
+
+ return s;
+}
+
+//==============================================================================
+
+void esd_track_add_param(AliEveTrack* track, const AliExternalTrackParam* tp)
{
- // Helper function
- Double_t pbuf[3], vbuf[3];
- Reve::RecTrack rt;
+ // Add additional track parameters as a path-mark to track.
- if(tp == 0) tp = at;
+ if (tp == 0)
+ return;
- rt.label = at->GetLabel();
- rt.index = index;
- rt.status = (Int_t) at->GetStatus();
- rt.sign = tp->GetSign();
+ Double_t pbuf[3], vbuf[3];
tp->GetXYZ(vbuf);
- rt.V.Set(vbuf);
tp->GetPxPyPz(pbuf);
- rt.P.Set(pbuf);
- Double_t ep = at->GetP(), mc = at->GetMass();
- rt.beta = ep/TMath::Sqrt(ep*ep + mc*mc);
-
- Reve::Track* track = new Reve::Track(&rt, rnrStyle);
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH track->SetName(Form("ESDTrack %d", rt.label));
- //PH track->SetTitle(Form("pT=%.3f, pZ=%.3f; V=(%.3f, %.3f, %.3f)",
- //PH rt.sign*TMath::Hypot(rt.P.x, rt.P.y), rt.P.z,
- //PH rt.V.x, rt.V.y, rt.V.z));
- char form[1000];
- sprintf(form,"Track %d", rt.index);
- track->SetName(form);
- sprintf(form,"Index=%d, Label=%d\npT=%.3f, pZ=%.3f\nV=(%.3f, %.3f, %.3f)",
- index, rt.label,
- rt.sign*TMath::Hypot(rt.P.x, rt.P.y), rt.P.z,
- rt.V.x, rt.V.y, rt.V.z);
- track->SetTitle(form);
- return track;
+
+ TEvePathMark pm(TEvePathMark::kReference);
+ pm.fV.Set(vbuf);
+ pm.fP.Set(pbuf);
+ track->AddPathMark(pm);
}
-Bool_t gkFixFailedITSExtr = kFALSE;
+//==============================================================================
-Reve::TrackList* esd_tracks(Double_t min_pt=0.1, Double_t max_pt=100)
+AliEveTrack* esd_make_track(AliESDtrack *at, TEveTrackList* cont)
{
- AliESDEvent* esd = Alieve::Event::AssertESD();
+ // Make a standard track representation and put it into given container.
- Double_t minptsq = min_pt*min_pt;
- Double_t maxptsq = max_pt*max_pt;
- Double_t ptsq;
+ // Choose which parameters to use a track's starting point.
+ // If gkFixFailedITSExtr is TRUE (FALSE by default) and
+ // if ITS refit failed, take track parameters at inner TPC radius.
- Reve::TrackList* cont = new Reve::TrackList("ESD Tracks");
- cont->SetMainColor(Color_t(6));
- Reve::TrackRnrStyle* rnrStyle = cont->GetRnrStyle();
- rnrStyle->SetMagField( esd->GetMagneticField() );
+ const AliExternalTrackParam* tp = at;
+
+ Bool_t innerTaken = kFALSE;
+ if ( ! at->IsOn(AliESDtrack::kITSrefit) && g_esd_tracks_use_ip_on_failed_its_refit)
+ {
+ tp = at->GetInnerParam();
+ innerTaken = kTRUE;
+ }
- gReve->AddRenderElement(cont);
+ AliEveTrack* track = new AliEveTrack(at, cont->GetPropagator());
+ track->SetAttLineAttMarker(cont);
+ track->SetName(Form("AliEveTrack %d", at->GetID()));
+ track->SetElementTitle(esd_track_title(at));
+ track->SetSourceObject(at);
- Int_t count = 0;
- Double_t pbuf[3];
- for (Int_t n=0; n<esd->GetNumberOfTracks(); n++)
+ // Add inner/outer track parameters as path-marks.
+ if (at->IsOn(AliESDtrack::kTPCrefit))
{
- AliESDtrack* at = esd->GetTrack(n);
+ if ( ! innerTaken)
+ {
+ esd_track_add_param(track, at->GetInnerParam());
+ }
+ esd_track_add_param(track, at->GetOuterParam());
+ }
- // Here would be sweet to have TObjectFormula.
- at->GetPxPyPz(pbuf);
- ptsq = pbuf[0]*pbuf[0] + pbuf[1]*pbuf[1];
- if(ptsq < minptsq || ptsq > maxptsq)
- continue;
+ return track;
+}
- ++count;
- // If gkFixFailedITSExtr is TRUE (FALSE by default) and
- // if ITS refit failed, take track parameters at inner TPC radius.
- AliExternalTrackParam* tp = at;
- if (gkFixFailedITSExtr && !at->IsOn(AliESDtrack::kITSrefit)) {
- tp = at->GetInnerParam();
- }
+//==============================================================================
+// esd_tracks()
+//==============================================================================
- Reve::Track* track = esd_make_track(rnrStyle, n, at, tp);
- track->SetAttLineAttMarker(cont);
- gReve->AddRenderElement(track, cont);
- }
+TEveTrackList* esd_tracks()
+{
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
+
+ TEveTrackList* cont = new TEveTrackList("ESD Tracks");
+ cont->SetMainColor(6);
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH const Text_t* tooltip = Form("pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
- char tooltip[1000];
- sprintf(tooltip,"pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
- cont->SetTitle(tooltip); // Not broadcasted automatically ...
- cont->UpdateItems();
+ esd_track_propagator_setup(cont->GetPropagator(),
+ 0.1*esd->GetMagneticField(), 520);
+ gEve->AddElement(cont);
+
+ Int_t count = 0;
+ for (Int_t n = 0; n < esd->GetNumberOfTracks(); ++n)
+ {
+ ++count;
+ AliEveTrack* track = esd_make_track(esd->GetTrack(n), cont);
+
+ cont->AddElement(track);
+ }
+ cont->SetTitle(Form("N=%d", count));
cont->MakeTracks();
- gReve->Redraw3D();
+ gEve->Redraw3D();
+
+ return cont;
+}
+
+TEveElementList* esd_tracks_MI()
+{
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
+
+ TEveElementList* cont = new TEveElementList("ESD Tracks MI");
+ gEve->AddElement(cont);
+
+ Int_t count = 0;
+ for (Int_t n = 0; n < esd->GetNumberOfTracks(); ++n)
+ {
+ ++count;
+ AliESDtrack* at = esd->GetTrack(n);
+ TEveLine* l = new TEveLine;
+ l->SetLineColor(5);
+ at->FillPolymarker(l, esd->GetMagneticField(), 0, 250, 5);
+ l->SetName(Form("ESDTrackMI %d", at->GetID()));
+ l->SetElementTitle(esd_track_title(at));
+ l->SetSourceObject(at);
+
+ cont->AddElement(l);
+ }
+ cont->SetTitle(Form("N=%d", count));
+
+ gEve->Redraw3D();
return cont;
}
-/**************************************************************************/
+
+//==============================================================================
// esd_tracks_from_array()
-/**************************************************************************/
+//==============================================================================
-Reve::TrackList* esd_tracks_from_array(TCollection* col, AliESDEvent* esd=0)
+TEveTrackList* esd_tracks_from_array(TCollection* col, AliESDEvent* esd=0)
{
// Retrieves AliESDTrack's from collection.
// See example usage with AliAnalysisTrackCuts in the next function.
- if(esd == 0) esd = Alieve::Event::AssertESD();
+ if (esd == 0) esd = AliEveEventManager::AssertESD();
- Reve::TrackList* cont = new Reve::TrackList("ESD Tracks");
- cont->SetMainColor(Color_t(6));
- Reve::TrackRnrStyle* rnrStyle = cont->GetRnrStyle();
- rnrStyle->SetMagField( esd->GetMagneticField() );
+ TEveTrackList* cont = new TEveTrackList("ESD Tracks");
+ cont->SetMainColor(6);
- gReve->AddRenderElement(cont);
+ esd_track_propagator_setup(cont->GetPropagator(),
+ 0.1*esd->GetMagneticField(), 520);
+
+ gEve->AddElement(cont);
Int_t count = 0;
TIter next(col);
TObject *obj;
- while((obj = next()) != 0)
+ while ((obj = next()) != 0)
{
- if(obj->IsA()->InheritsFrom("AliESDtrack") == kFALSE) {
- Warning("Object '%s', '%s' is not an AliESDtrack.",
+ if (obj->IsA()->InheritsFrom("AliESDtrack") == kFALSE)
+ {
+ Warning("esd_tracks_from_array", "Object '%s', '%s' is not an AliESDtrack.",
obj->GetName(), obj->GetTitle());
continue;
}
++count;
- AliESDtrack* at = (AliESDtrack*) obj;
+ AliESDtrack* at = reinterpret_cast<AliESDtrack*>(obj);
- Reve::Track* track = esd_make_track(rnrStyle, count, at);
- track->SetAttLineAttMarker(cont);
- gReve->AddRenderElement(track, cont);
+ AliEveTrack* track = esd_make_track(at, cont);
+ cont->AddElement(track);
}
-
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH const Text_t* tooltip = Form("N=%d", count);
- const tooltip[1000];
- sprintf(tooltip,"N=%d", count);
- cont->SetTitle(tooltip); // Not broadcasted automatically ...
- cont->UpdateItems();
-
+ cont->SetTitle(Form("N=%d", count));
cont->MakeTracks();
- gReve->Redraw3D();
+ gEve->Redraw3D();
return cont;
}
void esd_tracks_alianalcuts_demo()
{
- AliESDEvent* esd = Alieve::Event::AssertESD();
- gSystem->Load("libANALYSIS");
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
- AliAnalysisTrackCuts atc;
+ AliESDtrackCuts atc;
atc.SetPtRange(0.1, 5);
atc.SetRapRange(-1, 1);
- esd_tracks_from_array(atc.GetAcceptedParticles(esd), esd);
+ esd_tracks_from_array(atc.GetAcceptedTracks(esd), esd);
}
-/**************************************************************************/
-// esd_tracks_vertex_cut
-/**************************************************************************/
+
+//==============================================================================
+// esd_tracks_by_category
+//==============================================================================
Float_t get_sigma_to_vertex(AliESDtrack* esdTrack)
{
Float_t bRes[2];
Float_t bCov[3];
esdTrack->GetImpactParameters(b,bCov);
- if (bCov[0]<=0 || bCov[2]<=0) {
+ if (bCov[0] <= 0 || bCov[2] <= 0)
+ {
printf("Estimated b resolution lower or equal zero!\n");
- bCov[0]=0; bCov[2]=0;
+ bCov[0] = bCov[2] = 0;
}
bRes[0] = TMath::Sqrt(bCov[0]);
bRes[1] = TMath::Sqrt(bCov[2]);
// It means that for a 2-dim gauss: n_sigma(d) = Sqrt(2)*ErfInv(1 - Exp((-x**2)/2)
// Can this be expressed in a different way?
- if (bRes[0] == 0 || bRes[1] ==0)
+ if (bRes[0] == 0 || bRes[1] == 0)
return -1;
Float_t d = TMath::Sqrt(TMath::Power(b[0]/bRes[0],2) + TMath::Power(b[1]/bRes[1],2));
return d;
}
-Reve::RenderElementList* esd_tracks_vertex_cut()
+TEveElementList* esd_tracks_by_category()
{
- // Import ESD tracks, separate them into five containers according to
- // primary-vertex cut and ITS refit status.
+ // Import ESD tracks, separate them into several containers
+ // according to primary-vertex cut and ITS&TPC refit status.
+
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
- AliESDEvent* esd = Alieve::Event::AssertESD();
+ TEveElementList* cont = new TEveElementList("ESD Tracks by category");
+ gEve->AddElement(cont);
- Reve::RenderElementList* cont = new Reve::RenderElementList("ESD Tracks", 0, kTRUE);
- gReve->AddRenderElement(cont);
- Reve::TrackList *tl[5];
- Int_t tc[5];
- Int_t count = 0;
+ const Int_t nCont = 7;
+ const Float_t maxR = 520;
+ const Float_t magF = 0.1*esd->GetMagneticField();
- tl[0] = new Reve::TrackList("Sigma < 3");
+ TEveTrackList *tl[nCont];
+ Int_t tc[nCont];
+ Int_t count = 0;
+
+ tl[0] = new TEveTrackList("Sigma < 3");
tc[0] = 0;
- tl[0]->GetRnrStyle()->SetMagField( esd->GetMagneticField() );
- tl[0]->SetMainColor(Color_t(3));
- gReve->AddRenderElement(tl[0], cont);
+ esd_track_propagator_setup(tl[0]->GetPropagator(), magF, maxR);
+ tl[0]->SetMainColor(3);
+ cont->AddElement(tl[0]);
- tl[1] = new Reve::TrackList("3 < Sigma < 5");
+ tl[1] = new TEveTrackList("3 < Sigma < 5");
tc[1] = 0;
- tl[1]->GetRnrStyle()->SetMagField( esd->GetMagneticField() );
- tl[1]->SetMainColor(Color_t(7));
- gReve->AddRenderElement(tl[1], cont);
+ esd_track_propagator_setup(tl[1]->GetPropagator(), magF, maxR);
+ tl[1]->SetMainColor(7);
+ cont->AddElement(tl[1]);
- tl[2] = new Reve::TrackList("5 < Sigma");
+ tl[2] = new TEveTrackList("5 < Sigma");
tc[2] = 0;
- tl[2]->GetRnrStyle()->SetMagField( esd->GetMagneticField() );
- tl[2]->SetMainColor(Color_t(46));
- gReve->AddRenderElement(tl[2], cont);
+ esd_track_propagator_setup(tl[2]->GetPropagator(), magF, maxR);
+ tl[2]->SetMainColor(46);
+ cont->AddElement(tl[2]);
- tl[3] = new Reve::TrackList("no ITS refit; Sigma < 5");
+ tl[3] = new TEveTrackList("no ITS refit; Sigma < 5");
tc[3] = 0;
- tl[3]->GetRnrStyle()->SetMagField( esd->GetMagneticField() );
- tl[3]->SetMainColor(Color_t(41));
- gReve->AddRenderElement(tl[3], cont);
+ esd_track_propagator_setup(tl[3]->GetPropagator(), magF, maxR);
+ tl[3]->SetMainColor(41);
+ cont->AddElement(tl[3]);
- tl[4] = new Reve::TrackList("no ITS refit; Sigma > 5");
+ tl[4] = new TEveTrackList("no ITS refit; Sigma > 5");
tc[4] = 0;
- tl[4]->GetRnrStyle()->SetMagField( esd->GetMagneticField() );
- tl[4]->SetMainColor(Color_t(48));
- gReve->AddRenderElement(tl[4], cont);
-
- for (Int_t n=0; n<esd->GetNumberOfTracks(); n++)
+ esd_track_propagator_setup(tl[4]->GetPropagator(), magF, maxR);
+ tl[4]->SetMainColor(48);
+ cont->AddElement(tl[4]);
+
+ tl[5] = new TEveTrackList("no TPC refit");
+ tc[5] = 0;
+ esd_track_propagator_setup(tl[5]->GetPropagator(), magF, maxR);
+ tl[5]->SetMainColor(kRed);
+ cont->AddElement(tl[5]);
+
+ tl[6] = new TEveTrackList("ITS stand-alone");
+ tc[6] = 0;
+ esd_track_propagator_setup(tl[6]->GetPropagator(), magF, maxR);
+ tl[6]->SetMainColor(kMagenta - 9);
+ cont->AddElement(tl[6]);
+
+ for (Int_t n = 0; n < esd->GetNumberOfTracks(); ++n)
{
AliESDtrack* at = esd->GetTrack(n);
else if (s <= 5) ti = 1;
else ti = 2;
- AliExternalTrackParam* tp = at;
- // If ITS refit failed, optionally take track parameters at inner
- // TPC radius and put track in a special container.
- // This ignores state of gkFixFailedITSExtr (used in esd_tracks()).
- // Use BOTH functions to compare results.
- if (!at->IsOn(AliESDtrack::kITSrefit)) {
- // tp = at->GetInnerParam();
+ if (at->IsOn(AliESDtrack::kITSin) && ! at->IsOn(AliESDtrack::kTPCin))
+ {
+ ti = 6;
+ }
+ else if (at->IsOn(AliESDtrack::kTPCin) && ! at->IsOn(AliESDtrack::kTPCrefit))
+ {
+ ti = 5;
+ }
+ else if ( ! at->IsOn(AliESDtrack::kITSrefit))
+ {
ti = (ti == 2) ? 4 : 3;
}
- Reve::TrackList* tlist = tl[ti];
+ TEveTrackList* tlist = tl[ti];
+ ++tc[ti];
+ ++count;
+
+ AliEveTrack* track = esd_make_track(at, tlist);
+ track->SetName(Form("ESD Track idx=%d, sigma=%5.3f", at->GetID(), s));
+ tlist->AddElement(track);
+ }
+
+ for (Int_t ti = 0; ti < nCont; ++ti)
+ {
+ TEveTrackList* tlist = tl[ti];
+ tlist->SetName(Form("%s [%d]", tlist->GetName(), tlist->NumChildren()));
+ tlist->SetTitle(Form("N tracks=%d", tc[ti]));
+
+ tlist->MakeTracks();
+
+ Bool_t good_cont = ti <= 1;
+ if (AliEveTrackCounter::IsActive())
+ {
+ AliEveTrackCounter::fgInstance->RegisterTracks(tlist, good_cont);
+ }
+ else
+ {
+ if ( ! good_cont)
+ tlist->SetLineStyle(6);
+ }
+ }
+ cont->SetTitle(Form("N all tracks = %d", count));
+ // ??? The following does not always work:
+ cont->FindListTreeItem(gEve->GetListTree())->SetOpen(kTRUE);
+
+ gEve->Redraw3D();
+
+ return cont;
+}
+
+
+//==============================================================================
+// esd_tracks_by_anal_cuts
+//==============================================================================
+
+AliESDtrackCuts* g_esd_tracks_anal_cuts = 0;
+
+TEveElementList* esd_tracks_by_anal_cuts()
+{
+ AliESDEvent* esd = AliEveEventManager::AssertESD();
+
+ if (g_esd_tracks_anal_cuts == 0)
+ {
+ gSystem->Load("libPWG0base");
+ gROOT->ProcessLine(".L $ALICE_ROOT/PWG0/CreateStandardCuts.C");
+ AliPWG0Helper::AnalysisMode mode = AliPWG0Helper::kTPCITS;
+ if (TMath::Abs(esd->GetMagneticField()) > 0.01)
+ mode |= AliPWG0Helper::kFieldOn;
+ g_esd_tracks_anal_cuts = CreateTrackCuts(mode, kFALSE);
+ }
+
+ TEveElementList* cont = new TEveElementList("ESD Tracks by Analysis Cuts");
+ gEve->AddElement(cont);
+
+ const Int_t nCont = 2;
+ const Float_t maxR = 520;
+ const Float_t magF = 0.1*esd->GetMagneticField();
+
+ TEveTrackList *tl[nCont];
+ Int_t tc[nCont];
+ Int_t count = 0;
+
+ tl[0] = new TEveTrackList("Passed");
+ tc[0] = 0;
+ esd_track_propagator_setup(tl[0]->GetPropagator(), magF, maxR);
+ tl[0]->SetMainColor(3);
+ cont->AddElement(tl[0]);
+
+ tl[1] = new TEveTrackList("Rejected");
+ tc[1] = 0;
+ esd_track_propagator_setup(tl[1]->GetPropagator(), magF, maxR);
+ tl[1]->SetMainColor(kRed);
+ cont->AddElement(tl[1]);
+
+ for (Int_t n = 0; n < esd->GetNumberOfTracks(); ++n)
+ {
+ AliESDtrack* at = esd->GetTrack(n);
+
+ Float_t s = get_sigma_to_vertex(at);
+ Int_t ti = (g_esd_tracks_anal_cuts->AcceptTrack(at)) ? 0 : 1;
+
+ TEveTrackList* tlist = tl[ti];
++tc[ti];
++count;
- Reve::Track* track = esd_make_track(tlist->GetRnrStyle(), n, at, tp);
- track->SetAttLineAttMarker(tlist);
-
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH track->SetName(Form("track %d, sigma=%5.3f", at->GetLabel(), s));
- char form[1000];
- sprintf(form,"Track idx=%d, sigma=%5.3f", at->GetIndex(), s);
- track->SetName(form);
- gReve->AddRenderElement(track, tlist);
+ AliEveTrack* track = esd_make_track(at, tlist);
+ track->SetName(Form("ESD Track idx=%d, sigma=%5.3f", at->GetID(), s));
+ tlist->AddElement(track);
}
- for (Int_t ti=0; ti<5; ++ti) {
- Reve::TrackList* tlist = tl[ti];
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH const Text_t* tooltip = Form("N tracks=%d", tc[ti]);
- //MT Modified somewhat.
- char buff[1000];
- sprintf(buff, "%s [%d]", tlist->GetName(), tlist->GetNChildren());
- tlist->SetName(buff);
- sprintf(buff, "N tracks=%d", tc[ti]);
- tlist->SetTitle(buff); // Not broadcasted automatically ...
- tlist->UpdateItems();
+ for (Int_t ti = 0; ti < nCont; ++ti)
+ {
+ TEveTrackList* tlist = tl[ti];
+ tlist->SetName(Form("%s [%d]", tlist->GetName(), tlist->NumChildren()));
+ tlist->SetTitle(Form("N tracks=%d", tc[ti]));
tlist->MakeTracks();
+
+ Bool_t good_cont = ti < 1;
+ if (AliEveTrackCounter::IsActive())
+ {
+ AliEveTrackCounter::fgInstance->RegisterTracks(tlist, good_cont);
+ }
+ else
+ {
+ if ( ! good_cont)
+ tlist->SetLineStyle(6);
+ }
}
- //PH The line below is replaced waiting for a fix in Root
- //PH which permits to use variable siza arguments in CINT
- //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
- //PH cont->SetTitle(Form("N all tracks = %d", count));
- char form[1000];
- sprintf(form,"N all tracks = %d", count);
- cont->SetTitle(form);
- cont->UpdateItems();
- gReve->Redraw3D();
+ cont->SetTitle(Form("N all tracks = %d", count));
+ // ??? The following does not always work:
+ cont->FindListTreeItem(gEve->GetListTree())->SetOpen(kTRUE);
+
+ gEve->Redraw3D();
return cont;
}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff