// -*- mode: c++ -*-
//
// $Id: BrC1AdcGainCalModule.h,v 1.2 2001/10/23 20:51:23 ouerdane Exp $
// $Author: ouerdane $
// $Date: 2001/10/23 20:51:23 $
// $Copyright: (C) 2001 BRAHMS Collaboration <brahmlib@rcf.rhic.bnl.gov>
//
#ifndef BRAT_BrC1AdcGainCalModule
#define BRAT_BrC1fAdcGainCalModule

#ifndef ROOT_TF1
#include "TF1.h"
#endif
#ifndef ROOT_TH1
#include "TH1.h"
#endif
#ifndef BRAT_BrEventNode
#include "BrEventNode.h"
#endif
#ifndef BRAT_BrDetectorVolume
#include "BrDetectorVolume.h"
#endif
#ifndef BRAT_BrDetectorTrack
#include "BrDetectorTrack.h"
#endif
#ifndef BRAT_BrDetectorParamsTof
#include "BrDetectorParamsTof.h"
#endif
#ifndef BRAT_BrChkvCalModule
#include "BrChkvCalModule.h"
#endif
#ifndef BRAT_BrVector3D
#include "BrVector3D.h"
#endif

class BrC1AdcGainCalModule : public BrChkvCalModule {

private:
  // members
  BrDetectorVolume* fC1Volume;
  BrDetectorVolume* fT2Volume;
  BrPlane3D         fC1BackPlane;
  
  Float_t fMaxGainLimit;     // max limit for reasonable gain 
  Float_t fMinGainLimit;     // min limit for reasonable gain 
  Float_t fMomentumCutHigh;  // maximum momentum  
  Float_t fMomentumCutLow;   // minimum momentum  
  Float_t fWindowWidth;      // 

  TH1F**  fAdc;
  TH1F**  fAdcCal;
  TH1F*   fSum;
  TH1F*   fSumW;
  TH1F*   fChi2;

  // methods
  void     EvaluateGain(TH1F*, TF1*, Int_t);
  void     MakeAdcCalHist(Int_t tube);

protected:
  virtual void SaveAscii();
  virtual void ReadAscii();

public:

  BrC1AdcGainCalModule();
  BrC1AdcGainCalModule(Char_t *Name, Char_t *Title);
  virtual ~BrC1AdcGainCalModule();
  
  virtual void Init();
  virtual void Begin();
  virtual void End();
  virtual void Finish();

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

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

  BrVector3D PointToC1BackPlane(BrDetectorTrack* t2track);
  void       SetMaxGainLimit(Float_t mx = 400) { fMaxGainLimit = mx;  }
  void       SetMinGainLimit(Float_t mn = 20)  { fMinGainLimit = mn;  }
  void       SetMomentumCut(Float_t cutL = 4, Float_t cutH = 6)  
  { fMomentumCutLow = cutL; fMomentumCutHigh = cutH; }
  void       SetWindowWidth(Float_t w = 3.5) {fWindowWidth = w;}

  ClassDef(BrC1AdcGainCalModule,0) // C1 ADC gain calibration Module
    
};
#endif

//____________________________________________________________________
//
// $Log: BrC1AdcGainCalModule.h,v $
// Revision 1.2  2001/10/23 20:51:23  ouerdane
// Updated modules ala tof or bb, with Set[Save,Commit,Load]Ascii, etc
//
// Revision 1.1  2001/09/03 18:00:10  ekman
// Added BrC1AdcGainCalModule.
//
//