OltDataIntegrity.cxx

Go to the documentation of this file.

// $Id: OltDataIntegrity.cxx,v 1.1 2003/08/27 12:52:06 phil.adamson Exp $
//
// OltDataIntegrity.cxx
//
// Block Scan
//
// Author: Phil Adamson 5/2002

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "OfflineTrigger/OltDataIntegrity.h"
#include "MessageService/MsgService.h"
#include "RawData/RawRecord.h"
#include "RawData/RawHeader.h"
#include "RawData/RawDataBlock.h"
#include "RawData/RawDaqHeaderBlock.h"
#include "RawData/RawSnarlHeaderBlock.h"
#include "RawData/RawChannelId.h"
#include "RawData/RawVaDigit.h"
#include "RawData/RawQieDigit.h"
#include "RawData/RawDigitDataBlock.h"
#include "RawData/RawRunCommentBlock.h"
#include "RawData/RawRunConfigBlock.h"

#include "Persistency/PerOutputStreamManager.h"
#include "Persistency/PerOutputStream.h"
#include "Rotorooter/RotoObjectifier.h"
#include "Plex/PlexHandle.h"
#include "Plex/PlexSEIdAltL.h"
#include "Validity/VldContext.h"

#include "TString.h"
#include "TFile.h"
#include "TTree.h"
#include <string>

#include <algorithm>
#include <cstdio>
#include <cstring>
#include <iostream>
#include <fstream>
#ifndef __CINT__
#include "OnlineUtil/mdBlockDefinitions.h"
#endif
ClassImp(OltDataIntegrity)

  // ......

CVSID("$Id: OltDataIntegrity.cxx,v 1.1 2003/08/27 12:52:06 phil.adamson Exp $");
JOBMODULE(OltDataIntegrity,"OltDataIntegrity","The Offline Trigger Block Scanner");





OltDataIntegrity::OltDataIntegrity()
{
  fRecord = 0;
}

OltDataIntegrity::~OltDataIntegrity()
{
  TDirectory *save = gDirectory;
  fFile->cd();
  fTree->Write();
  fFile->Close();
  save->cd();
}

JobCResult OltDataIntegrity::Reco(MomNavigator* mom)
{
  //MSG("Olt",Msg::kVerbose) << "OltModule::Reco()\n";
  // Find RawRecord fragment in MOM.
  TIter iter = mom->FragmentIter();
  while (TObject *obj = iter.Next()) {
    RawRecord *rawrec = dynamic_cast<RawRecord *>(obj);
    if (rawrec) {
      ++fRecord;
      MSG("Olt",Msg::kVerbose) <<"---------------New Record---------------"<<endl;
      TIter recit = rawrec->GetRawBlockIter();
      while (TObject *obj = recit.Next()) {
        RawDataBlock *rawdata = dynamic_cast<RawDataBlock *>(obj);
        if (rawdata) {
          if (rawdata->InheritsFrom("RawSnarlHeaderBlock")) {
            RawSnarlHeaderBlock *rshb =dynamic_cast<RawSnarlHeaderBlock *>(rawdata);
            fTimeframe = rshb->GetTimeFrameNo();
            int run = rshb->GetRun();
            int subrun = rshb->GetSubRun();
            if ((run!=fRun)||(subrun!=fSubRun)) {
              MSG("Olt",Msg::kInfo)<<"Now starting run "<<run<<":"<<subrun<<endl;
              fRun = run;
              fSubRun = subrun;
            }
          } 
         
          MSG("Olt",Msg::kVerbose) <<fRecord<<": "<<rawdata->GetBlockId()<<endl;
          RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
          if (rddb) {
            TIter it = rddb->GetDatumIter();
            while (TObject *obj = it.Next()) {
              RawVaDigit *rvd = dynamic_cast<RawVaDigit *>(obj);
              if (rvd) {
                int tb = (rvd->GetTDC()/16000000);
                fTimeBlockHits[0][tb] = fTimeBlockHits[0][tb] + 1;
                fTimeBlockHits[0][40] = fTimeBlockHits[0][40] + 1;
              }
              RawQieDigit *rqd = dynamic_cast<RawQieDigit *>(obj);
              if (rqd) {
                int error = rqd->GetErrorCode();
                RawChannelId rcid = rqd->GetChannel();
                int crate = rcid.GetCrate();
                int tb = rqd->GetTDC()/int(1000000*77.0/58.0);
                fTimeBlockHits[crate][tb] = fTimeBlockHits[crate][tb] + 1;
                fTimeBlockHits[crate][40] = fTimeBlockHits[crate][40] + 1;
                if (error) {
                  MSG("Olt",Msg::kInfo)<<"Run "<<fRun<<":"<<fSubRun<<" Timeframe "<<fTimeframe<<" QIE error "<<error<<" in channel "<<rcid<<endl;               
                }
              }
            }
          }
        } 
      }  
    }
  }

  if (fTimeframe&&(fTimeframe%fCheckPeriod==0)) {
    // Check timeblocks...
    for (int crate=0;crate<3;crate++) {
      for (int tb = 0; tb<40; tb++) {
        if (float(fTimeBlockHits[crate][tb])/float(fTimeBlockHits[crate][40])
            <fTBThresh) 
          MSG("Olt",Msg::kInfo) << "Warning: at tf " <<fTimeframe<< " Crate "
                                <<crate<<" Timeblock "
                                <<tb<<" has only "
                                <<float(fTimeBlockHits[crate][tb])/float(fTimeBlockHits[crate][40])
                                <<" of the hits"<<endl;
        fTimeBlockHits[crate][tb]=0;
      }
      fTimeBlockHits[crate][40] = 0;
    }
  }
  return JobCResult::kAOK;
}


 
const Registry& OltDataIntegrity::DefaultConfig() const
{

  static Registry r; // Default configuration for module


  std::string name = this->JobCModule::GetName();
  name += ".config.default";
  r.SetName(name.c_str());

  // Set values in configuration
  r.UnLockValues();
  r.Set("DumpData", true);
  r.Set("NumDump",10);
  r.Set("DoTree", false);
  r.Set("StartTF",0);
  r.Set("OutFile","di.root");
  r.Set("CheckPeriod",10);
  r.Set("TBThresh",0.01);
  r.LockValues();

  return r;
}

//......................................................................

void OltDataIntegrity::Config(const Registry& r)
{
  int   tmpb;
  cout<<"Teapots\n";
  cout<<r<<endl;
  if (r.Get("DumpData",tmpb)) { cout<<"Hello "<<tmpb<<endl;fDump = bool(tmpb); }  
  if (r.Get("DoTree",tmpb)) fDoTree = bool(tmpb);
  if (r.Get("StartTF",tmpb)) fStartTF = tmpb;
  if (r.Get("NumDump",tmpb)) fNumDump = tmpb;

  const char *s;
  if (r.Get("OutFile",s)) fFileName = std::string(s);
   TDirectory *save = gDirectory;
  fFile = new TFile(fFileName.c_str(),"RECREATE");
  fTree = new TTree("ditree","Data Integrity Tree");
  save->cd();

  double f;
  if (r.Get("CheckPeriod",tmpb)) fCheckPeriod = tmpb;
 if (r.Get("TBThresh",f)) fTBThresh = f;
}

---