DetStatus.cc

Go to the documentation of this file.

#include "DPGAnalysis/Skims/interface/DetStatus.h"
#include "FWCore/Framework/interface/Event.h"
#include "DataFormats/Common/interface/Handle.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include <iostream>

using namespace std;

//
// -- Constructor
//
DetStatus::DetStatus(const edm::ParameterSet& pset) {
  verbose_ = pset.getUntrackedParameter<bool>("DebugOn", false);
  AndOr_ = pset.getParameter<bool>("AndOr");
  applyfilter_ = pset.getParameter<bool>("ApplyFilter");
  DetNames_ = pset.getParameter<std::vector<std::string> >("DetectorType");

  // build a map
  DetMap_ = 0;
  for (unsigned int detreq = 0; detreq < DetNames_.size(); detreq++) {
    for (unsigned int detlist = 0; detlist < DcsStatus::nPartitions; detlist++) {
      if (DetNames_[detreq] == DcsStatus::partitionName[detlist]) {
        DetMap_ |= (1 << DcsStatus::partitionList[detlist]);
        if (verbose_)
          std::cout << "DCSStatus filter: asked partition " << DcsStatus::partitionName[detlist] << " bit "
                    << DcsStatus::partitionList[detlist] << std::endl;
      }
    }
  }
  edm::InputTag scalersTag("scalersRawToDigi");
  scalersToken_ = consumes<DcsStatusCollection>(scalersTag);
}

//
// -- Destructor
//
DetStatus::~DetStatus() {}

bool DetStatus::filter(edm::Event& evt, edm::EventSetup const& es) {
  bool accepted = false;

  edm::Handle<DcsStatusCollection> dcsStatus;
  evt.getByToken(scalersToken_, dcsStatus);

  if (dcsStatus.isValid()) {
    if (dcsStatus->empty()) {
      // Only complain in case it's real data. Accept in any case.
      if (evt.eventAuxiliary().isRealData())
        edm::LogError("DetStatus") << "Error! dcsStatus has size 0, accept in any case";
      accepted = true;
    } else {
      unsigned int curr_dcs = (*dcsStatus)[0].ready();
      if (verbose_)
        std::cout << "curr_dcs = " << curr_dcs << std::endl;
      if (AndOr_)
        accepted = ((DetMap_ & curr_dcs) == DetMap_);
      else
        accepted = ((DetMap_ & curr_dcs) != 0);

      if (verbose_)

      {
        std::cout << "DCSStatus filter: requested map: " << DetMap_ << " dcs in event: " << curr_dcs
                  << " filter: " << accepted << std::endl;
        std::cout << "Partitions ON: ";
        for (unsigned int detlist = 0; detlist < DcsStatus::nPartitions; detlist++) {
          if ((*dcsStatus)[0].ready(DcsStatus::partitionList[detlist])) {
            std::cout << " " << DcsStatus::partitionName[detlist];
          }
        }
        std::cout << std::endl;
      }
    }
  } else {
    edm::LogError("DetStatus") << "Error! can't get the product: scalersRawToDigi, accept in any case";
    accepted = true;
  }

  if (!applyfilter_)
    accepted = true;
  return accepted;
}

#include "FWCore/Framework/interface/MakerMacros.h"
DEFINE_FWK_MODULE(DetStatus);