2 // Main authors: Paraskevi Ganoti: 2009
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 #if !defined(__CINT__) || defined(__MAKECINT__)
11 #include <TEveVector.h>
12 #include <TEveTrackPropagator.h>
13 #include <TEveVSDStructs.h>
14 #include <TEveManager.h>
16 #include <AliExternalTrackParam.h>
18 #include <AliESDEvent.h>
19 #include <AliESDtrack.h>
20 #include <AliESDkink.h>
21 #include <AliEveKink.h>
22 #include <AliEveEventManager.h>
25 void esd_kink_init_rectrackMother(TEveRecTrack& rt, const AliExternalTrackParam* tp)
27 Double_t pbuf[3], vbuf[3];
29 rt.fSign = tp->GetSign();
30 tp->GetXYZ(vbuf); rt.fV.Set(vbuf);
31 tp->GetPxPyPz(pbuf); rt.fP.Set(pbuf);
33 rt.fBeta = 1; // ep/TMath::Sqrt(ep*ep + mc*mc);
36 void esd_kink_init_rectrackDaughter(TEveRecTrack& rt, const AliExternalTrackParam* tp, TEveVector* svt,TEveVector* spt)
38 rt.fSign = tp->GetSign();
42 rt.fBeta = 1; // ep/TMath::Sqrt(ep*ep + mc*mc);
45 AliEveKink* esd_make_kink(TEveTrackPropagator* rnrStyleMoth,TEveTrackPropagator* rnrStyleDaugh, AliESDtrack* moth, AliESDtrack* daug, AliESDkink* kink, Int_t i)
51 const TVector3 p1(kink->GetMotherP());
52 rcKink.fPMother.Set(p1);
53 const TVector3 p2(kink->GetDaughterP());
54 rcKink.fPDaughter.Set(p2);
56 const TVector3 r1(kink->GetPosition());
57 rcKink.fVKink.Set(r1);
59 for (Int_t j=0; j<3; ++j) rcKink.fKinkAngle[j]=kink->GetAngle(j);
63 moth->GetTPCInnerParam()->GetXYZ(r2); rcKink.fVMother.Set(r2);
64 daug->GetOuterParam()->GetXYZ(r); rcKink.fVDaughter.Set(r);
66 esd_kink_init_rectrackMother(rcMoth, (moth->GetTPCInnerParam()));
67 rcMoth.fIndex = moth->GetID();
69 esd_kink_init_rectrackDaughter(rcDaug, daug->GetOuterParam(), &rcKink.fVKink, &rcKink.fPDaughter);
70 rcDaug.fIndex = daug->GetID();
72 AliEveKink* myKink = new AliEveKink(&rcMoth, &rcDaug, &rcKink, rnrStyleMoth,rnrStyleDaugh);
74 myKink->SetElementName(Form("ESDkink %d \n", i));
75 myKink->SetESDKinkIndex(i);
77 for (Int_t j=0; j<3; ++j) myKink->SetKinkAngle(j, kink->GetAngle(j));
78 Double_t daugProbability[10];
80 daug->GetESDpid(daugProbability);
83 // ****** Tentative particle type "concentrations"
84 Double_t c[5]={0.01, 0.01, 0.85, 0.10, 0.05};
87 AliPID daugPid(daugProbability);
89 Int_t daugMostProbPdg = 0;
91 switch (daugPid.GetMostProbable()){
93 daugMostProbPdg = 11; break;
95 daugMostProbPdg = 13; break;
97 daugMostProbPdg = 211; break;
99 daugMostProbPdg = 13; break;
102 Float_t daugMaxProbPid = daugPid.GetProbability(daugPid.GetMostProbable());
104 myKink->SetMaxProbPdgPid(daugMostProbPdg,daugMaxProbPid);//????????????
110 AliEveKinkList* esd_kink()
112 AliESDEvent* esd = AliEveEventManager::AssertESD();
113 AliEveKinkList* cont = new AliEveKinkList("ESD kink");
114 cont->SetMainColor(3); // green
115 TEveTrackPropagator* rnrStyleMoth = cont->GetPropagatorMoth();
116 rnrStyleMoth->SetMagField( 0.1*esd->GetMagneticField() );
117 TEveTrackPropagator* rnrStyleDaugh = cont->GetPropagatorDaugh();
118 rnrStyleDaugh->SetMagField( 0.1*esd->GetMagneticField() );
119 rnrStyleDaugh->SetMaxR(520);
120 gEve->AddElement(cont);
123 // for (Int_t n=0; n<esd->GetNumberOfKinks(); ++n)
125 // AliESDkink *kink = esd->GetKink(n); //???????????
126 // printf("kink number = %d, label of mother = %d , label of daug = %d --", n, kink->GetLabel(0), kink->GetLabel(1));
127 // } To be investigated...
128 for (Int_t n=0; n<esd->GetNumberOfTracks(); ++n)
130 AliESDtrack* mtrack = esd->GetTrack(n);
131 if(mtrack->GetKinkIndex(0)<0){
133 AliESDkink *kink = new AliESDkink;
135 kink=esd->GetKink(TMath::Abs(mtrack->GetKinkIndex(0))-1);
137 for (Int_t m=0; m<esd->GetNumberOfTracks(); ++m)
139 AliESDtrack * dtrack = esd->GetTrack(m);
141 if((dtrack->GetKinkIndex(0)>0)&&(dtrack->GetKinkIndex(0)==TMath::Abs(mtrack->GetKinkIndex(0)))) {
142 AliESDtrack* mothTr = esd->GetTrack(n);
143 AliESDtrack* daugTr = esd->GetTrack(m);
145 AliEveKink* myKink = esd_make_kink(rnrStyleMoth, rnrStyleDaugh, mothTr, daugTr, kink, (TMath::Abs(mtrack->GetKinkIndex(0))-1));
148 gEve->AddElement(myKink, cont);
152 } // inner track loop
155 } //mother kink index <0
156 } // Outer track loop
158 cont->SetTitle("test");