EVE/alice-macros/trd_qaRec.C

   1 void trd_qaRec()
   2 {
   3   gSystem->Load("libANALYSIS.so");
   4   gSystem->Load("libTRDqaRec.so");
   5
   6   if(!TFile::Open("TRD.DebugInfoGen.root")){
   7     printf("No debug file for InfoGen task.\n");
   8     return;
   9   }
  10   TTree *t = (TTree*)gFile->Get("trackInfo");
  11   AliTRDtrackInfo *fTrackInfo = 0x0;
  12   t->SetBranchAddress("TrackInfo.", &fTrackInfo);
  13   gROOT->cd();
  14
  15   AliCDBManager *cdb = AliCDBManager::Instance();
  16   cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  17   cdb->SetRun(0);
  18   AliEveEventManager::AssertGeometry();
  19   AliEveEventManager::AssertMagField();
  20
  21   AliTRDReconstructor *reco = new AliTRDReconstructor();
  22   reco->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
  23
  24   AliEveTRDTrackList *tracks = new AliEveTRDTrackList("TRD QA Tracks");
  25
  26
  27
  28   AliTRDtrackV1 *fTrack = 0x0, *track = 0x0;
  29   for (Int_t it=0; it<t->GetEntries(); it++){
  30     if(!t->GetEntry(it)) continue;
  31     if(!fTrackInfo) continue;
  32     if(!(fTrack = fTrackInfo->GetTrack())) continue;
  33
  34     track = new AliTRDtrackV1(*fTrack);
  35     track->SetOwner();
  36     track->SetReconstructor(reco);
  37     tracks->AddElement(new AliEveTRDTrack(track));
  38     printf("Trk[%3d] ESD[%d] Ncls[%d]\n", it, fTrackInfo->GetESDinfo()->GetId(), fTrack->GetNumberOfClusters());
  39     if(it>= 100) break;
  40   }
  41   gEve->AddElement(tracks);
  42   gEve->Redraw3D();
  43 }