Implemented restricting event subrange to read

diff --git a/HBTAN/AliHBTReader.cxx b/HBTAN/AliHBTReader.cxx
index fd4d25d213df324df7ca22f7258b070039f8a0d7..b6edaa3623744e480690238e0a69ec5575129463 100644 (file)
--- a/HBTAN/AliHBTReader.cxx
+++ b/HBTAN/AliHBTReader.cxx
@@ -26,7 +26,10 @@ AliHBTReader::AliHBTReader():
  fParticles(0x0),
  fTracks(0x0),
  fIsRead(kFALSE),
- fBufferEvents(kFALSE)
+ fBufferEvents(kFALSE),
+ fBlend(kFALSE),
+ fFirst(0),
+ fLast(0)
 {
 //constructor
 }
@@ -43,7 +46,10 @@ AliHBTReader::AliHBTReader(TObjArray* dirs):
  fParticles(0x0),
  fTracks(0x0),
  fIsRead(kFALSE),
- fBufferEvents(kFALSE)
+ fBufferEvents(kFALSE),
+ fBlend(kFALSE),
+ fFirst(0),
+ fLast(0)
 {
 //ctor with array of directories to read as parameter
 }
@@ -65,17 +71,22 @@ AliHBTReader::~AliHBTReader()
 Int_t AliHBTReader::Next()
 {
 //moves to next event
-  if ( ReadNext() == kTRUE)
-     return kTRUE;
+  if ((fNEventsRead > fLast) && (fLast > 0) )
+   
+  do
+   {
+    if ( ReadNext() == kTRUE)
+      return kTRUE;
+   }while (fNEventsRead < fFirst);
   
   if (fBlend) Blend();
   
   if (fBufferEvents)
    {
      if ( ReadsTracks() && fTracksEvent) 
-       fTracks->SetEvent(fNEventsRead-1,fTracksEvent);
+       fTracks->SetEvent(fNEventsRead-1-fFirst,fTracksEvent);
      if ( ReadsParticles() && fParticlesEvent)
-       fParticles->SetEvent(fNEventsRead-1,fParticlesEvent);
+       fParticles->SetEvent(fNEventsRead-1-fFirst,fParticlesEvent);
    }
   return kFALSE;
 }

diff --git a/HBTAN/AliHBTReader.h b/HBTAN/AliHBTReader.h
index e986b8d36f12562f4f525a2880313c594b116d90..f9c250dd327fe10eb9db75fcd5270f8641475e15 100644 (file)
--- a/HBTAN/AliHBTReader.h
+++ b/HBTAN/AliHBTReader.h
@@ -44,6 +44,7 @@ class AliHBTReader: public TNamed
     void                 SetEventBuffering(Bool_t flag){fBufferEvents = flag;}
     void          SetBlend(Bool_t flag = kTRUE){fBlend=flag;}
     virtual Int_t GetNumberOfDirs() const {return (fDirs)?fDirs->GetEntries():0;}
+    void          ReadEventsFromTo(Int_t first,Int_t last){fFirst = first; fLast = last;}
   protected:
     
     TObjArray*    fCuts;//array with particle cuts
@@ -64,6 +65,10 @@ class AliHBTReader: public TNamed
     Bool_t        fBufferEvents;//flag indicating if the data should be bufferred
     
     Bool_t        fBlend;// flag indicating if randomly change positions of the particles after reading
+    
+    Int_t         fFirst;//first event to return (all are before are skipped)
+    Int_t         fLast;//last
+    
     virtual Int_t ReadNext() = 0; //this methods reads next event and put result in fTracksEvent and/or fParticlesEvent
     Bool_t Pass(AliHBTParticle*);
     Bool_t Pass(Int_t pid);