HGCalMulticlusteringImpl Class Reference

#include <HGCalMulticlusteringImpl.h>

Public Member Functions
void	clusterizeDBSCAN (const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, l1t::HGCalMulticlusterBxCollection &multiclusters, const HGCalTriggerGeometryBase &triggerGeometry)
void	clusterizeDR (const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, l1t::HGCalMulticlusterBxCollection &multiclusters, const HGCalTriggerGeometryBase &triggerGeometry)
void	eventSetup (const edm::EventSetup &es)
	HGCalMulticlusteringImpl (const edm::ParameterSet &conf)
bool	isPertinent (const l1t::HGCalCluster &clu, const l1t::HGCalMulticluster &mclu, double dR) const

Private Member Functions
void	findNeighbor (const std::vector< std::pair< unsigned int, double >> &rankedList, unsigned int searchInd, const std::vector< edm::Ptr< l1t::HGCalCluster >> &clustersPtr, std::vector< unsigned int > &neigbors)

Private Attributes
double	distDbscan_ = 0.005
double	dr_
unsigned	minNDbscan_ = 3
std::string	multiclusterAlgoType_
double	ptC3dThreshold_
HGCalShowerShape	shape_
HGCalTriggerTools	triggerTools_

Detailed Description

Definition at line 14 of file HGCalMulticlusteringImpl.h.

Constructor & Destructor Documentation

HGCalMulticlusteringImpl::HGCalMulticlusteringImpl

(

const edm::ParameterSet &

conf

)

Definition at line 8 of file HGCalMulticlusteringImpl.cc.

References distDbscan_, dr_, minNDbscan_, multiclusterAlgoType_, and ptC3dThreshold_.

                                                                                :
    dr_(conf.getParameter<double>("dR_multicluster")),
    ptC3dThreshold_(conf.getParameter<double>("minPt_multicluster")),
    multiclusterAlgoType_(conf.getParameter<string>("type_multicluster")),
    distDbscan_(conf.getParameter<double>("dist_dbscan_multicluster")),
    minNDbscan_(conf.getParameter<unsigned>("minN_dbscan_multicluster"))
{    
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster dR for Near Neighbour search: " << dr_;  
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster minimum transverse-momentum: " << ptC3dThreshold_;
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster DBSCAN Clustering distance: " << distDbscan_;
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster clustering min number of subclusters: " << minNDbscan_;
    edm::LogInfo("HGCalMulticlusterParameters") << "Multicluster type of multiclustering algortihm: " << multiclusterAlgoType_;
}

Member Function Documentation

void HGCalMulticlusteringImpl::clusterizeDBSCAN	(	const std::vector< edm::Ptr< l1t::HGCalCluster >> &	clustersPtr,
		l1t::HGCalMulticlusterBxCollection &	multiclusters,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)

Definition at line 147 of file HGCalMulticlusteringImpl.cc.

References HGCalShowerShape::coreShowerLength(), HGCalShowerShape::eMax(), PVValHelper::eta, findNeighbor(), HGCalShowerShape::firstLayer(), mps_fire::i, HGCalShowerShape::maxLayer(), minNDbscan_, eostools::move(), phi, EnergyCorrector::pt, ptC3dThreshold_, BXVector< T >::push_back(), shape_, HGCalShowerShape::showerLength(), HGCalShowerShape::sigmaEtaEtaMax(), HGCalShowerShape::sigmaEtaEtaTot(), HGCalShowerShape::sigmaPhiPhiMax(), HGCalShowerShape::sigmaPhiPhiTot(), HGCalShowerShape::sigmaRRMax(), HGCalShowerShape::sigmaRRMean(), HGCalShowerShape::sigmaRRTot(), HGCalShowerShape::sigmaZZ(), TtFullHadEvtBuilder_cfi::sumPt, class-composition::visited, and HcalDigiParam_cfi::zside.

Referenced by eventSetup().

{

    std::vector<l1t::HGCalMulticluster> multiclustersTmp;
    l1t::HGCalMulticluster mcluTmp;
    std::vector<bool> visited(clustersPtrs.size(),false);
    std::vector<bool> merged (clustersPtrs.size(),false);
    std::vector<std::pair<unsigned int,double>>  rankedList;
    rankedList.reserve(clustersPtrs.size());
    std::vector<std::vector<unsigned int>> neighborList;
    neighborList.reserve(clustersPtrs.size());

    int iclu = 0, imclu = 0, neighNo;
    double dist = 0.;

    for(std::vector<edm::Ptr<l1t::HGCalCluster>>::const_iterator clu = clustersPtrs.begin(); clu != clustersPtrs.end(); ++clu, ++iclu){
        dist = (*clu)->centreProj().mag()*HGCalDetId((*clu)->detId()).zside();
        rankedList.push_back(std::make_pair(iclu,dist));
    }  
    iclu = 0;
    std::sort(rankedList.begin(), rankedList.end(), [](auto &left, auto &right) {
            return left.second < right.second;
            });

    for(const auto& cluRanked: rankedList){
        std::vector<unsigned int> neighbors;      

        if(!visited.at(iclu)){
            visited.at(iclu) = true;
            findNeighbor(rankedList, iclu, clustersPtrs, neighbors);
            neighborList.push_back(std::move(neighbors));

            if(neighborList.at(iclu).size() >= minNDbscan_) {
                multiclustersTmp.emplace_back( clustersPtrs[cluRanked.first] );
                merged.at(iclu) = true;
                /* dynamic range loop: range-based loop syntax cannot be employed */
                for(unsigned int neighInd = 0; neighInd < neighborList.at(iclu).size(); neighInd++){
                    neighNo = neighborList.at(iclu).at(neighInd);
                    /* This condition also ensures merging of clusters visited by other clusters but not merged. */
                    if(!merged.at(neighNo) ){
                        merged.at(neighNo) = true;          
                        multiclustersTmp.at(imclu).addConstituent( clustersPtrs[rankedList.at(neighNo).first] );

                        if(!visited.at(neighNo)){
                            visited.at(neighNo) = true;
                            std::vector<unsigned int> secNeighbors;
                            findNeighbor(rankedList, neighNo,clustersPtrs, secNeighbors);

                            if(secNeighbors.size() >= minNDbscan_){
                                neighborList.at(iclu).insert(neighborList.at(iclu).end(), secNeighbors.begin(), secNeighbors.end());
                            }
                        }
                    }
                }
                imclu++;
            }
        }
        else neighborList.push_back(std::move(neighbors));
        iclu++;    
    }
    /* making the collection of multiclusters */
    for( unsigned i(0); i<multiclustersTmp.size(); ++i ){

        // compute the eta, phi observables for multicluster starting from its barycenter x,y,z position + pT as scalar sum of pT of constituents
        double sumPt=0.;

        const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>> &pertinentClu = multiclustersTmp[i].constituents();
        for(const auto& id_cluster : pertinentClu) sumPt += id_cluster.second->pt();

        math::PtEtaPhiMLorentzVector multiclusterP4(  sumPt,
                multiclustersTmp[i].centre().eta(),
                multiclustersTmp[i].centre().phi(),
                0. );
        multiclustersTmp[i].setP4( multiclusterP4 );


        if( multiclustersTmp[i].pt() > ptC3dThreshold_ ){

            //compute shower shape
            multiclustersTmp[i].showerLength(shape_.showerLength(multiclustersTmp[i]));
            multiclustersTmp[i].coreShowerLength(shape_.coreShowerLength(multiclustersTmp[i], triggerGeometry));
            multiclustersTmp[i].firstLayer(shape_.firstLayer(multiclustersTmp[i]));
            multiclustersTmp[i].maxLayer(shape_.maxLayer(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaEtaEtaTot(shape_.sigmaEtaEtaTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaEtaEtaMax(shape_.sigmaEtaEtaMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaPhiPhiTot(shape_.sigmaPhiPhiTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaPhiPhiMax(shape_.sigmaPhiPhiMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaZZ(shape_.sigmaZZ(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRTot(shape_.sigmaRRTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRMax(shape_.sigmaRRMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRMean(shape_.sigmaRRMean(multiclustersTmp[i]));
            multiclustersTmp[i].eMax(shape_.eMax(multiclustersTmp[i]));

            multiclusters.push_back( 0, multiclustersTmp[i]);
        }
    }
}

void HGCalMulticlusteringImpl::clusterizeDR	(	const std::vector< edm::Ptr< l1t::HGCalCluster >> &	clustersPtr,
		l1t::HGCalMulticlusterBxCollection &	multiclusters,
		const HGCalTriggerGeometryBase &	triggerGeometry
	)

Definition at line 73 of file HGCalMulticlusteringImpl.cc.

References HGCalShowerShape::coreShowerLength(), edmIntegrityCheck::d, dr_, HGCalShowerShape::eMax(), PVValHelper::eta, HGCalShowerShape::firstLayer(), mps_fire::i, isPertinent(), mag(), HGCalShowerShape::maxLayer(), phi, EnergyCorrector::pt, ptC3dThreshold_, BXVector< T >::push_back(), shape_, HGCalShowerShape::showerLength(), HGCalShowerShape::sigmaEtaEtaMax(), HGCalShowerShape::sigmaEtaEtaTot(), HGCalShowerShape::sigmaPhiPhiMax(), HGCalShowerShape::sigmaPhiPhiTot(), HGCalShowerShape::sigmaRRMax(), HGCalShowerShape::sigmaRRMean(), HGCalShowerShape::sigmaRRTot(), HGCalShowerShape::sigmaZZ(), and TtFullHadEvtBuilder_cfi::sumPt.

Referenced by eventSetup().

{

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

    int iclu = 0;
    for(std::vector<edm::Ptr<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.empty() ){
            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 */
    for( unsigned i(0); i<multiclustersTmp.size(); ++i ){

        // compute the eta, phi observables for multicluster starting from its barycenter x,y,z position + pT as scalar sum of pT of constituents
        double sumPt=0.;

        const std::unordered_map<uint32_t, edm::Ptr<l1t::HGCalCluster>> &pertinentClu = multiclustersTmp[i].constituents();
        for(const auto& id_cluster : pertinentClu) sumPt += id_cluster.second->pt();
        math::PtEtaPhiMLorentzVector multiclusterP4(  sumPt,
                                                      multiclustersTmp[i].centre().eta(),
                                                      multiclustersTmp[i].centre().phi(),
                                                      0. );
        multiclustersTmp[i].setP4( multiclusterP4 );


        if( multiclustersTmp[i].pt() > ptC3dThreshold_ ){

            //compute shower shape
            multiclustersTmp[i].showerLength(shape_.showerLength(multiclustersTmp[i]));
            multiclustersTmp[i].coreShowerLength(shape_.coreShowerLength(multiclustersTmp[i], triggerGeometry));
            multiclustersTmp[i].firstLayer(shape_.firstLayer(multiclustersTmp[i]));
            multiclustersTmp[i].maxLayer(shape_.maxLayer(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaEtaEtaTot(shape_.sigmaEtaEtaTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaEtaEtaMax(shape_.sigmaEtaEtaMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaPhiPhiTot(shape_.sigmaPhiPhiTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaPhiPhiMax(shape_.sigmaPhiPhiMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaZZ(shape_.sigmaZZ(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRTot(shape_.sigmaRRTot(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRMax(shape_.sigmaRRMax(multiclustersTmp[i]));
            multiclustersTmp[i].sigmaRRMean(shape_.sigmaRRMean(multiclustersTmp[i]));
            multiclustersTmp[i].eMax(shape_.eMax(multiclustersTmp[i]));

            multiclusters.push_back( 0, multiclustersTmp[i]);
        }
    }
    
}

void HGCalMulticlusteringImpl::eventSetup ( const edm::EventSetup &  es )

inline

Definition at line 20 of file HGCalMulticlusteringImpl.h.

References clusterizeDBSCAN(), clusterizeDR(), PFRecoTauDiscriminationAgainstElectronDeadECAL_cfi::dR, HGCalShowerShape::eventSetup(), HGCalTriggerTools::eventSetup(), findNeighbor(), isPertinent(), shape_, and triggerTools_.

    {
        triggerTools_.eventSetup(es);
        shape_.eventSetup(es);
    }

void HGCalMulticlusteringImpl::findNeighbor ( const std::vector< std::pair< unsigned int, double >> &  rankedList, unsigned int  searchInd, const std::vector< edm::Ptr< l1t::HGCalCluster >> &  clustersPtr, std::vector< unsigned int > &  neigbors  )

private

Definition at line 41 of file HGCalMulticlusteringImpl.cc.

References distDbscan_, Exception, and mag().

Referenced by clusterizeDBSCAN(), and eventSetup().

                                             {

    if(clustersPtrs.size() <= searchInd || clustersPtrs.size() < rankedList.size()){
        throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");
    }

    for(unsigned int ind = searchInd+1; ind < rankedList.size() && fabs(rankedList.at(ind).second - rankedList.at(searchInd).second) < distDbscan_ ; ind++){

        if(clustersPtrs.size() <= rankedList.at(ind).first){
            throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");

        } else if(((*(clustersPtrs[rankedList.at(ind).first])).centreProj() - (*(clustersPtrs[rankedList.at(searchInd).first])).centreProj()).mag() < distDbscan_){
            neighbors.push_back(ind);
        }
    }

    for(unsigned int ind = 0; ind < searchInd && fabs(rankedList.at(searchInd).second - rankedList.at(ind).second) < distDbscan_ ; ind++){

        if(clustersPtrs.size() <= rankedList.at(ind).first){
            throw cms::Exception("IndexOutOfBound: clustersPtrs in 'findNeighbor'");

        } else if(((*(clustersPtrs[rankedList.at(ind).first])).centreProj() - (*(clustersPtrs[rankedList.at(searchInd).first])).centreProj()).mag() < distDbscan_){
            neighbors.push_back(ind);
        }
    }
}

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

Definition at line 23 of file HGCalMulticlusteringImpl.cc.

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

Referenced by clusterizeDR(), and eventSetup().

{
    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::distDbscan_ = 0.005

private

Definition at line 50 of file HGCalMulticlusteringImpl.h.

Referenced by findNeighbor(), and HGCalMulticlusteringImpl().

double HGCalMulticlusteringImpl::dr_

private

Definition at line 47 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDR(), and HGCalMulticlusteringImpl().

unsigned HGCalMulticlusteringImpl::minNDbscan_ = 3

private

Definition at line 51 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDBSCAN(), and HGCalMulticlusteringImpl().

std::string HGCalMulticlusteringImpl::multiclusterAlgoType_

private

Definition at line 49 of file HGCalMulticlusteringImpl.h.

Referenced by HGCalMulticlusteringImpl().

double HGCalMulticlusteringImpl::ptC3dThreshold_

private

Definition at line 48 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDBSCAN(), clusterizeDR(), and HGCalMulticlusteringImpl().

HGCalShowerShape HGCalMulticlusteringImpl::shape_

private

Definition at line 53 of file HGCalMulticlusteringImpl.h.

Referenced by clusterizeDBSCAN(), clusterizeDR(), and eventSetup().

HGCalTriggerTools HGCalMulticlusteringImpl::triggerTools_

private

Definition at line 54 of file HGCalMulticlusteringImpl.h.

Referenced by eventSetup().