]>
Commit | Line | Data |
---|---|---|
ba978640 | 1 | #if !defined(__CINT__) || defined(__MAKECINT__) |
2 | #include <TBranch.h> | |
3 | #include <TFile.h> | |
4 | #include <TSystem.h> | |
5 | #include <TTree.h> | |
6 | #include <TROOT.h> | |
7 | #include <TEveManager.h> | |
8 | ||
6c49a8e1 | 9 | #include <AliCDBManager.h> |
10 | #include <AliTRDarrayADC.h> | |
11 | #include <AliTRDReconstructor.h> | |
12 | #include <AliTRDrecoParam.h> | |
13 | #include <AliTRDtrackV1.h> | |
ba978640 | 14 | #include <PWGPP/TRD/info/AliTRDtrackInfo.h> |
6c49a8e1 | 15 | #include <AliEveEventManager.h> |
16 | #include <AliEveTRDTrackList.h> | |
17 | #include <AliEveTRDData.h> | |
ba978640 | 18 | #endif |
19 | ||
e0d63214 | 20 | void trd_qaRec() |
21 | { | |
230b7aff | 22 | gSystem->Load("libANALYSIS"); |
23 | gSystem->Load("libTender"); | |
24 | gSystem->Load("libPWGPP"); | |
e0d63214 | 25 | |
26 | if(!TFile::Open("TRD.DebugInfoGen.root")){ | |
27 | printf("No debug file for InfoGen task.\n"); | |
28 | return; | |
29 | } | |
30 | TTree *t = (TTree*)gFile->Get("trackInfo"); | |
e69303af | 31 | AliTRDtrackInfo *fTrackInfo = 0x0; |
e0d63214 | 32 | t->SetBranchAddress("TrackInfo.", &fTrackInfo); |
33 | gROOT->cd(); | |
34 | ||
35 | AliCDBManager *cdb = AliCDBManager::Instance(); | |
162637e4 | 36 | cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB"); |
e0d63214 | 37 | cdb->SetRun(0); |
38 | AliEveEventManager::AssertGeometry(); | |
8661a211 | 39 | AliEveEventManager::AssertMagField(); |
12896fe5 | 40 | |
41 | AliTRDReconstructor *reco = new AliTRDReconstructor(); | |
42 | reco->SetRecoParam(AliTRDrecoParam::GetLowFluxParam()); | |
43 | ||
e0d63214 | 44 | AliEveTRDTrackList *tracks = new AliEveTRDTrackList("TRD QA Tracks"); |
45 | ||
46 | ||
47 | ||
e69303af | 48 | AliTRDtrackV1 *fTrack = 0x0, *track = 0x0; |
e0d63214 | 49 | for (Int_t it=0; it<t->GetEntries(); it++){ |
50 | if(!t->GetEntry(it)) continue; | |
51 | if(!fTrackInfo) continue; | |
e69303af | 52 | if(!(fTrack = fTrackInfo->GetTrack())) continue; |
e0d63214 | 53 | |
e69303af | 54 | track = new AliTRDtrackV1(*fTrack); |
55 | track->SetOwner(); | |
56 | track->SetReconstructor(reco); | |
57 | tracks->AddElement(new AliEveTRDTrack(track)); | |
58 | printf("Trk[%3d] ESD[%d] Ncls[%d]\n", it, fTrackInfo->GetESDinfo()->GetId(), fTrack->GetNumberOfClusters()); | |
59 | if(it>= 100) break; | |
e0d63214 | 60 | } |
61 | gEve->AddElement(tracks); | |
62 | gEve->Redraw3D(); | |
f7a1cc68 | 63 | } |