OltBlockScan.cxx

Go to the documentation of this file.

// $Id: OltBlockScan.cxx,v 1.13 2003/10/03 13:32:58 phil.adamson Exp $
//
// OltBlockScan.cxx
//
// Block Scan
//
// Author: Phil Adamson 5/2002

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "OfflineTrigger/OltBlockScan.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(OltBlockScan)

  // ......

CVSID("$Id: OltBlockScan.cxx,v 1.13 2003/10/03 13:32:58 phil.adamson Exp $");
JOBMODULE(OltBlockScan,"OltBlockScan","The Offline Trigger Block Scanner");


OltBlockScan::OltBlockScan()
{
  fRecord = 0;
  TDirectory *save = gDirectory;
  fFile = new TFile("/home/pa/bs.root","RECREATE");
  fTree = new TTree("tree","Timing Summary");
  fTree->Branch("crate",&fTreeCrate,"crate/I");
  fTree->Branch("time",&fTreeTime,"time/I");
  fTree->Branch("sec",&fCrateSec,"sec/I");
  fTree->Branch("nsec",&fCrateNSec,"nsec/I");
  fTree->Branch("tf",&fTreeTF,"tf/I");
  fTree->Branch("index",&fTreeIndex,"index/I");
  save->cd();
}

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

JobCResult OltBlockScan::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::kInfo) <<"---------------New Record---------------"<<endl;
      TIter recit = rawrec->GetRawBlockIter();
      while (TObject *obj = recit.Next()) {
        RawDataBlock *rawdata = dynamic_cast<RawDataBlock *>(obj);
        if (rawdata) {
          if (rawdata->InheritsFrom("RawDaqHeaderBlock")) {
            RawDaqHeaderBlock *rshb =dynamic_cast<RawDaqHeaderBlock *>(rawdata);
            MSG("Olt",Msg::kInfo)<<rshb->GetRun()<<" "<<rshb->GetSubRun()<<" "<<rshb->GetTimeStamp()<<endl;
          }
          if (rawdata->InheritsFrom("RawSnarlHeaderBlock")) {
            RawSnarlHeaderBlock *rshb =dynamic_cast<RawSnarlHeaderBlock *>(rawdata);
            fTreeTF = rshb->GetTimeFrameNo();
            MSG("Olt",Msg::kInfo) <<"Timeframe  "<<rshb->GetTriggerTime()<<endl;
            MSG("Olt",Msg::kInfo) <<"Run  "<<rshb->GetRun()<<" subrun "
<< rshb->GetSubRun()<<endl;
            MSG("Olt",Msg::kInfo) <<"Timeframe  "<<fTreeTF<<endl;
          } 
          fTreeIndex = 0;
          MSG("Olt",Msg::kInfo) <<fRecord<<": "<<rawdata->GetBlockId()<<endl;
          if (fDump) {
            RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
            if (rddb) {
              if (fTreeTF>fStartTF) {
                TIter it = rddb->GetDatumIter();
                int vnumber = 0;
                int qnumber = 0;
                while (TObject *obj = it.Next()) {
                  RawVaDigit *rvd = dynamic_cast<RawVaDigit *>(obj);
                  if (rvd) {
                  ++vnumber;
                  if (fNumDump>0 && vnumber>fNumDump) continue;
                    MSG("Olt",Msg::kInfo)<<fTreeTF<<" "<<*rvd<<endl;
                    fTreeCrate = 0;
                    fTreeTime = rvd->GetTDC();
                    fCrateSec = rvd->GetCrateT0().GetSec();
                    fCrateNSec = rvd->GetCrateT0().GetNanoSec();
                    //fTreeTF = fRecord;
                    TDirectory *save = gDirectory;
                    if (fDoTree) {
                      fFile->cd();
                      fTree->Fill();
                      save->cd();
                    }
                    ++fTreeIndex;
                  }
                  RawQieDigit *rqd = dynamic_cast<RawQieDigit *>(obj);
                  if (rqd) {
                  ++qnumber;
                  if (fNumDump>0 && qnumber>fNumDump) break;
                    MSG("Olt",Msg::kInfo)<<fTreeTF<<" "<<*rqd<<endl;
                RawChannelId rcid = rqd->GetChannel();
                fTreeCrate = rcid.GetCrate();
//                  fTreeCrate = 1;
                    fTreeTime = rqd->GetTDC();
                    fCrateSec = rqd->GetCrateT0().GetSec();
                    fCrateNSec = rqd->GetCrateT0().GetNanoSec();
                    //fTreeTF = fRecord;
                    if (fDoTree) {
                      TDirectory *save = gDirectory;
                      fFile->cd();
                      fTree->Fill();
                      save->cd();
                    }
                    ++fTreeIndex;
                  }
                }
              }
            }   
          }  
          RawRunCommentBlock *rc = dynamic_cast<RawRunCommentBlock *>(rawdata);
          if (rc) cout<<rc->GetRunComment()<<endl;
          RawRunConfigBlock *rcc = dynamic_cast<RawRunConfigBlock *>(rawdata);
          if (rcc) cout<<rcc->GetRunConfig()<<endl;
        }
      }
      
    }
  }
  return JobCResult::kAOK;

}


 
const Registry& OltBlockScan::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.LockValues();

  return r;
}

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

void OltBlockScan::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;
}

---