StarParticleStack.h

#ifndef __StarParticleStack_h__
#define __StarParticleStack_h__

#include "TObjArray.h"
#include "StarCallf77.h"
#include "TVirtualMCStack.h"
#include "TDataSet.h"
#include "TMCProcess.h"
#include "TParticle.h"
#include "TClonesArray.h"
#include <vector>
#include <list>


class StarGenParticle;

class StarParticleStack : public TVirtualMCStack
{
 public:

  StarParticleStack( const Char_t *name = "PrimaryParticleStack" );
  ~StarParticleStack();

  virtual void  PushTrack(Int_t toBeDone, Int_t parent, Int_t pdg,
                          Double_t px, Double_t py, Double_t pz, Double_t e,
                          Double_t vx, Double_t vy, Double_t vz, Double_t tof,
                          Double_t polx, Double_t poly, Double_t polz,
                          TMCProcess mech, Int_t& ntr, Double_t weight,
                          Int_t is);



  virtual TParticle* PopNextTrack(Int_t& itrack);
  
  virtual TParticle* PopPrimaryForTracking(Int_t i);

  virtual void       SetCurrentTrack(Int_t trackNumber);
  
  virtual Int_t      GetNtrack()    const;
  
  virtual Int_t      GetNprimary()  const { return mNumPrimary; }
  
  virtual TParticle* GetCurrentTrack() const;

  virtual Int_t      GetCurrentTrackNumber() const;
  
  virtual Int_t      GetCurrentParentTrackNumber() const;

  virtual TParticle *GetParticle( const Int_t i ) const;

  virtual void Clear( const Option_t *opts="" );

 private:
 protected:

  Int_t               mNumPrimary;
  Int_t               mCurrent;

  Int_t               mArraySize;
  TClonesArray       *mArray; 

  Int_t               mStackSize;
  std::list  <TParticle *> mStack;
  std::list  <Int_t>       mStackIdx;


  ClassDef(StarParticleStack,1);
    
};

#endif