ADCMeasurement.cpp

Go to the documentation of this file.

// $Id: ADCMeasurement.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 "ADCMeasurement.h"

//-----------------------------------------------------------------------------
// Implementation file for class : ADCMeasurement
//
// 2004-07-23 : Fr�d�ric Machefert
//-----------------------------------------------------------------------------

//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
ADCMeasurement::ADCMeasurement(  ){
  setName ( "ADCMeasurement" );
  setType ( "NI6008" );
  setTitle( "NI6008 ADC measurement" );  
}

//=========================================================================
//  virtual function intialize
//=========================================================================
StatusCode ADCMeasurement::initialize ( ) {
  m_minRamp  = 0.; 
  m_maxRamp  = 1.2;
  m_stepRamp = 20;
  m_baseline = 1400;
  m_ramp = m_minRamp;
  time(&m_startTime); 
  m_device = dynamic_cast <NI6008*> (element());
  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,"Time/I");
    m_tree->Branch("Duration",&m_duration,"Duration/I");
    m_tree->Branch("ai0",&m_ai0,"ai0/F");
    m_tree->Branch("ai1",&m_ai1,"ai1/F");
    m_tree->Branch("ao0",&m_ao0,"ao0/F");
    //    m_tree->Branch("ao1",&m_ao1,"ao1/F");
    m_tree->Branch("di_0",&m_diA,"diA/I");
    m_tree->Branch("di_1",&m_diB,"diB/I");
  }
  info("Preparing dataStreams.");
  addDataStream("Duration","Duration");
  addDataStream("ai0", "Channel 0");
  addDataStream("ai1", "Channel 1");
  addDataStream("ao0", "Ramp Ampl");
  //  addDataStream("ao1", "ADC Chan Ctrl");
  addDataStream("diA","ADC A");
  addDataStream("diB","ADC B");
  addDataStream("ratioA","ratio amplitude A");
  addDataStream("ratioB","ratio amplitude B");
  m_runNumber=application()->options()->runNumber();
  return StatusCode::SUCCESS;
}

//=========================================================================
//  Virtual function execute
//=========================================================================
StatusCode ADCMeasurement::execute ( ) {
  // get current
  std::string val = m_device->cmd("Dev1", "ai");
  std::istringstream sv1(val);
  std::string v0, v1, tmp;
  sv1 >> v0;
  sv1 >> v1;
  //  m_ai0=toDouble(v0);
  //  m_ai1=toDouble(v1);
  m_ai0 = atof(v0.c_str());
  m_ai1 = atof(v1.c_str());
  std::cout << v0 << " " << v1 << " " << m_ai0 << " " << m_ai1 << std::endl;
  dataFill(1, m_ai0);
  dataFill(2, m_ai1);
  // set ramp value and channel A
  dataFill(3, m_ramp);
  m_ao0 = m_ramp;
  char cmd[80];
  sprintf(cmd,"ao %f %f", m_ramp, 0.);
  m_device->cmd("Dev1", cmd );
  val = m_device->cmd("Dev1", "di");
  std::istringstream sv2(val);
  sv2 >> tmp;
  sv2 >> tmp;
  sv2 >> tmp;
  sv2 >> v0 >> tmp >> v1 ;
  m_diA = atoi(v0.c_str());
  dataFill(4, m_diA);

  sprintf(cmd,"ao %f %f", m_ramp, 3.3);
  m_device->cmd("Dev1", cmd );
  val = m_device->cmd("Dev1", "di");
  std::istringstream sv3(val);
  sv3 >> tmp;
  sv3 >> tmp;
  sv3 >> tmp;
  sv3 >> v0 >> tmp >> v1 ;
  m_diB = atoi(v1.c_str());
  dataFill(5, m_diB);
  
  if (m_diA>m_baseline){ 
    dataFill(6, m_ramp/float(m_diA-m_baseline));
  }
  else {
    dataFill(6, 0.);
  }
  if (m_diB>m_baseline){ 
    dataFill(7, m_ramp/float(m_diB-m_baseline));
  }
  else {
    dataFill(7, 0.);
  }

  //  std::cout << "digit: " << m_ramp << " " << m_diA << " " << m_diB << std::endl;

  time_t t =time(0);
  m_duration=difftime(t, m_startTime);
  m_timestamp = t; 
  dataFill(0, m_duration);

  if (application()->options()->dataStorage()){
    m_evtNumber=application()->options()->nEvt()-1;
    m_tree->Fill();
  }
  wait(m_frequency);
  m_ramp += (m_maxRamp-m_minRamp)/float(m_stepRamp);
  if ( m_ramp > m_maxRamp ) { m_ramp = m_minRamp; }
  return StatusCode::SUCCESS;
}

//=========================================================================
//  Virtual function finalize
//=========================================================================
StatusCode ADCMeasurement::finalize ( ) {
  return StatusCode::SUCCESS;
}