CheckND.cxx

Go to the documentation of this file.

// $Id: CheckND.cxx,v 1.2 2007/03/01 17:02:25 rhatcher Exp $
//
// FILL_IN: [Document your code!!]
//
// n.tagg1@physics.ox.ac.uk
#include "CheckND.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 <vector>
#include <algorithm>
#include <string>
#include <TSystem.h>
#include <TFile.h>
#include <TH1.h>

JOBMODULE(CheckND, "CheckND",
          "CheckND");

CVSID("$Id: CheckND.cxx,v 1.2 2007/03/01 17:02:25 rhatcher Exp $");


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

CheckND::CheckND()
{
  fRawHist=fAdcHist=fHitHist=0;
}

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

CheckND::~CheckND()
{
  TFile f("checknd.root","RECREATE");
  if(fAdcHist) fAdcHist->Write();
  if(fHitHist) fHitHist->Write();
  f.Write();
  f.Close();
}

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

JobCResult CheckND::Ana(const MomNavigator* mom)
{

  static int latch53 = 0;

  const RawDaqSnarlHeader* header = 
    DataUtil::GetRawHeader<RawDaqSnarlHeader>(mom);
  if(header) {
    if((header->GetTrigSrc() & TRIGGER_BIT_SPILL)>0) {

      const RawDigitDataBlock* rddb = 
        DataUtil::GetRawBlock<RawDigitDataBlock>(mom);
      
      if(rddb) {
        
        rddb->SetDebugFlags(0);
        if(!fRawHist) fRawHist = new TH1F("hRaw","hRaw",5400000,0,54e6);
        if(!fAdcHist) fAdcHist = new TH1F("hAdc","hAdc",3000,-1000,2000);
        if(!fHitHist) fHitHist = new TH1F("hHit","hHit",3000,-1000,2000);
        
        fRawHist->Reset();

        double hits = 0;
        double ntdc = 0;
        double qtdc = 0;
        double q = 0;
        int first_tdc = 54000000;
        
        RawDigit *rawdigit;
        TIter iter = rddb->GetDatumIter();
        TObject* tobj;
        while( (tobj=iter.Next() ) ) {
          if( ( rawdigit = dynamic_cast<RawDigit*>(tobj) ) ) {;
            RawDigit* rawdigit = dynamic_cast<RawDigit*>(tobj);
            fRawHist->Fill(rawdigit->GetTDC(),rawdigit->GetADC());
            fAdcHist->Fill(rawdigit->GetTDC()-latch53,rawdigit->GetADC());
            fHitHist->Fill(rawdigit->GetTDC()-latch53);
            hits +=1;
            ntdc += rawdigit->GetTDC();
            qtdc += rawdigit->GetTDC()*rawdigit->GetADC();
            q    += rawdigit->GetADC();
            if(rawdigit->GetTDC() < first_tdc) first_tdc = rawdigit->GetTDC();
          }
        }
        cout << "Caught spill: "
             << "hits: " << hits
             << "\tq: " << q
             << "\tntdc: " << ntdc/hits
             << "\tqtdc: " << qtdc/q
             << "\tfirst: "  << first_tdc
             << "\tfirst-latch: "  << first_tdc-latch53
             << endl;
      }
    }
  }
  
  const std::vector<const RawVtmTimeInfoBlock*> vtmBlocks = 
    DataUtil::GetRawBlocks<RawVtmTimeInfoBlock>(mom);
  
  if(vtmBlocks.size()>0) {

    //for(UInt_t i=0;i<vtmBlocks.size(); i++) {
    //  vtmBlocks[i]->Print();
    // }
    
    // Find if a new value has been latched.
    static int sLast10MHz =0;
    if(vtmBlocks[0]->GetSpillTimeGPS() != sLast10MHz) {
      
      sLast10MHz = vtmBlocks[0]->GetSpillTimeGPS();
      
      // Get time to nearest spill.
      VldContext cx = vtmBlocks[0]->GetVldContext();
      VldTimeStamp spill = SpillTimeFinder::Instance().
        GetTimeOfRecentSpill(cx);
      
      int nsec_spill = spill.GetNanoSec();
      int nsec_10MHz = (int)(1e9*(double)vtmBlocks[0]->GetSpillTimeGPS()/10e6);
      int nsec_53MHz = (int)(1e9*(double)vtmBlocks[0]->GetSpillTimeRF()
                             /(double)vtmBlocks[0]->GetLastTimeStamp());
      
      static double offset_tot=0;
      static double nspill = 0;

      offset_tot +=  nsec_spill-nsec_10MHz;
      nspill += 1;
      
      std::cout << "Sec: " << spill.GetSec()
                << "\tGPS latch: " << nsec_10MHz
                << "\t RF latch: " << nsec_53MHz
                << "\t spill: " << nsec_spill
                << "\t spill-gps: " << nsec_spill-nsec_10MHz
                << "\t GPS-RF: " << nsec_10MHz - nsec_53MHz
                << std::endl;
      std::cout << "Avg spill-gps: " << offset_tot / nspill << endl;


      latch53 = vtmBlocks[0]->GetSpillTimeRF();
    }
  } 
  
  return JobCResult::kPassed;
}

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

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

  return r;
}

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

void CheckND::Config(const Registry& /*r*/)
{
//======================================================================
// Configure the module given the Registry r
//======================================================================

}

---