CurrentMeasurement.cpp

Go to the documentation of this file.

// $Id: CurrentMeasurement.cpp,v 1.4 2006/03/16 17:08:23 fmachefe Exp $

// Include files 
#include <ctime>
#include <iostream>
#include <sstream>
#include <string>

// ROOT include files
#include "TRandom.h"
#include "TH2D.h"
#include "TH1D.h"

// local
#include "CurrentMeasurement.h"

//-----------------------------------------------------------------------------
// Implementation file for class : CurrentMeasurement
//
// 2004-07-23 : Frédéric Machefert
//-----------------------------------------------------------------------------

//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
CurrentMeasurement::CurrentMeasurement(  ){
  setName ( "CurrentMeasurement" );
  setType ( "NI6008" );
  setTitle( "NI6008 current measurement" );  
  m_device = 0; 
}

//=========================================================================
//  virtual function intialize
//=========================================================================
StatusCode CurrentMeasurement::initialize ( ) {
  m_device = dynamic_cast <NI6008*> (element());
  m_numberOfDevices = m_device->numberOfDevices();
  time(&m_startTime); 
  if (application()->options()->dataStorage()){
    m_tree = new TTree(name().c_str(), title().c_str());
    m_tree->Branch("Run",&m_runNumber,"Run/I");
    m_tree->Branch("Event",&m_evtNumber,"Event/I");
    m_tree->Branch("Time",m_timestamp,(std::string("Time[")+
                       itos(m_device->numberOfDevices())+
                       std::string("]/I")).c_str());
    m_tree->Branch("Duration",m_duration,(std::string("Duration[")+
                       itos(m_device->numberOfDevices())+
                       std::string("]/I")).c_str());
    m_tree->Branch("ai0",m_ai0,
           (std::string("ai0[")+itos(m_device->numberOfDevices())+std::string("]/D")).c_str());
    m_tree->Branch("ai1",m_ai1,
           (std::string("ai1[")+itos(m_device->numberOfDevices())+std::string("]/D")).c_str());
    m_tree->Branch("ai2",m_ai2,
           (std::string("ai2[")+itos(m_device->numberOfDevices())+std::string("]/D")).c_str());
    m_tree->Branch("ai3",m_ai3,
           (std::string("ai3[")+itos(m_device->numberOfDevices())+std::string("]/D")).c_str());
  }
  for (int i = 0; i<m_device->numberOfDevices(); ++i){
    std::string dev = m_device->device(i);
    info("Preparing dataStream for device "+dev);
    addDataStream((dev+std::string(":Duration")).c_str(),"Duration");
    addDataStream((dev+std::string(":ai0")).c_str(),"Channel 0");
    addDataStream((dev+std::string(":ai1")).c_str(),"Channel 1");
    addDataStream((dev+std::string(":ai2")).c_str(),"Channel 2");
    addDataStream((dev+std::string(":ai3")).c_str(),"Channel 3");
  }
  m_runNumber=application()->options()->runNumber();
  return StatusCode::SUCCESS;
}

//=========================================================================
//  Virtual function execute
//=========================================================================
StatusCode CurrentMeasurement::execute ( ) {
  for (int i = 0; i<m_device->numberOfDevices(); ++i){
    std::string dev = m_device->device(i);
    std::string val = m_device->cmd(dev, "ai");
    std::istringstream splitval(val);
    std::string v0, v1, v2, v3;
    splitval >> v0;
    splitval >> v1;
    splitval >> v2;
    splitval >> v3;
    // m_ai0[i]=toDouble(v0);
    // m_ai1[i]=toDouble(v1);
    // m_ai2[i]=toDouble(v2);
    // m_ai3[i]=toDouble(v3);
    m_ai0[i]=atof(v0.c_str());
    m_ai1[i]=atof(v1.c_str());
    m_ai2[i]=atof(v2.c_str());
    m_ai3[i]=atof(v3.c_str());
    time_t t = time(0);
    m_timestamp[i]=t;
    m_duration[i]=difftime(t, m_startTime);
    dataFill(0+5*i, m_duration[i]);
    dataFill(1+5*i, m_ai0[i]);
    dataFill(2+5*i, m_ai1[i]);
    dataFill(3+5*i, m_ai2[i]);
    dataFill(4+5*i, m_ai3[i]);
  }
  if (application()->options()->dataStorage()){
    m_evtNumber=application()->options()->nEvt()-1;
    m_tree->Fill();
  }
  wait(m_frequency);
  return StatusCode::SUCCESS;
}

//=========================================================================
//  Virtual function finalize
//=========================================================================
StatusCode CurrentMeasurement::finalize ( ) {
  // if (application()->options()->dataStorage()){
  //   m_tree->Print();
  // }
  return StatusCode::SUCCESS;
}