CRAB added

[u/mrichter/AliRoot.git] / HBTAN / AliHBTReaderITSv1.cxx

diff --git a/HBTAN/AliHBTReaderITSv1.cxx b/HBTAN/AliHBTReaderITSv1.cxx
index 11dcfabf56d257cf937eef12396f392834d3945b..82c3e3360f106a0df221df2905c694fa28f1c704 100644 (file)
--- a/HBTAN/AliHBTReaderITSv1.cxx
+++ b/HBTAN/AliHBTReaderITSv1.cxx
@@ -20,6 +20,7 @@
 #include <AliMagF.h>
 #include <AliKalmanTrack.h>
 #include <AliITSIOTrack.h>
+#include "AliMC.h"
 
 ClassImp(AliHBTReaderITSv1)
 /********************************************************************/
@@ -170,7 +171,7 @@ Int_t AliHBTReaderITSv1::Read(AliHBTRun* particles, AliHBTRun *tracks)
       
       gAliceFile->cd();
       gAlice->GetEvent(currentEvent);
-      gAlice->Particles();
+      gAlice->GetMCApp()->Particles();
 
       Int_t nentr=(Int_t)tracktree->GetEntries();
       
@@ -188,7 +189,7 @@ Int_t AliHBTReaderITSv1::Read(AliHBTRun* particles, AliHBTRun *tracks)
            continue;
          }
 
-        TParticle *p = (TParticle*)gAlice->Particle(label);
+        TParticle *p = (TParticle*)gAlice->GetMCApp()->Particle(label);
         if(!p)
          {
            Warning("Read","Can not get particle with label &d",label);
@@ -197,7 +198,7 @@ Int_t AliHBTReaderITSv1::Read(AliHBTRun* particles, AliHBTRun *tracks)
         if(Pass(p->GetPdgCode())) continue; //check if we are intersted with particles of this type
                                            //if not take next partilce
 
-        AliHBTParticle* part = new AliHBTParticle(*p);
+        AliHBTParticle* part = new AliHBTParticle(*p,i);
         if(Pass(part)) { delete part; continue;}//check if meets all criteria of any of our cuts
                                                 //if it does not delete it and take next good track
         
@@ -211,7 +212,7 @@ Int_t AliHBTReaderITSv1::Read(AliHBTRun* particles, AliHBTRun *tracks)
         Double_t y= iotrack->GetY();
         Double_t z= iotrack->GetZ();
         
-        AliHBTParticle* track = new AliHBTParticle(p->GetPdgCode(), px, py , pz, tEtot, x, y, z, 0.);
+        AliHBTParticle* track = new AliHBTParticle(p->GetPdgCode(), i, px, py , pz, tEtot, x, y, z, 0.);
         if(Pass(track)) { delete  track;continue;}//check if meets all criteria of any of our cuts
                                                   //if it does not delete it and take next good track