FiltMBSpillModule.cxx

Go to the documentation of this file.

#include "FiltMBSpillModule.h"
#include "MessageService/MsgService.h"
#include "MinosObjectMap/MomNavigator.h"
#include "JobControl/JobCModuleRegistry.h" // For JOBMODULE macro
#include "JobControl/JobCommand.h"  
#include "Conventions/Munits.h"

#include "RawData/RawRecord.h"
#include "RawData/RawDaqSnarlHeader.h"

#include "OnlineUtil/mdTriggerCodes.h"

//Get MiniBooNE spill
#include "MiniBooNEAna/MBSpillAccessor.h"
#include "MiniBooNEAna/MBSpill.h"

// ROOT includes
#include "TSystem.h"
#include "TFile.h"
#include "TTree.h"
#include "TMath.h"

JOBMODULE(FiltMBSpillModule, "FiltMBSpillModule",
          "MBSpill");
CVSID("$Id: FiltMBSpillModule.cxx,v 1.5 2008/03/14 16:44:51 tjyang Exp $");
//......................................................................

FiltMBSpillModule::FiltMBSpillModule() : 
  fTWin(0.0),
  fTrigMask(0)
{

  // vector of timestamps to store miniboone data:
}
//......................................................................

FiltMBSpillModule::~FiltMBSpillModule()
{
}

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

void FiltMBSpillModule::BeginJob()
{
 

}

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

void FiltMBSpillModule::EndJob()
{
}

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

void FiltMBSpillModule::BeginFile()
{
}

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

void FiltMBSpillModule::EndFile()
{
}

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

void FiltMBSpillModule::BeginRun()
{
}

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

void FiltMBSpillModule::EndRun()
{
}

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

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

  JobCResult result = JobCResult::kFailed;
  TObject* recObj = 0;
  int    nrawrec=0;
  int    nrawblk=0;
  int    nsnarlblk=0;
  int trgsrc=0;
  VldTimeStamp vldts;  
  time_t trgsec;
  int trgnsec;

  int snarl;
  int run;

  double dt = -9e9;
  int mbsec;
  int mbnsec;

  // Loop over rawrecords
  TIter recItr = mom->FragmentIter();
  while ( ( recObj = recItr.Next() ) ) {
    RawRecord* rawrec = dynamic_cast<RawRecord*>(recObj);
    if (!rawrec) continue;  // skip things in mom that aren't RawRecords
    nrawrec++;

    // Get trigger source
    const RawDaqSnarlHeader* snarlHdr = 
      dynamic_cast<const RawDaqSnarlHeader*>(rawrec->GetHeader());
    if ( snarlHdr ) {
      trgsrc =  snarlHdr->GetTrigSrc();
      run = snarlHdr->GetRun();
      snarl = snarlHdr->GetSnarl();
    } // snarl header

    // Loop over rawblocks
    TObject* blkObj = 0;
    TIter blkItr = rawrec->GetRawBlockIter();
    while ( ( blkObj = blkItr.Next() ) ) {
      const RawDataBlock* rawBlk = dynamic_cast<const RawDataBlock*>(blkObj);
      if (!rawBlk) continue;  
      nrawblk++;
    
      // Get trigger time
      const RawSnarlHeaderBlock* snarlblk = 
        dynamic_cast<const RawSnarlHeaderBlock*>(rawBlk);
      if (!snarlblk) continue;
      nsnarlblk++;
      vldts = snarlblk->GetTriggerTime();
      trgsec = vldts.GetSec();
      trgnsec = vldts.GetNanoSec();
    
    
      VldTimeStamp tmpnd=vldts;
      
      //get MinoBooNE spill
      MBSpillAccessor& spillAccess = MBSpillAccessor::Get();
      
      const MBSpill *spillInfo = spillAccess.LoadSpill(vldts);
      if (spillInfo) {
        dt = double(spillInfo->GetTimeStamp()-vldts);
        mbsec = spillInfo->GetSec();
        mbnsec = spillInfo->GetNanoSec();
      }
    
    
      //Select "good" snarl here. Good means: inside MB spill window and allowed trigger source
      // Cosmic (M/N = 4 (hex=4); Activity = 16 (hex=10)) or SGATE (65536) trigger source:
      if (trgsrc & fTrigMask) {
        
        if (TMath::Abs(dt)<fTWin) {
          MSG("Filt",Msg::kDebug) << "-------------------------------------------------" << endl;
          MSG("Filt",Msg::kDebug) << "ND snarl within +-" << fTWin << " seconds of MiniBooNE spill. The trigger source is: " << trgsrc << endl;
          MSG("Filt",Msg::kDebug) << "Run, Snarl: " << run << ", " << snarl << endl;
          MSG("Filt",Msg::kDebug) << "dt(MB-MINOS snarl)= " << TMath::Abs(dt) <<" seconds"<< endl;
          MSG("Filt",Msg::kDebug) << "ND vldtimestamp in sec, nsec, double: " << vldts.GetSec() << ", " << vldts.GetNanoSec() << ", " << vldts.AsString() << endl;
          MSG("Filt",Msg::kDebug) << "MB vldtimestamp in sec, nsec, double: " << mbsec << ", " << mbnsec << ", " << (spillInfo->GetTimeStamp()).AsString() << endl;
          result = JobCResult::kPassed;
        }  
      }
    
    } // loop over blocks
  }  // loop over records

//  MSG("Filt",Msg::kDebug) << 
//    nrawrec << " rawrecords found " << 
//    nrawblk << " rawblocks found " << 
//    nsnarlblk << " snarlblocks found " << endl;

  return result; // kNoDecision, kFailed, etc.



}

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

JobCResult FiltMBSpillModule::Reco(MomNavigator* /* mom */)
{
  return JobCResult::kPassed; // kNoDecision, kFailed, etc.
}

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

const Registry& FiltMBSpillModule::DefaultConfig() const
{
  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("TWin",     1.0e-3);

  Int_t thisTrigMask = 0;

  // From OnlineUtil/mdTriggerCodes.h:
  // TRIGGER_BIT_NULL         0x1  /* Dump of Time Frame - ie singles */
  // TRIGGER_BIT_RAWLI        0x2  /* Raw LI events */
  // TRIGGER_BIT_PLANE        0x4  /* M/N Plane Trigger */
  // TRIGGER_BIT_TOF          0x8  /* Time Of Flight Trigger */
  // TRIGGER_BIT_ACTIVITY    0x10  /* Detector Activity Trigger */
  // TRIGGER_BIT_E4          0x20  /* E4 Trigger */
  // TRIGGER_BIT_MINBIAS   0x8000  /* Minimum Bias Sampling Trigger */
  // TRIGGER_BIT_SPILL    0x10000  /* Spill Trigger */
  // TRIGGER_BIT_SPILL_IP 0x20000  /* Remote Spill Trigger - ie spill over IP */
  thisTrigMask |= TRIGGER_BIT_PLANE;
  thisTrigMask |= TRIGGER_BIT_ACTIVITY;
  thisTrigMask |= TRIGGER_BIT_SPILL;
  //  cout << "Trigger mask set to: " << thisTrigMask << endl;
  
  r.Set("TrigMask", thisTrigMask);
  r.LockValues();

  return r;
}

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

void FiltMBSpillModule::Config(const Registry& r)
{
  fTWin = r.GetDouble("TWin");
  fTrigMask = r.GetInt("TrigMask");
}

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

void FiltMBSpillModule::Report()
{
}

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

void FiltMBSpillModule::Reset()
{
}

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

void FiltMBSpillModule::HandleCommand(JobCommand* /* c */)
{
}

---