OltRawLI.cxx

Go to the documentation of this file.

// $Id: OltRawLI.cxx,v 1.2 2003/09/04 13:10:55 phil.adamson Exp $
//
// OltBlockScan.cxx
//
// Block Scan
//
// Author: Phil Adamson 5/2002

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

  // ......

CVSID("$Id: OltRawLI.cxx,v 1.2 2003/09/04 13:10:55 phil.adamson Exp $");
JOBMODULE(OltRawLI,"OltRawLI","The Offline Trigger Time Scanner");


OltRawLI::OltRawLI()
{
  fRecord = 0;
  TDirectory *save = gDirectory;
  fFile = new TFile("/home/pa/rawli.root","RECREATE");
  fTree = new TTree("tree","Raw LI Timing Summary");
  fTree->Branch("crate",&fTreeCrate,"crate/I");
  fTree->Branch("time",&fTreeTime,"time/I");
  fTree->Branch("subtime",&fTreeSubTime,"subtime/I");
  fTree->Branch("elec",&fElec,"elec/I");
  fTree->Branch("tf",&fTreeTF,"tf/I");
  fTree->Branch("index",&fTreeIndex,"index/I");
  save->cd();
  fPlex = 0;
}

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

JobCResult OltRawLI::Reco(MomNavigator* mom)
{
  TIter iter = mom->FragmentIter();
  while (TObject *obj = iter.Next()) {
    RawRecord *rawrec = dynamic_cast<RawRecord *>(obj);
    if (rawrec) {
      ++fRecord;
      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();
            fTreeIndex = rshb->GetSnarl();
            MSG("Olt",Msg::kInfo)<<"Timeframe "<<fTreeTF<<"  "
                                 <<"Snarl "<<rshb->GetSnarl()<<endl;
            if (fDoPlex && !fPlex) fPlex = new PlexHandle(rshb->GetVldContext());
          } 
          //      fTreeIndex = 0;
          RawDigitDataBlock *rddb = dynamic_cast<RawDigitDataBlock *>(rawdata);
          if (rddb) {
            multimap<int,RawDigit *> hitmap;
            bool gottrig = false;
            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));
                RawChannelId rcid = rvd->GetChannel();
                if (rcid.GetCrate()==0 && CHAD_TRIGPMT==rcid.GetChAdd()) {
                  fTrigTime = rvd->GetTDC();
                  gottrig = true;
                  MSG("Olt",Msg::kVerbose)<<"Trigger PMT is at "<<fTrigTime<<endl;
                }
              }
              RawQieDigit *rqd = dynamic_cast<RawQieDigit *>(obj);
              if (rqd) {
                hitmap.insert(pair<int,RawDigit *>(int(rqd->GetTDC()*16*58.0/77.0 + 133),rqd));
              }
            }
            if (!gottrig) MSG("Olt",Msg::kWarning)<<"Missing Trigger PMT in tf "<<fTreeTF<<" snarl "<<fTreeIndex<<endl;
            for (multimap<int,RawDigit *>::iterator it = hitmap.begin();
                 it!=hitmap.end();++it) {
              RawQieDigit *q = dynamic_cast<RawQieDigit *>((*it).second);
              if (q) {
                fElec = 2;
                fTreeTime = int(q->GetTDC()*16*58.0/77.0 + 133);
              } else {
                fElec = 1;
                fTreeTime = (*it).second->GetTDC();
              }
              RawChannelId rcid = (*it).second->GetChannel();
              fTreeCrate = rcid.GetCrate();
              fTreeSubTime = fTreeTime - fTrigTime;
              MSG("Olt",Msg::kVerbose)<<"Time "<<fTreeTime<<" Trig "<<fTrigTime<<endl;
              fTree->Fill();

        
            }
          }
          
          
          
        }
      }
      
    }
  }
  return JobCResult::kAOK;
  
}


 
const Registry& OltRawLI::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 OltRawLI::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); }
}

---