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
Public Member Functions inherited from edm::stream::EDProducerBase
	EDProducerBase ()
ModuleDescription const &	moduleDescription () const
	~EDProducerBase () override
Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)
std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const
	ProducerBase ()
std::vector< edm::ProductResolverIndex > const &	putTokenIndexToProductResolverIndex () const
void	registerProducts (ProducerBase *, ProductRegistry *, ModuleDescription const &)
std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...
void	resolvePutIndicies (BranchType iBranchType, ModuleToResolverIndicies const &iIndicies, std::string const &moduleLabel)
virtual	~ProducerBase () noexcept(false)
Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const
void	convertCurrentProcessAlias (std::string const &processName)
	Convert "@currentProcess" in InputTag process names to the actual current process name. More...
	EDConsumerBase ()
	EDConsumerBase (EDConsumerBase const &)=delete
	EDConsumerBase (EDConsumerBase &&)=default
ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const
void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const
std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const
void	labelsForToken (EDGetToken iToken, Labels &oLabels) const
void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const
EDConsumerBase const &	operator= (EDConsumerBase const &)=delete
EDConsumerBase &	operator= (EDConsumerBase &&)=default
bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const
bool	registeredToConsumeMany (TypeID const &, BranchType) const
void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)
virtual	~EDConsumerBase () noexcept(false)

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

Additional Inherited Members
Public Types inherited from edm::stream::EDProducer<>
typedef CacheContexts< T... >	CacheTypes
typedef CacheTypes::GlobalCache	GlobalCache
typedef AbilityChecker< T... >	HasAbility
typedef CacheTypes::LuminosityBlockCache	LuminosityBlockCache
typedef LuminosityBlockContextT< LuminosityBlockCache, RunCache, GlobalCache >	LuminosityBlockContext
typedef CacheTypes::LuminosityBlockSummaryCache	LuminosityBlockSummaryCache
typedef CacheTypes::RunCache	RunCache
typedef RunContextT< RunCache, GlobalCache >	RunContext
typedef CacheTypes::RunSummaryCache	RunSummaryCache
Public Types inherited from edm::stream::EDProducerBase
typedef EDProducerAdaptorBase	ModuleType
Public Types inherited from edm::ProducerBase
using	ModuleToResolverIndicies = std::unordered_multimap< std::string, std::tuple< edm::TypeID const *, const char *, edm::ProductResolverIndex >>
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList
Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels
Static Public Member Functions inherited from edm::stream::EDProducerBase
static const std::string &	baseType ()
static void	fillDescriptions (ConfigurationDescriptions &descriptions)
static void	prevalidate (ConfigurationDescriptions &descriptions)
Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)
EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)
ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...
template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()
void	consumesMany (const TypeToGet &id)
template<BranchType B>
void	consumesMany (const TypeToGet &id)
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)
template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 32 of file Phase2TrackerClusterizer.cc.

Constructor & Destructor Documentation

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

explicit

Definition at line 52 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 60 of file Phase2TrackerClusterizer.cc.

{ }

Member Function Documentation

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

override

Definition at line 66 of file Phase2TrackerClusterizer.cc.

References edmNew::DetSetVector< T >::FastFiller::abort(), patPFMETCorrections_cff::algo, GetRecoTauVFromDQM_MC_cff::cl, fastPrimaryVertexProducer_cfi::clusters, gather_cfg::cout, SimDataFormats::CaloAnalysis::cp, DEFINE_FWK_MODULE, edm::EventSetup::get(), mps_fire::i, TrackerGeometry::idToDetUnit(), createfilelist::int, eostools::move(), token_, and mitigatedMETSequence_cff::U.

                                                                                           {

        // 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 (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 43 of file Phase2TrackerClusterizer.cc.

Referenced by produce().