Phase2TrackerClusterizer Class Reference

Inheritance diagram for Phase2TrackerClusterizer:

Public Member Functions
	Phase2TrackerClusterizer (const edm::ParameterSet &conf)
void	produce (edm::Event &event, const edm::EventSetup &eventSetup) override
	~Phase2TrackerClusterizer () override
Public Member Functions inherited from edm::stream::EDProducer<>
	EDProducer ()=default
	EDProducer (const EDProducer &)=delete
bool	hasAbilityToProduceInBeginLumis () const final
bool	hasAbilityToProduceInBeginProcessBlocks () const final
bool	hasAbilityToProduceInBeginRuns () const final
bool	hasAbilityToProduceInEndLumis () const final
bool	hasAbilityToProduceInEndProcessBlocks () const final
bool	hasAbilityToProduceInEndRuns () const final
const EDProducer &	operator= (const EDProducer &)=delete

Private Attributes
edm::EDGetTokenT < edm::DetSetVector < Phase2TrackerDigi > >	token_

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
using	CacheTypes = CacheContexts< T...>
using	GlobalCache = typename CacheTypes::GlobalCache
using	HasAbility = AbilityChecker< T...>
using	InputProcessBlockCache = typename CacheTypes::InputProcessBlockCache
using	LuminosityBlockCache = typename CacheTypes::LuminosityBlockCache
using	LuminosityBlockContext = LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >
using	LuminosityBlockSummaryCache = typename CacheTypes::LuminosityBlockSummaryCache
using	RunCache = typename CacheTypes::RunCache
using	RunContext = RunContextT< RunCache, GlobalCache >
using	RunSummaryCache = typename CacheTypes::RunSummaryCache

Detailed Description

Definition at line 30 of file Phase2TrackerClusterizer.cc.

Constructor & Destructor Documentation

Phase2TrackerClusterizer::Phase2TrackerClusterizer ( const edm::ParameterSet &  conf )

explicit

Definition at line 47 of file Phase2TrackerClusterizer.cc.

    :
#ifdef VERIFY_PH2_TK_CLUS
      clusterizer_(new Phase2TrackerClusterizerAlgorithm(conf.getParameter<unsigned int>("maxClusterSize"),
                                                         conf.getParameter<unsigned int>("maxNumberClusters"))),
#endif
      token_(consumes<edm::DetSetVector<Phase2TrackerDigi> >(conf.getParameter<edm::InputTag>("src"))) {
  produces<Phase2TrackerCluster1DCollectionNew>();
}

Phase2TrackerClusterizer::~Phase2TrackerClusterizer ( )

override

Definition at line 57 of file Phase2TrackerClusterizer.cc.

{}

Member Function Documentation

void Phase2TrackerClusterizer::produce ( edm::Event &  event, const edm::EventSetup &  eventSetup  )

override

Definition at line 63 of file Phase2TrackerClusterizer.cc.

References edmNew::DetSetVector< T >::FastFiller::abort(), cms::cuda::assert(), haddnano::cl, HLT_FULL_cff::clusters, gather_cfg::cout, CommonMethods::cp(), edm::EventSetup::get(), mps_fire::i, TrackerGeometry::idToDetUnit(), eostools::move(), slimmedTrackExtras_cff::outputClusters, and token_.

                                                                                       {
  // Get the Digis
  edm::Handle<edm::DetSetVector<Phase2TrackerDigi> > digis;
  event.getByToken(token_, digis);

#ifdef VERIFY_PH2_TK_CLUS
  // Get the geometry
  edm::ESHandle<TrackerGeometry> geomHandle;
  eventSetup.get<TrackerDigiGeometryRecord>().get(geomHandle);
  const TrackerGeometry* tkGeom(&(*geomHandle));
  // Global container for the clusters of each modules
  auto outputClustersOld = std::make_unique<Phase2TrackerCluster1DCollectionNew>();
#endif
  auto outputClusters = std::make_unique<Phase2TrackerCluster1DCollectionNew>();

  // Go over all the modules
  for (const auto& DSViter : *digis) {
    DetId detId(DSViter.detId());

    Phase2TrackerCluster1DCollectionNew::FastFiller clusters(*outputClusters, DSViter.detId());
    Phase2TrackerClusterizerSequentialAlgorithm algo;
    algo.clusterizeDetUnit(DSViter, clusters);
    if (clusters.empty())
      clusters.abort();

#ifdef VERIFY_PH2_TK_CLUS
    if (!clusters.empty()) {
      auto cp = clusters[0].column();
      auto sp = clusters[0].firstStrip();
      for (auto const& cl : clusters) {
        if (cl.column() < cp)
          std::cout << "column not in order! " << std::endl;
        if (cl.column() == cp && cl.firstStrip() < sp)
          std::cout << "strip not in order! " << std::endl;
        cp = cl.column();
        sp = cl.firstStrip();
      }
    }
#endif

#ifdef VERIFY_PH2_TK_CLUS
    // Geometry
    const GeomDetUnit* geomDetUnit(tkGeom->idToDetUnit(detId));
    const PixelGeomDetUnit* pixDet = dynamic_cast<const PixelGeomDetUnit*>(geomDetUnit);
    if (!pixDet)
      assert(0);

    // Container for the clusters that will be produced for this modules
    Phase2TrackerCluster1DCollectionNew::FastFiller clustersOld(*outputClustersOld, DSViter.detId());

    // Setup the clusterizer algorithm for this detector (see ClusterizerAlgorithm for more details)
    clusterizer_->setup(pixDet);

    // Pass the list of Digis to the main algorithm
    // This function will store the clusters in the previously created container
    clusterizer_->clusterizeDetUnit(DSViter, clustersOld);
    if (clustersOld.empty())
      clustersOld.abort();

    if (clusters.size() != clustersOld.size()) {
      std::cout << "SIZEs " << int(detId) << ' ' << clusters.size() << ' ' << clustersOld.size() << std::endl;
      for (auto const& cl : clusters)
        std::cout << cl.size() << ' ' << cl.threshold() << ' ' << cl.firstRow() << ' ' << cl.column() << std::endl;
      std::cout << "Old " << std::endl;
      for (auto const& cl : clustersOld)
        std::cout << cl.size() << ' ' << cl.threshold() << ' ' << cl.firstRow() << ' ' << cl.column() << std::endl;
    }
#endif
  }

#ifdef VERIFY_PH2_TK_CLUS
  // std::cout << "SIZEs " << outputClusters->dataSize() << ' ' << outputClustersOld->dataSize() << std::endl;
  assert(outputClusters->dataSize() == outputClustersOld->dataSize());
  for (auto i = 0U; i < outputClusters->dataSize(); ++i) {
    assert(outputClusters->data()[i].size() == outputClustersOld->data()[i].size());
    assert(outputClusters->data()[i].threshold() == outputClustersOld->data()[i].threshold());
    assert(outputClusters->data()[i].firstRow() == outputClustersOld->data()[i].firstRow());
    assert(outputClusters->data()[i].column() == outputClustersOld->data()[i].column());
  }
#endif

  // Add the data to the output
  outputClusters->shrink_to_fit();
  event.put(std::move(outputClusters));
}

Member Data Documentation

edm::EDGetTokenT<edm::DetSetVector<Phase2TrackerDigi> > Phase2TrackerClusterizer::token_

private

Definition at line 40 of file Phase2TrackerClusterizer.cc.

Referenced by produce().