Correct GetPrimary when reading events

diff --git a/STEER/AliRun.cxx b/STEER/AliRun.cxx
index c90c14e9e5756eb6cf62b027b716cc07c5134765..39b670fcf1f296fab53edc3c861becf4d6d86bd8 100644 (file)
--- a/STEER/AliRun.cxx
+++ b/STEER/AliRun.cxx
@@ -15,6 +15,9 @@
 
 /*
 $Log$
+Revision 1.62  2001/04/06 11:12:33  morsch
+Clear fParticles after each event. (Ivana Hrivnacova)
+
 Revision 1.61  2001/03/30 07:04:10  morsch
 Call fGenerator->FinishRun() for final print-outs, cross-section and weight calculations.
 
@@ -1023,7 +1026,7 @@ void AliRun::GetNextTrack(Int_t &mtrack, Int_t &ipart, Float_t *pmom,
 }
 
 //_____________________________________________________________________________
-Int_t AliRun::GetPrimary(Int_t track) const
+Int_t AliRun::GetPrimary(Int_t track) 
 {
   //
   // return number of primary that has generated track
@@ -1035,6 +1038,7 @@ Int_t AliRun::GetPrimary(Int_t track) const
   while (1) {
     current=parent;
     part = (TParticle *)fParticleMap->At(current);
+    if(!part) part = Particle(current);
     parent=part->GetFirstMother();
     if(parent<0) return current;
   }

diff --git a/STEER/AliRun.h b/STEER/AliRun.h
index dd43ade7ee7170332d6bfd05a81891e46b40cb4b..e3807990b3ac650488c317869acf2550fea096e3 100644 (file)
--- a/STEER/AliRun.h
+++ b/STEER/AliRun.h
@@ -87,7 +87,7 @@ public:
                               Float_t &e, Float_t *vpos, Float_t *polar, 
                               Float_t &tof);
    Int_t          GetNtrack() const {return fNtrack;}
-   virtual  Int_t GetPrimary(Int_t track) const;
+   virtual  Int_t GetPrimary(Int_t track);
    virtual  void  Hits2Digits(const char *detector=0); 
    virtual  void  Hits2SDigits(const char *detector=0)   {Tree2Tree("S",detector);}
    virtual  void  SDigits2Digits(const char *detector=0) {Tree2Tree("D",detector);}