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

//#############################################################
//#                                                           # 
//#         Class PairAnalysisVarCuts                           #
//#         Provide cuts for all variables handled in         #
//#           PairAnalysisVarManager                            #
//#                                                           #
//#  Authors:                                                 #
//#   Julian    Book,     Uni Ffm / Julian.Book@cern.ch       #
//#                                                           #
//#############################################################

#include <Rtypes.h>
#include <TBits.h>

#include "AnalysisCuts.h"
#include "PairAnalysisVarManager.h"

class THnBase;
class PairAnalysisVarCuts : public AnalysisCuts {
public:
  // Whether all cut criteria have to be fulfilled of just any
  enum CutType { kAll=0, kAny };

  PairAnalysisVarCuts();
  PairAnalysisVarCuts(const char* name, const char* title);
  virtual ~PairAnalysisVarCuts();
  //TODO: make copy constructor and assignment operator public
  void AddCut(PairAnalysisVarManager::ValueTypes type, Double_t min, Double_t max, Bool_t excludeRange=kFALSE);
  void AddCut(const char *formula, Double_t min, Double_t max, Bool_t excludeRange=kFALSE);
  void AddCut(PairAnalysisVarManager::ValueTypes type, Double_t value, Bool_t excludeRange=kFALSE);
  void AddBitCut(PairAnalysisVarManager::ValueTypes type, UInt_t bit, Bool_t excludeRange=kFALSE);
  void AddCut(PairAnalysisVarManager::ValueTypes type, Double_t min, THnBase * const max,  Bool_t excludeRange=kFALSE);

  // setters
  void    SetCutOnMCtruth(Bool_t mc=kTRUE) { fCutOnMCtruth=mc; }
  void    SetCutType(CutType type)         { fCutType=type;    }

  // getters
  Bool_t  GetCutOnMCtruth() const { return fCutOnMCtruth; }
  CutType GetCutType()      const { return fCutType;      }

  Int_t GetNCuts() { return fNActiveCuts; }

  //
  //Analysis cuts interface
  //
  virtual Bool_t IsSelected(TObject* track);
  virtual Bool_t IsSelected(TList*   /* list */ ) {return kFALSE;}

  //   virtual Bool_t IsSelected(TObject* track, TObject */*event*/=0);
  //   virtual Long64_t Merge(TCollection* /* list */)      { return 0; }

  //
  // Cut information
  //
  virtual UInt_t GetSelectedCutsMask() const { return fSelectedCutsMask; }

  virtual void Print(const Option_t* option = "") const;


 private:

  TBits     *fUsedVars;            // list of used variables
  UShort_t  fActiveCuts[PairAnalysisVarManager::kNMaxValuesMC];       // list of activated cuts
  UShort_t  fNActiveCuts;                      // number of acive cuts
  UInt_t    fActiveCutsMask;                   // mask of active cuts

  UInt_t   fSelectedCutsMask;                 // Maks of selected cuts, is available after calling IsSelected

  Bool_t   fCutOnMCtruth;                     // whether to cut on the MC truth of the particle

  CutType  fCutType;                          // type of the cut: any, all

  Double_t fCutMin[PairAnalysisVarManager::kNMaxValuesMC];           // minimum values for the cuts
  Double_t fCutMax[PairAnalysisVarManager::kNMaxValuesMC];           // maximum values for the cuts
  Bool_t fCutExclude[PairAnalysisVarManager::kNMaxValuesMC];         // inverse cut logic?
  Bool_t fBitCut[PairAnalysisVarManager::kNMaxValuesMC];             // bit cut
  TObject    *fUpperCut[PairAnalysisVarManager::kNMaxValuesMC];          // use object as upper cut

  PairAnalysisVarCuts(const PairAnalysisVarCuts &c);
  PairAnalysisVarCuts &operator=(const PairAnalysisVarCuts &c);

  ClassDef(PairAnalysisVarCuts,8)         //Cut class providing cuts to all infomation available for the AliVParticle interface
};


//
//Inline functions
//
inline void PairAnalysisVarCuts::AddCut(PairAnalysisVarManager::ValueTypes type, Double_t value, Bool_t excludeRange)
{
  //
  // Set cut in a small delta around value
  //
  const Double_t kDelta=1e-20;
  AddCut(type,value-kDelta,value+kDelta, excludeRange);
}

#endif