1 // Import tracks from kinematics-tree / particle-stack.
2 // Preliminary/minimal solution.
4 Reve::TrackList* kine_tracks(Double_t min_pt=0.5, Double_t max_pt=100)
6 AliRunLoader* rl = Alieve::Event::AssertRunLoader();
8 AliStack* stack = rl->Stack();
10 Error("kine_tracks.C", "can not get kinematics.");
14 Reve::TrackList* cont = new Reve::TrackList("Kine Tracks");
15 cont->SetMainColor(Color_t(6));
16 Reve::TrackRnrStyle* rnrStyle = cont->GetRnrStyle();
18 // !!! Watch the '-', apparently different sign convention then for ESD.
19 rnrStyle->SetMagField( - gAlice->Field()->SolenoidField() );
21 gReve->AddRenderElement(cont);
24 Int_t N = stack->GetNtrack();
25 for (Int_t i=0; i<N; ++i)
27 if(stack->IsPhysicalPrimary(i))
29 TParticle* p = stack->Particle(i);
30 Double_t pT = p->Pt();
31 if (pT<min_pt || pT>max_pt) continue;
34 Reve::Track* track = new Reve::Track(p, i, rnrStyle);
36 //PH The line below is replaced waiting for a fix in Root
37 //PH which permits to use variable siza arguments in CINT
38 //PH on some platforms (alphalinuxgcc, solariscc5, etc.)
39 //PH track->SetName(Form("%s [%d]", p->GetName(), i));
41 sprintf(form,"%s [%d]", p->GetName(), i);
43 track->SetRnrStyle(rnrStyle);
44 gReve->AddRenderElement(cont, track);
51 kt.SetPathMarks(cont,stack, rl->TreeTR());
52 cont->SetEditPathMarks(kTRUE);
55 //PH const Text_t* tooltip = Form("pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
57 sprintf(tooltip,"pT ~ (%.2lf, %.2lf), N=%d", min_pt, max_pt, count);
58 cont->SetTitle(tooltip); // Not broadcasted automatically ...