PhaserRampExec.cpp

Go to the documentation of this file.

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

// Include files 
#include <stdio.h>
#include <stdlib.h>
#include <time.h>
#include <ctime>

// ROOT include files
#include "TROOT.h"
#include "TObject.h"
#include "TFile.h"
#include "TH1F.h"


// CATCore include files
#include "Element.h"
#include "Phaser.h"

// CATProc include files
#include "Histo1D.h"
#include "Plot.h"

// local
#include "PhaserRampExec.h"

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

//=============================================================================
// Standard constructor, initializes variables
//=============================================================================
PhaserRampExec::PhaserRampExec(  ) {
  setName ( "PhaserRampExec" );
  setType ( "Phaser" );
  setTitle( "Phaser I2C Configuration Test" );
}
//=============================================================================
// Destructor
//=============================================================================
PhaserRampExec::~PhaserRampExec () {}

//=========================================================================
//  virtual function intialize
//=========================================================================
StatusCode PhaserRampExec::initialize ( ) {

  info("" );
  info("____________________________");
  info("Phaser Configuration Control");
  info(" =>" + element()->path() + " " + 
       element()->name() + " " + element()->type() );
  
  m_phaser=dynamic_cast<Phaser*>( element() );
  
  srand( (unsigned)time( NULL ) );

  m_phaser->read();

  m_initPhase[0]=m_phaser->phase(0);
  m_initPhase[1]=m_phaser->phase(1);
  m_initPhase[2]=m_phaser->phase(2);
  m_initPhase[3]=m_phaser->phase(3);
  
  openRootFile ();

  hv1=new TH1D( "Delay1" , "Delay Channel 1" , 25 , 0. , 25. );
  hv2=new TH1D( "Delay2" , "Delay Channel 2" , 25 , 0. , 25. );
  hv3=new TH1D( "Delay3" , "Delay Channel 3" , 25 , 0. , 25. );
  hv4=new TH1D( "Delay4" , "Delay Channel 4" , 25 , 0. , 25. );

  h1=new TH1D( "Error1" , "Error Channel 1" , 25 , 0. , 25. );
  h2=new TH1D( "Error2" , "Error Channel 2" , 25 , 0. , 25. );
  h3=new TH1D( "Error3" , "Error Channel 3" , 25 , 0. , 25. );
  h4=new TH1D( "Error4" , "Error Channel 4" , 25 , 0. , 25. );
  
  m_nProcessedEvents=0;
  errorCount=0;
  startChrono();
  
  return StatusCode::SUCCESS;  
}

//=========================================================================
//  Virtual function execute
//=========================================================================
StatusCode PhaserRampExec::execute ( ) { 

  m_nProcessedEvents++;
  
  int v1 =  (int) ( (float)rand()/RAND_MAX * 25. );
  int v2 =  (int) ( (float)rand()/RAND_MAX * 25. );
  int v3 =  (int) ( (float)rand()/RAND_MAX * 25. );
  int v4 =  (int) ( (float)rand()/RAND_MAX * 25. );
  
  m_phaser->setPhase(v1,v2,v3,v4);
  m_phaser->write();
  m_phaser->read();

  int d0=m_phaser->phase(0);
  int d1=m_phaser->phase(1);
  int d2=m_phaser->phase(2);
  int d3=m_phaser->phase(3);
  
  hv1->Fill(v1);
  hv2->Fill(v2);
  hv3->Fill(v3);
  hv4->Fill(v4);

  
  if( v1!=d0 || v2!=d1 || v3!=d2 || v4!=d3 ){

    if ( d0 != v1 ) h1->Fill(int(d0));
    if ( d1 != v2 ) h2->Fill(int(d1));
    if ( d2 != v3 ) h3->Fill(int(d2));
    if ( d3 != v4 ) h4->Fill(int(d3));

    errorCount++;
    info("WRITE/READ INCOMPATIBILTY ["+
         itos(errorCount)+"] " +
         itos(v1)+"/"+itos(d0)+" "+
         itos(v2)+"/"+itos(d1)+" "+
         itos(v3)+"/"+itos(d2)+" "+
         itos(v4)+"/"+itos(d3),
         "Phaser::ramp");
    //     return false;
  }
  
  return StatusCode::SUCCESS;
}


//=========================================================================
//  Virtual function finalize
//=========================================================================
StatusCode PhaserRampExec::finalize ( ) {
  info("__" );
  info("Errors                             "+itos(errorCount));
  info("Processed Events                   "+itos(m_nProcessedEvents));
  info("Elapsed Time                       "+ftos(elapsedTime())+" s.");
  info("****************************" );
 
  m_phaser->setPhase(m_initPhase[0],
                     m_initPhase[1],
                     m_initPhase[2],
                     m_initPhase[3]);
  
  m_phaser->write();

  closeRootFile();

  hv1->Delete();
  hv2->Delete();
  hv3->Delete();
  hv4->Delete();
  h1->Delete();
  h2->Delete();
  h3->Delete();
  h4->Delete();
  
  return StatusCode::SUCCESS;  
}