]> git.uio.no Git - u/mrichter/AliRoot.git/blame - HLT/EVE/AliHLTEveHLT.cxx
git://git.uio.no / u / mrichter / AliRoot.git / blame
? search:
summary | shortlog | log | commit | commitdiff | tree
blob | blame (incremental) | history | HEAD
Update master to aliroot
[u/mrichter/AliRoot.git] / HLT / EVE / AliHLTEveHLT.cxx
CommitLineData
35e470ad 1// $Id$
33791895 2/**************************************************************************
3 * This file is property of and copyright by the ALICE HLT Project *
4 * ALICE Experiment at CERN, All rights reserved. *
5 * *
6 * Primary Authors: Svein Lindal <slindal@fys.uio.no > *
7 * for The ALICE HLT Project. *
8 * *
9 * Permission to use, copy, modify and distribute this software and its *
10 * documentation strictly for non-commercial purposes is hereby granted *
11 * without fee, provided that the above copyright notice appears in all *
12 * copies and that both the copyright notice and this permission notice *
13 * appear in the supporting documentation. The authors make no claims *
14 * about the suitability of this software for any purpose. It is *
15 * provided "as is" without express or implied warranty. *
16 **************************************************************************/
17
18/// @file AliHLTEvePhos.cxx
19/// @author Svein Lindal <slindal@fys.uio.no>
20/// @brief HLT class for the HLT EVE display
21
22#include "AliHLTEveHLT.h"
23#include "AliHLTHOMERBlockDesc.h"
24#include "AliHLTEveBase.h"
fd2adb88 25#include "AliEveHLTEventManager.h"
c8fc4838 26#include "AliHLTGlobalTriggerDecision.h"
33791895 27#include "TEveManager.h"
28#include "TEvePointSet.h"
29#include "TEveTrack.h"
73ac9862 30#include "TEveElement.h"
33791895 31#include "TCanvas.h"
32#include "AliESDEvent.h"
4883ef27 33#include "AliESDtrack.h"
33791895 34#include "TEveTrackPropagator.h"
35#include "AliEveTrack.h"
36#include "TEveVSDStructs.h"
37#include "TString.h"
aac9c09b 38#include "tracking-ca/AliHLTTPCCATrackParam.h"
39#include "tracking-ca/AliHLTTPCCATrackConvertor.h"
33791895 40#include "AliEveMagField.h"
f1a55ef4 41#include "TH1.h"
33791895 42#include "TH1F.h"
43#include "TH2F.h"
fd2adb88 44#include "TThread.h"
33791895 45
46ClassImp(AliHLTEveHLT)
47
48AliHLTEveHLT::AliHLTEveHLT() :
fd2adb88 49 AliHLTEveBase("TPC tracks"),
33791895 50 fTrueField(kFALSE),
51 fUseIpOnFailedITS(kFALSE),
52 fUseRkStepper(kFALSE),
353f10b3 53 fTrackList(NULL),
73ac9862 54 fTrackLists(NULL),
fd2adb88 55 fOldTrackList(NULL),
bc1a13f1 56 fPointSetVertex(NULL),
353f10b3 57 fTrCanvas(NULL),
f1a55ef4 58 fVertexCanvas(NULL),
353f10b3 59 fHistEta(NULL),
353f10b3 60 fHistPhi(NULL),
61 fHistnClusters(NULL),
f1a55ef4 62 fHistMult(NULL),
4883ef27 63 fHistDCAr(NULL),
f1a55ef4 64 fTrCount(0),
73ac9862 65 fVCount(0),
66 fNTrackBins(10)
33791895 67{
68 // Constructor.
4883ef27 69 CreateHistograms();
70
33791895 71}
c8fc4838 72///___________________________________________________________________
33791895 73AliHLTEveHLT::~AliHLTEveHLT()
74{
75 //Destructor, not implemented
76 if(fTrackList)
77 delete fTrackList;
78 fTrackList = NULL;
73ac9862 79
80 if(fTrackLists)
81 delete fTrackLists;
82 fTrackLists = NULL;
83
84
353f10b3 85}
4883ef27 86
73ac9862 87
4883ef27 88///________________________________________________________________________
89void AliHLTEveHLT::CreateHistograms(){
90 //See header file for documentation
91
92 //fHistPt = new TH1F("fHistPt", "transverse momentum", 100, 0, 10); // KK
93 //fHistP = new TH1F("fHistP", "signed momentum", 100,-7, 7);
94
95 //fHistPt ->SetXTitle("p_{t} (GeV/c)"); // KK
96 //fHistP ->SetXTitle("P*charge (GeV/c)");
97
98 // fHistTheta = new TH1F("fHistTheta", "polar angle", 180, 0,180);
99 // fHistTheta->SetXTitle("#theta (degrees)");
100
101 fHistEta = new TH1F("fHistEta", "pseudorapidity", 100,-2, 2);
102 fHistEta ->SetXTitle("#eta");
103
104 fHistPhi = new TH1F("fHistPhi", "azimuthal angle", 180, 0,360);
105 fHistPhi ->SetXTitle("#phi (degrees)");
106
107 fHistnClusters = new TH1F("fHistnClusters","TPC clusters per track", 160, 0,160);
108
437159d3 109 fHistMult = new TH1F("fHistMult", "event track multiplicity",150, 0, 15000);
4883ef27 110
111 fHistDCAr = new TH1F("fHistDCAr", "DCA r", 200, -100, 100);
112
113}
114
115
fd2adb88 116///_____________________________________________________________________
33791895 117void AliHLTEveHLT::ProcessBlock(AliHLTHOMERBlockDesc * block) {
118 //See header file for documentation
119 if ( ! block->GetDataType().CompareTo("ALIESDV0") ) {
73ac9862 120 ProcessEsdBlock(block);
33791895 121 }
122
123 else if ( ! block->GetDataType().CompareTo("ROOTTOBJ") ) {
124 //processROOTTOBJ( block, gHLTText );
125 }
126
127 else if ( ! block->GetDataType().CompareTo("HLTRDLST") ) {
c8fc4838 128 cout << "ignoring hlt rdlst"<<endl;
33791895 129 //processHLTRDLST( block );
130 }
131
c8fc4838 132 else if ( ! block->GetDataType().CompareTo("GLOBTRIG") ) {
133 ProcessGlobalTrigger( block );
134 }
135
33791895 136 else if ( !block->GetDataType().CompareTo("ROOTHIST") ) {
137 if( !fCanvas ) {
4883ef27 138 fCanvas = CreateCanvas("PVtx/Tr QA", "PrimV");
139 fCanvas->Divide(3, 3);
33791895 140 }
141 ProcessHistograms( block , fCanvas);
fd2adb88 142 }
33791895 143}
144
fd2adb88 145///____________________________________________________________________________
33791895 146void AliHLTEveHLT::UpdateElements() {
147 //See header file for documentation
73ac9862 148
4883ef27 149 DrawHistograms();
73ac9862 150
151 if(fTrackLists) {
152 for(Int_t il = 0; il < fNTrackBins; il++) {
153 TEveTrackList * trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FindChild(Form("Tracks_%d", il)));
040245b5 154 if(trackList) trackList->ElementChanged();
73ac9862 155 }
156 }
157
bc1a13f1 158 if(fPointSetVertex) fPointSetVertex->ResetBBox();
f1a55ef4 159 if(fTrCanvas) fTrCanvas->Update();
160 if(fVertexCanvas) fVertexCanvas->Update();
161 if(fCanvas) fCanvas->Update();
162
33791895 163}
73ac9862 164
fd2adb88 165///_________________________________________________________________________________
33791895 166void AliHLTEveHLT::ResetElements(){
fd2adb88 167 //See header file for documentation
168
169 cout << "destroy"<<endl;
73ac9862 170
171 if(fTrackLists) {
172 for(Int_t il = 0; il < fNTrackBins; il++) {
173 TEveTrackList * trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FindChild(Form("Tracks_%d", il)));
040245b5 174 if(trackList) trackList->DestroyElements();
73ac9862 175 }
fd2adb88 176 }
73ac9862 177
fd2adb88 178
2d91aa68 179 fTrCount = 0;
180 fVCount = 0;
181
bc1a13f1 182 if(fPointSetVertex) fPointSetVertex->Reset();
fd2adb88 183 cout<< "reset done"<<endl;
33791895 184 fHistoCount = 0;
185
186}
187
fd2adb88 188///_____________________________________________________________________________________
189void * AliHLTEveHLT::DestroyGarbage(void * arg) {
190 AliHLTEveHLT * hlt = reinterpret_cast<AliHLTEveHLT*>(arg);
191 if(hlt) hlt->DestroyOldTrackList();
192 return (void*)0;
193}
194///_____________________________________________________________________________________
195void AliHLTEveHLT::DestroyOldTrackList() {
196 cout << "Destroying the old tracklist's elements"<<endl;
197 fOldTrackList->DestroyElements();
198 cout << "Destroying the old tracklist itself"<<endl;
199 fOldTrackList->Destroy();
200}
201
202///_____________________________________________________________________________________
33791895 203void AliHLTEveHLT::ProcessHistograms(AliHLTHOMERBlockDesc * block, TCanvas * canvas) {
204 //See header file for documentation
f1a55ef4 205 if (!fTrCanvas) {
4883ef27 206 fTrCanvas = CreateCanvas("ESD", "ESD");
f1a55ef4 207 fTrCanvas->Divide(4, 2);
208 }
4883ef27 209
f1a55ef4 210 if(!fVertexCanvas) {
4883ef27 211 fVertexCanvas = CreateCanvas("Vtx", "Vtx");
f1a55ef4 212 fVertexCanvas->Divide(4, 2);
213 }
214
215
216
33791895 217 if ( ! block->GetClassName().CompareTo("TH1F")) {
218 TH1F* histo = reinterpret_cast<TH1F*>(block->GetTObject());
219 if( histo ){
220 TString name(histo->GetName());
f1a55ef4 221 cout << "TH1F " <<name << endl;
33791895 222 if( !name.CompareTo("primVertexZ") ){
223 canvas->cd(2);
224 histo->Draw();
225 }else if( !name.CompareTo("primVertexX") ){
226 canvas->cd(3);
227 histo->Draw();
228 }else if( !name.CompareTo("primVertexY") ){
229 canvas->cd(4);
230 histo->Draw();
f1a55ef4 231 } else if ( name.Contains("Track")) {
232 AddHistogramToCanvas(histo, fTrCanvas, fTrCount);
233 } else {
234 AddHistogramToCanvas(histo, fVertexCanvas, fVCount);
33791895 235 }
236 }
237 } else if ( ! block->GetClassName().CompareTo("TH2F")) {
238 TH2F *hista = reinterpret_cast<TH2F*>(block->GetTObject());
239 if (hista ){
240 TString name(hista->GetName());
f1a55ef4 241 cout << "TH2F " << name << endl;
33791895 242 if( !name.CompareTo("primVertexXY")) {
243 canvas->cd(1);
244 hista->Draw();
f1a55ef4 245 } else if ( name.Contains("Track")) {
246 AddHistogramToCanvas(hista, fTrCanvas, fTrCount);
247 } else {
248 AddHistogramToCanvas(hista, fVertexCanvas, fVCount);
249 }
33791895 250 }
251 }
252 canvas->cd();
33791895 253}
254
f1a55ef4 255///_____________________________________________________________________________________
256void AliHLTEveHLT::AddHistogramToCanvas(TH1* histogram, TCanvas * canvas, Int_t &cdCount) {
257 canvas->cd(++cdCount);
258 histogram->Draw();
259}
260
261
fd7a3071 262///________________________________________________________________________________________
33791895 263void AliHLTEveHLT::CreateTrackList() {
264 //See header file for documentation
265 fTrackList = new TEveTrackList("ESD Tracks");
73ac9862 266 fTrackList->SetMainColor(kOrange);
fd2adb88 267 AddElement(fTrackList);
33791895 268}
73ac9862 269
270///________________________________________________________________________________________
271void AliHLTEveHLT::CreateTrackLists() {
272 //See header file for documentation
273 fTrackLists = new TEveElementList("ESD Track lists");
274 AddElement(fTrackLists);
275 for(Int_t i = 0; i < 10; i++) {
276 TEveTrackList * trackList = new TEveTrackList(Form("Tracks_%d", i));
277 trackList->SetMainColor(kOrange-i);
278 fTrackLists->AddElement(trackList);
279 }
280}
281
33791895 282
fd7a3071 283///_________________________________________________________________________________________
bc1a13f1 284void AliHLTEveHLT::CreateVertexPointSet() {
285 //See header file for documentation
286 fPointSetVertex = new TEvePointSet("Primary Vertex");
287 fPointSetVertex->SetMainColor(6);
288 fPointSetVertex->SetMarkerStyle((Style_t)kFullStar);
289
fd2adb88 290 AddElement(fPointSetVertex);
bc1a13f1 291}
292
c8fc4838 293///________________________________________________________________________
294void AliHLTEveHLT::ProcessGlobalTrigger( AliHLTHOMERBlockDesc * block ) {
295 //See header file for documentation
296 AliHLTGlobalTriggerDecision * decision = dynamic_cast<AliHLTGlobalTriggerDecision*>(block->GetTObject());
040245b5 297 if(decision) decision->Print();
c8fc4838 298
299}
300
33791895 301
c8fc4838 302///______________________________________________________________________
73ac9862 303void AliHLTEveHLT::ProcessEsdBlock( AliHLTHOMERBlockDesc * block) {
33791895 304 //See header file for documentation
305
306 AliESDEvent* esd = (AliESDEvent *) (block->GetTObject());
fd2adb88 307 if (!esd) return;
308
73ac9862 309 ProcessEsdEvent(esd);
310}
311
312
313
314///_________________________________________________________________________________
315Color_t AliHLTEveHLT::GetColor(Float_t pt) {
316 //See header file
317 Color_t baseColor = kOrange;
318
319 Float_t binlimit = 0.1;
320 for(Int_t i = 0; i< 10; i++) {
321
322 if (pt < binlimit)
323 return baseColor - i;
324
325 binlimit +=0.1;
326 }
327
328 return baseColor - 9;
329
fd2adb88 330}
73ac9862 331
332///_________________________________________________________________________________
333Int_t AliHLTEveHLT::GetColorBin(Float_t pt) {
334 //See header file
335
336 Float_t binlimit = 0.1;
337 for(Int_t i = 0; i< 10; i++) {
338
339 if (pt < binlimit)
340 return i;
341
342 binlimit +=0.1;
343 }
344
345 return 9;
346
347}
348
349///____________________________________________________________________
350void AliHLTEveHLT::ProcessEsdEvent( AliESDEvent * esd ) {
351 //See header file for documentation
352 if(!fPointSetVertex) CreateVertexPointSet();
353 if(!fTrackLists) CreateTrackLists();
fd2adb88 354
fd7a3071 355 cout << "ProcessESDEvent() :"<< esd->GetEventNumberInFile()<< " " << esd->GetNumberOfCaloClusters() << " tracks : " << esd->GetNumberOfTracks() << endl;
fd2adb88 356
357 //fEventManager->SetRunNumber(esd->GetRunNumber());
358
bc1a13f1 359 Double_t vertex[3];
cb9a310a 360 const AliESDVertex * esdVertex = esd->GetPrimaryVertex();
bc1a13f1 361
cb9a310a 362 if(esdVertex) {
363 esdVertex->GetXYZ(vertex);
364 fPointSetVertex->SetNextPoint(vertex[0], vertex[1], vertex[2]);
365 }
33791895 366
040245b5 367 TEveTrackList * trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FirstChild());
368 if(trackList) SetUpTrackPropagator(trackList->GetPropagator(),-0.1*esd->GetMagneticField(), 520);
33791895 369 for (Int_t iter = 0; iter < esd->GetNumberOfTracks(); ++iter) {
4883ef27 370
73ac9862 371 AliESDtrack * esdTrack = dynamic_cast<AliESDtrack*>(esd->GetTrack(iter));
372 FillTrackList(esdTrack);
373 FillHistograms(esdTrack);
353f10b3 374
33791895 375 }
4883ef27 376 fHistMult->Fill(esd->GetNumberOfTracks()); // KK
33791895 377
73ac9862 378 //BALLE hardcoded size
379 for(Int_t il = 0; il < fNTrackBins; il++) {
040245b5 380 trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FindChild(Form("Tracks_%d", il)));
73ac9862 381 trackList->MakeTracks();
382 }
383}
384///__________________________________________________________________________
385void AliHLTEveHLT::FillTrackList(AliESDtrack * esdTrack) {
386 //See header file for documentation
040245b5 387 TEveTrackList * trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FirstChild());
388 if (!trackList) return;
389
390 AliEveTrack* track = dynamic_cast<AliEveTrack*>(MakeEsdTrack(esdTrack, trackList));
73ac9862 391 Int_t bin = GetColorBin(esdTrack->Pt());
040245b5 392 trackList = dynamic_cast<TEveTrackList*>(fTrackLists->FindChild(Form("Tracks_%d", bin)));
73ac9862 393 if(trackList) {
394 track->SetAttLineAttMarker(trackList);
395 trackList->AddElement(track);
396 } else cout << "BALLE"<<endl;
397
398
33791895 399}
400
d0cbdfc0 401
402///____________________________________________________________________________________
403void AliHLTEveHLT::FillHistograms(AliESDtrack * esdTrack) {
404
405 if(esdTrack->GetTPCNcls() == 0) return;
406
407 fHistEta->Fill(esdTrack->Eta());
408 // fHistTheta->Fill(esdTrack->Theta()*TMath::RadToDeg());
409 fHistPhi->Fill(esdTrack->Phi()*TMath::RadToDeg());
410
411
412 Float_t DCAr, DCAz = -99;
413 esdTrack->GetImpactParametersTPC(DCAr, DCAz);
414 fHistDCAr->Fill(DCAr);
415
416
417 if(esdTrack->GetStatus()&AliESDtrack::kTPCin || (esdTrack->GetStatus()&AliESDtrack::kTPCin && esdTrack->GetStatus()&AliESDtrack::kITSin)){
418 fHistnClusters->Fill(esdTrack->GetTPCNcls());
419 }
420}
421
353f10b3 422void AliHLTEveHLT::DrawHistograms(){
423 //See header file for documentation
4883ef27 424 if(!fCanvas) {
425 fCanvas = CreateCanvas("PVtx/Tr QA", "Primary vertex, Track QA");
426 fCanvas->Divide(3, 3);
427 }
428
429 fCanvas->cd(5);
430 fHistEta->Draw();
431
432 fCanvas->cd(6);
433 fHistPhi->Draw();
434
435 fCanvas->cd(7);
436 fHistnClusters->Draw();
437
438 fCanvas->cd(8);
439 fHistMult->Draw();
353f10b3 440
4883ef27 441 fCanvas->cd(9);
442 fHistDCAr->Draw();
353f10b3 443
4883ef27 444 fCanvas->cd();
353f10b3 445
446}
447
33791895 448AliEveTrack* AliHLTEveHLT::MakeEsdTrack (AliESDtrack *at, TEveTrackList* cont) {
449 //See header file for documentation
73ac9862 450
451
452
33791895 453
454
455 const double kCLight = 0.000299792458;
456 double bz = - kCLight*10.*( cont->GetPropagator()->GetMagField(0,0,0).fZ);
457
458 Bool_t innerTaken = kFALSE;
459 if ( ! at->IsOn(AliESDtrack::kITSrefit) && fUseIpOnFailedITS)
460 {
461 //tp = at->GetInnerParam();
462 innerTaken = kTRUE;
463 }
464
465 // Add inner/outer track parameters as path-marks.
466
467 Double_t pbuf[3], vbuf[3];
468
469 AliExternalTrackParam trackParam = *at;
470
471 // take parameters constrained to vertex (if they are)
472
473 if( at->GetConstrainedParam() ){
474 trackParam = *at->GetConstrainedParam();
475 }
476 else if( at->GetInnerParam() ){
db742b8e 477 float x = at->GetX();
478 float y = at->GetY();
479 if(sqrt(x*x+y*y)>.5 ) trackParam = *(at->GetInnerParam());
33791895 480 }
481 if( at->GetStatus()&AliESDtrack::kTRDin ){
482 // transport to TRD in
483 trackParam = *at;
484 trackParam.PropagateTo( 290.45, -10.*( cont->GetPropagator()->GetMagField(0,0,0).fZ) );
485 }
486
487 TEveRecTrack rt;
488 {
489 rt.fLabel = at->GetLabel();
490 rt.fIndex = (Int_t) at->GetID();
491 rt.fStatus = (Int_t) at->GetStatus();
492 rt.fSign = (Int_t) trackParam.GetSign();
493 trackParam.GetXYZ(vbuf);
494 trackParam.GetPxPyPz(pbuf);
495 rt.fV.Set(vbuf);
496 rt.fP.Set(pbuf);
497 Double_t ep = at->GetP(), mc = at->GetMass();
498 rt.fBeta = ep/TMath::Sqrt(ep*ep + mc*mc);
499 }
500
501 AliEveTrack* track = new AliEveTrack(&rt, cont->GetPropagator());
33791895 502 track->SetName(Form("AliEveTrack %d", at->GetID()));
503 track->SetElementTitle(CreateTrackTitle(at));
504 track->SetSourceObject(at);
505
506
507 // Set reference points along the trajectory
508 // and the last point
509
510 {
511 TEvePathMark startPoint(TEvePathMark::kReference);
512 trackParam.GetXYZ(vbuf);
513 trackParam.GetPxPyPz(pbuf);
514 startPoint.fV.Set(vbuf);
515 startPoint.fP.Set(pbuf);
516 rt.fV.Set(vbuf);
517 rt.fP.Set(pbuf);
518 Double_t ep = at->GetP(), mc = at->GetMass();
519 rt.fBeta = ep/TMath::Sqrt(ep*ep + mc*mc);
520
521 track->AddPathMark( startPoint );
522 }
523
2d91aa68 524 bool ok = 1;
7036e5c1 525 if( at->GetOuterParam() && at->GetTPCPoints(2)>80 ){
2d91aa68 526
33791895 527 //
528 // use AliHLTTPCCATrackParam propagator
529 // since AliExternalTrackParam:PropagateTo()
530 // has an offset at big distances
531 //
7036e5c1 532 double rot = at->GetOuterParam()->GetAlpha() - trackParam.GetAlpha();
533 double crot = cos(rot), srot = sin(rot);
2d91aa68 534 double xEnd = at->GetTPCPoints(2)*crot - at->GetTPCPoints(3)*srot;
7036e5c1 535 // take parameters constrained to vertex (if they are)
536
537
33791895 538 AliHLTTPCCATrackParam t;
539 AliHLTTPCCATrackConvertor::SetExtParam( t, trackParam );
540
541 Double_t x0 = trackParam.GetX();
7036e5c1 542 Double_t dx = xEnd - x0;
33791895 543
2d91aa68 544 if( dx<0 ) ok = 0;
33791895 545 //
546 // set a reference at the half of trajectory for better drawing
547 //
548
e090fd32 549 if( ok ){
550 double dxx=dx/2;
551 if( TMath::Abs(dxx)>=1. ){
552
553 if( !t.TransportToX(x0+dxx, bz, .999 ) ){
554 ok = 0;
555 } else {
556 AliExternalTrackParam tt;
557 AliHLTTPCCATrackConvertor::GetExtParam( t, tt, trackParam.GetAlpha() );
558 tt.GetXYZ(vbuf);
559 tt.GetPxPyPz(pbuf);
560 TEvePathMark midPoint(TEvePathMark::kReference);
561 midPoint.fV.Set(vbuf);
562 midPoint.fP.Set(pbuf);
563 track->AddPathMark( midPoint );
564 }
2d91aa68 565 }
33791895 566 }
7036e5c1 567
33791895 568 //
569 // Set a reference at the end of the trajectory
570 // and a "decay point", to let the event display know where the track ends
571 //
572
2d91aa68 573
da0437de 574 ///Removing meaningless loop, please have a look Sergey!!
575
576 // FIXME this loop does not make sense anymore
577 // Matthias 2010-05-02
578 if( !t.TransportToX(x0+dx, bz, .999 ) ){
579 ok = 0;
580 }
581
2d91aa68 582 {
583 if( !ok ){
584 AliHLTTPCCATrackConvertor::SetExtParam( t, trackParam );
585 }
586 AliExternalTrackParam tt;
7036e5c1 587 AliHLTTPCCATrackConvertor::GetExtParam( t, tt, trackParam.GetAlpha() );
588 tt.GetXYZ(vbuf);
589 tt.GetPxPyPz(pbuf);
33791895 590 TEvePathMark endPoint(TEvePathMark::kReference);
591 TEvePathMark decPoint(TEvePathMark::kDecay);
592 endPoint.fV.Set(vbuf);
593 endPoint.fP.Set(pbuf);
594 decPoint.fV.Set(vbuf);
595 decPoint.fP.Set(pbuf);
596 track->AddPathMark( endPoint );
2d91aa68 597 track->AddPathMark( decPoint );
33791895 598 }
599 }
600
2d91aa68 601 if ( ok && at->IsOn(AliESDtrack::kTPCrefit))
33791895 602 {
603 if ( ! innerTaken)
604 {
605 AddTrackParamToTrack(track, at->GetInnerParam());
606 }
607 AddTrackParamToTrack(track, at->GetOuterParam());
608 }
2d91aa68 609
33791895 610 return track;
611}
612
613void AliHLTEveHLT::SetUpTrackPropagator(TEveTrackPropagator* trkProp, Float_t magF, Float_t maxR) {
614 //See header file for documentation
615
616 if (fTrueField) {
617 trkProp->SetMagFieldObj(new AliEveMagField);
618
619 } else {
620 trkProp->SetMagField(magF);
621 }
622
623 if (fUseRkStepper) {
624 trkProp->SetStepper(TEveTrackPropagator::kRungeKutta);
625 }
626
627 trkProp->SetMaxR(maxR);
628}
629
630
631void AliHLTEveHLT::AddTrackParamToTrack(AliEveTrack* track, const AliExternalTrackParam* tp) {
632 //See header file for documentation
633
634 if (tp == 0)
635 return;
636
637 Double_t pbuf[3], vbuf[3];
638 tp->GetXYZ(vbuf);
639 tp->GetPxPyPz(pbuf);
640
641 TEvePathMark pm(TEvePathMark::kReference);
642 pm.fV.Set(vbuf);
643 pm.fP.Set(pbuf);
644 track->AddPathMark(pm);
645}
646
647
648
649TString AliHLTEveHLT::CreateTrackTitle(AliESDtrack* t) {
650 // Add additional track parameters as a path-mark to track.
651
652 TString s;
653
654 Int_t label = t->GetLabel(), index = t->GetID();
655 TString idx(index == kMinInt ? "<undef>" : Form("%d", index));
656 TString lbl(label == kMinInt ? "<undef>" : Form("%d", label));
657
658 Double_t p[3], v[3];
659 t->GetXYZ(v);
660 t->GetPxPyPz(p);
661 Double_t pt = t->Pt();
662 Double_t ptsig = TMath::Sqrt(t->GetSigma1Pt2());
663 Double_t ptsq = pt*pt;
664 Double_t ptm = pt / (1.0 + pt*ptsig);
665 Double_t ptM = pt / (1.0 - pt*ptsig);
666
667 s = Form("Index=%s, Label=%s\nChg=%d, Pdg=%d\n"
668 "pT = %.3f + %.3f - %.3f [%.3f]\n"
669 "P = (%.3f, %.3f, %.3f)\n"
670 "V = (%.3f, %.3f, %.3f)\n",
671 idx.Data(), lbl.Data(), t->Charge(), 0,
672 pt, ptM - pt, pt - ptm, ptsig*ptsq,
673 p[0], p[1], p[2],
674 v[0], v[1], v[2]);
675
676 Int_t o;
677 s += "Det (in,out,refit,pid):\n";
678 o = AliESDtrack::kITSin;
679 s += Form("ITS (%d,%d,%d,%d) ", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
680 o = AliESDtrack::kTPCin;
681 s += Form("TPC(%d,%d,%d,%d)\n", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
682 o = AliESDtrack::kTRDin;
683 s += Form("TRD(%d,%d,%d,%d) ", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
684 o = AliESDtrack::kTOFin;
685 s += Form("TOF(%d,%d,%d,%d)\n", t->IsOn(o), t->IsOn(o<<1), t->IsOn(o<<2), t->IsOn(o<<3));
686 o = AliESDtrack::kHMPIDout;
687 s += Form("HMPID(out=%d,pid=%d)\n", t->IsOn(o), t->IsOn(o<<1));
688 s += Form("ESD pid=%d", t->IsOn(AliESDtrack::kESDpid));
689
690 if (t->IsOn(AliESDtrack::kESDpid))
691 {
692 Double_t pid[5];
693 t->GetESDpid(pid);
694 s += Form("\n[%.2f %.2f %.2f %.2f %.2f]", pid[0], pid[1], pid[2], pid[3], pid[4]);
695 }
696
697 return s;
698}
699
fd2adb88 700