// -*- mode: c++ -*-
//
// $Id: BrMultCentCalModule.h,v 1.2 2001/11/21 19:13:35 cholm Exp $
// $Author: cholm $
// $Date: 2001/11/21 19:13:35 $
// $Copyright: (C) 2001 BRAHMS Collaboration <brahmlib@rhic.bnl.gov>
//
#ifndef BRAT_BrMultCentCalModule
#define BRAT_BrMultCentCalModule
#ifndef BRAT_BrModule
#include "BrModule.h"
#endif
#ifndef BRAT_BrEventNode
#include "BrEventNode.h"
#endif
#ifndef ROOT_TH2
#include <TH2.h>
#endif
#ifndef ROOT_TH1
#include <TH1.h>
#endif
#ifndef ROOT_TGraph
#include <TGraph.h>
#endif
#ifndef ROOT_TString
#include <TString.h>
#endif
#ifndef ROOT_TDirectory
#include "TDirectory.h"
#endif

class BrMultCentCalModule : public BrModule
{
protected:
  TString  fOutputName;       // 

  TDirectory* fHistDir;          // Histogram directory  
  TH2F*       fVtxVsEHisto;      // Two dim scatter plot to fill. 
  TH1D**      fCutHistos;        // 
  TH1D*       fVtxHisto;         // 

  Float_t     fMax;              // Max Energy for single histograms
  Int_t       fSingleNBins;      // Number of single bins  
  Float_t     fVertexMax;        // Max absolut Z for vertex. 
  Int_t       fVertexNBins;      // Number of vertex bins.

  Int_t       fFunctionOrder;    // order of polynomials
  Int_t       fNCuts;            // number of cut functions 
  Ssiz_t      fSCuts;            // Size of cut array
  Float_t*    fCuts;             // Cut values for functions. 

  Float_t FindVertex(BrEventNode* node);
public:
  BrMultCentCalModule();
  BrMultCentCalModule(const Char_t* name, const Char_t* title);
  virtual ~BrMultCentCalModule();

  virtual void SetOutputName(const Char_t* name=""); // *MENU*
  virtual void SetMax(Float_t max=4);                // *MENU*
  virtual void SetSingleNBins(Int_t n=100);          // *MENU*  
  virtual void SetVertexMax(Float_t z=50);           // *MENU*
  virtual void SetVertexNBins(Int_t n=8);            // *MENU*
  virtual void SetFunctionOrder(Int_t order=3);      // *MENU*
  virtual void AddCut(Float_t cut);                  // *MENU*

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

  ClassDef(BrMultCentCalModule,0) // DOCUMENT ME
};

#endif
//____________________________________________________________________
//
// $Log: BrMultCentCalModule.h,v $
// Revision 1.2  2001/11/21 19:13:35  cholm
// Some fixes, and removed old unused and redundant code.
//
// Revision 1.1  2001/09/04 16:53:24  cholm
// New centrality calibration module BrMultCentCalModule, based on cuts
// in energy deposited in the single detectors SMA and TMA, as outline in
// BAN 34.
//
//