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());
64 Double_t negProbability[10], posProbability[10];
65 Double_t negP = 0.0, posP = 0.0;
66 neg->GetESDpid(negProbability);
67 pos->GetESDpid(posProbability);
71 // ****** Tentative particle type "concentrations"
72 Double_t c[5]={0.01, 0.01, 0.85, 0.10, 0.05};
75 AliPID negPid(negProbability);
76 AliPID posPid(posProbability);
78 Int_t negMostProbPdg = 0;
79 Int_t posMostProbPdg = 0;
81 switch (negPid.GetMostProbable()){
83 negMostProbPdg = 11; break;
85 negMostProbPdg = 13; break;
87 negMostProbPdg = 211; break;
89 negMostProbPdg = 321; break;
91 negMostProbPdg = 2212; break;
93 negMostProbPdg = 211; break;
96 switch (posPid.GetMostProbable()){
98 posMostProbPdg = 11; break;
100 posMostProbPdg = 13; break;
102 posMostProbPdg = 211; break;
104 posMostProbPdg = 321; break;
106 posMostProbPdg = 2212; break;
108 posMostProbPdg = 211; break;
111 Float_t negMaxProbPid = negPid.GetProbability(negPid.GetMostProbable());
112 Float_t posMaxProbPid = posPid.GetProbability(posPid.GetMostProbable());
114 myV0->SetMaxProbPdgPid(0,negMostProbPdg,negMaxProbPid);
115 myV0->SetMaxProbPdgPid(1,posMostProbPdg,posMaxProbPid);
121 AliEveV0List* esd_V0(Bool_t onFly=kFALSE)
123 AliESDEvent* esd = AliEveEventManager::AssertESD();
125 AliESDVertex* primVertex = (AliESDVertex*) esd->GetPrimaryVertex();
127 AliEveV0List* cont = new AliEveV0List("ESD v0");
128 cont->SetMainColor(3); // green
129 TEveTrackPropagator* rnrStyle = cont->GetPropagator();
130 rnrStyle->SetMagField( 0.1*esd->GetMagneticField() );
132 gEve->AddElement(cont);
135 for (Int_t n=0; n<esd->GetNumberOfV0s(); ++n)
137 AliESDv0 *v0 = esd->GetV0(n);
139 if (v0->GetOnFlyStatus() != onFly) continue;
141 Int_t negInd = v0->GetNindex();
142 Int_t posInd = v0->GetPindex();
143 AliESDtrack* negTr = esd->GetTrack(negInd);
144 AliESDtrack* posTr = esd->GetTrack(posInd);
146 AliEveV0* myV0 = esd_make_v0(rnrStyle, primVertex, negTr,posTr, v0, n);
149 gEve->AddElement(myV0, cont);
154 cont->SetTitle("test");