mcr/prean/ZebraAccess/AccessZebra_Appl_1.C

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

#include "AccessZebra_MultiRun.C"


class ZebraAcc_Appl1: public TZebraAccess{
  
private:
  TTree *fTofTree;

public:
  ZebraAcc_Appl1();
  ~ZebraAcc_Appl1();

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

  void Pre_Run_Operation();
  void Post_Run_Operation();

  void UseProcessRec_GES();
  void UseProcessRec_DST();

  ClassDef(ZebraAcc_Appl1,1)
    };

ClassImp(ZebraAcc_Appl1)

//_______________________________
ZebraAcc_Appl1::ZebraAcc_Appl1():
  fTofTree()
{
  
  cout << "ZebraAcc_Appl1 constructor called" << endl;
 
}

//_______________________________
ZebraAcc_Appl1::~ZebraAcc_Appl1(){
  
}

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

//_______________________________
void ZebraAcc_Appl1::SaveTrees(const TString &fnam){
  
  TFile fout(fnam,"RECREATE");
  fout.cd();
  fTofTree->Write();

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

//_______________________________
void ZebraAcc_Appl1::Pre_Run_Operation(){

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

}

//_______________________________
void ZebraAcc_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 ZebraAcc_Appl1::UseProcessRec_GES(){

  /*
  Int_t JFGES = fPrec->GetJFGES();
  Int_t trigbit = fPrec->GetIQ(JFGES+3);
  Int_t rnum = fPrec->GetIQ(JFGES+1);
  Int_t evnum = fPrec->GetIQ(JFGES+2);
  cout << "JFGES = " << JFGES 
       << " ---> trig: " << trigbit
       << "; rnum: " << rnum
       << "; evnum: " << evnum << endl;
  
  Int_t L_DCH1 = fPrec->GetLQ(JFGES-4);
  while(L_DCH1){
    Int_t dch_num = fPrec->GetIQ(L_DCH1+1);
    cout << "dch_num = " << dch_num << endl;
    
    L_DCH1 = fPrec->GetLQ(L_DCH1);
  }
  */
}

//_______________________________
void ZebraAcc_Appl1::UseProcessRec_DST(){

  Int_t JFDST = fPrec->GetJFDST();
  if(!JFDST) return;
  Int_t trigbit = fPrec->GetIQ(JFDST+3);
  Int_t rnum = fPrec->GetIQ(JFDST+1);
  Int_t evnum = fPrec->GetIQ(JFDST+2);
  cout << "JFDST = " << JFDST 
       << " ---> trig: " << trigbit
       << "; rnum: " << rnum
       << "; evnum: " << evnum << endl;
  
  // --- eval

  // ---

  fTofTree->Fill();
}


//_______________________________
//_______________________________
//_______________________________
void AccessZebra_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){
  
  ZebraAcc_Appl1 *zac = new ZebraAcc_Appl1();
  Int_t debug_lev = 0;
  //
  zac->InitTrees();
  zac->SetFillDst(1);
  zac->Init("$MSQL_DB_HOST","$RDT",run_type,first_run,last_run,n_evs_run,n_evs_tot,debug_lev);

  zac->RunLoop();

  zac->SaveTrees("test.root");
  delete zac;
}

---