NDSgateTimeLooter.cxx

Go to the documentation of this file.

// $Id: NDSgateTimeLooter.cxx,v 1.4 2007/11/11 05:51:05 rhatcher Exp $
//
// FILL_IN: [Document your code!!]
//
// n.tagg1@physics.ox.ac.uk
#include "NDSgateTimeLooter.h"
#include "DatabaseInterface/DbiWriter.h"
#include "SpillTiming/SpillTimeND.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro
#include "TTree.h"
#include "TFile.h"
#include "RawData/RawDigitDataBlock.h"
#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"
#include "DataUtil/GetRawBlock.h"
#include "DataUtil/GetRawHeader.h"
#include "SpillTiming/SpillTimeFinder.h"
#include "Calibrator/Calibrator.h"
#include "RawData/RawDigit.h"
#include "RawData/RawVtmTimeInfoBlock.h"
#include "OnlineUtil/mdTriggerCodes.h"
#include "Conventions/Munits.h"
#include <vector>
#include <algorithm>
#include <string>
#include <TSystem.h>
#include <TFile.h>
#include <TH1.h>

JOBMODULE(NDSgateTimeLooter, "NDSgateTimeLooter",
          "NDSgateTimeLooter");

ClassImp(NdSgateRec)

CVSID("$Id: NDSgateTimeLooter.cxx,v 1.4 2007/11/11 05:51:05 rhatcher Exp $");


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

NDSgateTimeLooter::NDSgateTimeLooter()
{
  fFile = 0;
  fTree = 0;
  fInitialized = false;
}

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

NDSgateTimeLooter::~NDSgateTimeLooter()
{
  Write();
  if(fFile) {
    if(fTree) fTree->Write();
    fFile->Write();
    fFile->Close();
    delete fFile;
  }
}

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

JobCResult NDSgateTimeLooter::Ana(const MomNavigator* mom)
{
  const RawDaqHeader* header = 
    DataUtil::GetRawHeader<RawDaqHeader>(mom);
  
  if(header) {
    if(!fInitialized) {
      fInitialized = true;
      fStart = header->GetVldContext().GetTimeStamp();
      fEnd = fStart;
      fStamps.clear();                         
      fRunStart = fRunStop = header->GetRun();
    }
    
    if(header->GetVldContext().GetTimeStamp()>fEnd) 
      fEnd = header->GetVldContext().GetTimeStamp();
    if(header->GetRun() > fRunStop) 
      fRunStop =  header->GetRun();

  }

  const std::vector<const RawVtmTimeInfoBlock*> vtmBlocks = 
    DataUtil::GetRawBlocks<RawVtmTimeInfoBlock>(mom);
  
  if(vtmBlocks.size()>0) {
    
    static int last_latch = -1;
    if(vtmBlocks[0]->GetSpillTimeGPS() != last_latch) {
      last_latch = vtmBlocks[0]->GetSpillTimeGPS();

    if(vtmBlocks[0]->GetTimeFrameNum()!=0) 
    {
    
      // Find the 10MHz latch time. See how it varies from crate to crate.
      int latch10_tot = 0;
      int latch10_low  = 100000000;
      int latch10_high = 0;
      int nlatch  = 0;

      for(UInt_t i=0;i<vtmBlocks.size();i++) {
        int l = vtmBlocks[i]->GetSpillTimeGPS();
        if(l<latch10_low) latch10_low = l;
        if(l>latch10_high) latch10_high = l;
        latch10_tot += l;         
        nlatch++;
      }
      
      double latch10_avg = (double)latch10_tot/(double)nlatch;

      // 53 Mhz quantites:
      int nsec_53MHz = (int)(1e9*(double)vtmBlocks[0]->GetSpillTimeRF()
                             /(double)vtmBlocks[0]->GetLastTimeStamp());
        

      if((latch10_high-latch10_low)>1) {
        MSG("SpillTime",Msg::kWarning) << "Not all crates agree on 10MHz time." << endl;
        MSG("SpillTime",Msg::kWarning) << "Earliest: " << latch10_low 
                                         << "  Latest: " << latch10_high 
                                         << "  Mean: " << latch10_avg << endl;
      }
          
      // Get time of nearest second.
      VldContext cx = vtmBlocks[0]->GetVldContext();

      // Get second.
      int second = cx.GetTimeStamp().GetSec()                   // Raw second
        + (int)TMath::Nint(cx.GetTimeStamp().GetNanoSec()/1e9)  // Round up or down.
        -1;                                                     // Start of prev second.

      VldTimeStamp sgate_time(second, (int)(latch10_avg * 100.));
      VldTimeStamp sgate_53  (second, nsec_53MHz);
      
      VldContext check_context(Detector::kNear,SimFlag::kData,sgate_time);
      VldTimeStamp t_database = SpillTimeFinder::Instance().
        GetTimeOfNearestSpill(cx);
      
      MSG("SpillTime",Msg::kDebug) << "SGATE: " << sgate_time.AsString() << endl;
      MSG("SpillTime",Msg::kDebug) << "   DB: " << t_database.AsString() << endl;
      MSG("SpillTime",Msg::kDebug) << "   dt: " << (double)(sgate_time - t_database)*1e9 << "ns"
                                    << "\t tf: " << vtmBlocks[0]->GetTimeFrameNum() << endl;
      if(fWriteToFile) {
        static NdSgateRec* rec = 0;
        if(fFile==0) {
          rec = new NdSgateRec;
          fFile = new TFile(GetConfig().GetCharString("FileName"),"RECREATE");
          fTree = new TTree("vtm","vtm");
          fTree->Branch("ndrec","NdSgateRec",&rec,32000,1);
        }

        rec->run       = vtmBlocks[0]->GetRun();
        rec->subrun    = vtmBlocks[0]->GetSubRun();
        rec->monitorts = cx.GetTimeStamp();
        rec->sgate    = sgate_time;
        rec->sgate_53 = sgate_53;
        rec->spill    = t_database;
        rec->dspill   = (double)(sgate_time-t_database);
        fTree->Fill();
      }

      if(fWriteToDB) 
        fStamps.push_back(sgate_time);
    }
  } 
  }
  
  if((double)(fEnd-fStart) > fTableDuration) {
    Write();
    fStart = fEnd;
    fEnd = fStart;
  }
  
  return JobCResult::kPassed;
}

void NDSgateTimeLooter::Write() 
{
  if(fWriteToDB) {
    
    MSG("SpillTime",Msg::kInfo) << "Time start: " << fStart.AsString() << endl;
    MSG("SpillTime",Msg::kInfo) << "Time end: "   << fEnd.AsString() << endl;  
    
    VldTimeStamp vldCreate = fStart;
    
    VldRange range( (int)Detector::kNear,
                    (int)SimFlag::kData,
                    fStart,
                    fEnd,
                    Form("Timestamps looted from ND VTM data, runs %d-%d",fRunStart,fRunStop)
                    );
    
    MSG("SpillTime",Msg::kInfo) << "Starting write: " << range.AsString() << endl;
    
    DbiWriter<SpillTimeND> writer (range,
                                   -1, //Aggregate
                                   SpillTimeND::kTask_Vtm, // task
                                   vldCreate, // Creation date
                                   0,   // db no
                                   Form("Timestamps looted from ND VTM data, runs %d-%d",fRunStart,fRunStop)
                                   );
    
    for(UInt_t i = 0; i<fStamps.size(); i++) {
      SpillTimeND datum(fStamps[i]);
      writer << datum;
    }
    
    writer.Close();
  }
  
  
  fStamps.clear();
  fInitialized = false;
}


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

const Registry& NDSgateTimeLooter::DefaultConfig() const
{
//======================================================================
// Supply the default configuration for the module
//======================================================================
  static Registry r; // Default configuration for module

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

  // Set values in configuration
  r.UnLockValues();

  r.Set("TableDuration",8.0*Munits::hour);
  r.Set("WriteToDB",1);
  r.Set("WriteToFile",0);
  r.Set("FileName","nd_sgates.root");

  r.LockValues();

  return r;
}

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

void NDSgateTimeLooter::Config(const Registry& r)
{
  r.Get("TableDuration",fTableDuration);
  r.Get("WriteToDB",fWriteToDB);
  r.Get("WriteToFile",fWriteToFile);
}

---