2 // Main authors: Matevz Tadel & Alja Mrak-Tadel: 2006, 2007
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 **************************************************************************/
10 void esd_v0_init_rectrack(TEveRecTrack& rt, AliExternalTrackParam* tp)
12 Double_t pbuf[3], vbuf[3];
14 rt.fSign = tp->GetSign();
15 tp->GetXYZ(vbuf); rt.fV.Set(vbuf);
16 tp->GetPxPyPz(pbuf); rt.fP.Set(pbuf);
17 // Double_t ep = at->GetP(), mc = at->GetMass();
18 rt.fBeta = 1; // ep/TMath::Sqrt(ep*ep + mc*mc);
21 AliEveV0* esd_make_v0(TEveTrackPropagator* rnrStyle, AliESDVertex* primVtx,
22 AliESDtrack* neg, AliESDtrack* pos, AliESDv0* v0, Int_t i)
29 v0->GetNPxPyPz(p[0], p[1], p[2]);
31 v0->GetPPxPyPz(p[0], p[1], p[2]);
34 v0->GetPxPyPz(p[0], p[1], p[2]);
38 v0->GetXYZ(v[0], v[1], v[2]);
41 v0->GetParamN()->GetXYZ(v); rcV0.fVNeg.Set(v);
42 v0->GetParamP()->GetXYZ(v); rcV0.fVPos.Set(v);
44 rcV0.fV0Birth.Set(primVtx->GetXv(), primVtx->GetYv(), primVtx->GetZv());
46 // Simulation data not directly available in AliESDv0
47 //rcV0.fDLabel[0] = v0->GetNindex();
48 //rcV0.fDLabel[1] = v0->GetPindex();
50 esd_v0_init_rectrack(rcNeg, v0->GetParamN());
51 rcNeg.fIndex = v0->GetNindex();
52 esd_v0_init_rectrack(rcPos, v0->GetParamP());
53 rcPos.fIndex = v0->GetPindex();
55 AliEveV0* myV0 = new AliEveV0(&rcNeg, &rcPos, &rcV0, rnrStyle);
56 myV0->SetElementName(Form("ESDv0 %d", i));
57 myV0->SetElementTitle(Form("OnFly: %d\nDCA %f",
59 v0->GetDcaV0Daughters()));
61 myV0->SetOnFlyStatus(v0->GetOnFlyStatus());
62 myV0->SetDaughterDCA(v0->GetDcaV0Daughters());
68 AliEveV0List* esd_V0(Bool_t onFly=kFALSE)
70 AliESDEvent* esd = AliEveEventManager::AssertESD();
72 AliESDVertex* primVertex = (AliESDVertex*) esd->GetPrimaryVertex();
74 AliEveV0List* cont = new AliEveV0List("ESD v0");
75 cont->SetMainColor(3); // green
76 TEveTrackPropagator* rnrStyle = cont->GetPropagator();
77 rnrStyle->SetMagField( 0.1*esd->GetMagneticField() );
79 gEve->AddElement(cont);
82 for (Int_t n=0; n<esd->GetNumberOfV0s(); ++n)
84 AliESDv0 *v0 = esd->GetV0(n);
86 if (v0->GetOnFlyStatus() != onFly) continue;
88 Int_t negInd = v0->GetNindex();
89 Int_t posInd = v0->GetPindex();
90 AliESDtrack* negTr = esd->GetTrack(negInd);
91 AliESDtrack* posTr = esd->GetTrack(posInd);
93 AliEveV0* myV0 = esd_make_v0(rnrStyle, primVertex, negTr,posTr, v0, n);
96 gEve->AddElement(myV0, cont);
101 cont->SetTitle("test");