kine_tracks.C
[u/mrichter/AliRoot.git] / 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