CalibrationTask Class Reference

#include <CalibrationTask.h>

Inheritance diagram for CalibrationTask:

Public Member Functions
	CalibrationTask (DQMStore *, const FedChannelConnection &, const sistrip::RunType &, const char *filename, uint32_t run, const edm::EventSetup &setup)
virtual	~CalibrationTask ()
Public Member Functions inherited from CommissioningTask
void	bookHistograms ()
	CommissioningTask (DQMStore *, const FedChannelConnection &, const std::string &my_name)
void	eventSetup (const edm::EventSetup *)
const uint32_t &	fillCntr () const
void	fillHistograms (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &)
void	fillHistograms (const SiStripEventSummary &, const uint16_t &fed_id, const std::map< uint16_t, float > &fed_ch)
const std::string &	myName () const
const uint32_t &	updateFreq () const
void	updateFreq (const uint32_t &)
void	updateHistograms ()
virtual	~CommissioningTask ()

Private Member Functions
virtual void	book ()
void	checkAndSave (const uint16_t &)
void	directory (std::stringstream &, uint32_t run_number=0)
virtual void	fill (const SiStripEventSummary &, const edm::DetSet< SiStripRawDigi > &)
virtual void	update ()

Private Attributes
MonitorElement *	calchanElement_
std::vector< HistoSet >	calib_
std::string	filename_
uint16_t	lastCalChan_
uint16_t	nBins_
std::vector< uint16_t >	ped
uint32_t	run_
sistrip::RunType	runType_

Additional Inherited Members
Protected Member Functions inherited from CommissioningTask
const FedChannelConnection &	connection () const
DQMStore *const	dqm () const
const edm::EventSetup *const	eventSetup () const
const uint32_t &	fecKey () const
const uint32_t &	fedKey () const
void	updateHistoSet (HistoSet &, const uint32_t &bin, const float &value)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin, const short &value)
void	updateHistoSet (HistoSet &, const uint32_t &bin)
void	updateHistoSet (CompactHistoSet &, const uint32_t &bin)
void	updateHistoSet (HistoSet &, const float &value)
void	updateHistoSet (CompactHistoSet &)
void	updateHistoSet (HistoSet &)

Detailed Description

Definition at line 10 of file CalibrationTask.h.

Constructor & Destructor Documentation

CalibrationTask::CalibrationTask	(	DQMStore *	dqm,
		const FedChannelConnection &	conn,
		const sistrip::RunType &	rtype,
		const char *	filename,
		uint32_t	run,
		const edm::EventSetup &	setup
	)

Definition at line 20 of file CalibrationTask.cc.

References FedChannelConnection::apvPairNumber(), FedChannelConnection::detId(), edm::EventSetup::get(), SiStripPedestals::getPed(), SiStripPedestals::getRange(), createfilelist::int, LogDebug, ped, and digi_MixPreMix_cfi::strip.

                                              :
  CommissioningTask( dqm, conn, "CalibrationTask" ),
  runType_(rtype),
  nBins_(65),lastCalChan_(8),
  filename_(filename),
  run_(run)
{
  LogDebug("Commissioning") << "[CalibrationTask::CalibrationTask] Constructing object...";
  // load the pedestals
  edm::ESHandle<SiStripPedestals> pedestalsHandle;
  setup.get<SiStripPedestalsRcd>().get(pedestalsHandle);
  SiStripPedestals::Range detPedRange = pedestalsHandle->getRange(conn.detId());
  int start = conn.apvPairNumber()*256; 
  int stop  = start + 256;
  int value = 0;
  ped.reserve(256);
  for(int strip = start; strip < stop; ++strip) {
    value = int(pedestalsHandle->getPed(strip,detPedRange));
    if(value>895) value -= 1024;
    ped.push_back(value);
  }
}

CalibrationTask::~CalibrationTask ( )

virtual

Definition at line 50 of file CalibrationTask.cc.

References LogDebug.

                                  {
  LogDebug("Commissioning") << "[CalibrationTask::CalibrationTask] Destructing object...";
}

Member Function Documentation

void CalibrationTask::book ( )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 56 of file CalibrationTask.cc.

References sistrip::APV, DQMStore::bookInt(), calchanElement_, calib_, CommissioningTask::connection(), sistrip::DET_KEY, CommissioningTask::dqm(), sistrip::EXPERT_HISTO, spr::find(), DQMStore::get(), DQMStore::getMEs(), mps_fire::i, LogDebug, nBins_, alignCSCRings::pwd, DQMStore::pwd(), sistrip::root_, runType_, DQMStore::setCurrentFolder(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, and SiStripHistoTitle::title().

                           {
  LogDebug("Commissioning") << "[CalibrationTask::book]";

  // book 32 histograms, one for each strip in a calibration group on 2 APVs
  calib_.resize(32);
  for(int apv=0;apv<2;++apv) {
    for(int i=0;i<16;++i) {
      std::string title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                                             runType_, 
                                             sistrip::DET_KEY, 
                                             connection().detId(),
                                             sistrip::APV, 
                                             connection().i2cAddr(apv) ).title(); 
      std::stringstream complement;
      complement << "_" << i ;
      title += complement.str();            
      calib_[apv*16+i].histo( dqm()->book1D( title, title, nBins_, 0, 203.125) );
      calib_[apv*16+i].isProfile_=false;
      calib_[apv*16+i].vNumOfEntries_.resize(nBins_,0);
    }
  }

  // book the calchan values
  std::string pwd = dqm()->pwd();
  std::string rootDir = pwd.substr(0,pwd.find(std::string(sistrip::root_) + "/")+(sizeof(sistrip::root_) - 1));
  dqm()->setCurrentFolder( rootDir );
  std::vector<std::string> existingMEs = dqm()->getMEs();
  if(find(existingMEs.begin(),existingMEs.end(),"calchan")!=existingMEs.end()) {
    calchanElement_ = dqm()->get(dqm()->pwd()+"/calchan");
  } else {
    calchanElement_ = dqm()->bookInt("calchan");
  }
  dqm()->setCurrentFolder(pwd);
  
  LogDebug("Commissioning") << "[CalibrationTask::book] done";
  
}

void CalibrationTask::checkAndSave ( const uint16_t &  calchan )

private

Definition at line 127 of file CalibrationTask.cc.

References sistrip::APV, calchanElement_, calib_, CommissioningTask::connection(), sistrip::DET_KEY, dir, directory(), CommissioningTask::dqm(), sistrip::EXPERT_HISTO, cmsPerfPublish::fail(), filename_, MonitorElement::Fill(), mps_fire::i, recoMuon::in, lastCalChan_, LogTrace, dataset::name, nBins_, NULL, run_, runType_, DQMStore::save(), AlCaHLTBitMon_QueryRunRegistry::string, fftjetcommon_cfi::title, and update().

Referenced by fill().

                                                          {
  //for the first time
  if(lastCalChan_==8) {
    // set new parameter value
    lastCalChan_ = calchan;
    calchanElement_->Fill(lastCalChan_);
    return;
  }

  //check if CDRV register (ical) has changed
  // in that case, save the histograms, reset them, change the title and continue
  if(calchan!=lastCalChan_) {
    
    // set histograms before saving
    update(); //TODO: we must do this for all tasks, otherwise only the first is updated.

    // change the title
    for(int apv=0;apv<2;++apv) {
      for(int i=0;i<16;++i) {
        std::stringstream complement;
        complement << "STRIP_" << (apv*128+lastCalChan_+i*8) ;
        std::string title = SiStripHistoTitle( sistrip::EXPERT_HISTO, 
                                               runType_, 
                                               sistrip::DET_KEY, 
                                               connection().detId(),
                                               sistrip::APV, 
                                               connection().i2cAddr(apv),
                                               complement.str() ).title(); 
        calib_[apv*16+i].histo()->setTitle(title);
    calib_[apv*16+i].histo()->setAxisTitle(title);
      }
    }

    // Strip filename of ".root" extension
    std::string name;
    if ( filename_.find(".root",0) == std::string::npos ) { name = filename_; }
    else { name = filename_.substr( 0, filename_.find(".root",0) ); }
  
    // Retrieve SCRATCH directory
    std::string scratch = "SCRATCH"; //@@ remove trailing slash!!!
    std::string dir = "";
    if ( getenv(scratch.c_str()) != NULL ) {
      dir = getenv(scratch.c_str());
    }
  
    // Save file with appropriate filename 
    std::stringstream ss;
    if ( !dir.empty() ) { ss << dir << "/"; }
    else { ss << "/tmp/"; }
    ss << name << "_";
    directory(ss,run_); // Add filename with run number, host ip, pid and .root extension
    ss << "_CALCHAN" << lastCalChan_; // Add CalChan value
    ss << "_000"; // Add FU instance number (fake)
    ss << ".root"; // Append ".root" extension
    if ( !filename_.empty() ) {
      if(std::ifstream(ss.str().c_str(),std::ifstream::in).fail()) { //save only once. Skip if the file already exist
        dqm()->save( ss.str() ); 
        LogTrace("DQMsource")
          << "[SiStripCommissioningSource::" << __func__ << "]"
          << " Saved all histograms to file \""
          << ss.str() << "\"";
      } else {
        LogTrace("DQMsource")
          << "[SiStripCommissioningSource::" << __func__ << "]"
          << " Skipping creation of file \""
          << ss.str() << "\" that already exists" ;
      }
    } else {
      edm::LogWarning("DQMsource")
        << "[SiStripCommissioningSource::" << __func__ << "]"
        << " NULL value for filename! No root file saved!";
    }
  
    // reset
    for(std::vector<HistoSet>::iterator it=calib_.begin();it<calib_.end();++it) {
      it->vNumOfEntries_.clear();
      it->vNumOfEntries_.resize(nBins_,0);
    }
    // set new parameter values
    lastCalChan_ = calchan;
    calchanElement_->Fill(lastCalChan_);
  }
}

void CalibrationTask::directory ( std::stringstream &  dir, uint32_t  run_number = 0  )

private

Definition at line 213 of file CalibrationTask.cc.

References hcalSimParameters_cfi::he, sysUtil::pid, AlCaHLTBitMon_QueryRunRegistry::string, and tmp.

Referenced by checkAndSave().

                                              {

  // Get details about host
  char hn[256];
  gethostname( hn, sizeof(hn) );
  struct hostent* he;
  he = gethostbyname(hn);

  // Extract host name and ip
  std::string host_name;
  std::string host_ip;
  if ( he ) { 
    host_name = std::string(he->h_name);
    host_ip = std::string( inet_ntoa( *(struct in_addr*)(he->h_addr) ) );
  } else {
    host_name = "unknown.cern.ch";
    host_ip = "255.255.255.255";
  }

  // Reformat IP address
  std::string::size_type pos = 0;
  std::stringstream ip;
  //for ( uint16_t ii = 0; ii < 4; ++ii ) {
  while ( pos != std::string::npos ) {
    std::string::size_type tmp = host_ip.find(".",pos);
    if ( tmp != std::string::npos ) {
      ip << std::setw(3)
     << std::setfill('0')
     << host_ip.substr( pos, tmp-pos ) 
     << ".";
      pos = tmp+1; // skip the delimiter "."
    } else {
      ip << std::setw(3)
     << std::setfill('0')
     << host_ip.substr( pos );
      pos = std::string::npos;
    }
  }
  
  // Get pid
  pid_t pid = getpid();

  // Construct string
  dir << std::setw(8) 
      << std::setfill('0') 
      << run_number
      << "_" 
      << ip.str()
      << "_"
      << std::setw(5) 
      << std::setfill('0') 
      << pid;
}

void CalibrationTask::fill ( const SiStripEventSummary &  summary, const edm::DetSet< SiStripRawDigi > &  digis  )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 96 of file CalibrationTask.cc.

References stringResolutionProvider_cfi::bin, SiStripEventSummary::calChan(), calib_, SiStripEventSummary::calSel(), checkAndSave(), edm::DetSet< T >::data, gen::k, SiStripEventSummary::latency(), LogDebug, ped, update(), and CommissioningTask::updateHistoSet().

                                                         {
  LogDebug("Commissioning") << "[CalibrationTask::fill]";
  // Check if CalChan changed. In that case, save, reset histo, change title, and continue
  int isub,ical = summary.calChan();
  isub = ical<4 ? ical+4 : ical-4;
  checkAndSave(ical);
  // retrieve the delay from the EventSummary
  int bin = (100-summary.latency())*8+(7-summary.calSel());
  // Fill the histograms.
  // data-ped -(data-ped)_isub
  // the second term corresponds to the common mode substraction, looking at a strip far away.
  for (int apv=0; apv<2; ++apv) {
    for (int k=0;k<16;++k)  {
      updateHistoSet( calib_[apv*16+k],bin,digis.data[apv*128+ical+k*8].adc()-ped[apv*128+ical+k*8]-(digis.data[apv*128+isub+k*8].adc()-ped[apv*128+isub+k*8]));
    }
  }
  update(); //TODO: temporary: find a better solution later
}

void CalibrationTask::update ( )

privatevirtual

Reimplemented from CommissioningTask.

Definition at line 118 of file CalibrationTask.cc.

References calib_, and CommissioningTask::updateHistoSet().

Referenced by progressbar.ProgressBar::__next__(), MatrixUtil.Matrix::__setitem__(), MatrixUtil.Steps::__setitem__(), Vispa.Gui.VispaWidget.VispaWidget::autosize(), checkAndSave(), Vispa.Views.LineDecayView.LineDecayContainer::createObject(), Vispa.Views.LineDecayView.LineDecayContainer::deselectAllObjects(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::deselectAllWidgets(), Vispa.Gui.VispaWidget.VispaWidget::enableAutosizing(), fill(), dqm-mbProfile.Profile::finish(), progressbar.ProgressBar::finish(), Vispa.Gui.MenuWidget.MenuWidget::leaveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseMoveEvent(), Vispa.Gui.MenuWidget.MenuWidget::mouseMoveEvent(), Vispa.Views.LineDecayView.LineDecayContainer::mouseMoveEvent(), Vispa.Gui.VispaWidgetOwner.VispaWidgetOwner::mouseReleaseEvent(), Vispa.Views.LineDecayView.LineDecayContainer::objectMoved(), MatrixUtil.Steps::overwrite(), Vispa.Views.LineDecayView.LineDecayContainer::removeObject(), Vispa.Gui.ConnectableWidget.ConnectableWidget::removePorts(), Vispa.Gui.FindDialog.FindDialog::reset(), Vispa.Gui.PortConnection.PointToPointConnection::select(), Vispa.Gui.VispaWidget.VispaWidget::select(), Vispa.Views.LineDecayView.LineDecayContainer::select(), Vispa.Gui.VispaWidget.VispaWidget::setText(), Vispa.Gui.VispaWidget.VispaWidget::setTitle(), Vispa.Gui.ZoomableWidget.ZoomableWidget::setZoom(), Vispa.Views.LineDecayView.LineDecayContainer::setZoom(), and Vispa.Gui.PortConnection.PointToPointConnection::updateConnection().

                             {
  // LogDebug("Commissioning") << "[CalibrationTask::update]"; // huge output
  for(std::vector<HistoSet>::iterator it=calib_.begin();it<calib_.end();++it) {
    updateHistoSet( *it );
  }
}

Member Data Documentation

MonitorElement* CalibrationTask::calchanElement_

private

Definition at line 37 of file CalibrationTask.h.

Referenced by book(), and checkAndSave().

std::vector<HistoSet> CalibrationTask::calib_

private

Definition at line 30 of file CalibrationTask.h.

Referenced by book(), checkAndSave(), fill(), and update().

std::string CalibrationTask::filename_

private

Definition at line 34 of file CalibrationTask.h.

Referenced by checkAndSave().

uint16_t CalibrationTask::lastCalChan_

private

Definition at line 33 of file CalibrationTask.h.

Referenced by checkAndSave().

uint16_t CalibrationTask::nBins_

private

Definition at line 32 of file CalibrationTask.h.

Referenced by book(), and checkAndSave().

std::vector<uint16_t> CalibrationTask::ped

private

Definition at line 35 of file CalibrationTask.h.

Referenced by CalibrationTask(), and fill().

uint32_t CalibrationTask::run_

private

Definition at line 36 of file CalibrationTask.h.

Referenced by checkAndSave().

sistrip::RunType CalibrationTask::runType_

private

Definition at line 28 of file CalibrationTask.h.

Referenced by book(), and checkAndSave().