mcr/hdt_read/AccessHDT_Appl_1.C

#include "TFile.h"
#include "TTree.h"

#include "AccessHDT_MultiRun.C"

class HdtAcc_Appl1: public THdtAccess{

private:
  TTree *fTofTree;

public:
  HdtAcc_Appl1();
  ~HdtAcc_Appl1();

  void InitTrees();
  void SaveTrees(const TString &fnam);

  void Pre_Run_Operation();
  void Post_Run_Operation();

  void UseHdtEvent();

  ClassDef(HdtAcc_Appl1,1)
    };

ClassImp(HdtAcc_Appl1)

//_______________________________
HdtAcc_Appl1::HdtAcc_Appl1():
  fTofTree()
{

  //  cout << "HdtAcc_Appl1 constructor called" << endl;

}

//_______________________________
HdtAcc_Appl1::~HdtAcc_Appl1(){

}

//_______________________________
void HdtAcc_Appl1::InitTrees(){
  fTofTree = new TTree("TofTree","TOF tree");
}

//_______________________________
void HdtAcc_Appl1::SaveTrees(const TString &fnam){

  TFile fout(fnam,"RECREATE");
  fout.cd();
  fTofTree->Write();

  fout.Purge();
  fout.Close();
}

//_______________________________
void HdtAcc_Appl1::Pre_Run_Operation(){

  //  cout << "Performing pre-run operation for run " << fCurrentRun << endl;
  fTofTree->Reset();

}

//_______________________________
void HdtAcc_Appl1::Post_Run_Operation(){

  //  cout << "Performing post-run operation for run " << fCurrentRun << endl;

  TString fnam;
  fnam.Form("%s-%d.test.root",fRunType.Data(),fCurrentRun);

  SaveTrees(fnam);
}

//_______________________________
void HdtAcc_Appl1::UseHdtEvent(){
  
  //  cout << "...............UseHdtEvent.................... " << endl;
  //  fCurHdtEv->PrintHits();
}

  //_______________________________________________
  //_______________________________________________
  //_______________________________________________
void AccessHDT_Appl_1(TString run_type="FINU",Int_t first_run = -1,Int_t last_run = -1,Int_t n_evs_run=-1,Int_t n_evs_tot=-1){

  HdtAcc_Appl1 *hdtac = new HdtAcc_Appl1();
  Int_t debug_lev = 0;
  //
  hdtac->SetVerboseLevel(3);
  hdtac->SetUpdateEvGap(3000);
  //
  hdtac->InitTrees();
  hdtac->Init("$MSQL_DB_HOST","$RDT",run_type,first_run,last_run,n_evs_run,n_evs_tot,debug_lev);

  hdtac->RunLoop();

  hdtac->SaveTrees("test.root");
  delete hdtac;


}

---