// -*- mode: c++ -*-
//
// $Id: BrCentFilter.h,v 1.3 2001/11/22 17:37:07 cholm Exp $
//
#ifndef BRAT_BrCentFilter
#define BRAT_BrCentFilter

#ifndef BRAT_BrModule
#include "BrModule.h"
#endif
#ifndef BRAT_BrEventNode
#include "BrEventNode.h"
#endif
#ifndef ROOT_TH1
#include "TH1.h"
#endif

class BrCentFilter : public BrModule
{
private:
  Float_t fTileLowLimit;  // Minimum centrality to except (most central)
  Float_t fTileHighLimit; // Maximum centrality to except (least central)
  Float_t fSiLowLimit;    // Minimum centrality to except (most central)
  Float_t fSiHighLimit;   // Maximum centrality to except (least central)
  Float_t fMultLowLimit;  // Minimum centrality to except (most central)
  Float_t fMultHighLimit; // Maximum centrality to except (least central)

  TH1F*   fTileAccCentHisto; 
  TH1F*   fTileAllCentHisto; 
  TH1F*   fSiAccCentHisto; 
  TH1F*   fSiAllCentHisto; 
  TH1F*   fMultAccCentHisto; 
  TH1F*   fMultAllCentHisto; 

  Int_t  fNoTimesCalled;   // Counter of events 
  Int_t  fNoAcceptedTile;  // Counter of "good" TMA centralities
  Int_t  fNoAcceptedSi;    // Counter of "good" SMA centralities
  Int_t  fNoAcceptedMult;  // Counter of "good" TMA+SMA centralities

  Bool_t CheckTile(BrEventNode* inNode);  
  Bool_t CheckSi(BrEventNode* inNode);  
  Bool_t CheckMult(BrEventNode* inNode);
public:
  BrCentFilter();
  BrCentFilter(const Char_t* name, const Char_t* title); 

  void         SetLowLimit(Float_t low = 0);         //*MENU*
  void         SetHighLimit(Float_t high = 100);     //*MENU*
  void         SetTileLowLimit(Float_t low = 0);     //*MENU* 
  void         SetTileHighLimit(Float_t high = 100); //*MENU*
  void         SetSiLowLimit(Float_t low = 0);       //*MENU*
  void         SetSiHighLimit(Float_t high = 100);   //*MENU*
  void         SetMultLowLimit(Float_t low = 0);     //*MENU*
  void         SetMultHighLimit(Float_t high = 100); //*MENU*

  virtual void DefineHistograms();
  virtual void Init();
  virtual void Event(BrEventNode* inNode, BrEventNode* outNode);
  virtual void Finish();
  virtual void Print(Option_t* option="B") const; //*MENU*

  ClassDef(BrCentFilter,0) // Centrality filter 
};

#endif
//
// $Log: BrCentFilter.h,v $
// Revision 1.3  2001/11/22 17:37:07  cholm
// Updates to use new SDE Centrality methods
//
// Revision 1.2  2001/08/10 14:28:57  cholm
// Changed a few things, so that centrality data not present can be
// ignored, and also added a statistics output (conditional on verbosity
// level ofcourse).
//
// Revision 1.1.1.1  2001/06/21 14:55:14  hagel
// Initial revision of brat2
//
// Revision 1.3  2001/04/19 15:46:22  cholm
// Updated BrCentFilter for Si and Mult cent
//
// Revision 1.2  2001/03/07 12:13:12  cholm
// Changed Info to Print.
//
// Revision 1.1  2001/02/13 09:37:17  cholm
// Added the class BrCentFilter to ease making cuts in centrality.
//
//