In muon-related macros in EVE:
[u/mrichter/AliRoot.git] / EVE / alice-macros / kine_tracks.C
index 902cc68..e9e3618 100644 (file)
 // Import tracks from kinematics-tree / particle-stack.
 // Preliminary/minimal solution.
 
-#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) {
+  if (!stack)
+  {
     Error("kine_tracks.C", "can not get kinematics.");
     return 0;
   }
@@ -29,22 +99,12 @@ kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
 
   TEveTrackList* cont = new TEveTrackList("Kine Tracks");
   cont->SetMainColor(3);
-  TEveTrackPropagator* rnrStyle = cont->GetPropagator();
-  Float_t mag_field = 0;
-  if (gAlice && gAlice->Field())
-  {
-    // !!! Watch the '-', apparently different sign convention then for ESD.
-    mag_field = -0.1*gAlice->Field()->SolenoidField();
-  }
-  else
-  {
-    Warning("kine_tracks.C", "Could not determine magnetic field from gAlice - using zero.");
-  }
-  rnrStyle->SetMagField(mag_field);
+  TEveTrackPropagator* trkProp = cont->GetPropagator();
+
+  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)
   {
@@ -54,7 +114,7 @@ kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
       if (p->Pt() < min_pt && p->P() < min_p) continue;
 
       ++count;
-      TEveTrack* track = new TEveTrack(p, i, rnrStyle);
+      AliEveTrack* track = new AliEveTrack(p, i, trkProp);
 
       //PH The line below is replaced waiting for a fix in Root
       //PH which permits to use variable siza arguments in CINT
@@ -86,7 +146,7 @@ kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
   if (use_track_refs && rl->LoadTrackRefs() == 0)
   {
     kt.SetTrackReferences(cont, rl->TreeTR(), recurse);
-    rnrStyle->SetEditPathMarks(kTRUE);
+    trkProp->SetEditPathMarks(kTRUE);
   }
   kt.SortPathMarks(cont, recurse);
 
@@ -102,7 +162,7 @@ kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
   return cont;
 }
 
-void kine_daughters(TEveTrack* parent,  AliStack* stack,
+void kine_daughters(AliEveTrack* parent,  AliStack* stack,
                    Double_t     min_pt,  Double_t  min_p,
                    Bool_t       pdg_col, Bool_t    recurse)
 {
@@ -115,7 +175,7 @@ void kine_daughters(TEveTrack* parent,  AliStack* stack,
       TParticle* dp = stack->Particle(d);
       if (dp->Pt() < min_pt && dp->P() < min_p) continue;
 
-      TEveTrack* dtrack = new TEveTrack(dp, d, rs);
+      AliEveTrack* dtrack = new AliEveTrack(dp, d, rs);
       char form[1000];
       sprintf(form,"%s [%d]", dp->GetName(), d);
       dtrack->SetName(form);
@@ -130,7 +190,9 @@ void kine_daughters(TEveTrack* parent,  AliStack* stack,
   }
 }
 
-Color_t set_track_color(TEveTrack* 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()));
@@ -166,7 +228,8 @@ Color_t get_pdg_color(Int_t pdg)
   // 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;
@@ -186,14 +249,13 @@ Color_t get_pdg_color(Int_t pdg)
   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
@@ -217,8 +279,6 @@ kine_track(Int_t  label,
 
   if (import_mother || (import_daughters && p->GetNDaughters()))
   {
-    TEveTrack* toptrack = 0;
-    TEveTrackList* tracklist = 0;
     TEveTrackPropagator* rs = 0;
 
     if (cont == 0)
@@ -227,22 +287,23 @@ kine_track(Int_t  label,
        (Form("Kinematics of %d", label, p->GetNDaughters()));
       cont = tlist;
 
-      TEveTrackPropagator* rnrStyle = tlist->GetPropagator();
-      // !!! Watch the '-', apparently different sign convention then for ESD.
-      rnrStyle->SetMagField( -0.1*gAlice->Field()->SolenoidField() );
+      TEveTrackPropagator* trkProp = tlist->GetPropagator();
+
+      kine_track_propagator_setup(trkProp);
+
       char tooltip[1000];
       sprintf(tooltip,"Ndaughters=%d", p->GetNDaughters());
       tlist->SetTitle(tooltip);
-      rnrStyle->fMaxOrbs = 2;
-      rnrStyle->SetEditPathMarks(kTRUE);
+      trkProp->SetMaxOrbs(2);
+      trkProp->SetEditPathMarks(kTRUE);
 
       gEve->AddElement(cont);
       rs = tlist->GetPropagator();
     }
     else
     {
-      // check if container is TEveTrackList or TEveTrack (has rnr-style)
-      TEveTrack* t = dynamic_cast<TEveTrack*>(cont);
+      // check if container is TEveTrackList or AliEveTrack (has rnr-style)
+      AliEveTrack* t = dynamic_cast<AliEveTrack*>(cont);
       if (t) {
        rs = t->GetPropagator();
       } else {
@@ -256,7 +317,7 @@ kine_track(Int_t  label,
 
     if (import_mother)
     {
-      TEveTrack* track = new TEveTrack(p, label, rs);
+      AliEveTrack* track = new AliEveTrack(p, label, rs);
       char form[1000];
       sprintf(form,"%s [%d]", p->GetName(), label);
       track->SetName(form);
@@ -273,7 +334,7 @@ kine_track(Int_t  label,
       for (int d=p->GetFirstDaughter(); d>0 && d<=p->GetLastDaughter(); ++d)
       {
        TParticle* dp = stack->Particle(d);
-       TEveTrack* track = new TEveTrack(dp, d, rs);
+       AliEveTrack* track = new AliEveTrack(dp, d, rs);
        char form[1000];
        sprintf(form,"%s [%d]", dp->GetName(), d);
        track->SetName(form);
@@ -293,3 +354,26 @@ kine_track(Int_t  label,
   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();
+}