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");
  17   cdb->SetRun(0);
  18   AliEveEventManager::AssertGeometry();
  19   AliMagFMaps *field = new AliMagFMaps("Maps","Maps", 2, 1., 10., AliMagFMaps::k5kG);
  20   AliTracker::SetFieldMap(field, kTRUE);
  21
  22   AliTRDReconstructor *reco = new AliTRDReconstructor();
  23   reco->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
  24
  25   AliEveTRDTrackList *tracks = new AliEveTRDTrackList("TRD QA Tracks");
  26
  27
  28
  29   AliTRDtrackV1 *fTrack = 0x0, *track = 0x0;
  30   for (Int_t it=0; it<t->GetEntries(); it++){
  31     if(!t->GetEntry(it)) continue;
  32     if(!fTrackInfo) continue;
  33     if(!(fTrack = fTrackInfo->GetTrack())) continue;
  34
  35     track = new AliTRDtrackV1(*fTrack);
  36     track->SetOwner();
  37     track->SetReconstructor(reco);
  38     tracks->AddElement(new AliEveTRDTrack(track));
  39     printf("Trk[%3d] ESD[%d] Ncls[%d]\n", it, fTrackInfo->GetESDinfo()->GetId(), fTrack->GetNumberOfClusters());
  40     if(it>= 100) break;
  41   }
  42   gEve->AddElement(tracks);
  43   gEve->Redraw3D();
  44 }