SiPixelClusterShapeCacheProducer.cc

Go to the documentation of this file.

#include "FWCore/Framework/interface/global/EDProducer.h"
#include "FWCore/Framework/interface/MakerMacros.h"
#include "FWCore/Framework/interface/Event.h"
#include "FWCore/Framework/interface/EventSetup.h"
#include "FWCore/Utilities/interface/EDGetToken.h"
 
#include "FWCore/ParameterSet/interface/ParameterSet.h"
#include "FWCore/Utilities/interface/InputTag.h"
 
#include "FWCore/Framework/interface/ESHandle.h"
#include "DataFormats/Common/interface/Handle.h"
#include "DataFormats/SiPixelCluster/interface/SiPixelClusterShapeCache.h"
 
#include "Geometry/Records/interface/TrackerDigiGeometryRecord.h"
#include "Geometry/TrackerGeometryBuilder/interface/TrackerGeometry.h"
#include "Geometry/CommonDetUnit/interface/PixelGeomDetUnit.h"
 
#include "FWCore/ParameterSet/interface/ConfigurationDescriptions.h"
#include "FWCore/ParameterSet/interface/ParameterSetDescription.h"
 
#include "RecoPixelVertexing/PixelLowPtUtilities/interface/ClusterShape.h"
#include "RecoPixelVertexing/PixelLowPtUtilities/interface/ClusterData.h"
 
#include <cassert>
 
class SiPixelClusterShapeCacheProducer : public edm::global::EDProducer<> {
public:
  explicit SiPixelClusterShapeCacheProducer(const edm::ParameterSet& iConfig);
  ~SiPixelClusterShapeCacheProducer() override;
 
  static void fillDescriptions(edm::ConfigurationDescriptions& descriptions);
 
  void produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup& iSetup) const override;
 
private:
  using InputCollection = edmNew::DetSetVector<SiPixelCluster>;
 
  const edm::EDGetTokenT<InputCollection> token_;
  const edm::ESGetToken<TrackerGeometry, TrackerDigiGeometryRecord> geomToken_;
};
 
SiPixelClusterShapeCacheProducer::SiPixelClusterShapeCacheProducer(const edm::ParameterSet& iConfig)
    : token_(consumes<InputCollection>(iConfig.getParameter<edm::InputTag>("src"))), geomToken_(esConsumes()) {
  if (iConfig.getParameter<bool>("onDemand")) {
    throw cms::Exception("OnDemandNotAllowed")
        << "Use of the `onDemand` feature of SiPixelClusterShapeCacheProducer is no longer supported";
  }
  produces<SiPixelClusterShapeCache>();
}
 
SiPixelClusterShapeCacheProducer::~SiPixelClusterShapeCacheProducer() {}
 
void SiPixelClusterShapeCacheProducer::fillDescriptions(edm::ConfigurationDescriptions& descriptions) {
  edm::ParameterSetDescription desc;
  desc.add<edm::InputTag>("src", edm::InputTag("siPixelClusters"));
  desc.add<bool>("onDemand", false)->setComment("The on demand feature is no longer supported");
  descriptions.add("siPixelClusterShapeCache", desc);
}
 
void SiPixelClusterShapeCacheProducer::produce(edm::StreamID, edm::Event& iEvent, const edm::EventSetup& iSetup) const {
  edm::Handle<InputCollection> input;
  iEvent.getByToken(token_, input);
 
  const auto& geom = &iSetup.getData(geomToken_);
 
  auto output = std::make_unique<SiPixelClusterShapeCache>(input);
  output->resize(input->data().size());
 
  ClusterData data;  // reused
  ClusterShape clusterShape;
 
  for (const auto& detSet : *input) {
    const GeomDetUnit* genericDet = geom->idToDetUnit(detSet.detId());
    const PixelGeomDetUnit* pixDet = dynamic_cast<const PixelGeomDetUnit*>(genericDet);
    assert(pixDet);
 
    edmNew::DetSet<SiPixelCluster>::const_iterator iCluster = detSet.begin(), endCluster = detSet.end();
    for (; iCluster != endCluster; ++iCluster) {
      SiPixelClusterShapeCache::ClusterRef clusterRef = edmNew::makeRefTo(input, iCluster);
      if (not output->isFilled(clusterRef)) {
        data.size.clear();
        clusterShape.determineShape(*pixDet, *iCluster, data);
        output->insert(clusterRef, data);
      }
    }
  }
  output->shrink_to_fit();
 
  iEvent.put(std::move(output));
}
 
DEFINE_FWK_MODULE(SiPixelClusterShapeCacheProducer);