HGCClusterAlgo< FECODEC, DATA > Class Template Reference

Inheritance diagram for HGCClusterAlgo< FECODEC, DATA >:

Public Member Functions
	HGCClusterAlgo (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
void	putInEvent (edm::Event &evt) final
void	reset () final
void	run (const l1t::HGCFETriggerDigiCollection &coll, const edm::EventSetup &es, edm::Event &evt) final
void	setProduces (edm::stream::EDProducer<> &prod) const final
Public Member Functions inherited from HGCalTriggerBackend::Algorithm< FECODEC >
	Algorithm (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
void	setGeometry (const HGCalTriggerGeometryBase *const geom) final
Public Member Functions inherited from HGCalTriggerBackendAlgorithmBase
	HGCalTriggerBackendAlgorithmBase (const edm::ParameterSet &conf, edm::ConsumesCollector &cc)
const std::string &	name () const
virtual	~HGCalTriggerBackendAlgorithmBase ()

Private Types
enum	ClusterType { dRC2d, NNC2d, dRNNC2d }
enum	MulticlusterType { dRC3d, DBSCANC3d }

Private Attributes
HGCalTriggerCellCalibration	calibration_
std::unique_ptr< l1t::HGCalClusterBxCollection >	cluster_product_
HGCalClusteringImpl	clustering_
ClusterType	clusteringAlgorithmType_
std::unique_ptr< l1t::HGCalMulticlusterBxCollection >	multicluster_product_
HGCalMulticlusteringImpl	multiclustering_
MulticlusterType	multiclusteringAlgoType_
std::unique_ptr< l1t::HGCalTriggerCellBxCollection >	trgcell_product_
double	triggercell_threshold_scintillator_
double	triggercell_threshold_silicon_
edm::ESHandle< HGCalTriggerGeometryBase >	triggerGeometry_

Additional Inherited Members
Protected Attributes inherited from HGCalTriggerBackend::Algorithm< FECODEC >
FECODEC	codec_
Protected Attributes inherited from HGCalTriggerBackendAlgorithmBase
const HGCalTriggerGeometryBase *	geometry_

Detailed Description

template<typename FECODEC, typename DATA>
class HGCClusterAlgo< FECODEC, DATA >

Definition at line 18 of file HGCClusterAlgo.cc.

Member Enumeration Documentation

template<typename FECODEC , typename DATA >

enum HGCClusterAlgo::ClusterType

private

Enumerator
dRC2d
NNC2d
dRNNC2d

Definition at line 27 of file HGCClusterAlgo.cc.

                        {
            dRC2d,
            NNC2d,
            dRNNC2d
        };

template<typename FECODEC , typename DATA >

enum HGCClusterAlgo::MulticlusterType

private

Enumerator
dRC3d
DBSCANC3d

Definition at line 32 of file HGCClusterAlgo.cc.

                             {
            dRC3d,
            DBSCANC3d
        };

Constructor & Destructor Documentation

template<typename FECODEC , typename DATA >

HGCClusterAlgo< FECODEC, DATA >::HGCClusterAlgo ( const edm::ParameterSet &  conf, edm::ConsumesCollector &  cc  )

inline

Definition at line 39 of file HGCClusterAlgo.cc.

References Exception, edm::ParameterSet::getParameter(), edm::ParameterSet::getParameterSet(), and AlCaHLTBitMon_QueryRunRegistry::string.

                                                                              :
        Algorithm<FECODEC>(conf, cc),
        trgcell_product_( new l1t::HGCalTriggerCellBxCollection ),
        cluster_product_( new l1t::HGCalClusterBxCollection ),
        multicluster_product_( new l1t::HGCalMulticlusterBxCollection ),
        calibration_( conf.getParameterSet("calib_parameters") ),
        clustering_( conf.getParameterSet("C2d_parameters") ),
        multiclustering_( conf.getParameterSet("C3d_parameters" ) ),
        triggercell_threshold_silicon_( conf.getParameter<double>("triggercell_threshold_silicon") ),
        triggercell_threshold_scintillator_( conf.getParameter<double>("triggercell_threshold_scintillator") )
        {
            std::string typeCluster(conf.getParameterSet("C2d_parameters").getParameter<std::string>("clusterType"));
            if(typeCluster=="dRC2d"){
                clusteringAlgorithmType_ = dRC2d;
            }else if(typeCluster=="NNC2d"){
                clusteringAlgorithmType_ = NNC2d;
            }else if(typeCluster=="dRNNC2d"){
                clusteringAlgorithmType_ = dRNNC2d;
            }else {
                throw cms::Exception("HGCTriggerParameterError")
                    << "Unknown clustering type '" << typeCluster;
            }
            std::string typeMulticluster(conf.getParameterSet("C3d_parameters").getParameter<std::string>("type_multicluster"));
            if(typeMulticluster=="dRC3d"){
                multiclusteringAlgoType_ = dRC3d;
            }else if(typeMulticluster=="DBSCANC3d"){
                multiclusteringAlgoType_ = DBSCANC3d;
            }else {
                throw cms::Exception("HGCTriggerParameterError")
                    << "Unknown Multiclustering type '" << typeMulticluster;
            }

        }

Member Function Documentation

template<typename FECODEC , typename DATA >

void HGCClusterAlgo< FECODEC, DATA >::putInEvent ( edm::Event &  evt )

inlinefinalvirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 83 of file HGCClusterAlgo.cc.

        {

        }

template<typename FECODEC , typename DATA >

void HGCClusterAlgo< FECODEC, DATA >::reset ( void  )

inlinefinalvirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 89 of file HGCClusterAlgo.cc.

        {
            trgcell_product_.reset( new l1t::HGCalTriggerCellBxCollection );            
            cluster_product_.reset( new l1t::HGCalClusterBxCollection );
            multicluster_product_.reset( new l1t::HGCalMulticlusterBxCollection );
        }

template<typename FECODEC , typename DATA >

void HGCClusterAlgo< FECODEC, DATA >::run ( const l1t::HGCFETriggerDigiCollection &  coll, const edm::EventSetup &  es, edm::Event &  evt  )

finalvirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 120 of file HGCClusterAlgo.cc.

References cscdqm::DATA, data, edm::EventSetup::get(), HGCHEB, mps_fire::i, l1t::HGCalTriggerCell::mipPt(), eostools::move(), edm::Event::put(), and BXVector< T >::size().

{
    es.get<CaloGeometryRecord>().get("", triggerGeometry_);
    calibration_.eventSetup(es);
    clustering_.eventSetup(es);
    multiclustering_.eventSetup(es);

    for( const auto& digi : coll ){
        
        HGCalDetId module_id( digi.id() );
        
        DATA data;
        data.reset();
        digi.decode(codec_, data);
        
        for(const auto& triggercell : data.payload)
        {
            
            if( triggercell.hwPt() > 0 )
            {
                l1t::HGCalTriggerCell calibratedtriggercell( triggercell );
                calibration_.calibrateInGeV( calibratedtriggercell); 
                double triggercell_threshold = (triggercell.subdetId()==HGCHEB ? triggercell_threshold_scintillator_ : triggercell_threshold_silicon_);
                if(calibratedtriggercell.mipPt()<triggercell_threshold) continue;
                trgcell_product_->push_back( 0, calibratedtriggercell );
            }           
        
        }
    
    }

    /* orphan handles to the collections of trigger-cells, clusters and multiclusters */
    edm::OrphanHandle<l1t::HGCalTriggerCellBxCollection> triggerCellsHandle;
    edm::OrphanHandle<l1t::HGCalClusterBxCollection> clustersHandle;
    edm::OrphanHandle<l1t::HGCalMulticlusterBxCollection> multiclustersHandle;
    
    /* retrieve the orphan handle to the trigger-cells collection and put the collection in the event */
    triggerCellsHandle = evt.put( std::move( trgcell_product_ ), "calibratedTriggerCells");

    /* create a persistent vector of pointers to the trigger-cells */
    std::vector<edm::Ptr<l1t::HGCalTriggerCell>> triggerCellsPtrs;
    for( unsigned i = 0; i < triggerCellsHandle->size(); ++i ) {
        edm::Ptr<l1t::HGCalTriggerCell> ptr(triggerCellsHandle,i);
        triggerCellsPtrs.push_back(ptr);
    }
    
    /* call to C2d clustering */
    switch(clusteringAlgorithmType_){
        case dRC2d : 
            clustering_.clusterizeDR( triggerCellsPtrs, *cluster_product_);
            break;
        case NNC2d:
            clustering_.clusterizeNN( triggerCellsPtrs, *cluster_product_, *triggerGeometry_ );
            break;
        case dRNNC2d:
            clustering_.clusterizeDRNN( triggerCellsPtrs, *cluster_product_, *triggerGeometry_ );
            break;
        default:
            // Should not happen, clustering type checked in constructor
            break;
    }

    /* retrieve the orphan handle to the clusters collection and put the collection in the event */
    clustersHandle = evt.put( std::move( cluster_product_ ), "cluster2D");

    /* create a persistent vector of pointers to the trigger-cells */
    std::vector<edm::Ptr<l1t::HGCalCluster>> clustersPtrs;
    for( unsigned i = 0; i < clustersHandle->size(); ++i ) {
        edm::Ptr<l1t::HGCalCluster> ptr(clustersHandle,i);
        clustersPtrs.push_back(ptr);
    }
    
    /* call to multiclustering and compute shower shape*/
    switch(multiclusteringAlgoType_){
        case dRC3d : 
            multiclustering_.clusterizeDR( clustersPtrs, *multicluster_product_, *triggerGeometry_);
            break;
        case DBSCANC3d:
            multiclustering_.clusterizeDBSCAN( clustersPtrs, *multicluster_product_, *triggerGeometry_);
            break;
        default:
            // Should not happen, clustering type checked in constructor
            break;
    }

    /* retrieve the orphan handle to the multiclusters collection and put the collection in the event */
    multiclustersHandle = evt.put( std::move( multicluster_product_ ), "cluster3D");



}

template<typename FECODEC , typename DATA >

void HGCClusterAlgo< FECODEC, DATA >::setProduces ( edm::stream::EDProducer<> &  prod ) const

inlinefinalvirtual

Implements HGCalTriggerBackendAlgorithmBase.

Definition at line 73 of file HGCClusterAlgo.cc.

References coll, parseEventContent::prod, and findQualityFiles::run.

        {
            prod.produces<l1t::HGCalTriggerCellBxCollection>( "calibratedTriggerCells" );            
            prod.produces<l1t::HGCalClusterBxCollection>( "cluster2D" );
            prod.produces<l1t::HGCalMulticlusterBxCollection>( "cluster3D" );   
        }

Member Data Documentation

template<typename FECODEC , typename DATA >

HGCalTriggerCellCalibration HGCClusterAlgo< FECODEC, DATA >::calibration_

private

Definition at line 107 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

std::unique_ptr<l1t::HGCalClusterBxCollection> HGCClusterAlgo< FECODEC, DATA >::cluster_product_

private

Definition at line 101 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

HGCalClusteringImpl HGCClusterAlgo< FECODEC, DATA >::clustering_

private

Definition at line 108 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

ClusterType HGCClusterAlgo< FECODEC, DATA >::clusteringAlgorithmType_

private

Definition at line 112 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

std::unique_ptr<l1t::HGCalMulticlusterBxCollection> HGCClusterAlgo< FECODEC, DATA >::multicluster_product_

private

Definition at line 102 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

HGCalMulticlusteringImpl HGCClusterAlgo< FECODEC, DATA >::multiclustering_

private

Definition at line 109 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

MulticlusterType HGCClusterAlgo< FECODEC, DATA >::multiclusteringAlgoType_

private

Definition at line 115 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

std::unique_ptr<l1t::HGCalTriggerCellBxCollection> HGCClusterAlgo< FECODEC, DATA >::trgcell_product_

private

Definition at line 100 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

double HGCClusterAlgo< FECODEC, DATA >::triggercell_threshold_scintillator_

private

Definition at line 114 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

double HGCClusterAlgo< FECODEC, DATA >::triggercell_threshold_silicon_

private

Definition at line 113 of file HGCClusterAlgo.cc.

template<typename FECODEC , typename DATA >

edm::ESHandle<HGCalTriggerGeometryBase> HGCClusterAlgo< FECODEC, DATA >::triggerGeometry_

private

Definition at line 104 of file HGCClusterAlgo.cc.