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
bool	hasAbilityToProduceInLumis () const final
bool	hasAbilityToProduceInRuns () const final

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

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, haddnano::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(), slimmedTrackExtras_cff::outputClusters, 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().