PixTestDacScanCurrent.cc

// -- author: Daniel Pitzl
// scan a DAC, measure currents

#include <stdlib.h>  // atof, atoi
#include <algorithm> // std::find

#include <TStopwatch.h>

#include "PixTestDacScanCurrent.hh"
#include "log.h"

using namespace std;
using namespace pxar;

ClassImp(PixTestDacScanCurrent)

//------------------------------------------------------------------------------
PixTestDacScanCurrent::PixTestDacScanCurrent( PixSetup *a, std::string name )
: PixTest(a, name), fParDAC("nada")
{
  PixTest::init();
  init();
}

//------------------------------------------------------------------------------
PixTestDacScanCurrent::PixTestDacScanCurrent() : PixTest()
{
}

//------------------------------------------------------------------------------
bool PixTestDacScanCurrent::setParameter( string parName, string sval )
{
  bool found(false);

  for( uint32_t i = 0; i < fParameters.size(); ++i ) {

    if( fParameters[i].first == parName ) {

      found = true;

      sval.erase( remove( sval.begin(), sval.end(), ' ' ), sval.end() );

      if( !parName.compare( "dac" ) ) {
    fParDAC = sval;
    LOG(logDEBUG) << "PixTestDacScanCurrent setting fParDAC  ->" << fParDAC
             << "<- from sval = " << sval;
      }
      break;
    }
  }
  return found;
}

//------------------------------------------------------------------------------
void PixTestDacScanCurrent::init()
{
  fDirectory = gFile->GetDirectory( fName.c_str() );
  if( !fDirectory )
    fDirectory = gFile->mkdir( fName.c_str() );
  fDirectory->cd();
}

// ----------------------------------------------------------------------
void PixTestDacScanCurrent::setToolTips()
{
  fTestTip = string( "measure ROC currents vs DAC");
  fSummaryTip = string("summary plot to be implemented");
}

//------------------------------------------------------------------------------
void PixTestDacScanCurrent::bookHist(string name) // general booking routine
{
  fDirectory->cd();

  // Improvement: change binning for 4-bit DACs...

  TH1D *h1(0);
  fHistList.clear();
  for( uint32_t i = 0; i < fPixSetup->getConfigParameters()->getNrocs(); ++i ) {

    h1 = new TH1D( Form( "Ia_vs_%s_C%d", name.c_str(), i),
           Form( "Ia vs %s C%d", name.c_str(), i),
           256, 0, 256 );
    h1->SetMinimum(0);
    h1->SetStats(0);
    setTitles( h1, Form( "%s [DAC]", name.c_str() ), "analog current [mA]" );
    fHistList.push_back(h1);

    h1 = new TH1D( Form( "Id_vs_%s_C%d", name.c_str(), i),
           Form( "Id vs %s C%d", name.c_str(), i),
           256, 0, 256 );
    h1->SetMinimum(0);
    h1->SetStats(0);
    setTitles( h1, Form( "%s [DAC]", name.c_str() ), "digital current [mA]" );
    fHistList.push_back(h1);
  }
}

//------------------------------------------------------------------------------
PixTestDacScanCurrent::~PixTestDacScanCurrent()
{
  LOG(logDEBUG) << "PixTestDacScanCurrent dtor";
  std::list<TH1*>::iterator il;
  fDirectory->cd();
  for( il = fHistList.begin(); il != fHistList.end(); ++il) {
    LOG(logDEBUG) << "Write out " << (*il)->GetName();
    (*il)->SetDirectory(fDirectory);
    (*il)->Write();
  }
}

//------------------------------------------------------------------------------
void PixTestDacScanCurrent::doTest()
{
  fDirectory->cd();
  PixTest::update();

  LOG(logINFO) << "PixTestDacScanCurrent::doTest() DAC = " << fParDAC;

  uint8_t maxDac = fApi->getDACRange( fParDAC ); // 15 or 255

  if( maxDac == 0 ) {
    LOG(logINFO) << "ERROR: " << fParDAC << " is not a ROC register";
    return;
  }

  fApi->_dut->testAllPixels(false);

  bookHist( fParDAC );

  TH1D *hia(0);
  TH1D *hid(0);
  TStopwatch sw;

  size_t nRocs = fPixSetup->getConfigParameters()->getNrocs();

  for( uint32_t roc = 0; roc < nRocs; ++roc ) {

    hia = (TH1D*)fDirectory->Get( Form( "Ia_vs_%s_C%d", fParDAC.c_str(), roc ) );
    hid = (TH1D*)fDirectory->Get( Form( "Id_vs_%s_C%d", fParDAC.c_str(), roc ) );

    if( hia && hid ) {

      // remember DAC

      uint8_t dacval = fApi->_dut->getDAC( roc, fParDAC );

      fApi->setDAC( fParDAC, 0, roc ); // start at zero

      // delay:

      sw.Start(kTRUE); // reset stopwatch
      do {
    sw.Start(kFALSE); // continue
    fApi->getTBia();
      }
      while( sw.RealTime() < 0.1 ); // [s]

      // loop over DAC:

      for( uint32_t idac = 0; idac <= maxDac; ++idac ) {

    fApi->setDAC( fParDAC, idac, roc );
    // delay?
    hia->SetBinContent( idac+1, fApi->getTBia()*1E3 );
    hid->SetBinContent( idac+1, fApi->getTBid()*1E3 );

      }

      fApi->setDAC( fParDAC, dacval, roc ); // restore
    }
    else {
      LOG(logINFO) << "XX did not find "
           << Form( "Ia_vs_%s_C%d", fParDAC.c_str(), roc );
    }
    fDisplayedHist = find( fHistList.begin(), fHistList.end(), hid );
    if( hid ) hid->Draw();
    PixTest::update();

  } // rocs

}