* full copyright notice. *
**************************************************************************/
-// Import tracks from kinematics-tree / particle-stack.
-// Preliminary/minimal solution.
+/// \ingroup evemacros
+/// \file kine_tracks.C
+/// \brief Import tracks from kinematics-tree / particle-stack.
+///
+/// Preliminary/minimal solution.
+///
+/// \author Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
-#include "TParticlePDG.h"
+#if !defined(__CINT__) || defined(__MAKECINT__)
+
+#include <TParticle.h>
+#include <TParticlePDG.h>
+
+#include <TEveManager.h>
+#include <TEveTrackPropagator.h>
+
+#include <EveBase/AliEveEventManager.h>
+#include <EveBase/AliEveMagField.h>
+#include <EveBase/AliEveTrack.h>
+#include <EveBase/AliEveKineTools.h>
+
+#include <AliRunLoader.h>
+#include <AliStack.h>
+#include <AliMagF.h>
+
+#endif
+
+// Use magnetic-field as retrieved from GRP.
+Bool_t g_kine_tracks_true_field = kTRUE;
+
+// Use Runge-Kutta track stepper.
+Bool_t g_kine_tracks_rk_stepper = kFALSE;
+
+//==============================================================================
+
+void kine_track_propagator_setup(TEveTrackPropagator* trkProp);
TEveTrackList*
kine_tracks(Double_t min_pt = 0, Double_t min_p = 0,
Bool_t pdg_col = kTRUE, Bool_t recurse = kTRUE,
- Bool_t use_track_refs = kTRUE)
+ Bool_t use_track_refs = kTRUE);
+
+void kine_daughters(AliEveTrack* parent, AliStack* stack,
+ Double_t min_pt, Double_t min_p,
+ Bool_t pdg_col, Bool_t recurse);
+
+void set_track_color(AliEveTrack* t, Bool_t pdg_col);
+Color_t get_pdg_color(Int_t pdg);
+
+TEveElement*
+kine_track(Int_t label,
+ Bool_t import_mother = kTRUE, Bool_t import_daughters = kTRUE,
+ Bool_t pdg_col = kTRUE, Bool_t recurse = kTRUE,
+ TEveElement* cont = 0);
+
+void kine_hide_neutrals(TEveElement* el=0, Int_t level=0);
+
+//==============================================================================
+
+void kine_track_propagator_setup(TEveTrackPropagator* trkProp)
+{
+ AliMagF* fld = AliEveEventManager::AssertMagField();
+
+ if (g_kine_tracks_true_field)
+ {
+ trkProp->SetMagFieldObj(new AliEveMagField(fld));
+ }
+ else
+ {
+ trkProp->SetMagField(-0.1*fld->SolenoidField());
+ }
+ if (g_kine_tracks_rk_stepper)
+ {
+ trkProp->SetStepper(TEveTrackPropagator::kRungeKutta);
+ }
+}
+
+//==============================================================================
+
+TEveTrackList*
+kine_tracks(Double_t min_pt, Double_t min_p,
+ Bool_t pdg_col, Bool_t recurse,
+ Bool_t use_track_refs)
{
AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
rl->LoadKinematics();
AliStack* stack = rl->Stack();
if (!stack)
{
- Error("kine_tracks.C", "can not get kinematics.");
+ Error("kine_tracks", "can not get kinematics.");
return 0;
}
cont->SetMainColor(3);
TEveTrackPropagator* trkProp = cont->GetPropagator();
- AliMagF* fld = AliEveEventManager::AssertMagField();
- // !!! Watch the '-', apparently different sign convention then for ESD.
- trkProp->SetMagField(fld ? -0.1*fld->SolenoidField() : 0);
+ kine_track_propagator_setup(trkProp);
gEve->AddElement(cont);
Int_t count = 0;
- Int_t N = stack->GetNtrack();
Int_t Np = stack->GetNprimary();
for (Int_t i = 0; i < Np; ++i)
{
}
}
-Color_t set_track_color(AliEveTrack* t, Bool_t pdg_col)
+//==============================================================================
+
+void set_track_color(AliEveTrack* t, Bool_t pdg_col)
{
if (pdg_col)
t->SetMainColor(get_pdg_color(t->GetPdg()));
// mesons and barions
else if (pdga < 100000) {
Int_t i = pdga;
- Int_t i0 = i%10; i /= 10;
+ // Int_t i0 = i%10; // Not used at the moment.
+ i /= 10;
Int_t i1 = i%10; i /= 10;
Int_t i2 = i%10; i /= 10;
Int_t i3 = i%10; i /= 10;
return col;
}
-/******************************************************************************/
+//==============================================================================
TEveElement*
kine_track(Int_t label,
- Bool_t import_mother = kTRUE, Bool_t import_daughters = kTRUE,
- Bool_t pdg_col = kTRUE, Bool_t recurse = kTRUE,
- TEveElement* cont = 0)
-
+ Bool_t import_mother, Bool_t import_daughters,
+ Bool_t pdg_col, Bool_t recurse,
+ TEveElement* cont)
{
// Create mother and daughters tracks with given label.
// mother -> particle with label
AliRunLoader* rl = AliEveEventManager::AssertRunLoader();
rl->LoadKinematics();
AliStack* stack = rl->Stack();
+ if (!stack)
+ {
+ Warning("kine_track", "can not get kinematics.");
+ return 0;
+ }
if (label >= stack->GetNtrack())
{
Warning("kine_track", "label out of range.");
if (import_mother || (import_daughters && p->GetNDaughters()))
{
- AliEveTrack* toptrack = 0;
- TEveTrackList* tracklist = 0;
TEveTrackPropagator* rs = 0;
if (cont == 0)
{
TEveTrackList* tlist = new TEveTrackList
- (Form("Kinematics of %d", label, p->GetNDaughters()));
+ (Form("Kinematics of %d %d", label, p->GetNDaughters()));
cont = tlist;
TEveTrackPropagator* trkProp = tlist->GetPropagator();
- AliMagF* fld = AliEveEventManager::AssertMagField();
- trkProp->SetMagField(fld ? -0.1*fld->SolenoidField() : 0);
+ kine_track_propagator_setup(trkProp);
char tooltip[1000];
sprintf(tooltip,"Ndaughters=%d", p->GetNDaughters());
tlist->SetTitle(tooltip);
- trkProp->fMaxOrbs = 2;
+ trkProp->SetMaxOrbs(2);
trkProp->SetEditPathMarks(kTRUE);
gEve->AddElement(cont);
return cont;
}
+//==============================================================================
+
+void kine_hide_neutrals(TEveElement* el, Int_t level)
+{
+ if (el == 0)
+ {
+ el = gEve->GetCurrentEvent()->FindChild("Kine Tracks");
+ if (!el)
+ return;
+ }
+
+ TEveTrack* t = dynamic_cast<TEveTrack*>(el);
+ if (t && t->GetCharge() == 0)
+ t->SetRnrSelf(kFALSE);
+
+ for (TEveElement::List_i i = el->BeginChildren(); i != el->EndChildren(); ++i)
+ {
+ kine_hide_neutrals(*i, level + 1);
+ }
+
+ if (level == 0)
+ gEve->Redraw3D();
+}
git://git.uio.no / u / mrichter / AliRoot.git / blobdiff