// -*- mode: c++ -*- 
//
// $Id: BrChkvRdoModule.h,v 1.4 2001/10/17 16:33:58 ekman Exp $
// $Author: ekman $
// $Date: 2001/10/17 16:33:58 $
// $Copyright: (C) 2001 BRAHMS Collaboration <brahmlib@rhic.bnl.gov>
//
//
// Chkv reconstruction Module
//
///////////////////////////////////////////////////////////////////////
//                                                                   //
//    BrChkvRdoModule                                                //
//                                                                   //
//    BRAHMS Chkv analysis class                                     //
//                                                                   //
//    Take digit data and give rdo object with hits -                //
//    a hit is defined as a adc signal (in a tube) above pedestal    //     
//                                                                   //
//                                                                   //
///////////////////////////////////////////////////////////////////////
#ifndef BRAT_BrChkvRdoModule
#define BRAT_BrChkvRdoModule
#ifndef ROOT_TNtuple
#include "TNtuple.h"
#endif
#ifndef ROOT_TH1
#include "TH1.h"
#endif
#ifndef ROOT_TProfile	
#include "TProfile.h"
#endif
#ifndef ROOT_TH2
#include "TH2.h"
#endif
#ifndef BRAT_BrModule
#include "BrModule.h"
#endif
#ifndef BRAT_BrEventNode
#include <BrEventNode.h>
#endif
#ifndef BRAT_BrChkvCalibration
#include <BrChkvCalibration.h>
#endif
#ifndef BRAT_BrDetectorVolume
#include <BrDetectorVolume.h>
#endif
#ifndef BRAT_BrMagnetVolume
#include "BrMagnetVolume.h"
#endif
#ifndef BRAT_BrChkvParameters
#include "BrChkvParameters.h"
#endif
#ifndef BRAT_BrSpectrometerTracks
#include "BrSpectrometerTracks.h"
#endif

class BrChkvRdoModule : public BrModule 
{

private:
  
  // needed stuff
  BrDetectorVolume*  fChkvVol;
  BrChkvParameters*  fChkvParams;
  BrChkvCalibration* fCalParams;

  Double_t fThreshold;         // cut on cal adc for hit selection
                               // default = 0 for C1, 0 for RICH 
  Double_t fNoPedWidthCut;     // requires adc > ped+ n* pedW 
                               // default = 0.5 for C1 and 2 for RICH

  // histograms
  TH1F*     fNHits;
  TH2F*     fAdcVsTube;
  TH2F*     fCalAdcVsTube;
  TH2F*     fCalAdcVsTubeAcc;

  TH1F**    fTubeAdc;
  TH1F**    fTubeCalAdc;


public:
  //  constructors and destructors
  //
  BrChkvRdoModule();
  BrChkvRdoModule(Text_t *Name, Char_t *Title);
  virtual ~BrChkvRdoModule();

  //
  //  Methods
  //

  virtual void Init();
  virtual void Begin();
  virtual void DefineHistograms();
  virtual void Event(BrEventNode*, BrEventNode*);

  void SetThreshold (Double_t e)    { fThreshold = e; }
  void SetNoPedWidthCut(Double_t n)  { fNoPedWidthCut = n; }
  void Print(Option_t* option="B") const;

  ClassDef(BrChkvRdoModule, 0) // BRAHMS Chkv hit recontruction module


};
    
#endif
    
// $Log: BrChkvRdoModule.h,v $
// Revision 1.4  2001/10/17 16:33:58  ekman
// Introduced histograms for each tube and deleted the ntuple.
//
// Revision 1.3  2001/10/17 00:23:00  ekman
// Introduced a new cut method to select hits - now one can select by cutting
// relative to cal adc = 1 (threshold) or requiring adc > ped + n* pedW. Updated the histograms.
//
// Revision 1.2  2001/08/29 10:46:41  ekman
// Changed energythreshold to threshold.
//
// Revision 1.1  2001/08/09 16:01:06  ekman
// Added the classes BrChkvRdo and BrChkvRcoModule.
//