DetectorStateFilter.cc

Go to the documentation of this file.

#include "DQM/TrackerCommon/plugins/DetectorStateFilter.h" 
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/Scalers/interface/DcsStatus.h"

#include "DataFormats/FEDRawData/interface/FEDRawDataCollection.h"
#include "DataFormats/FEDRawData/interface/FEDRawData.h"
#include "DataFormats/FEDRawData/interface/FEDNumbering.h"

#include <iostream>
 
//
// -- Constructor
//
DetectorStateFilter::DetectorStateFilter( const edm::ParameterSet & pset ) {
   verbose_        = pset.getUntrackedParameter<bool>( "DebugOn", false );
   detectorType_   = pset.getUntrackedParameter<std::string>( "DetectorType", "sistrip");
   dcsStatusLabel_ = pset.getUntrackedParameter<edm::InputTag>( "DcsStatusLabel", edm::InputTag("scalersRawToDigi") );

   nEvents_         = 0;
   nSelectedEvents_ = 0;
   detectorOn_  = false;
}
//
// -- Destructor
//
DetectorStateFilter::~DetectorStateFilter() {
}
 
bool DetectorStateFilter::filter( edm::Event & evt, edm::EventSetup const& es) {
  
  nEvents_++;
  // Check Detector state Only for Real Data and return true for MC
  if (evt.isRealData()) {
    edm::Handle<DcsStatusCollection> dcsStatus;
    evt.getByLabel(dcsStatusLabel_, dcsStatus);
    if (dcsStatus.isValid()) {
      if (detectorType_ == "pixel" && dcsStatus->size() > 0 ) {
      if ((*dcsStatus)[0].ready(DcsStatus::BPIX) && 
          (*dcsStatus)[0].ready(DcsStatus::FPIX)) {
    detectorOn_ = true;
    nSelectedEvents_++;
      } else detectorOn_ = false;
      if ( verbose_ ) std::cout << " Total Events " << nEvents_ 
                << " Selected Events " << nSelectedEvents_ 
                << " DCS States : " << " BPix " << (*dcsStatus)[0].ready(DcsStatus::BPIX) 
                << " FPix " << (*dcsStatus)[0].ready(DcsStatus::FPIX)
                << " Detector State " << detectorOn_<<  std::endl;           
      } else if (detectorType_ == "sistrip" && dcsStatus->size() > 0) {  
    if ((*dcsStatus)[0].ready(DcsStatus::TIBTID) &&
        (*dcsStatus)[0].ready(DcsStatus::TOB) &&   
        (*dcsStatus)[0].ready(DcsStatus::TECp) &&  
        (*dcsStatus)[0].ready(DcsStatus::TECm)) {
      detectorOn_ = true;             
      nSelectedEvents_++;
    } else detectorOn_ = false;
    if ( verbose_ ) std::cout << " Total Events " << nEvents_ 
                  << " Selected Events " << nSelectedEvents_ 
                  << " DCS States : " << " TEC- " << (*dcsStatus)[0].ready(DcsStatus::TECm) 
                  << " TEC+ " << (*dcsStatus)[0].ready(DcsStatus::TECp)
                  << " TIB/TID " << (*dcsStatus)[0].ready(DcsStatus::TIBTID) 
                  << " TOB " << (*dcsStatus)[0].ready(DcsStatus::TOB)   
                  << " Detector States " << detectorOn_<<  std::endl;      
      }
    } else {
      edm::LogError("DetectorStatusFilter")<<"ERROR: DcsStatusCollection not found !";
    }
  } else {
    detectorOn_ = true;
    nSelectedEvents_++;
    if ( verbose_ ) std::cout << "Total MC Events " << nEvents_
                              << " Selected Events " << nSelectedEvents_
                              << " Detector States " << detectorOn_<<  std::endl;
  }
  return detectorOn_;
}

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