OltTimeScan.cxx

Go to the documentation of this file.

// $Id: OltTimeScan.cxx,v 1.7 2004/03/18 17:10:16 rhatcher Exp $
//
// OltBlockScan.cxx
//
// Block Scan
//
// Author: Phil Adamson 5/2002

#include "JobControl/JobCModuleRegistry.h"
#include "MinosObjectMap/MomNavigator.h"
#include "OfflineTrigger/OltTimeScan.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/RawDeadChipBlock.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>
#include <map>
#ifndef __CINT__
#include "OnlineUtil/mdBlockDefinitions.h"
#endif
ClassImp(OltTimeScan)

  // ......

CVSID("$Id: OltTimeScan.cxx,v 1.7 2004/03/18 17:10:16 rhatcher Exp $");
JOBMODULE(OltTimeScan,"OltTimeScan","The Offline Trigger Time Scanner");


OltTimeScan::OltTimeScan()
{
  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();
  fPlex = 0;
}

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

JobCResult OltTimeScan::Reco(MomNavigator* mom)
{
  int plane[2][61] = {{0}, {0}};
  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("RawSnarlHeaderBlock")) {
            RawSnarlHeaderBlock *rshb =dynamic_cast<RawSnarlHeaderBlock *>(rawdata);
            fTreeTF = rshb->GetTimeFrameNo();
            MSG("Olt",Msg::kInfo)<<"Timeframe "<<fTreeTF<<"  "
                                 <<"Snarl "<<rshb->GetSnarl()<<endl;
            if (fDoPlex && !fPlex) fPlex = new PlexHandle(rshb->GetVldContext());
          } 
          fTreeIndex = 0;
          MSG("Olt",Msg::kInfo) <<fRecord<<": "<<rawdata->GetBlockId()<<endl;
          if (fDump) {
          RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
          if (rddb) {
            multimap<int,RawDigit *> hitmap;
            TIter it = rddb->GetDatumIter();
            while (TObject *obj = it.Next()) {
              RawVaDigit *rvd = dynamic_cast<RawVaDigit *>(obj);
              if (rvd) {
                hitmap.insert(pair<int,RawDigit *>(rvd->GetTDC(),rvd));
              }
              RawQieDigit *rqd = dynamic_cast<RawQieDigit *>(obj);
              if (rqd) {
                hitmap.insert(pair<int,RawDigit *>(int(rqd->GetTDC()*16*58.0/77.0 + 133),rqd));
              }
             }

            for (multimap<int,RawDigit *>::iterator it = hitmap.begin();
                 it!=hitmap.end();++it) {
              if (!fDoPlex) {         MSG("Olt",Msg::kInfo)<<(*it).first<<" "<<*((*it).second)<<endl; }
              else {
                // Use the plex
                RawChannelId rcid = (*it).second->GetChannel();
                PlexSEIdAltL altl= fPlex->GetSEIdAltL(rcid);
                MSG("Olt",Msg::kInfo)<<"Time: "<<(*it).first<<" Plane "<<altl.GetPlane()
                                     <<" End "<<altl.GetEnd()
                                     <<"   "<<rcid<<endl;
                int p = altl.GetPlane();
                if (p>-1 && p < 61) ++plane[altl.GetEnd()-1][p];
              }
            }
            for (int q=0;q<2;q++) {
            for (int i=0;i<61;i++) {
              if (plane[q][i]==0) {cout<<"-";}
              else if (plane[q][i]>9) {cout<<"9";}
              else {cout<<plane[q][i];}

            }
            cout<<endl;
            }
           }
          }       
          RawRunCommentBlock *rc = dynamic_cast<RawRunCommentBlock *>(rawdata);
          if (rc) cout<<rc->GetRunComment()<<endl;

          RawDeadChipBlock *dead = dynamic_cast<RawDeadChipBlock *>(rawdata);
          if (dead) cout<<"Dead Chips: "<<dead->GetNumberOfEntries()<<endl;
        }
      }
      
    }
  }
  return JobCResult::kAOK;

}


 
const Registry& OltTimeScan::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("Plex", true);
  r.LockValues();

  return r;
}

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

void OltTimeScan::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("Plex",tmpb)) { cout<<"Plex "<<tmpb<<endl;fDoPlex = bool(tmpb); }
}

---