EVE/alice-macros/kine_tracks.C

   1 // $Id$
   2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
   3
   4 /**************************************************************************
   5  * Copyright(c) 1998-2008, ALICE Experiment at CERN, all rights reserved. *
   6  * See http://aliceinfo.cern.ch/Offline/AliRoot/License.html for          *
   7  * full copyright notice.                                                 *
   8  **************************************************************************/
   9
  10 // Import tracks from kinematics-tree / particle-stack.
  11 // Preliminary/minimal solution.
  12
  13 #include "TParticlePDG.h"
  14
  15 TEveTrackList*
  16 kine_tracks(Double_t min_pt  = 0,     Double_t min_p   = 0,
  17             Bool_t   pdg_col = kTRUE, Bool_t   recurse = kTRUE,
  18             Bool_t   use_track_refs = kTRUE)
  19 {
  20   AliRunLoader* rl =  AliEveEventManager::AssertRunLoader();
  21   rl->LoadKinematics();
  22   AliStack* stack = rl->Stack();
  23   if (!stack) {
  24     Error("kine_tracks.C", "can not get kinematics.");
  25     return 0;
  26   }
  27
  28   gEve->DisableRedraw();
  29
  30   TEveTrackList* cont = new TEveTrackList("Kine Tracks");
  31   cont->SetMainColor(3);
  32   TEveTrackPropagator* rnrStyle = cont->GetPropagator();
  33   // !!! Watch the '-', apparently different sign convention then for ESD.
  34   rnrStyle->SetMagField( -0.1*gAlice->Field()->SolenoidField() );
  35
  36   gEve->AddElement(cont);
  37   Int_t count = 0;
  38   Int_t N = stack->GetNtrack();
  39   for (Int_t i=0; i<N; ++i)
  40   {
  41     if (stack->IsPhysicalPrimary(i))
  42     {
  43       TParticle* p = stack->Particle(i);
  44       if (p->Pt() < min_pt && p->P() < min_p) continue;
  45
  46       ++count;
  47       TEveTrack* track = new TEveTrack(p, i, rnrStyle);
  48
  49       //PH The line below is replaced waiting for a fix in Root
  50       //PH which permits to use variable siza arguments in CINT
  51       //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
  52       //PH    track->SetName(Form("%s [%d]", p->GetName(), i));
  53       char form[1000];
  54       sprintf(form,"%s [%d]", p->GetName(), i);
  55       track->SetName(form);
  56       track->SetStdTitle();
  57       Int_t ml = p->GetMother(0);
  58       if (ml != -1)
  59       {
  60         track->SetTitle(Form("%s\nMother label=%d\nMother Pdg=%d",
  61                              track->GetElementTitle(),
  62                              ml, stack->Particle(ml)->GetPdgCode()));
  63       }
  64       set_track_color(track, pdg_col);
  65
  66       gEve->AddElement(track, cont);
  67
  68       if (recurse)
  69         kine_daughters(track, stack, min_pt, min_p, pdg_col, recurse);
  70     }
  71   }
  72
  73   // set path marks
  74   AliEveKineTools kt;
  75   kt.SetDaughterPathMarks(cont, stack, recurse);
  76   if (use_track_refs && rl->LoadTrackRefs() == 0)
  77   {
  78     kt.SetTrackReferences(cont, rl->TreeTR(), recurse);
  79     rnrStyle->SetEditPathMarks(kTRUE);
  80   }
  81   kt.SortPathMarks(cont, recurse);
  82
  83   //PH  const Text_t* tooltip = Form("min pT=%.2lf, min P=%.2lf), N=%d", min_pt, min_p, count);
  84   char tooltip[1000];
  85   sprintf(tooltip,"min pT=%.2lf, min P=%.2lf), N=%d", min_pt, min_p, count);
  86   cont->SetTitle(tooltip); // Not broadcasted automatically ...
  87
  88   cont->MakeTracks(recurse);
  89   gEve->EnableRedraw();
  90   gEve->Redraw3D();
  91
  92   return cont;
  93 }
  94
  95 void kine_daughters(TEveTrack* parent,  AliStack* stack,
  96                     Double_t     min_pt,  Double_t  min_p,
  97                     Bool_t       pdg_col, Bool_t    recurse)
  98 {
  99   TParticle *p = stack->Particle(parent->GetLabel());
 100   if (p->GetNDaughters() > 0)
 101   {
 102     TEveTrackPropagator* rs = parent->GetPropagator();
 103     for (int d=p->GetFirstDaughter(); d>0 && d<=p->GetLastDaughter(); ++d)
 104     {
 105       TParticle* dp = stack->Particle(d);
 106       if (dp->Pt() < min_pt && dp->P() < min_p) continue;
 107
 108       TEveTrack* dtrack = new TEveTrack(dp, d, rs);
 109       char form[1000];
 110       sprintf(form,"%s [%d]", dp->GetName(), d);
 111       dtrack->SetName(form);
 112       dtrack->SetStdTitle();
 113       set_track_color(dtrack, pdg_col);
 114
 115       gEve->AddElement(dtrack, parent);
 116
 117       if (recurse)
 118         kine_daughters(dtrack, stack, min_pt, min_p, pdg_col, recurse);
 119     }
 120   }
 121 }
 122
 123 Color_t set_track_color(TEveTrack* t, Bool_t pdg_col)
 124 {
 125   if (pdg_col)
 126     t->SetMainColor(get_pdg_color(t->GetPdg()));
 127   else
 128     t->SetMainColor(30);
 129 }
 130
 131 Color_t get_pdg_color(Int_t pdg)
 132 {
 133   // PDG color indices
 134   static const Color_t DefCol   = 30;
 135   static const Color_t ECol     = 5;
 136   static const Color_t MuCol    = 6;
 137   static const Color_t GammaCol = 7;
 138   static const Color_t MesCol1  = 3;
 139   static const Color_t MesCol2  = 38;
 140   static const Color_t BarCol   = 10;
 141
 142   Int_t pdga  = TMath::Abs(pdg);
 143   Color_t col = DefCol;
 144
 145   // elementary  particles
 146   if (pdga < 100) {
 147     switch (pdga) {
 148       case 11:
 149         col = ECol; break;
 150       case 12:
 151         col = MuCol; break;
 152       case 22:
 153         col = GammaCol; break;
 154     }
 155   }
 156   // mesons and barions
 157   else if (pdga < 100000) {
 158     Int_t i  = pdga;
 159     Int_t i0 = i%10; i /= 10;
 160     Int_t i1 = i%10; i /= 10;
 161     Int_t i2 = i%10; i /= 10;
 162     Int_t i3 = i%10; i /= 10;
 163     Int_t i4 = i%10;
 164     //printf("pdg(%d) quark indices (%d,%d,%d,%d,%d) \n",pdg, i4,i3,i2, i1, i0);
 165     // meson
 166     if ((i3 == 0) && ( i4 < 2)){
 167       col = MesCol1; // quarks: i1,i2 (spin = i0)
 168       if(i1 == 3 || i2 == 3)
 169         col = MesCol2;
 170     } // barion
 171     else if ( i2 >= i1 && i3 >= i2 ) {
 172       col = BarCol; // quarks: i1,i2, i3 (spin = i0))
 173     }
 174   }
 175
 176   return col;
 177 }
 178
 179 /******************************************************************************/
 180
 181 TEveElement*
 182 kine_track(Int_t  label,
 183            Bool_t import_mother    = kTRUE, Bool_t import_daughters = kTRUE,
 184            Bool_t pdg_col          = kTRUE, Bool_t recurse          = kTRUE,
 185            TEveElement* cont = 0)
 186
 187 {
 188   // Create mother and daughters tracks with given label.
 189   // mother     -> particle with label
 190   // daughters  -> daughters of label
 191
 192   if (label < 0) {
 193     Warning("kine_track", "label not set.");
 194     return 0;
 195   }
 196
 197   AliRunLoader* rl =  AliEveEventManager::AssertRunLoader();
 198   rl->LoadKinematics();
 199   AliStack* stack = rl->Stack();
 200   if (label >= stack->GetNtrack())
 201   {
 202     Warning("kine_track", "label out of range.");
 203     return 0;
 204   }
 205
 206   TParticle* p = stack->Particle(label);
 207
 208   if (import_mother || (import_daughters && p->GetNDaughters()))
 209   {
 210     TEveTrack* toptrack = 0;
 211     TEveTrackList* tracklist = 0;
 212     TEveTrackPropagator* rs = 0;
 213
 214     if (cont == 0)
 215     {
 216       TEveTrackList* tlist = new TEveTrackList
 217         (Form("Kinematics of %d", label, p->GetNDaughters()));
 218       cont = tlist;
 219
 220       TEveTrackPropagator* rnrStyle = tlist->GetPropagator();
 221       // !!! Watch the '-', apparently different sign convention then for ESD.
 222       rnrStyle->SetMagField( -0.1*gAlice->Field()->SolenoidField() );
 223       char tooltip[1000];
 224       sprintf(tooltip,"Ndaughters=%d", p->GetNDaughters());
 225       tlist->SetTitle(tooltip);
 226       rnrStyle->fMaxOrbs = 2;
 227       rnrStyle->SetEditPathMarks(kTRUE);
 228
 229       gEve->AddElement(cont);
 230       rs = tlist->GetPropagator();
 231     }
 232     else
 233     {
 234       // check if container is TEveTrackList or TEveTrack (has rnr-style)
 235       TEveTrack* t = dynamic_cast<TEveTrack*>(cont);
 236       if (t) {
 237         rs = t->GetPropagator();
 238       } else {
 239         TEveTrackList* l = dynamic_cast<TEveTrackList*>(cont);
 240         if (l)
 241           rs = l->GetPropagator();
 242         else
 243           Error("kine_tracks.C", "TrackRenderStyle not set.");
 244       }
 245     }
 246
 247     if (import_mother)
 248     {
 249       TEveTrack* track = new TEveTrack(p, label, rs);
 250       char form[1000];
 251       sprintf(form,"%s [%d]", p->GetName(), label);
 252       track->SetName(form);
 253       track->SetStdTitle();
 254       set_track_color(track, pdg_col);
 255
 256       track->MakeTrack();
 257       gEve->AddElement(track, cont);
 258       cont = track;
 259     }
 260
 261     if (import_daughters && p->GetNDaughters())
 262     {
 263       for (int d=p->GetFirstDaughter(); d>0 && d<=p->GetLastDaughter(); ++d)
 264       {
 265         TParticle* dp = stack->Particle(d);
 266         TEveTrack* track = new TEveTrack(dp, d, rs);
 267         char form[1000];
 268         sprintf(form,"%s [%d]", dp->GetName(), d);
 269         track->SetName(form);
 270         track->SetStdTitle();
 271         set_track_color(track, pdg_col);
 272
 273         track->MakeTrack();
 274         gEve->AddElement(track, cont);
 275
 276         if (recurse)
 277           kine_daughters(track, stack, 0, 0, pdg_col, recurse);
 278       }
 279     }
 280   }
 281
 282   gEve->Redraw3D();
 283   return cont;
 284 }
 285