fixed the tainted variables
[u/mrichter/AliRoot.git] / PMD / AliPMDTrackRefRead.C
1 //
2 // Reading macro for the track reference tree
3 // minimum stuff
4 //
5 void AliPMDTrackRefRead()
6 {
7
8   AliRunLoader*    runLoader =  AliRunLoader::Open("galice.root");
9
10   if (!runLoader)
11     { 
12       printf("Could not open galice.root");
13     }
14
15   runLoader->LoadTrackRefs("READ");
16
17   Int_t nEvents = runLoader->GetNumberOfEvents();
18
19   printf("Total number of Events = %d\n",nEvents);
20
21
22   FILE *fpw = fopen("trackref.dat","w");
23
24
25
26   TClonesArray *arrayTR = 0x0;
27   arrayTR = new TClonesArray("AliTrackReference");
28
29   TTree *treeTR;
30
31
32   AliStack *stack = runLoader->Stack();
33
34
35   // Starts the event loop
36
37   for (Int_t iev = 0; iev < nEvents; iev++)
38     {
39       runLoader->GetEvent(iev);
40       treeTR = runLoader->TreeTR();
41       treeTR->SetBranchAddress("TrackReferences",  &arrayTR);
42       
43       Int_t nEntries = treeTR->GetEntries();
44       
45       printf("Total number of Entries = %d \n",nEntries);
46       
47       for (Int_t i = 0; i < nEntries; i++)
48         {
49           Int_t trRead  = treeTR->GetEntry(i);
50           if (trRead <= 0) continue;
51           Int_t nTrackRefs = arrayTR->GetEntries();
52           
53           for (Int_t j = 0; j < nTrackRefs; j++)
54             {
55               AliTrackReference* trackRef = static_cast<AliTrackReference*>(arrayTR->At(j));
56               if (!trackRef) continue;
57               
58               
59               if (trackRef->DetectorId() != AliTrackReference::kPMD) continue;
60               
61               //TParticle* mpart = 0;
62               
63               //Int_t trackno = trackRef->GetTrack();
64               Int_t trackno = trackRef->Label();
65               
66               //printf("Track number = %d \n",trackno);
67               if (iev == 0) fprintf(fpw,"Track number = %d \n",trackno);
68               
69               
70               //mpart = stack->Particle(trackno);
71
72               //Int_t   MPID   = mpart->GetFirstMother();
73               //Int_t   M2PID  = mpart->GetSecondMother();
74
75             }
76           
77         }
78     }
79   
80   
81   
82   
83   
84 }