AliClusterContainer.h

Go to the documentation of this file.

#ifndef AliClusterContainer_H
#define AliClusterContainer_H
/* Copyright(c) 1998-1999, ALICE Experiment at CERN, All rights reserved. *
 * See cxx source for full Copyright notice                               */

class TLorentzVector;

class AliVEvent;

#include <iterator>
#include <TArrayI.h>
#include <AliVCluster.h>

#include "AliEmcalContainer.h"

class AliClusterContainer : public AliEmcalContainer {
 public:
  typedef enum AliVCluster::VCluUserDefEnergy_t VCluUserDefEnergy_t;
  
  class accept_iterator : public std::iterator<std::bidirectional_iterator_tag,
                                               AliVCluster,std::ptrdiff_t,
                                               AliVCluster **, AliVCluster *&>
  {
  public:
    accept_iterator(const AliClusterContainer *cont, int startpos, bool forward  = true);
    accept_iterator(const accept_iterator &other);
    ~accept_iterator() {}
    accept_iterator &operator=(const accept_iterator &other);
    bool operator!=(const accept_iterator &other) const;

    accept_iterator &operator++();
    accept_iterator operator++(int);
    accept_iterator &operator--();
    accept_iterator operator--(int);

    AliVCluster *operator*() const;

  private:
    accept_iterator();

    const AliClusterContainer       *fkContainer;       
    TArrayI                          fAcceptIndices;    
    int                              fCurrentPos;       
    bool                             fForward;          
  };

  class all_iterator : public std::iterator<std::bidirectional_iterator_tag,
                                            AliVCluster,std::ptrdiff_t,
                                            AliVCluster **, AliVCluster *&>
  {
  public:
    all_iterator(const AliClusterContainer *cont, int startpos, bool forward  = true);
    all_iterator(const all_iterator &other);
    ~all_iterator() {}
    all_iterator &operator=(const all_iterator &other);
    bool operator!=(const all_iterator &other) const;

    all_iterator &operator++();
    all_iterator operator++(int);
    all_iterator &operator--();
    all_iterator operator--(int);

    AliVCluster *operator*() const;

  private:
    all_iterator();

    const AliClusterContainer       *fkContainer;       
    int                              fCurrentPos;       
    bool                             fForward;          
  };

  AliClusterContainer();
  AliClusterContainer(const char *name); 
  virtual ~AliClusterContainer(){;}

  virtual Bool_t              AcceptObject(Int_t i, UInt_t &rejectionReason) const              { return AcceptCluster(i, rejectionReason);}
  virtual Bool_t              AcceptObject(const TObject* obj, UInt_t &rejectionReason) const   { return AcceptCluster(dynamic_cast<const AliVCluster*>(obj), rejectionReason);}
  virtual Bool_t              AcceptCluster(Int_t i, UInt_t &rejectionReason)                 const;
  virtual Bool_t              AcceptCluster(const AliVCluster* vp, UInt_t &rejectionReason)   const;
  virtual Bool_t              ApplyClusterCuts(const AliVCluster* clus, UInt_t &rejectionReason) const;
  AliVCluster                *GetAcceptCluster(Int_t i)              const;
  AliVCluster                *GetAcceptClusterWithLabel(Int_t lab)   const;
  void                        SetClusECut(Double_t cut)                    { SetMinE(cut)     ; }
  void                        SetClusPtCut(Double_t cut)                   { SetMinPt(cut)    ; }
  Double_t                    GetClusPtCut()                         const { return GetMinPt(); }
  AliVCluster                *GetCluster(Int_t i)                    const;
  AliVCluster                *GetClusterWithLabel(Int_t lab)         const;
  AliVCluster                *GetLeadingCluster(const char* opt="")       ;
  Bool_t                      GetMomentum(TLorentzVector &mom, const AliVCluster* vc, Double_t mass) const;
  Bool_t                      GetMomentum(TLorentzVector &mom, const AliVCluster* clus) const;
  Bool_t                      GetMomentum(TLorentzVector &mom, Int_t i) const;
  Bool_t                      GetAcceptMomentum(TLorentzVector &mom, Int_t i) const;
  Bool_t                      GetNextMomentum(TLorentzVector &mom);
  Bool_t                      GetNextAcceptMomentum(TLorentzVector &mom);
  AliVCluster                *GetNextAcceptCluster();
  AliVCluster                *GetNextCluster();
  Int_t                       GetNClusters()                         const { return GetNEntries();   }
  Int_t                       GetNAcceptedClusters()                 const;
  void                        SetClusTimeCut(Double_t min, Double_t max)   { fClusTimeCutLow  = min ; fClusTimeCutUp = max ; }
  void                        SetMinMCLabel(Int_t s)                       { fMinMCLabel      = s   ; }
  void                        SetMaxMCLabel(Int_t s)                       { fMaxMCLabel      = s   ; }
  void                        SetMCLabelRange(Int_t min, Int_t max)        { SetMinMCLabel(min)     ; SetMaxMCLabel(max)    ; }
  void                        SetExoticCut(Bool_t e)                       { fExoticCut       = e   ; }

  void                        SetClusUserDefEnergyCut(Int_t t, Double_t cut);
  Double_t                    GetClusUserDefEnergyCut(Int_t t) const;

  void                        SetClusNonLinCorrEnergyCut(Double_t cut)                     { SetClusUserDefEnergyCut(AliVCluster::kNonLinCorr, cut); }
  void                        SetClusHadCorrEnergyCut(Double_t cut)                        { SetClusUserDefEnergyCut(AliVCluster::kHadCorr, cut)   ; }
  void                        SetDefaultClusterEnergy(Int_t d)                             { fDefaultClusterEnergy = d                             ; }

  Int_t                       GetDefaultClusterEnergy() const                              { return fDefaultClusterEnergy                          ; }

  const char*                 GetTitle() const;

  accept_iterator             accept_begin() const;
  accept_iterator             accept_end() const;
  accept_iterator             accept_rbegin() const;
  accept_iterator             accept_rend() const;

  all_iterator                begin() const;
  all_iterator                end() const;
  all_iterator                rbegin() const;
  all_iterator                rend() const;

 protected:
  
  Double_t         fClusTimeCutLow;             
  Double_t         fClusTimeCutUp;              
  Bool_t           fExoticCut;                  
  Double_t         fUserDefEnergyCut[AliVCluster::kLastUserDefEnergy+1]; 
  Int_t            fDefaultClusterEnergy;       

 private:
  AliClusterContainer(const AliClusterContainer& obj); // copy constructor
  AliClusterContainer& operator=(const AliClusterContainer& other); // assignment

  ClassDef(AliClusterContainer,5);
};

int TestClusterContainerIterator(const AliClusterContainer *const cont, int iteratorType = 0, bool verbose = false);

#endif