CMS 3D CMS Logo

Public Types | Public Member Functions | Private Member Functions | Private Attributes

SiStripCommissioningSource Class Reference

#include <SiStripCommissioningSource.h>

Inheritance diagram for SiStripCommissioningSource:
edm::EDAnalyzer

List of all members.

Public Types

typedef std::map< unsigned int,
CommissioningTask * > 
TaskMap
typedef std::vector
< CommissioningTask * > 
VecOfTasks
typedef std::vector< VecOfTasksVecOfVecOfTasks

Public Member Functions

void analyze (const edm::Event &, const edm::EventSetup &)
void beginRun (edm::Run const &, const edm::EventSetup &)
void endJob ()
 SiStripCommissioningSource (const edm::ParameterSet &)
 ~SiStripCommissioningSource ()

Private Member Functions

void clearCablingTasks ()
void clearTasks ()
void createCablingTasks ()
void createRunNumber ()
void createTask (const SiStripEventSummary *const, const edm::EventSetup &)
void createTasks (sistrip::RunType, const edm::EventSetup &)
void directory (std::stringstream &, uint32_t run_number=0)
DQMStore *const dqm (std::string method="") const
void fillCablingHistos (const SiStripEventSummary *const, const edm::DetSetVector< SiStripRawDigi > &)
void fillHistos (const SiStripEventSummary *const, const edm::DetSetVector< SiStripRawDigi > &)
void remove ()
 SiStripCommissioningSource ()

Private Attributes

std::string base_
bool cablingTask_
TaskMap cablingTasks_
DQMStoredqm_
SiStripFecCablingfecCabling_
SiStripFedCablingfedCabling_
std::string filename_
std::string inputModuleLabel_
std::string inputModuleLabelSummary_
edm::ParameterSet parameters_
uint32_t run_
sistrip::RunType task_
std::string taskConfigurable_
VecOfVecOfTasks tasks_
bool tasksExist_
int32_t time_
int updateFreq_
std::string view_

Detailed Description

Definition at line 26 of file SiStripCommissioningSource.h.


Member Typedef Documentation

typedef std::map<unsigned int, CommissioningTask*> SiStripCommissioningSource::TaskMap

Map of task objects, identified through FedChanelId

Definition at line 31 of file SiStripCommissioningSource.h.

Definition at line 32 of file SiStripCommissioningSource.h.

Definition at line 33 of file SiStripCommissioningSource.h.


Constructor & Destructor Documentation

SiStripCommissioningSource::SiStripCommissioningSource ( const edm::ParameterSet pset)

Definition at line 55 of file SiStripCommissioningSource.cc.

References LogTrace, sistrip::mlDqmSource_, and tasks_.

                                                                                    :
  dqm_(0),
  fedCabling_(0),
  fecCabling_(0),
  inputModuleLabel_( pset.getParameter<std::string>( "InputModuleLabel" ) ),
  inputModuleLabelSummary_( pset.getParameter<std::string>( "SummaryInputModuleLabel" ) ),
  filename_( pset.getUntrackedParameter<std::string>("RootFileName",sistrip::dqmSourceFileName_) ),
  run_(0),
  time_(0),
  taskConfigurable_( pset.getUntrackedParameter<std::string>("CommissioningTask","UNDEFINED") ),
  task_(sistrip::UNDEFINED_RUN_TYPE),
  tasks_(),
  cablingTasks_(),
  tasksExist_(false),
  cablingTask_(false),
  updateFreq_( pset.getUntrackedParameter<int>("HistoUpdateFreq",1) ),
  base_(""),
  view_( pset.getUntrackedParameter<std::string>("View", "Default") ),
  parameters_(pset)
{
  LogTrace(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Constructing object...";
  tasks_.clear();
  tasks_.resize( 1024, VecOfTasks(96,static_cast<CommissioningTask*>(0)) ); 
}
SiStripCommissioningSource::~SiStripCommissioningSource ( )

Definition at line 84 of file SiStripCommissioningSource.cc.

References LogTrace, and sistrip::mlDqmSource_.

                                                        {
  LogTrace(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Destructing object...";
}
SiStripCommissioningSource::SiStripCommissioningSource ( ) [private]

Private default constructor.


Member Function Documentation

void SiStripCommissioningSource::analyze ( const edm::Event event,
const edm::EventSetup setup 
) [virtual]

Implements edm::EDAnalyzer.

Definition at line 258 of file SiStripCommissioningSource.cc.

References sistrip::APV_LATENCY, sistrip::APV_TIMING, cablingTask_, sistrip::CALIBRATION, sistrip::CALIBRATION_DECO, sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, createRunNumber(), createTask(), sistrip::DAQ_SCOPE_MODE, edm::EventID::event(), sistrip::FAST_CABLING, sistrip::FED_CABLING, sistrip::FED_SCOPE_MODE, sistrip::FED_TIMING, sistrip::FED_VIRGIN_RAW, SiStripEnumsAndStrings::fedReadoutMode(), fillCablingHistos(), fillHistos(), sistrip::FINE_DELAY, edm::EventBase::id(), inputModuleLabel_, inputModuleLabelSummary_, sistrip::mlDqmSource_, sistrip::NOISE, NULL, sistrip::OPTO_SCAN, sistrip::PEDESTALS, sistrip::PEDS_FULL_NOISE, sistrip::PEDS_ONLY, edm::Handle< T >::product(), edm::EventID::run(), run_, environment_file_cfi::runType, HcalObjRepresent::setup(), summarizeEdmComparisonLogfiles::summary, task_, tasksExist_, cond::rpcobgas::time, time_, sistrip::UNDEFINED_RUN_TYPE, sistrip::UNKNOWN_RUN_TYPE, updateFreq_, and sistrip::VPSP_SCAN.

                                                                       {
  // Retrieve commissioning information from "event summary" 
  edm::Handle<SiStripEventSummary> summary;
  event.getByLabel( inputModuleLabelSummary_, summary );

  // Check if EventSummary has info attached
  if ( ( summary->runType() == sistrip::UNDEFINED_RUN_TYPE ||
         summary->runType() == sistrip::UNKNOWN_RUN_TYPE ) &&
       summary->nullParams() ) {
    edm::LogWarning(mlDqmSource_)
      << "[SiStripCommissioningSource::" << __func__ << "]"
      << " Unknown/undefined RunType and NULL parameter values!"
      << " It may be that the 'trigger FED' object was not found!"; 
  }

  // Check if need to rebuild FED/FEC cabling objects for connection run
  //cablingForConnectionRun( summary->runType() ); //@@ do not use!
  
  // Extract run number and forward to client
  if ( event.id().run() != run_ ) { 
    run_ = event.id().run(); 
    createRunNumber();
  }
  
  // Coarse event rate counter
  if ( !(event.id().event()%updateFreq_) ) {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " The last " << updateFreq_ 
       << " events were processed at a rate of ";
    if ( time(NULL) == time_ ) { ss << ">" << updateFreq_ << " Hz"; }
    else { ss << (updateFreq_/(time(NULL)-time_)) << " Hz"; }
    edm::LogVerbatim(mlDqmSource_) << ss.str();
    time_ = time(NULL);
  }
  
  // Create commissioning task objects 
  if ( !tasksExist_ ) { createTask( summary.product(), setup ); }

  // Retrieve raw digis with mode appropriate to task 
  edm::Handle< edm::DetSetVector<SiStripRawDigi> > raw;
  if ( task_ == sistrip::DAQ_SCOPE_MODE ) { 
    if ( summary->fedReadoutMode() == FED_VIRGIN_RAW ) {
      event.getByLabel( inputModuleLabel_, "VirginRaw", raw );
    } else if ( summary->fedReadoutMode() == FED_SCOPE_MODE ) {
      event.getByLabel( inputModuleLabel_, "ScopeMode", raw );
    } else {
      std::stringstream ss;
      ss << "[SiStripCommissioningSource::" << __func__ << "]"
         << " Requested DAQ_SCOPE_MODE but unknown FED"
         << " readout mode retrieved from SiStripEventSummary: " 
         << SiStripEnumsAndStrings::fedReadoutMode( summary->fedReadoutMode() );
      edm::LogWarning(mlDqmSource_) << ss.str();
    }
  } else if ( task_ == sistrip::FAST_CABLING ||
              task_ == sistrip::FED_CABLING ||
              task_ == sistrip::APV_TIMING ||
              task_ == sistrip::FED_TIMING ||
              task_ == sistrip::OPTO_SCAN ) { 
    event.getByLabel( inputModuleLabel_, "ScopeMode", raw );
  } else if ( task_ == sistrip::VPSP_SCAN ||
              task_ == sistrip::CALIBRATION ||
              task_ == sistrip::CALIBRATION_DECO ||
              task_ == sistrip::CALIBRATION_SCAN ||
              task_ == sistrip::CALIBRATION_SCAN_DECO ||
              task_ == sistrip::PEDESTALS ||
              task_ == sistrip::PEDS_ONLY ||
              task_ == sistrip::NOISE ||
              task_ == sistrip::PEDS_FULL_NOISE ) {
    event.getByLabel( inputModuleLabel_, "VirginRaw", raw );
  } else if ( task_ == sistrip::APV_LATENCY ||
              task_ == sistrip::FINE_DELAY ) {
    event.getByLabel( inputModuleLabel_, "FineDelaySelection", raw );
  } else {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Unknown CommissioningTask: " 
       << SiStripEnumsAndStrings::runType( task_ )
       << " Unable to establish FED readout mode and retrieve digi container!"
       << " Check if SiStripEventSummary object is found/present in Event";
    edm::LogWarning(mlDqmSource_) << ss.str();
    return;
  }

  // Check for NULL pointer to digi container
  if ( &(*raw) == 0 ) {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]" << std::endl
       << " NULL pointer to DetSetVector!" << std::endl
       << " Unable to fill histograms!";
    edm::LogWarning(mlDqmSource_) << ss.str();
    return;
  }
  
  if ( !cablingTask_ ) { fillHistos( summary.product(), *raw );  }
  else { fillCablingHistos( summary.product(), *raw ); }
  
}
void SiStripCommissioningSource::beginRun ( edm::Run const &  run,
const edm::EventSetup setup 
) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 107 of file SiStripCommissioningSource.cc.

References base_, cablingTask_, clearCablingTasks(), clearTasks(), SiStripFecCabling::crates(), dir, dqm(), dqm_, fecCabling_, fedCabling_, SiStripFedCabling::feds(), edm::EventSetup::get(), LogDebug, LogTrace, sistrip::mlDqmSource_, cmsCodeRules::cppFunctionSkipper::operator, edm::ESHandle< T >::product(), DQMStore::setVerbose(), task_, tasksExist_, and sistrip::UNDEFINED_RUN_TYPE.

                                                                                           {
  LogTrace(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Configuring..." << std::endl;
  
  // ---------- DQM back-end interface ----------
  
  dqm_ = edm::Service<DQMStore>().operator->();
  edm::LogInfo(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " DQMStore service: " 
    << dqm_;
  dqm(__func__);
  dqm()->setVerbose(0);
  
  // ---------- Base directory ----------

  std::stringstream dir(""); 
  base_ = dir.str();
  
  // ---------- FED and FEC cabling ----------
  
  edm::ESHandle<SiStripFedCabling> fed_cabling;
  setup.get<SiStripFedCablingRcd>().get( fed_cabling ); 
  fedCabling_ = const_cast<SiStripFedCabling*>( fed_cabling.product() ); 
  LogDebug(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << "Initialized FED cabling. Number of FEDs is " << fedCabling_->feds().size();
  fecCabling_ = new SiStripFecCabling( *fed_cabling );
  if ( fecCabling_->crates().empty() ) {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Empty std::vector returned by FEC cabling object!" 
       << " Check if database connection failed...";
    edm::LogWarning(mlDqmSource_) << ss.str();
  }

  // ---------- Reset ---------- 

  tasksExist_ = false;
  task_ = sistrip::UNDEFINED_RUN_TYPE;
  cablingTask_ = false;
  
  remove();

  clearCablingTasks();
  clearTasks();
  
}
void SiStripCommissioningSource::clearCablingTasks ( ) [private]

Definition at line 1100 of file SiStripCommissioningSource.cc.

References cablingTasks_.

Referenced by beginRun().

                                                   {
  if ( cablingTasks_.empty() ) { return; }
  for ( TaskMap::iterator itask = cablingTasks_.begin(); itask != cablingTasks_.end(); itask++ ) { 
    if ( itask->second ) { delete itask->second; }
  }
  cablingTasks_.clear();
}
void SiStripCommissioningSource::clearTasks ( ) [private]

Definition at line 1110 of file SiStripCommissioningSource.cc.

References tasks_.

Referenced by beginRun().

                                            {
  if ( tasks_.empty() ) { return; }
  VecOfVecOfTasks::iterator ifed = tasks_.begin();
  for ( ; ifed !=  tasks_.end(); ifed++ ) { 
    VecOfTasks::iterator ichan = ifed->begin();
    for ( ; ichan != ifed->end(); ichan++ ) { 
      if ( *ichan ) { delete *ichan; *ichan = 0; }
    }
    ifed->resize(96,0);
  }
  tasks_.resize(1024);
}
void SiStripCommissioningSource::createCablingTasks ( ) [private]

Definition at line 789 of file SiStripCommissioningSource.cc.

References base_, cablingTasks_, SiStripFecCabling::crates(), dir, dqm(), fecCabling_, sistrip::FED_CABLING, combine::key, sistrip::mlDqmSource_, getHLTPrescaleColumns::path, environment_file_cfi::runType, DQMStore::setCurrentFolder(), task_, sistrip::UNDEFINED_RUN_TYPE, and sistrip::UNKNOWN_RUN_TYPE.

Referenced by createTask().

                                                    {
  
  // Iterate through FEC cabling and create commissioning task objects
  uint16_t booked = 0;
  for ( std::vector<SiStripFecCrate>::const_iterator icrate = fecCabling_->crates().begin(); icrate != fecCabling_->crates().end(); icrate++ ) {
    for ( std::vector<SiStripFec>::const_iterator ifec = icrate->fecs().begin(); ifec != icrate->fecs().end(); ifec++ ) {
      for ( std::vector<SiStripRing>::const_iterator iring = ifec->rings().begin(); iring != ifec->rings().end(); iring++ ) {
        for ( std::vector<SiStripCcu>::const_iterator iccu = iring->ccus().begin(); iccu != iring->ccus().end(); iccu++ ) {
          for ( std::vector<SiStripModule>::const_iterator imodule = iccu->modules().begin(); imodule != iccu->modules().end(); imodule++ ) {
              
            // Build FEC key
            SiStripFecKey path( icrate->fecCrate(), 
                                ifec->fecSlot(), 
                                iring->fecRing(), 
                                iccu->ccuAddr(), 
                                imodule->ccuChan() );
            
            // Check if FEC key is invalid
            if ( !path.isValid() ) { continue; }
            
            // Set working directory prior to booking histograms
            std::string dir = base_ + path.path();
            dqm()->setCurrentFolder( dir );
            
            // Iterate through all APV pairs for this module
            for ( uint16_t ipair = 0; ipair < imodule->nApvPairs(); ipair++ ) {
                
              // Retrieve active APV devices
              SiStripModule::PairOfU16 apvs = imodule->activeApvPair( imodule->lldChannel(ipair) );
              
              // Create connection object to hold all relevant info
              FedChannelConnection conn( icrate->fecCrate(), 
                                         ifec->fecSlot(), 
                                         iring->fecRing(), 
                                         iccu->ccuAddr(), 
                                         imodule->ccuChan(),
                                         apvs.first,
                                         apvs.second,
                                         imodule->dcuId(),
                                         imodule->detId(),
                                         imodule->nApvPairs() );
              
              // Define key encoding control path 
              uint32_t key = SiStripFecKey( icrate->fecCrate(), 
                                            ifec->fecSlot(), 
                                            iring->fecRing(), 
                                            iccu->ccuAddr(), 
                                            imodule->ccuChan(),
                                            imodule->lldChannel(ipair) ).key();
                
              // Check key is non zero
              if ( !key ) { 
                edm::LogWarning(mlDqmSource_)
                  << "[SiStripCommissioningSource::" << __func__ << "]"
                  << " Unexpected NULL value for FEC key!";
                continue; 
              }
                
              // Create cabling task objects if not already existing
              if ( cablingTasks_.find( key ) == cablingTasks_.end() ) {
                
                if ( task_ == sistrip::FED_CABLING ) { 
                  cablingTasks_[key] = new FedCablingTask( dqm(), conn ); 
                } else if ( task_ == sistrip::UNDEFINED_RUN_TYPE ) { 
                  edm::LogWarning(mlDqmSource_)
                    << "[SiStripCommissioningSource::" << __func__ << "]"
                    << " Undefined CommissioningTask" 
                    << " Unable to create FedCablingTask object!";
                } else if ( task_ == sistrip::UNKNOWN_RUN_TYPE ) { 
                  edm::LogWarning(mlDqmSource_)
                    << "[SiStripCommissioningSource::" << __func__ << "]"
                    << " Unknown CommissioningTask" 
                    << " Unable to create FedCablingTask object!";
                } else { 
                  edm::LogWarning(mlDqmSource_)
                    << "[SiStripCommissioningSource::" << __func__ << "]"
                    << " Unexpected CommissioningTask: " 
                    << SiStripEnumsAndStrings::runType( task_ )
                    << " Unable to create FedCablingTask object!";
                }
                
                // Check if key is found and, if so, book histos and set update freq
                if ( cablingTasks_.find( key ) != cablingTasks_.end() ) {
                  if ( cablingTasks_[key] ) {
                    cablingTasks_[key]->bookHistograms(); 
                    cablingTasks_[key]->updateFreq(1); //@@ hardwired to update every event!!! 
                    booked++;
                    //std::stringstream ss;
                    //ss << "[SiStripCommissioningSource::" << __func__ << "]"
                    //<< " Booking histograms for '" << cablingTasks_[key]->myName()
                    //<< "' object with key 0x" << std::hex << std::setfill('0') << std::setw(8) << key << std::dec
                    //<< " in directory \"" << dir << "\"";
                    //LogTrace(mlDqmSource_) << ss.str();
                  } else {
                    std::stringstream ss;
                    ss << "[SiStripCommissioningSource::" << __func__ << "]"
                       << " NULL pointer to CommissioningTask for key 0x"
                       << std::hex << std::setfill('0') << std::setw(8) << key << std::dec
                       << " in directory " << dir 
                       << " Unable to book histograms!";
                    edm::LogWarning(mlDqmSource_) << ss.str();
                  }
                } else {
                  std::stringstream ss;
                  ss << "[SiStripCommissioningSource::" << __func__ << "]"
                     << " Unable to find CommissioningTask for key 0x"
                     << std::hex << std::setfill('0') << std::setw(8) << key << std::dec
                     << " in directory " << dir
                     << " Unable to book histograms!";
                  edm::LogWarning(mlDqmSource_) << ss.str();
                }
              
              } else {
                std::stringstream ss;
                ss << "[SiStripCommissioningSource::" << __func__ << "]"
                   << " CommissioningTask object already exists for key 0x"
                   << std::hex << std::setfill('0') << std::setw(8) << key << std::dec
                   << " in directory " << dir 
                   << " Unable to create FedCablingTask object!";
                edm::LogWarning(mlDqmSource_) << ss.str();
              }
              
            } // loop through apv pairs
          } // loop through modules
        } // loop through ccus
      } // loop through rings
    } // loop through fecs
  } // loop through crates

  edm::LogVerbatim(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Created " << booked 
    << " CommissioningTask objects and booked histograms";
  
}
void SiStripCommissioningSource::createRunNumber ( ) [private]

Definition at line 680 of file SiStripCommissioningSource.cc.

References base_, DQMStore::bookString(), dqm(), sistrip::mlDqmSource_, sistrip::root_, DTTTrigCorrFirst::run, run_, sistrip::runNumber_, sistrip::sep_, and DQMStore::setCurrentFolder().

Referenced by analyze().

                                                 {
  
  // Set commissioning task to default ("undefined") value
  if ( !run_ ) {
    edm::LogWarning(mlDqmSource_) 
      << "[SiStripCommissioningSource::" << __func__ << "]"
      << " NULL run number!";
    return;
  }
  
  // Create MonitorElement that identifies run number
  dqm()->setCurrentFolder( base_ + sistrip::root_ );
  std::stringstream run;
  run << run_;
  dqm()->bookString( std::string(sistrip::runNumber_) + sistrip::sep_ + run.str(), run.str() ); 
  
}
void SiStripCommissioningSource::createTask ( const SiStripEventSummary * const  summary,
const edm::EventSetup setup 
) [private]

Definition at line 700 of file SiStripCommissioningSource.cc.

References base_, DQMStore::bookString(), cablingTask_, createCablingTasks(), createTasks(), dqm(), sistrip::FED_CABLING, LogTrace, sistrip::mlDqmSource_, sistrip::root_, environment_file_cfi::runType, SiStripEventSummary::runType(), SiStripEnumsAndStrings::runType(), sistrip::sep_, DQMStore::setCurrentFolder(), task_, taskConfigurable_, sistrip::taskId_, tasksExist_, sistrip::UNDEFINED_RUN_TYPE, and sistrip::UNKNOWN_RUN_TYPE.

Referenced by analyze().

                                                                                                                  {
  
  // Set commissioning task to default ("undefined") value
  task_ = sistrip::UNDEFINED_RUN_TYPE;
  
  // Retrieve commissioning task from EventSummary
  if ( summary ) { 
    task_ = summary->runType(); 
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Identified CommissioningTask from EventSummary to be \"" 
       << SiStripEnumsAndStrings::runType( task_ )
       << "\"";
    LogTrace(mlDqmSource_) << ss.str();
  } else { 
    task_ = sistrip::UNKNOWN_RUN_TYPE; 
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " NULL pointer to SiStripEventSummary!" 
       << " Check SiStripEventSummary is found/present in Event";
    edm::LogWarning(mlDqmSource_) << ss.str();
  } 
  
  // Override task with ParameterSet configurable (if defined)
  sistrip::RunType configurable = SiStripEnumsAndStrings::runType( taskConfigurable_ );
  if ( configurable != sistrip::UNDEFINED_RUN_TYPE &&
       configurable != sistrip::UNKNOWN_RUN_TYPE ) { 
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Overriding CommissioningTask from EventSummary (\"" 
       << SiStripEnumsAndStrings::runType( task_ )
       << "\") with value retrieved from .cfg file (\""
       << SiStripEnumsAndStrings::runType( configurable )
       << "\")!";
    LogTrace(mlDqmSource_) << ss.str();
    task_ = configurable; 
  }
  
  // Create ME (std::string) that identifies commissioning task
  dqm()->setCurrentFolder( base_ + sistrip::root_ );
  std::string task_str = SiStripEnumsAndStrings::runType( task_ );
  dqm()->bookString( std::string(sistrip::taskId_) + sistrip::sep_ + task_str, task_str ); 
  
  // Check commissioning task is known / defined
  if ( task_ == sistrip::UNKNOWN_RUN_TYPE ||
       task_ == sistrip::UNDEFINED_RUN_TYPE ) {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Unexpected CommissioningTask found (" 
       << static_cast<uint16_t>(task_) << ") \""
       << SiStripEnumsAndStrings::runType( task_ ) << "\""
       << " Unexpected value found in SiStripEventSummary and/or cfg file"
       << " If SiStripEventSummary is not present in Event,"
       << " check 'CommissioningTask' configurable in cfg file";
    edm::LogWarning(mlDqmSource_) << ss.str();
    return; 
  } else {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]"
       << " Identified CommissioningTask to be \"" 
       << SiStripEnumsAndStrings::runType( task_ )
       << "\"";
    LogTrace(mlDqmSource_) << ss.str();
  }
  
  // Check if commissioning task is FED cabling 
  if ( task_ == sistrip::FED_CABLING ) { cablingTask_ = true; }
  else { cablingTask_ = false; }

  std::stringstream ss;
  ss << "[SiStripCommissioningSource::" << __func__ << "]"
     << " CommissioningTask: "
     << SiStripEnumsAndStrings::runType( summary->runType() );
  LogTrace(mlDqmSource_) << ss.str();

  edm::LogVerbatim(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Creating CommissioningTask objects and booking histograms...";
  if ( cablingTask_ ) { createCablingTasks(); }
  else { createTasks( task_ , setup ); }
  edm::LogVerbatim(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Finished booking histograms!";
  tasksExist_ = true;
  
}
void SiStripCommissioningSource::createTasks ( sistrip::RunType  run_type,
const edm::EventSetup setup 
) [private]

Definition at line 927 of file SiStripCommissioningSource.cc.

References sistrip::APV_LATENCY, sistrip::APV_TIMING, base_, sistrip::CALIBRATION, sistrip::CALIBRATION_DECO, sistrip::CALIBRATION_SCAN, sistrip::CALIBRATION_SCAN_DECO, SiStripFedCabling::connections(), sistrip::DAQ_SCOPE_MODE, sistrip::detectorView_, dir, sistrip::dir_, dqm(), sistrip::FAST_CABLING, SiStripFedKey::feChan(), sistrip::FED_TIMING, fedCabling_, SiStripFedCabling::feds(), SiStripFedKey::feUnit(), filename_, sistrip::FINE_DELAY, sistrip::invalid_, sistrip::mlDqmSource_, sistrip::NOISE, sistrip::OPTO_SCAN, parameters_, sistrip::PEDESTALS, sistrip::PEDS_FULL_NOISE, sistrip::PEDS_ONLY, sistrip::root_, run_, DQMStore::setCurrentFolder(), HcalObjRepresent::setup(), task_, tasks_, sistrip::UNDEFINED_RUN_TYPE, updateFreq_, view_, and sistrip::VPSP_SCAN.

Referenced by createTask().

                                                                                                  {

  uint16_t booked = 0;

  // latency has partition-level histos, so treat separately
  if (task_ == sistrip::APV_LATENCY) {

    for (uint16_t partition = 0; partition < 4; ++partition) {
      // make a task for every partition; tracker-wide histo is shared
      tasks_[0][partition] = new LatencyTask( dqm(), FedChannelConnection(partition+1,sistrip::invalid_,sistrip::invalid_,sistrip::invalid_,sistrip::invalid_) );
      tasks_[0][partition]->eventSetup( &setup );
      tasks_[0][partition]->bookHistograms();
      tasks_[0][partition]->updateFreq( updateFreq_ );
      booked++;
    }

  // fine-delay has 1 histo for the whole tracker, so treat separately
  } else if (task_ == sistrip::FINE_DELAY) {

    tasks_[0][0] = new FineDelayTask( dqm(), FedChannelConnection() );
    tasks_[0][0]->eventSetup( &setup );
    tasks_[0][0]->bookHistograms();
    tasks_[0][0]->updateFreq( updateFreq_ );
    booked++;

  } else { // now do any other task

    // Iterate through FED ids and channels 
    std::vector<uint16_t>::const_iterator ifed = fedCabling_->feds().begin();
    for ( ; ifed != fedCabling_->feds().end(); ifed++ ) {

      // Iterate through connected FED channels
      // S.L.: currently copy the vector, because it changes later when
      // reading in peds for calibration run -> not understood memory corruption!
      std::vector<FedChannelConnection> conns = fedCabling_->connections(*ifed);
      std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
      for ( ; iconn != conns.end(); iconn++ ) {

        // Create FEC key
        SiStripFecKey fec_key( iconn->fecCrate(), 
                               iconn->fecSlot(), 
                               iconn->fecRing(), 
                               iconn->ccuAddr(), 
                               iconn->ccuChan() );
        // Create FED key and check if non-zero
        SiStripFedKey fed_key( iconn->fedId(), 
                               SiStripFedKey::feUnit(iconn->fedCh()),
                               SiStripFedKey::feChan(iconn->fedCh()) );
        if ( !iconn->isConnected() ) { continue; }

        // define the view in which to work and paths for histograms
        //   currently FecView (control view) and FedView (readout view)
        //   DetView (detector view) implementation has started
        // Set working directory prior to booking histograms 
        std::stringstream dir;
        dir << base_;
        if (view_ == "Default") { // default
          if ( run_type == sistrip::FAST_CABLING ) {
            dir << fed_key.path(); // cabling in fed view
          } else {
            dir << fec_key.path(); // all other runs in control view
          }
        } else if (view_ == "FecView") {
          dir << fec_key.path();
        } else if (view_ == "FedView") {
          dir << fed_key.path();
        } else if (view_ == "DetView") {
          // currently just by detid from the connection, which is empty...
          dir << sistrip::root_ << sistrip::dir_
              << sistrip::detectorView_ << sistrip::dir_
              << iconn->detId();
        } else {
          edm::LogWarning(mlDqmSource_)
            << "[SiStripCommissioningSource::" << __func__ << "]"
            << " Invalid view " << view_ << std::endl
            << " Histograms will end up all in the top directory.";
        } // end if view_ == ...
        dqm()->setCurrentFolder( dir.str() );

        // Create commissioning task objects
        if ( !tasks_[iconn->fedId()][iconn->fedCh()] ) { 
          if ( task_ == sistrip::FAST_CABLING ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new FastFedCablingTask( dqm(), *iconn ); 
          } else if ( task_ == sistrip::APV_TIMING ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new ApvTimingTask( dqm(), *iconn ); 
          } else if ( task_ == sistrip::FED_TIMING ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new FedTimingTask( dqm(), *iconn );
          } else if ( task_ == sistrip::OPTO_SCAN ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new OptoScanTask( dqm(), *iconn );
          } else if ( task_ == sistrip::VPSP_SCAN ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new VpspScanTask( dqm(), *iconn );
          } else if ( task_ == sistrip::PEDESTALS ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new PedestalsTask( dqm(), *iconn );
          } else if ( task_ == sistrip::PEDS_ONLY ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new PedsOnlyTask( dqm(), *iconn );
          } else if ( task_ == sistrip::NOISE ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new NoiseTask( dqm(), *iconn );
          } else if ( task_ == sistrip::PEDS_FULL_NOISE ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new PedsFullNoiseTask( dqm(), *iconn, parameters_ );
          } else if ( task_ == sistrip::DAQ_SCOPE_MODE ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new DaqScopeModeTask( dqm(), *iconn );
          } else if ( task_ == sistrip::CALIBRATION_SCAN || 
                      task_ == sistrip::CALIBRATION_SCAN_DECO ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new CalibrationScanTask( dqm(), 
                                                                              *iconn, 
                                                                              task_, 
                                                                              filename_.c_str(), 
                                                                              run_, 
                                                                              setup );
          } else if ( task_ == sistrip::CALIBRATION || 
                      task_ == sistrip::CALIBRATION_DECO ) { 
            tasks_[iconn->fedId()][iconn->fedCh()] = new CalibrationTask( dqm(), 
                                                                          *iconn, 
                                                                          task_, 
                                                                          filename_.c_str(), 
                                                                          run_, 
                                                                          setup );
          } else if ( task_ == sistrip::UNDEFINED_RUN_TYPE ) { 
            edm::LogWarning(mlDqmSource_)  
              << "[SiStripCommissioningSource::" << __func__ << "]"
              << " Undefined CommissioningTask" 
              << " Unable to create CommissioningTask object!";
          } else { 
            edm::LogWarning(mlDqmSource_)
              << "[SiStripCommissioningSource::" << __func__ << "]"
              << " Unknown CommissioningTask" 
              << " Unable to create CommissioningTask object!";
          }

          // Check if fed_key is found and, if so, book histos and set update freq
          if ( tasks_[iconn->fedId()][iconn->fedCh()] ) {
            tasks_[iconn->fedId()][iconn->fedCh()]->eventSetup( &setup );
            tasks_[iconn->fedId()][iconn->fedCh()]->bookHistograms(); 
            tasks_[iconn->fedId()][iconn->fedCh()]->updateFreq( updateFreq_ ); 
            booked++;
            //std::stringstream ss;
            //ss << "[SiStripCommissioningSource::" << __func__ << "]"
            //<< " Booking histograms for '" << tasks_[iconn->fedId()][iconn->fedCh()]->myName()
            //<< "' object with key 0x" << std::hex << std::setfill('0') << std::setw(8) << fed_key.key() << std::dec
            //<< " in directory " << dir.str();
            //LogTrace(mlDqmSource_) << ss.str();
          } else {
            std::stringstream ss;
            ss << "[SiStripCommissioningSource::" << __func__ << "]"
               << " NULL pointer to CommissioningTask for key 0x"
               << std::hex << std::setfill('0') << std::setw(8) << fed_key.key() << std::dec
               << " in directory " << dir.str()
               << " Unable to book histograms!";
            edm::LogWarning(mlDqmSource_) << ss.str();
          }

        } else {
          std::stringstream ss;
          ss << "[SiStripCommissioningSource::" << __func__ << "]"
             << " CommissioningTask object already exists for key 0x"
             << std::hex << std::setfill('0') << std::setw(8) << fed_key.key() << std::dec
             << " in directory " << dir.str()
             << " Unable to create CommissioningTask object!";
          edm::LogWarning(mlDqmSource_) << ss.str();
        }

      } // loop over fed channels
    } // loop over feds
  } // end other tasks

  edm::LogVerbatim(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Created " << booked 
    << " CommissioningTask objects and booked histograms";
}
void SiStripCommissioningSource::directory ( std::stringstream &  dir,
uint32_t  run_number = 0 
) [private]

Definition at line 1136 of file SiStripCommissioningSource.cc.

References evf::utils::pid, pos, and tmp.

Referenced by endJob().

                                                                  {

  // 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
  if ( run_number ) {
    dir << std::setw(8) 
        << std::setfill('0') 
        << run_number
        << "_";
  }
  dir << ip.str()
      << "_"
      << std::setw(5) 
      << std::setfill('0') 
      << pid;
  
}
DQMStore *const SiStripCommissioningSource::dqm ( std::string  method = "") const [private]

Definition at line 92 of file SiStripCommissioningSource.cc.

References dqm_, and sistrip::mlDqmSource_.

Referenced by beginRun(), createCablingTasks(), createRunNumber(), createTask(), createTasks(), endJob(), and remove().

                                                                        {
  if ( !dqm_ ) { 
    std::stringstream ss;
    if ( method != "" ) { ss << "[SiStripCommissioningSource::" << method << "]" << std::endl; }
    else { ss << "[SiStripCommissioningSource]" << std::endl; }
    ss << " NULL pointer to DQMStore";
    edm::LogWarning(mlDqmSource_) << ss.str();
    return 0;
  } else { return dqm_; }
}
void SiStripCommissioningSource::endJob ( void  ) [virtual]

Reimplemented from edm::EDAnalyzer.

Definition at line 159 of file SiStripCommissioningSource.cc.

References sistrip::APV_LATENCY, cablingTasks_, SiStripFedCabling::connections(), dir, directory(), dqm(), fecCabling_, fedCabling_, SiStripFedCabling::feds(), filename_, sistrip::FINE_DELAY, LogTrace, sistrip::mlDqmSource_, mergeVDriftHistosByStation::name, NULL, run_, DQMStore::save(), task_, and tasks_.

                                        {
  LogTrace(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Halting..." << std::endl;
  
  // ---------- Update histograms ----------
  
  // Cabling task
  for ( TaskMap::iterator itask = cablingTasks_.begin(); itask != cablingTasks_.end(); itask++ ) { 
    if ( itask->second ) { itask->second->updateHistograms(); }
  }
  
  if (task_ == sistrip::APV_LATENCY) {
    for (uint16_t partition = 0; partition < 4; ++partition) {
      tasks_[0][partition]->updateHistograms();
    }
  } else if (task_ == sistrip::FINE_DELAY ) {
    tasks_[0][0]->updateHistograms();
  } else {
    // All tasks except cabling 
    uint16_t fed_id = 0;
    uint16_t fed_ch = 0;
    std::vector<uint16_t>::const_iterator ifed = fedCabling_->feds().begin(); 
    for ( ; ifed != fedCabling_->feds().end(); ifed++ ) { 
      const std::vector<FedChannelConnection>& conns = fedCabling_->connections(*ifed);
      std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
      for ( ; iconn != conns.end(); iconn++ ) {
        if ( !iconn->isConnected() ) { continue; }
        fed_id = iconn->fedId();
        fed_ch = iconn->fedCh();
        if ( tasks_[fed_id][fed_ch] ) { 
          tasks_[fed_id][fed_ch]->updateHistograms();
          delete tasks_[fed_id][fed_ch];
        }
      }
    }
  }
  
  // ---------- Save histos to root file ----------

  // 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()); 
  }

  // Add directory path 
  std::stringstream ss; 
  if ( !dir.empty() ) { ss << dir << "/"; }
  else { ss << "/tmp/"; }

  // Add filename with run number, host ip, pid and .root extension
  ss << name << "_";
  directory(ss,run_);
  ss << ".root";

  // Save file with appropriate filename (if run number is known)
  if ( !filename_.empty() ) { 
    if ( run_ != 0 ) { dqm()->save( ss.str() ); }
    else {
      edm::LogWarning(mlDqmSource_)
        << "[SiStripCommissioningSource::" << __func__ << "]"
        << " NULL value for RunNumber! No root file saved!";
    }
  } else {
    edm::LogWarning(mlDqmSource_)
      << "[SiStripCommissioningSource::" << __func__ << "]"
      << " NULL value for filename! No root file saved!";
  }

  LogTrace(mlDqmSource_)
    << "[SiStripCommissioningSource::" << __func__ << "]"
    << " Saved all histograms to file \""
    << ss.str() << "\"";
  
  // ---------- Delete histograms ----------
  
  // Remove all MonitorElements in "SiStrip" dir and below
  // remove();
  
  // Delete histogram objects
  // clearCablingTasks();
  // clearTasks();
  
  // ---------- Delete cabling ----------

  if ( fedCabling_ ) { fedCabling_ = 0; }
  if ( fecCabling_ ) { delete fecCabling_; fecCabling_ = 0; }
  
}
void SiStripCommissioningSource::fillCablingHistos ( const SiStripEventSummary * const  summary,
const edm::DetSetVector< SiStripRawDigi > &  raw 
) [private]

Definition at line 360 of file SiStripCommissioningSource.cc.

References Averages::add(), PDRates::average, cablingTasks_, Averages::calc(), sistrip::CCU_CHAN, SiStripFecKey::ccuAddr(), SiStripFecKey::ccuChan(), SiStripKey::channel(), SiStripEventSummary::dcuId(), SiStripEventSummary::deviceId(), edm::DetSetVector< T >::end(), fecCabling_, SiStripFecKey::fecCrate(), SiStripFecKey::fecRing(), SiStripFecKey::fecSlot(), fedCabling_, SiStripFedCabling::feds(), edm::DetSetVector< T >::find(), sistrip::invalid_, findQualityFiles::jj, SiStripKey::key(), SiStripModule::key(), LogTrace, Averages::Params::median_, sistrip::mlDqmSource_, sistrip::mlTest_, SiStripFecCabling::module(), getHLTPrescaleColumns::path, Averages::Params::rms_, and tmp.

Referenced by analyze().

                                                                                                 {
  
  // Create FEC key using DCU id and LLD channel from SiStripEventSummary
  const SiStripModule& module = fecCabling_->module( summary->dcuId() );
  uint16_t lld_channel = ( summary->deviceId() & 0x3 ) + 1;
  SiStripFecKey key_object( module.key().fecCrate(),
                            module.key().fecSlot(),
                            module.key().fecRing(),
                            module.key().ccuAddr(),
                            module.key().ccuChan(),
                            lld_channel );
  uint32_t fec_key = key_object.key();
  std::stringstream sss;
  sss << "[SiStripCommissioningSource::" << __func__ << "]" 
      << " Found DcuId 0x"
      << std::hex << std::setw(8) << std::setfill('0') << summary->dcuId() << std::dec 
      << " with Crate/FEC/Ring/CCU/Module/LLD: "
      << module.key().fecCrate() << "/"
      << module.key().fecSlot() << "/"
      << module.key().fecRing() << "/"
      << module.key().ccuAddr() << "/"
      << module.key().ccuChan() << "/"
      << lld_channel;
  edm::LogWarning(mlDqmSource_) << sss.str();
  
  //LogTrace(mlTest_) << "TEST : " << key_object;
  
  // Check on whether DCU id is found
  if ( key_object.isInvalid( sistrip::CCU_CHAN ) ) {
    std::stringstream ss;
    ss << "[SiStripCommissioningSource::" << __func__ << "]" 
       << " DcuId 0x"
       << std::hex << std::setw(8) << std::setfill('0') << summary->dcuId() << std::dec 
       << " in 'DAQ register' field not found in cabling map!"
       << " (NULL values returned for FEC path)";
    edm::LogWarning(mlDqmSource_) << ss.str();
    return;
  }
  
  // Iterate through FED ids
  std::vector<uint16_t>::const_iterator ifed = fedCabling_->feds().begin(); 
  for ( ; ifed != fedCabling_->feds().end(); ifed++ ) {

    // Check if FedId is non-zero
    if ( *ifed == sistrip::invalid_ ) { continue; }
    
    // Container to hold median signal level for FED cabling task
    std::map<uint16_t,float> medians; medians.clear(); 
    std::map<uint16_t,float> medians1; medians1.clear(); 
    
    // Iterate through FED channels
    for ( uint16_t ichan = 0; ichan < 96; ichan++ ) {
      
      //       // Retrieve digis for given FED key
      //       uint32_t fed_key = SiStripFedKey( *ifed, 
      //                                        SiStripFedKey::feUnit(ichan),
      //                                        SiStripFedKey::feChan(ichan) ).key();
      
      // Retrieve digis for given FED key
      uint32_t fed_key = ( ( *ifed & sistrip::invalid_ ) << 16 ) | ( ichan & sistrip::invalid_ );
      
      std::vector< edm::DetSet<SiStripRawDigi> >::const_iterator digis = raw.find( fed_key );
      if ( digis != raw.end() ) { 
        if ( digis->data.empty() ) { continue; }
        
        //      if ( digis->data[0].adc() > 500 ) {
        //        std::stringstream ss;
        //        ss << " HIGH SIGNAL " << digis->data[0].adc() << " FOR"
        //           << " FedKey: 0x" << std::hex << std::setw(8) << std::setfill('0') << fed_key << std::dec
        //           << " FedId/Ch: " << *ifed << "/" << ichan;
        //        LogTrace(mlDqmSource_) << ss.str();
        //      }
        
        Averages ave;
        for ( uint16_t idigi = 0; idigi < digis->data.size(); idigi++ ) { 
          ave.add( static_cast<uint32_t>(digis->data[idigi].adc()) ); 
        }
        Averages::Params params;
        ave.calc(params);
        medians[ichan] = params.median_; // Store median signal level
        medians1[ichan] = digis->data[0].adc(); 

        //      if ( !digis->data.empty() ) { medians[ichan] = digis->data[0].adc(); }
        //      else { 
        //        edm::LogWarning(mlTest_) << "TEST : NO DIGIS!";
        //      }
        
        //              std::stringstream ss;
        //              ss << "Channel Averages:" << std::endl
        //                 << "  nDigis: " << digis->data.size() << std::endl
        //                 << "  num/mean/MEDIAN/rms/max/min: "
        //                 << params.num_ << "/"
        //                 << params.mean_ << "/"
        //                 << params.median_ << "/"
        //                 << params.rms_ << "/"
        //                 << params.max_ << "/"
        //                 << params.min_ << std::endl;
        //              LogTrace(mlDqmSource_) << ss.str();

      }
      
    } // fed channel loop

    // Calculate mean and spread on all (median) signal levels
    Averages average;
    std::map<uint16_t,float>::const_iterator ii = medians.begin();
    for ( ; ii != medians.end(); ii++ ) { average.add( ii->second ); }
    Averages::Params tmp;
    average.calc(tmp);
      
    //     std::stringstream ss;
    //     ss << "FED Averages:" << std::endl
    //        << "  nChans: " << medians.size() << std::endl
    //        << "  num/mean/median/rms/max/min: "
    //        << tmp.num_ << "/"
    //        << tmp.mean_ << "/"
    //        << tmp.median_ << "/"
    //        << tmp.rms_ << "/"
    //        << tmp.max_ << "/"
    //        << tmp.min_ << std::endl;
    //     LogTrace(mlDqmSource_) << ss.str();
      
    // Calculate mean and spread on "filtered" data
    Averages truncated;
    std::map<uint16_t,float>::const_iterator jj = medians.begin();
    for ( ; jj != medians.end(); jj++ ) { 
      if ( jj->second < tmp.median_+tmp.rms_ ) { 
        truncated.add( jj->second ); 
      }
    }
    Averages::Params params;
    truncated.calc(params);
      
    //     std::stringstream ss1;
    //     ss1 << "Truncated Averages:" << std::endl
    //  << "  nChans: " << medians.size() << std::endl
    //  << "  num/mean/median/rms/max/min: "
    //  << params.num_ << "/"
    //  << params.mean_ << "/"
    //  << params.median_ << "/"
    //  << params.rms_ << "/"
    //  << params.max_ << "/"
    //  << params.min_ << std::endl;
    //     LogTrace(mlDqmSource_) << ss1.str();

    // Identify channels with signal
    std::stringstream ss2;
    std::stringstream ss3;
    //     ss2 << "Number of possible connections: " << medians.size()
    //  << " channel/signal: ";
    std::map<uint16_t,float> channels;
    std::map<uint16_t,float>::const_iterator ichan = medians.begin();
    for ( ; ichan != medians.end(); ichan++ ) { 
      //             cout << " mean: " << params.mean_
      //           << " rms: " << params.rms_
      //           << " thresh: " << params.mean_ + 5.*params.rms_
      //           << " value: " << ichan->second
      //           << " strip: " << ichan->first << std::endl;
      //if ( ichan->second > params.mean_ + 5.*params.rms_ ) { 
      if ( ichan->second > 200. ) { 
        LogTrace(mlTest_) << "TEST FOUND SIGNAL HIGH: " << *ifed << " " << ichan->first << " " << ichan->second;
        channels[ichan->first] = ichan->second;
      }
      ss2 //<< ichan->first << "/" 
        << ichan->second << " ";
      ss3 //<< ichan->first << "/" 
        << medians1[ichan->first] << " ";
    }

    ss2 << std::endl;
    ss3 << std::endl;
    LogTrace(mlTest_) << "DUMP for FED  " << *ifed << ": " << ss2.str();
    LogTrace(mlTest_) << "FIRST ADC VAL " << *ifed << ": " << ss3.str();

    //     LogTrace(mlDqmSource_)
    //       << "[FedCablingTask::" << __func__ << "]"
    //       << " Found candidate connection between device: 0x"
    //       << std::setfill('0') << std::setw(8) << std::hex << summary.deviceId() << std::dec
    //       << " with Crate/FEC/Ring/CCU/Module/LLDchannel: " 
    //       << connection().fecCrate() << "/"
    //       << connection().fecSlot() << "/"
    //       << connection().fecRing() << "/"
    //       << connection().ccuAddr() << "/"
    //       << connection().ccuChan() << "/"
    //       << connection().lldChannel()
    //       << " and FedId/Ch: "
    //       << fed_id << "/" << ichan->first
    //       << " with signal " << ichan->second 
    //       << " [adc] over background " << "XXX +/- YYY [adc]" 
    //       << " (S/N = " << "ZZZ" << ")";

    
    // Fill cabling histograms
    if ( cablingTasks_.find(fec_key) != cablingTasks_.end() ) {
      if ( !channels.empty() ) { 
        cablingTasks_[fec_key]->fillHistograms( *summary, *ifed, channels );
        SiStripFecKey path( fec_key );
        std::stringstream ss;
        ss << "[SiStripCommissioningSource::" << __func__ << "]"
           << " Filled histogram for '" << cablingTasks_[fec_key]->myName()
           << "' object with FecKey: 0x" 
           << std::hex << std::setfill('0') << std::setw(8) << fec_key << std::dec
           << " and Crate/FEC/ring/CCU/module/LLDchan: " 
           << path.fecCrate() << "/"
           << path.fecSlot() << "/"
           << path.fecRing() << "/"
           << path.ccuAddr() << "/"
           << path.ccuChan() << "/"
           << path.channel();
        LogTrace(mlDqmSource_) << ss.str();
      }
    } else {
      SiStripFecKey path( fec_key );
      std::stringstream ss;
      ss << "[SiStripCommissioningSource::" << __func__ << "]"
         << " Unable to find CommissioningTask object with FecKey: 0x" 
         << std::hex << std::setfill('0') << std::setw(8) << fec_key << std::dec
         << " and Crate/FEC/ring/CCU/module/LLDchan: " 
         << path.fecCrate() << "/"
         << path.fecSlot() << "/"
         << path.fecRing() << "/"
         << path.ccuAddr() << "/"
         << path.ccuChan() << "/"
         << path.channel();
      edm::LogWarning(mlDqmSource_) << ss.str();
    }
  
  } // fed id loop
  
}
void SiStripCommissioningSource::fillHistos ( const SiStripEventSummary * const  summary,
const edm::DetSetVector< SiStripRawDigi > &  raw 
) [private]

Definition at line 594 of file SiStripCommissioningSource.cc.

References sistrip::APV_LATENCY, SiStripFedCabling::connections(), edm::DetSetVector< T >::end(), fedCabling_, SiStripFedCabling::feds(), edm::DetSetVector< T >::find(), sistrip::FINE_DELAY, sistrip::invalid_, sistrip::mlDqmSource_, task_, tasks_, and sistrip::valid_.

Referenced by analyze().

                                                                                          {

    // Iterate through FED ids and channels
    std::vector<uint16_t>::const_iterator ifed = fedCabling_->feds().begin();
    for ( ; ifed != fedCabling_->feds().end(); ifed++ ) {

      // Iterate through connected FED channels
      const std::vector<FedChannelConnection>& conns = fedCabling_->connections(*ifed);
      std::vector<FedChannelConnection>::const_iterator iconn = conns.begin();
      for ( ; iconn != conns.end(); iconn++ ) {

        if ( !(iconn->fedId()) || iconn->fedId() > sistrip::valid_ ) { continue; }
     
        //       // Create FED key and check if non-zero
        //       uint32_t fed_key = SiStripFedKey( iconn->fedId(), 
        //                                      SiStripFedKey::feUnit(iconn->fedCh()),
        //                                      SiStripFedKey::feChan(iconn->fedCh()) ).key();

        // Create FED key and check if non-zero
        // note: the key is not computed using the same formula as in commissioning histograms.
        // beware that changes here must match changes in raw2digi and in SiStripFineDelayHit
        uint32_t fed_key = ( ( iconn->fedId() & sistrip::invalid_ ) << 16 ) | ( iconn->fedCh() & sistrip::invalid_ );

        // Retrieve digis for given FED key and check if found
        std::vector< edm::DetSet<SiStripRawDigi> >::const_iterator digis = raw.find( fed_key ); 

        if ( digis != raw.end() ) {
          // tasks involving tracking have partition-level histos, so treat separately
          if (task_ == sistrip::APV_LATENCY) {
            if ( tasks_[0][iconn->fecCrate()-1] ) {
              tasks_[0][iconn->fecCrate()-1]->fillHistograms( *summary, *digis );
            } else {
              std::stringstream ss;
              ss << "[SiStripCommissioningSource::" << __func__ << "]"
                 << " Unable to find CommissioningTask for FEC crate " 
                 << iconn->fecCrate() << ". Unable to fill histograms!";
              edm::LogWarning(mlDqmSource_) << ss.str();
            }
          } else if (task_ == sistrip::FINE_DELAY) {
            if ( tasks_[0][0] ) {
              tasks_[0][0]->fillHistograms( *summary, *digis );
            } else {
              std::stringstream ss;
              ss << "[SiStripCommissioningSource::" << __func__ << "]"
                 << " Unable to find global CommissioningTask for FineDelay. Unable to fill histograms!";
              edm::LogWarning(mlDqmSource_) << ss.str();
            }
          // now do any other task
          } else {
            if ( tasks_[iconn->fedId()][iconn->fedCh()] ) {
              tasks_[iconn->fedId()][iconn->fedCh()]->fillHistograms( *summary, *digis );
            } else {
              std::stringstream ss;
              ss << "[SiStripCommissioningSource::" << __func__ << "]"
                 << " Unable to find CommissioningTask object with FED key " 
                 << std::hex << std::setfill('0') << std::setw(8) << fed_key << std::dec
                 << " and FED id/ch " 
                 << iconn->fedId() << "/"
                 << iconn->fedCh()
                 << " Unable to fill histograms!";
              edm::LogWarning(mlDqmSource_) << ss.str();
            }
          }
        } else {
          // issue a warning only for standard runs, as latency and fine delay only deliver 
          // pseudo zero-suppressed data
          if ( task_ != sistrip::APV_LATENCY &&
               task_ != sistrip::FINE_DELAY ) {
            std::stringstream ss;
            ss << "[SiStripCommissioningSource::" << __func__ << "]"
               << " Unable to find any DetSet containing digis for FED key " 
               << std::hex << std::setfill('0') << std::setw(8) << fed_key << std::dec
               << " and FED id/ch " 
               << iconn->fedId() << "/"
               << iconn->fedCh();
            edm::LogWarning(mlDqmSource_) << ss.str();
          }
        }
      } // fed channel loop
    } // fed id loop

}
void SiStripCommissioningSource::remove ( ) [private]

Member Data Documentation

std::string SiStripCommissioningSource::base_ [private]

Flag to indicate whether task is FED cabling or not.

Definition at line 131 of file SiStripCommissioningSource.h.

Referenced by analyze(), beginRun(), and createTask().

Map of cabling task objects (indexed using FEC key).

Definition at line 125 of file SiStripCommissioningSource.h.

Referenced by clearCablingTasks(), createCablingTasks(), endJob(), and fillCablingHistos().

Interface to Data Quality Monitoring framework.

Definition at line 90 of file SiStripCommissioningSource.h.

Referenced by beginRun(), and dqm().

Filename of output root file containing source histos.

Definition at line 105 of file SiStripCommissioningSource.h.

Referenced by createTasks(), and endJob().

Name of digi input module.

Definition at line 101 of file SiStripCommissioningSource.h.

Referenced by analyze().

Definition at line 102 of file SiStripCommissioningSource.h.

Referenced by analyze().

parameters to pass to the tasks

Definition at line 143 of file SiStripCommissioningSource.h.

Referenced by createTasks().

Run number used for naming of root file.

Definition at line 108 of file SiStripCommissioningSource.h.

Referenced by analyze(), createRunNumber(), createTasks(), and endJob().

Identifies commissioning task.

Definition at line 119 of file SiStripCommissioningSource.h.

Referenced by analyze(), beginRun(), createCablingTasks(), createTask(), createTasks(), endJob(), and fillHistos().

Identifies commissioning task read from cfg file.

Definition at line 116 of file SiStripCommissioningSource.h.

Referenced by createTask().

Vector of vector of task objects (indexed using FED id.ch.

Definition at line 122 of file SiStripCommissioningSource.h.

Referenced by clearTasks(), createTasks(), endJob(), fillHistos(), and SiStripCommissioningSource().

Flag to indicate whether histo objects exist or not.

Definition at line 128 of file SiStripCommissioningSource.h.

Referenced by analyze(), beginRun(), and createTask().

Record of time used to calculate event rate.

Definition at line 111 of file SiStripCommissioningSource.h.

Referenced by analyze().

Update frequency for histograms (ignored for cabling).

Definition at line 134 of file SiStripCommissioningSource.h.

Referenced by analyze(), and createTasks().

std::string SiStripCommissioningSource::view_ [private]

flag for choosing the organizational 'view' the DQM histogram tree

Definition at line 140 of file SiStripCommissioningSource.h.

Referenced by createTasks().