HGCalMulticlusteringImpl Class Reference

#include <HGCalMulticlusteringImpl.h>

Public Member Functions
void	clusterize (const edm::PtrVector< l1t::HGCalCluster > &clustersPtr, l1t::HGCalMulticlusterBxCollection &multiclusters)
	HGCalMulticlusteringImpl (const edm::ParameterSet &conf)
bool	isPertinent (const l1t::HGCalCluster &clu, const l1t::HGCalMulticluster &mclu, double dR) const

Private Attributes
double	dr_

Detailed Description

Definition at line 9 of file HGCalMulticlusteringImpl.h.

Constructor & Destructor Documentation

HGCalMulticlusteringImpl::HGCalMulticlusteringImpl

(

const edm::ParameterSet &

conf

)

Definition at line 5 of file HGCalMulticlusteringImpl.cc.

References dr_.

                                                                                :
    dr_(conf.getParameter<double>("dR_multicluster"))
{    
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster dR for Near Neighbour search: " << dr_;  
}

Member Function Documentation

void HGCalMulticlusteringImpl::clusterize	(	const edm::PtrVector< l1t::HGCalCluster > &	clustersPtr,
		l1t::HGCalMulticlusterBxCollection &	multiclusters
	)

Definition at line 30 of file HGCalMulticlusteringImpl.cc.

References edm::PtrVector< T >::begin(), edmIntegrityCheck::d, dr_, edm::PtrVector< T >::end(), mps_fire::i, isPertinent(), mag(), BXVector< T >::resize(), and BXVector< T >::set().

{
           
    std::vector<l1t::HGCalMulticluster> multiclustersTmp;

    int iclu = 0;
    for(edm::PtrVector<l1t::HGCalCluster>::const_iterator clu = clustersPtrs.begin(); clu != clustersPtrs.end(); ++clu, ++iclu){
        
        int imclu=0;
        vector<int> tcPertinentMulticlusters;
        for( const auto& mclu : multiclustersTmp ){
            if( this->isPertinent(**clu, mclu, dr_) ){
                tcPertinentMulticlusters.push_back(imclu);
            }
            ++imclu;
        }
        if( tcPertinentMulticlusters.size() == 0 ){
            multiclustersTmp.emplace_back( *clu );
        }
        else{
            unsigned minDist = 1;
            unsigned targetMulticlu = 0; 
            for( int imclu : tcPertinentMulticlusters ){
                double d = ( multiclustersTmp.at(imclu).centreProj() - (*clu)->centreProj() ).mag() ;
                if( d < minDist ){
                    minDist = d;
                    targetMulticlu = imclu;
                }
            } 

            multiclustersTmp.at( targetMulticlu ).addConstituent( *clu );
            
        }        
    }

    /* making the collection of multiclusters */
    multiclusters.resize(0, multiclustersTmp.size());
    for( unsigned i(0); i<multiclustersTmp.size(); ++i ){
        multiclusters.set( 0, i, multiclustersTmp.at(i) );
    }
    
}

bool HGCalMulticlusteringImpl::isPertinent	(	const l1t::HGCalCluster &	clu,
		const l1t::HGCalMulticluster &	mclu,
		double	dR
	)		const

Definition at line 12 of file HGCalMulticlusteringImpl.cc.

References l1t::HGCalClusterT< C >::centreProj(), l1t::HGCalClusterT< C >::detId(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, and mag().

Referenced by clusterize().

{
    HGCalDetId cluDetId( clu.detId() );
    HGCalDetId firstClusterDetId( mclu.detId() );
    
    if( cluDetId.zside() != firstClusterDetId.zside() ){
        return false;
    }
    if( ( mclu.centreProj() - clu.centreProj() ).mag() < dR ){
        return true;
    }
    return false;

}

Member Data Documentation

double HGCalMulticlusteringImpl::dr_

private

Definition at line 24 of file HGCalMulticlusteringImpl.h.

Referenced by clusterize(), and HGCalMulticlusteringImpl().