L1TFED.cc

Go to the documentation of this file.

/*
 * \file L1TFED.cc
 *
 * \author J. Berryhill
 *
 */

#include "DQM/L1TMonitor/interface/L1TFED.h"

using namespace std;
using namespace edm;

L1TFED::L1TFED(const ParameterSet& ps)
{
  // verbosity switch
  verbose_        = ps.getUntrackedParameter<bool>("verbose", false);
  l1feds_         = ps.getParameter<std::vector<int> >("L1FEDS");
  directory_      = ps.getUntrackedParameter<std::string>("FEDDirName","L1T/FEDIntegrity");
  stableROConfig_ = ps.getUntrackedParameter<bool>("stableROConfig", true);
  rawl_           = consumes<FEDRawDataCollection>(ps.getParameter< InputTag >("rawTag"));

  if(verbose_) cout << "L1TFED: constructor...." << endl;
  nev_ = 0;
}

L1TFED::~L1TFED()
{
}

void L1TFED::bookHistograms(DQMStore::IBooker & ibooker, edm::Run const & iRun, edm::EventSetup const & iSetup) 
{
  ibooker.setCurrentFolder(directory_);
  
  fedentries = ibooker.book1D("FEDEntries", "Fed ID occupancy", l1feds_.size(), 0.,l1feds_.size() );      
  fedfatal = ibooker.book1D("FEDFatal", "Fed ID non present ", l1feds_.size(), 0., l1feds_.size());   
  fednonfatal = ibooker.book1D("FEDNonFatal", "Fed corrupted data ", l1feds_.size(), 0.,l1feds_.size() );
  hfedprof = ibooker.bookProfile("fedprofile","FED Size by ID", l1feds_.size(), 0., l1feds_.size(),0,0.,5000.);
  for(unsigned int i=0;i<l1feds_.size();i++){
    ostringstream sfed;
    sfed << l1feds_[i];
    fedentries->setBinLabel(i+1,"FED "+ sfed.str());
    fedfatal->setBinLabel(i+1,"FED "+ sfed.str());
    fednonfatal->setBinLabel(i+1,"FED "+ sfed.str());
  }
  
  hfedsize = ibooker.book1D("fedsize","FED Size Distribution",100,0.,10000.);
}


void L1TFED::analyze(const Event& e, const EventSetup& c)
{
  nev_++; 
  if(verbose_) cout << "L1T FED Integrity: analyze...." << endl;

  edm::Handle<FEDRawDataCollection> rawdata;
  bool t = e.getByToken(rawl_,rawdata);
  
  if ( ! t ) {
    if(verbose_) cout << "can't find FEDRawDataCollection "<< endl;
  }
  
  else {

     if(verbose_) cout << "fedlist size = " << l1feds_.size() << endl;

     for (unsigned int i = 0; i<l1feds_.size(); i++){
        int fedId = l1feds_[i];
        if(verbose_) cout << "fedId = " << fedId << endl;
       
        const FEDRawData & data = rawdata->FEDData(fedId);
        
    if(size_t size=data.size()){
      
      fedentries->Fill(i);
      hfedsize->Fill(float(size));
      hfedprof->Fill(float(i),float(size));
      if(verbose_) cout << "header check = " << FEDHeader(data.data()).check() << endl;
      if(verbose_) cout << "trailer check = " << FEDTrailer(data.data()).check() << endl;
      
      if(!FEDHeader(data.data()).check()) fedfatal->Fill(i);
      
        } else {
      
      if(verbose_) cout << "empty fed " << i << endl;
      if(stableROConfig_) fedfatal->Fill(i);
      
    }
     }  
  } 
}