Disable redraw while adding Tracks/RenderElements.
[u/mrichter/AliRoot.git] / EVE / alice-macros / kine_tracks.C
1 // Import tracks from kinematics-tree / particle-stack.
2 // Preliminary/minimal solution.
3 #include "TParticlePDG.h"
4
5 // PDG color indices
6 static Color_t DefCol  = 30;
7 static Color_t ECol    = 5;
8 static Color_t MuCol   = 6;
9 static Color_t GamaCol = 7; 
10 static Color_t MesCol1 = 3;
11 static Color_t MesCol2 = 38;
12 static Color_t BarCol  = 10;
13
14
15 Reve::TrackList*
16 kine_tracks(Double_t min_pt=0.5, Double_t max_pt=100, Bool_t pdg_col= kFALSE)
17 {
18   AliRunLoader* rl =  Alieve::Event::AssertRunLoader();
19   rl->LoadKinematics();
20   AliStack* stack = rl->Stack();
21   if (!stack) {
22     Error("kine_tracks.C", "can not get kinematics.");
23     return 0;
24   }
25
26   gReve->DisableRedraw();
27  
28   Reve::TrackList* cont = new Reve::TrackList("Kine Tracks"); 
29   cont->SetMainColor(Color_t(6));
30   Reve::TrackRnrStyle* rnrStyle = cont->GetRnrStyle();
31   rnrStyle->fColor = 8;
32   // !!! Watch the '-', apparently different sign convention then for ESD.
33   rnrStyle->SetMagField( - gAlice->Field()->SolenoidField() );
34
35   gReve->AddRenderElement(cont);
36   Int_t count = 0;
37   Int_t N = stack->GetNtrack();
38   for (Int_t i=0; i<N; ++i) 
39   {
40     if(stack->IsPhysicalPrimary(i)) 
41     {
42       TParticle* p = stack->Particle(i);
43       Double_t  pT = p->Pt();
44       if (pT<min_pt || pT>max_pt) continue;
45
46       ++count;
47       Reve::Track* track = new Reve::Track(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       TParticlePDG* pdgp = p->GetPDG();
57       track->SetMainColor(get_pdg_color(pdgp->PdgCode()));
58       gReve->AddRenderElement(cont, track);
59     }
60   }
61   // set path marks
62   Alieve::KineTools kt; 
63   kt.SetDaughterPathMarks(cont, stack);
64   rl->LoadTrackRefs();
65   kt.SetTrackReferences(cont, rl->TreeTR());
66   cont->SetEditPathMarks(kTRUE);
67
68   //PH  const Text_t* tooltip = Form("pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
69   char tooltip[1000];
70   sprintf(tooltip,"pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
71   cont->SetTitle(tooltip); // Not broadcasted automatically ...
72   cont->UpdateItems();
73
74   cont->MakeTracks();
75   cont->MakeMarkers();
76
77   gReve->EnableRedraw();
78   gReve->Redraw3D();
79
80   return cont;
81 }
82
83
84 Color_t get_pdg_color(Int_t pdg){
85   Int_t pdga = TMath::Abs(pdg);
86   Color_t col = Reve::DefCol;
87
88   // elementary  particles
89   if (pdga < 100) {
90     switch (pdga) {
91       case 11:  
92         col = ECol; break; 
93       case 12:
94         col = MuCol; break;
95       case 22:
96         col = GammaCol; break;
97     }
98   }
99   else if (pdga < 100000){ 
100     Int_t i  = pdga;
101     Int_t i0 = i%10; i /= 10;
102     Int_t i1 = i%10; i /= 10; 
103     Int_t i2 = i%10; i /= 10; 
104     Int_t i3 = i%10; i /= 10; 
105     Int_t i4 = i%10;
106     //printf("pdg(%d) quark indices (%d,%d,%d,%d,%d) \n",pdg, i4,i3,i2, i1, i0);
107     // meson
108     if ((i3 == 0) && ( i4 < 2)){
109       col = MesCol1; // quarks: i1,i2 (spin = i0)
110       if(i1 == 3 || i2 == 3)
111         col = MesCol2;
112     } // barion
113     else if ( i2 >= i1 && i3 >= i2 ) {
114       col = BarCol; // quarks: i1,i2, i3 (spin = i0))
115     }
116   }
117   return col;
118 }
119
120
121 // Create mother and daughters tracks with given label.
122
123 Reve::RenderElement*
124 kine_track(Int_t  label,
125            Bool_t import_mother    = kTRUE,
126            Bool_t import_daughters = kTRUE,
127            Reve::RenderElement*    cont = 0)
128
129 {
130   if (label < 0) {
131     Warning("kine_track", "label not set.");
132     return 0;
133   }
134  
135   AliRunLoader* rl =  Alieve::Event::AssertRunLoader();
136   rl->LoadKinematics();
137   AliStack* stack = rl->Stack();
138   TParticle* p = stack->Particle(label);
139
140   if (import_mother || (import_daughters && p->GetNDaughters()))
141   {
142     Track* toptrack = 0;
143     TrackList* tracklist = 0;  
144     TrackRnrStyle* rs = 0;
145
146     if (cont == 0)
147     {
148       cont = new TrackList(Form("Kinematics of %d", label, p->GetNDaughters()));
149
150       Reve::TrackRnrStyle* rnrStyle = cont->GetRnrStyle();
151       // !!! Watch the '-', apparently different sign convention then for ESD.
152       rnrStyle->SetMagField( - gAlice->Field()->SolenoidField() );
153       char tooltip[1000];
154       sprintf(tooltip,"Ndaughters=%d", p->GetNDaughters());
155       cont->SetTitle(tooltip);
156       cont->SelectByPt(0.2, 100);
157       rnrStyle->fColor   = 8;
158       rnrStyle->fMaxOrbs = 8;
159       cont->SetEditPathMarks(kTRUE);
160       gReve->AddRenderElement(cont);
161       rs = cont->GetRnrStyle();
162     }
163     else {
164       // check if argument is TrackList
165       Reve::Track* t = dynamic_cast<Reve::Track*>(cont);
166       if(t) 
167       {
168         rs = t->GetRnrStyle();
169       }
170       else {
171         Reve::TrackList* l = dynamic_cast<Reve::TrackList*>(cont);
172         if(l)
173         {
174           rs = l->GetRnrStyle();
175         }
176         else {
177           Error("kine_tracks.C", "TrackRenderStyle not set.");
178         }
179       }
180     }
181
182     if (import_mother)
183     {
184       Track* track = new Reve::Track(p, label, rs);  
185       char form[1000];
186       sprintf(form,"%s [%d]", p->GetName(), label);
187       track->SetName(form);
188       gReve->AddRenderElement(cont, track);
189
190     }
191
192     if (import_daughters && p->GetNDaughters()) 
193     {
194       for (int d=p->GetFirstDaughter(); d>0 && d<=p->GetLastDaughter(); ++d) 
195       { 
196         TParticle* dp = stack->Particle(d);
197         Track* track = new Reve::Track(dp, d, rs);  
198         char form[1000];
199         sprintf(form,"%s [%d]", dp->GetName(), d);
200         track->SetName(form);
201         track->MakeTrack();
202         gReve->AddRenderElement(cont, track);
203       }
204     }
205   }
206
207   cont->UpdateItems();
208   gReve->Redraw3D();
209   return cont;
210 }
211