test/TestTFndPreanMan_MultiRun.C

#include <Riostream.h>
#include "TROOT.h"
#include "TSystem.h"
#include "TApplication.h"

#include "TString.h"

#include "TFndHdt.h"
#include "TFndPreanMan.h"
#include "TFndPrean.h"

TFndProcessRec *fPrec;

TFndPreanMan *fPreanMan;

TString fDBhost;
TString fDataPath;
TString fRunType;

Int_t fFirstRun;
Int_t fLastRun;
Bool_t fIsLastRun;

Int_t fDebugLevel;

Int_t fTotEvCount;
Int_t fTotEvTgt;
Int_t fEvPerRunTgt;
Int_t fCurEv;
Int_t fIncrement;
Int_t fCurrentRun;

//___________________________________
void 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 UseProcessRec_DST(){

  Int_t JFDST = fPrec->GetJFDST();
  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;
  


}

//___________________________________
void EventLoop(){
  
  fPrec = fPreanMan->GetProcessRec();
  if(!fPrec) TerminateFroot(0,"fPrec not available!");

  while(1){
    if( fCurEv % 300 == 0 ) gROOT->Info("FinudaPreanalysis","Processing run %s-%d ===> event %d",fRunType.Data(),fCurrentRun,fCurEv);
    if(fTotEvTgt != -1 && fTotEvCount > fTotEvTgt){
      Printf("PreanalyzeRun ---> maximum number of events reached for current preanalysis process");
      return;
    }
    if(fEvPerRunTgt != -1 && fCurEv > fEvPerRunTgt){
      Printf("PreanalyzeRun ---> maximum number of events reached for current run");
      return;
    }
    if( fPreanMan->GetNextFidaEvent() != 0) return; // Fill FGES
    UseProcessRec_GES();
    //    fPreanMan->AnalyzeJFGES(); // used by pre-analysis
    if( fPreanMan->ReconstructFidaEvent() !=0) return; // Fill FDST (exit loop in case of problems)
    UseProcessRec_DST();
    //    fPreanMan->AnalyzeJFDST(); // used by pre-analysis
    fCurEv++;
    fTotEvCount++;
  } // event loop completed
  
}

//___________________________________
Int_t ProcessSingleRun(){
  // return value:
  //              0: ok
  //             -1: run initialization not possible 

  Int_t initval = fPreanMan->InitNewPreanRun(fRunType,fCurrentRun); // offline mode
  if(initval!=0){
    //      Printf("initval: %d",initval);
    fPreanMan->FinishPreanRun(kFALSE);
    fCurrentRun += fIncrement;
    Printf("\n --- --- --- --- \n"); // simply nice ;)
    return -1;
  }
  
  Printf("New run started (%s %d)",fRunType.Data(),fCurrentRun);
  
  fCurEv = 1; // reset single run event counter
  
  // ---
  EventLoop();
  // ---
  
  fPreanMan->FinishPreanRun(fIsLastRun); // last run completed
  
  Printf("\n --- --- --- --- \n"); // simply nice ;)
  
  return 0;
}

//___________________________________
Int_t RunLoop(){
  
  if(!fPreanMan){
    gROOT->Error("RunLoop","fPreanMan not defined");
    TerminateFroot();
  }

  gROOT->Info("RunLoop","Starting loop on runs (%d -> %d)",fFirstRun,fLastRun);

  fCurrentRun = fFirstRun;
  fIncrement = (fLastRun >= fFirstRun) ? +1 : -1;
  
  while(1){ // run loop
    if(fLastRun >= fFirstRun && fCurrentRun > fLastRun) break;    
    else if(fLastRun < fFirstRun && fCurrentRun < fLastRun) break; // allow reverse processing
    
    if(fTotEvTgt != -1 && fTotEvCount >= fTotEvTgt){
      Printf("RunLoop ---> maximum number of events reached for current preanalysis process");
      break;
    }
    
    fIsLastRun = (fCurrentRun == fLastRun)? kTRUE : kFALSE;
    if(ProcessSingleRun() == -1) continue;

    //    
    if(fCurrentRun == -1) break; // one initialization only in online-mode
    //

    //
    fCurrentRun += fIncrement;
  } // run loop completed
  
  fPreanMan->CloseInterface(); // close used reconstruction (dynamic) library
  
  return 0;
  
}

//___________________________________
//___________________________________
//___________________________________
void TestTFndPreanMan_MultiRun(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){
  
  // --- INIT ---
  Int_t debug_lev = 0;
  //
  fDBhost = "$MSQL_DB_HOST";
  fDataPath = "$RDT";
  fRunType = run_type;
  fEvPerRunTgt = n_evs_run;
  fTotEvTgt = n_evs_tot;
  fDebugLevel = debug_lev;
  //
  fFirstRun = first_run;
  fLastRun = last_run;
  //
  fPreanMan = new TFndPreanMan(fDBhost,fDataPath,fDebugLevel);

  RunLoop();
}

---