mcr/hdt_read/Tofino_PID_Check.C

#include "TFndGenInfo.h"
#include "TFndHdt.h"
#include "Gui.h"

#include "AccessHDT_MultiRun.C"

#include "FROOT.h"
using namespace FROOT;

static const Float_t K_EN_THR = 1000; // ADC calib. channel

class HdtAppl_TofiPID: public THdtAccess{

private:

  TFndHTof *fHtof;  
  

public:
  HdtAppl_TofiPID();
  ~HdtAppl_TofiPID();
  
  void InitProcess();
  void FinishProcess();
  
  void Pre_Run_Operation();
  void Post_Run_Operation();
  
  void UseHdtEvent();
  
  // --- added for current application
  static Bool_t CheckBtb(const Int_t &sl_ref,const Int_t &sl_opp);

  // ---

  ClassDef(HdtAppl_TofiPID,1)
    };

ClassImp(HdtAppl_TofiPID)

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

//_______________________________
HdtAppl_TofiPID::~HdtAppl_TofiPID(){

  FinishRun();

  fHtof = 0; // owned by 

}

//_______________________________
void HdtAppl_TofiPID::InitProcess(){

  cout << "Initializing process" << endl;
  
}

//_______________________________
void HdtAppl_TofiPID::FinishProcess(){

  cout << "Finishing process" << endl;

}

//_______________________________
void HdtAppl_TofiPID::Pre_Run_Operation(){

  cout << "Performing pre-run operation for run " << fCurrentRun << endl;

}

//_______________________________
void HdtAppl_TofiPID::Post_Run_Operation(){

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

}

//_______________________________
void HdtAppl_TofiPID::UseHdtEvent(){
  
  Int_t evnum = fCurHdtEv->GetEventNumber();
  cout << "...............Using Hdt-Event (" << evnum << ")" << endl;
  //  fCurHdtEv->PrintHits();
  Int_t nhits = fCurHdtEv->GetNumberOfTofHits();

  const Int_t nsl = K_N_TOFINO_SLABS;

  Bool_t IsHighTdc[nsl] = { kFALSE };
  Bool_t IsHighAdc[nsl] = { kFALSE };
  Bool_t IsTrigOk[nsl] = { kFALSE };


  for(Int_t nh=0; nh<nhits; nh++){ // loop on TOF hits
    fHtof = fCurHdtEv->GetTofHit(nh);
    Int_t slab = fHtof->GetSlab();
    if(slab < 1 || slab > nsl) continue;
    // tofino hits have been selected
    // ---
    cout << endl<< "nh: " << nh << "--------------" << endl;
    fHtof->PrintHit(2);
    // ---

    if(fHtof->GetHighMeanTimer() > -1) IsHighTdc[slab-1] = kTRUE;
    
    Float_t chrg = -1;
    chrg = fHtof->GetCalibAdcE()  * fHtof->GetCalibAdcP();
    if(chrg > 0){
      chrg = TMath::Sqrt(chrg);
      if( chrg > K_EN_THR) IsHighAdc[slab-1] = kTRUE;
    }
  }
  
  // --- check of obtained results
  Bool_t chk_res1 = 1; // hard-coded
  if(chk_res1){
    cout << "Slabs over TDC high threshold: ";
    for(Int_t i=0; i<nsl; i++) if( IsHighTdc[i] ) cout << " " << i+1;
    cout << endl;
    //
    cout << "Slabs over charge high threshold: ";
    for(Int_t i=0; i<nsl; i++) if( IsHighAdc[i] ) cout << " " << i+1;
    cout << endl << endl;
  }
  // ---
  








}



Bool_t HdtAppl_TofiPID::CheckBtb(const Int_t &sl_ref,const Int_t &sl_opp){
  
  for(Int_t of = 5 ; of <= 7 ; of++){
    Int_t check = ( (sl_ref+of) <= 12 ) ? ( sl_ref + of ) : ( sl_ref - (12-of) ) ;
    if(check == sl_opp) return kTRUE;
  }

  return kFALSE;
}



//_______________________________________________
//_______________________________________________
//_______________________________________________
void Tofino_PID_Check(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){
  // this is an application of the official method for reading HDT-files
  //  (see $FROOTSYS/mcr_included/AccessHDT_MultiRun.C


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

  hdt_appl->InitProcess();

  Int_t res_loop =  hdt_appl->RunLoop();
  cout << "Number of missing runs: " << res_loop << endl;

  // uncomment the following line if you want to
  // exclude processes in which at least one run is missing
  //  if( res_loop != 0 ) return;

  hdt_appl->FinishProcess();  
  // --- 
  delete hdt_appl; // need to keep histos
  
}

---