ShallowDigisProducer.cc

Go to the documentation of this file.

#include "CalibTracker/SiStripCommon/interface/ShallowDigisProducer.h"

#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/MessageLogger/interface/MessageLogger.h"
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "DataFormats/Common/interface/DetSetVector.h"
#include "DataFormats/Common/interface/DetSetVectorNew.h"
#include "DataFormats/SiStripDigi/interface/SiStripDigi.h"
#include "CondFormats/DataRecord/interface/SiStripNoisesRcd.h"
#include "CondFormats/SiStripObjects/interface/SiStripNoises.h"
#include "boost/foreach.hpp"

ShallowDigisProducer::ShallowDigisProducer(const edm::ParameterSet& conf)
  : inputTags(conf.getParameter<std::vector<edm::InputTag> >("DigiProducersList")) 
{
  produces<std::vector<unsigned> >("id");
  produces<std::vector<unsigned> >("subdet");
  produces<std::vector<unsigned> >("strip");
  produces<std::vector<unsigned> >("adc");
  produces<std::vector<float> >("noise");
}  

void ShallowDigisProducer::
insert(products& p, edm::Event& e) {
  e.put(std::move(p.id),     "id");
  e.put(std::move(p.subdet), "subdet");
  e.put(std::move(p.strip),  "strip");
  e.put(std::move(p.adc),    "adc");
  e.put(std::move(p.noise),  "noise");
}

template<class T>
inline
void ShallowDigisProducer::
recordDigis(const T& digiCollection, products& p) {
  BOOST_FOREACH(const typename T::value_type set, digiCollection) {
    SiStripNoises::Range detNoiseRange = noiseHandle->getRange(set.detId());
    BOOST_FOREACH(const SiStripDigi digi, set) {
      p.id->push_back(set.detId());
      p.subdet->push_back((set.detId()>>25)&0x7);
      p.strip->push_back(digi.strip());
      p.adc->push_back(digi.adc());
      p.noise->push_back(noiseHandle->getNoise( digi.strip(), detNoiseRange));
    }
  }
}

void ShallowDigisProducer::
produce(edm::Event& e, const edm::EventSetup& es) {  
  products p;
  edm::Handle< edm::DetSetVector<SiStripDigi> >     inputOld;  
  edm::Handle< edmNew::DetSetVector<SiStripDigi> >  inputNew;  
  es.get<SiStripNoisesRcd>().get(noiseHandle);
  if( findInput(inputOld, e) ) recordDigis(*inputOld, p); else 
    if( findInput(inputNew, e) ) recordDigis(*inputNew, p); else
      edm::LogWarning("Input Not Found");
  insert(p,e);
}

template<class T>
inline
bool ShallowDigisProducer::
findInput(edm::Handle<T>& handle, const edm::Event& e) {
  BOOST_FOREACH( const edm::InputTag inputTag, inputTags) {
    e.getByLabel(inputTag, handle);
    if( handle.isValid() && !handle->empty() ) {
      LogDebug("Input") << inputTag;
      return true;
    }
  }
  return false;
}