In This Package:

RollingGain.h

Go to the documentation of this file.

/*****************************************

   Rolling gain calibration and dark noise monitor
   Details: DocDB 5323, 5361

   Author: Qing He 
   Date: Sep., 2010

   Oct. 7, 2010 
   Zhe Wang
       Update to accommodate all PMT detectors 
       Fit range fine tuning
       Add all fit and error status
   Oct. 12, 2010
       First stable version. Result for AD2 dry run is in DocDB 5527
   Jan. 26, 2011
       Update to work with DaqFormat

   Mar 08, 2011
   Qing He 
       Update to fit (Adc-preAdc) instead of fit Adc, preAdc separately. 

*****************************************/

#ifndef _ROLLINGGAIN_H_
#define _ROLLINGGAIN_H_

#include "GaudiAlg/GaudiAlgorithm.h"
#include "Conventions/Detectors.h"
#include "Event/ReadoutHeader.h"
#include "Event/CalibReadoutHeader.h"
#include "DataSvc/ICableSvc.h"
#include "DataSvc/IRunDataSvc.h"
#include "TH1F.h"
#include "TF1.h"
#include "TFile.h"
#include "TTree.h"
#include <map>

using namespace std;
using namespace DayaBay;

class RollingGain : public GaudiAlgorithm 
{

 public:
  
  RollingGain(const std::string& name, ISvcLocator* pSvcLocator);
  virtual ~RollingGain();
  
  virtual StatusCode initialize();
  virtual StatusCode execute();
  virtual StatusCode finalize();
  
  static double poisson(double u, int n);
  static double gn(double x,double n, double q0, double sigma0,
                   double q1, double sigma1);
  static double NIMmodel(double* xpar,double* par);
  static double Pedestal(double* xpar,double* par);
  void Fit();

  int InitPmtProp(const ServiceMode& svc);
  int ResetPmtProp();

 private:
  ICableSvc* m_cableSvc;
  IRunDataSvc* m_runDataSvc;

  TF1* f1;
  TF1* f2;
  TFile* m_rootfile;
  TTree* m_tree;
  FILE* m_masterfile;
  string m_fileName;

  unsigned long m_NTrigger;

  struct PmtProp {
    TH1F  *h_Adc;
    TH1F  *h_Tdc;
    TH1F  *h_PreAdc;
    double Pedestal;
    double PedestalErr;
    double DarkRate;
    double DarkRateErr;
    double Gain;
    double GainErr;
    double Sigma;
    double SigmaErr;
    int    Status;
    int    NHits;
    int    FitStatus;
    double Chi2;
    string PmtLabel;
    int    PmtId;
  };

  map<int/*PmtId*/,PmtProp> m_PmtPropMap;

  double m_StartTime;
  double m_EndTime;

  int m_Fix;
  int m_NStop;
  int m_FirstTrigSecond;
  int m_FirstTrigNanoSec;
  int m_LastTrigSecond;
  int m_LastTrigNanoSec;
  int m_startrun;
  int m_currentrun;

  int m_site;                   // Site determined by the first trigger 
  int m_simFlag;                // SimFlag determined by the first trigger
  double m_currentTime;         // Time of current trigger
  double m_prevTrigTime;        // Time of previous trigger
  double m_lastFitTime;         // Time when last fit occurs
  double m_fitPeriod;           // Fitting period in second
};

#endif  // _ROLLINGGAIN_H_

---