KVINDRADB.h

/***************************************************************************
                          kvINDRADB.h  -  description
                             -------------------
    begin                : 9/12 2003
    copyright            : (C) 2003 by J.D. Frankland
    email                : frankland@ganil.fr
 ***************************************************************************/
 
/***************************************************************************
 *                                                                         *
 *   This program is free software; you can redistribute it and/or modify  *
 *   it under the terms of the GNU General Public License as published by  *
 *   the Free Software Foundation; either version 2 of the License, or     *
 *   (at your option) any later version.                                   *
 *                                                                         *
 ***************************************************************************/
 
#ifndef KV_INDRADB_H
#define KV_INDRADB_H
 
#include "KVExpDB.h"
#include "KVDataSet.h"
#include "KVINDRADBRun.h"
#include "KVDBSystem.h"
#include "KVDBTape.h"
#include "KVDetector.h"
#include "Riostream.h"
#include "TEnv.h"
#include "KVINDRARunListReader.h"
#include "KVINDRAPulserDataTree.h"
#include "KVNumberList.h"
 
#define MAX_NUM_RUN_RANGES 100
 
class KVNumberList;
 
class KVINDRADB: public KVExpDB, public KVINDRARunListReader {
 
private:
 
    std::ifstream __ifpeaks;
 
protected:
 
 
    KVDBTable* fChIoPressures;
    KVDBTable* fGains;
    KVDBTable* fTapes;
    KVDBTable* fCsILumCorr;
    KVDBTable* fPedestals;
    KVDBTable* fChanVolt;
    KVDBTable* fVoltMeVChIoSi;
    KVDBTable* fLitEnerCsI;
    KVDBTable* fLitEnerCsIZ1;
    KVDBTable* fAbsentDet;
    KVDBTable* fOoODet;
    KVDBTable* fOoOACQPar;
 
    KVINDRAPulserDataTree* fPulserData;
 
   virtual void ReadGainList();
   virtual void ReadChIoPressures();
   virtual void ReadCsITotalLightGainCorrections();
   virtual void ReadChannelVolt();
   virtual void ReadVoltEnergyChIoSi();
   virtual void ReadLightEnergyCsI(const Char_t*, KVDBTable*);
   virtual void ReadCalibCsI();
   virtual void ReadPedestalList();
 
   virtual void ReadAbsentDetectors();
   virtual void ReadOoODetectors();
   virtual void ReadOoOACQParams();
 
 
   Bool_t OpenCalibrationPeakFile();
   void CloseCalibrationPeakFile();
   std::ifstream& GetPeakFileStream()
   {
      return __ifpeaks;
   };
 
    mutable TF1* fMeanGDir;
   void _init_MeanGDir_func() const;
 
   void init();
 
   void make_minimal_database() override;
   
   void extend_minimal_database() override;
   
public:
   KVINDRADB();
   KVINDRADB(const Char_t* name);
   virtual ~ KVINDRADB();
 
   void cd() override;
 
   void AddTape(KVDBTape* r)
   {
      fTapes->AddRecord(r);
   }
   virtual KVDBTape* GetTape(Int_t tape) const
   {
      return (KVDBTape*) fTapes->GetRecord(tape);
   }
   virtual KVSeqCollection* GetTapes() const
   {
      return fTapes->GetRecords();
   }
 
   void GoodRunLine() override;
   KVINDRADBRun* GetRun(Int_t run) const
   {
      return (KVINDRADBRun*) GetDBRun(run);
   }
   void ReadNewRunList();
 
   KVList* GetCalibrationPeaks(Int_t run, KVDetector* detector = 0,
                               Int_t peak_type = -1, Int_t signal_type = 0,
                               Double_t peak_energy = -1.0);
 
   Double_t GetEventCrossSection(Int_t run, Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetTotalCrossSection(Int_t run, Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetEventCrossSection(Int_t run1, Int_t run2,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetTotalCrossSection(Int_t run1, Int_t run2,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetEventCrossSection(KVNumberList runs,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetTotalCrossSection(KVNumberList runs,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetEventCrossSection(const Char_t* system, Int_t Mult_trigger,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetTotalCrossSection(const Char_t* system, Int_t Mult_trigger,
                                 Double_t Q_apres_cible,
                                 Double_t Coul_par_top = 1.e-10) const;
   Double_t GetTotalCrossSection(TH1* events_histo, Double_t Q_apres_cible, Double_t Coul_par_top = 1.e-10);
 
   Double_t GetMeanGDirFreq(Int_t run) const;
 
   KVINDRAPulserDataTree* GetPulserData()
   {
      return fPulserData;
   }
   void WriteObjects(TFile*) override;
   void ReadObjects(TFile*) override;
 
   virtual Double_t GetRunLengthFromGeneDirect(int run) const;
   Float_t GetMeanPulser(const Char_t*, int);
 
   ClassDefOverride(KVINDRADB, 5)       //DataBase of parameters for an INDRA campaign
};
 
R__EXTERN KVINDRADB* gIndraDB;
 
#endif