EVE/alice-macros/trd_qaRec.C

   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
   9 #include <AliCDBManager.h>
  10 #include <AliTRDarrayADC.h>
  11 #include <AliTRDReconstructor.h>
  12 #include <AliTRDrecoParam.h>
  13 #include <AliTRDtrackV1.h>
  14 #include <PWGPP/TRD/info/AliTRDtrackInfo.h>
  15 #include <AliEveEventManager.h>
  16 #include <AliEveTRDTrackList.h>
  17 #include <AliEveTRDData.h>
  18 #endif
  19
  20 void trd_qaRec()
  21 {
  22   gSystem->Load("libANALYSIS");
  23   gSystem->Load("libTender");
  24   gSystem->Load("libPWGPP");
  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");
  31   AliTRDtrackInfo *fTrackInfo = 0x0;
  32   t->SetBranchAddress("TrackInfo.", &fTrackInfo);
  33   gROOT->cd();
  34
  35   AliCDBManager *cdb = AliCDBManager::Instance();
  36   cdb->SetDefaultStorage("local://$ALICE_ROOT/OCDB");
  37   cdb->SetRun(0);
  38   AliEveEventManager::AssertGeometry();
  39   AliEveEventManager::AssertMagField();
  40
  41   AliTRDReconstructor *reco = new AliTRDReconstructor();
  42   reco->SetRecoParam(AliTRDrecoParam::GetLowFluxParam());
  43
  44   AliEveTRDTrackList *tracks = new AliEveTRDTrackList("TRD QA Tracks");
  45
  46
  47
  48   AliTRDtrackV1 *fTrack = 0x0, *track = 0x0;
  49   for (Int_t it=0; it<t->GetEntries(); it++){
  50     if(!t->GetEntry(it)) continue;
  51     if(!fTrackInfo) continue;
  52     if(!(fTrack = fTrackInfo->GetTrack())) continue;
  53
  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;
  60   }
  61   gEve->AddElement(tracks);
  62   gEve->Redraw3D();
  63 }