SpillTime_module.cc

Go to the documentation of this file.

////////////////////////////////////////////////////////////////////////
/// \brief Print the time of each NuMI spill, or another specified
/// trigger type. This is intended for use by EliminateBeamSpills.
///
/// The time printed is the NOvA timestamp, which is a monotonically
/// increasing integer free of leap seconds. The only leap second
/// trouble is that since Accelerator Division systems don't handle leap
/// seconds well, we might be taking triggers that are 1s (or other
/// amounts) off of the physical spills near leap seconds.
///
/// This method of getting timestamps relies on the spill server having
/// been up and running at the time the data being processed was taken.
/// It is possible that we have taken some data with other triggers
/// while the beam was running, but the spill server was down, in which
/// case we do not learn the spill times from this method.
///
/// \author M. Strait <mstrait@fnal.gov>
////////////////////////////////////////////////////////////////////////

#include "art/Framework/Core/ModuleMacros.h"
#include "art/Framework/Core/EDAnalyzer.h"

#include "RawData/RawTrigger.h"
#include "DAQDataFormats/TriggerDefines.h"

class spilltime : public art::EDAnalyzer {
  public:
  explicit spilltime(fhicl::ParameterSet const& pset): 
    EDAnalyzer(pset), trigger2spill(pset.get<double>("trigger2spill")),
                      trigid(pset.get<int>("trigid")) { }
  virtual ~spilltime() { } // compiles, but does not run, without this
  void analyze(const art::Event& evt);

  private:
  double trigger2spill;
  int trigid;
};

void spilltime::analyze(const art::Event & evt)
{
  art::Handle< std::vector<rawdata::RawTrigger> > rawtrigger;
  evt.getByLabel("daq", rawtrigger);
  if(rawtrigger.failedToGet()){
    fprintf(stderr, "Unexpected lack of 'daq' product!\n");
    return;
  }
  if(rawtrigger->empty()) return;

  if(((*rawtrigger)[0].fTriggerMask_TriggerType)+1 != trigid) return;

  const double tdc_per_second = 64e6;
  const uint64_t trigger2spills_tdc = uint64_t(trigger2spill * tdc_per_second);

  // This is a monotonically increasing timestamp with no worries about leap
  // seconds.  
  const uint64_t triggertime = (*rawtrigger)[0].fTriggerTimingMarker_TimeStart;
  
  const uint64_t spilltime = triggertime + trigger2spills_tdc;

  printf("Spilltime: %lu\n", spilltime);
}

DEFINE_ART_MODULE(spilltime)