FarDetVarcErrorsDump.cxx

Go to the documentation of this file.

#include "FarDetVarcErrorsDump.h"

#include "MessageService/MsgService.h"     
#include "MinosObjectMap/MomNavigator.h"   
#include "JobControl/JobCommand.h"         
#include "JobControl/JobCModuleRegistry.h" 

#include "RawData/RawRecord.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawDaqHeaderBlock.h"
#include "RawData/RawVarcErrorInTfBlock.h"

#include "Validity/VldContext.h"
#include "Validity/VldTimeStamp.h"

#include "TDirectory.h"

CVSID("$Id: FarDetVarcErrorsDump.cxx,v 1.1 2006/04/21 17:40:11 blake Exp $");

JOBMODULE(FarDetVarcErrorsDump,"FarDetVarcErrorsDump","FarDetVarcErrorsDump");

FarDetVarcErrorsDump::FarDetVarcErrorsDump() :
  fVarcErrorsFile(0), 
  fVarcErrorsTree(0)
{

}

FarDetVarcErrorsDump::~FarDetVarcErrorsDump()
{

}

void FarDetVarcErrorsDump::BeginJob()
{
  MSG("FarDetVarcErrorsDump",Msg::kInfo) << " *** FarDetVarcErrorsDump::BeginJob() *** " << endl;

}

JobCResult FarDetVarcErrorsDump::Ana(const MomNavigator *mom) 
{
  MSG("FarDetVarcErrorsDump",Msg::kDebug) << " *** FarDetVarcErrorsDump::Ana(...) *** " << endl;

  UShort_t timeBlock,crateId,varcId,etcId,errCode;

  JobCResult result(JobCResult::kPassed);

  TIter momitr(mom->FragmentIter());
  while(const RawRecord* rawrec = dynamic_cast<const RawRecord*>(momitr())){
  
    MSG("FarDetVarcErrorsDump",Msg::kDebug) << "  *** NEXT RECORD ***  " << endl;  

    TIter rawrecitr = rawrec->GetRawBlockIter();
    TObject *tob;
    while((tob = rawrecitr())){
      MSG("FarDetVarcErrorsDump",Msg::kDebug) << tob->GetName() << endl;

      if(tob->InheritsFrom("RawDaqHeaderBlock")){
        MSG("FarDetVarcErrorsDump",Msg::kDebug) << " ... found RawDaqHeaderBlock " << endl;
        RawDaqHeaderBlock* rdb = (RawDaqHeaderBlock*)(tob);
        run = rdb->GetRun(); 
        timeframe = rdb->GetTimeFrameNum(); 
        date = (((VldTimeStamp)(rdb->GetTimeStamp())).GetSec()-1059696000)/(3600*24);
        time = (((VldTimeStamp)(rdb->GetTimeStamp())).GetSec()-1059696000)%(3600*24);
        MSG("FarDetVarcErrorsDump",Msg::kDebug) << "   ... timeframe=" << timeframe << endl; 
      }

      if(tob->InheritsFrom("RawVarcErrorInTfBlock")){
        MSG("FarDetDataQualityModule",Msg::kDebug) << " ... found RawVarcErrorInTfBlock" << endl;
        RawVarcErrorInTfBlock* rdb = (RawVarcErrorInTfBlock*)(tob);
    
        if( !fVarcErrorsFile ){
          //TString mystring("./results/fardet.rate");
          //mystring.Append("."); mystring+=run; 
          //mystring.Append(".root");
          TString mystring("fardet.varcerrors.root");
          TDirectory* tmpd = gDirectory;
          fVarcErrorsFile = new TFile(mystring.Data(),"RECREATE");
          fVarcErrorsTree = new TTree("VarcErrorsTree","VarcErrorsTree");
          fVarcErrorsTree->SetAutoSave(100);
          fVarcErrorsTree->Branch("date",&date,"date/I");
          fVarcErrorsTree->Branch("time",&time,"time/I");
          fVarcErrorsTree->Branch("run",&run,"run/I");
          fVarcErrorsTree->Branch("timeframe",&timeframe,"timeframe/I");
          fVarcErrorsTree->Branch("timeblock",&timeblock,"timeblock/I");
          fVarcErrorsTree->Branch("crate",&crate,"crate/I");
          fVarcErrorsTree->Branch("varc",&varc,"varc/I");
          fVarcErrorsTree->Branch("etc",&etc,"etc/I");
          fVarcErrorsTree->Branch("sparserr",&sparserr,"sparserr/I");
          fVarcErrorsTree->Branch("etcerr",&etcerr,"etcerr/I");
          fVarcErrorsTree->Branch("errcode",&errcode,"errcode/I");
          gDirectory = tmpd; 
        }

        crateId = rdb->GetCrateId();
        crateId-=1792;

        for(Int_t k=0;k<rdb->GetNumberOfSparsErrors();k++){
          rdb->UnpackSparsError(k,timeBlock,varcId,errCode);
          timeblock = timeBlock;
          crate = crateId;
          varc = varcId;
          etc = -1;
          sparserr = 1;
          etcerr = 0;
          errcode = errCode;
          TDirectory* tmpddd = gDirectory;
          fVarcErrorsFile->cd();
          fVarcErrorsTree->Fill();
          gDirectory = tmpddd;
        }

        for(Int_t k=0;k<rdb->GetNumberOfEtcErrors();k++){
          rdb->UnpackEtcError(k,timeBlock,varcId,etcId,errCode);
          timeblock = timeBlock;
          crate = crateId;
          varc = varcId;
          etc = etcId;
          sparserr = 0;
          etcerr = 1;
          errcode = errCode;
          TDirectory* tmpddd = gDirectory;
          fVarcErrorsFile->cd();
          fVarcErrorsTree->Fill();
          gDirectory = tmpddd;
        }

      }

    }
  }
  
  MSG("FarDetVarcErrorsDump",Msg::kDebug) << " *** FarDetVarcErrorsDump::Ana( ) FINISHED *** " << endl;

  return result;

}

void FarDetVarcErrorsDump::EndJob() 
{
  MSG("FarDetVarcErrorsDump",Msg::kInfo) << " *** FarDetVarcErrorsDump::EndJob() *** " << endl;
  if(fVarcErrorsFile){
    MSG("FarDetVarcErrorsDump",Msg::kInfo) << " *** saving data to file *** " << endl;
    TDirectory* tmpddd = gDirectory;
    fVarcErrorsFile->cd();
    fVarcErrorsTree->Write();
    fVarcErrorsFile->Close();
    gDirectory = tmpddd;
    MSG("FarDetVarcErrorsDump",Msg::kInfo) << " *** data saved to file *** " << endl;
  }
}

---