ClusterClusterMapping Class Reference

#include <ClusterClusterMapping.h>

Public Member Functions
	ClusterClusterMapping ()
	~ClusterClusterMapping ()

Static Public Member Functions
static int	checkOverlap (const reco::PFCluster &pfc, const std::vector< const reco::SuperCluster * > &sc, float minfrac=0.01, bool debug=false)
static int	checkOverlap (const reco::PFCluster &pfc, const std::vector< reco::SuperClusterRef > &sc, float minfrac=0.01, bool debug=false)
static int	checkOverlap (const reco::PFClusterRef &pfc, const std::vector< reco::SuperClusterRef > &sc, const edm::ValueMap< reco::CaloClusterPtr > &pfclusassoc)
static bool	overlap (const reco::CaloCluster &sc1, const reco::CaloCluster &sc, float minfrac=0.01, bool debug=false)
static bool	overlap (const reco::PFClusterRef &pfclustest, const reco::SuperCluster &sc, const edm::ValueMap< reco::CaloClusterPtr > &pfclusassoc)

Detailed Description

Definition at line 12 of file ClusterClusterMapping.h.

Constructor & Destructor Documentation

ClusterClusterMapping::ClusterClusterMapping ( )

inline

Definition at line 14 of file ClusterClusterMapping.h.

{;}

ClusterClusterMapping::~ClusterClusterMapping ( )

inline

Definition at line 15 of file ClusterClusterMapping.h.

{;}

Member Function Documentation

int ClusterClusterMapping::checkOverlap ( const reco::PFCluster &  pfc, const std::vector< const reco::SuperCluster * > &  sc, float  minfrac = 0.01, bool  debug = false  )

static

Definition at line 48 of file ClusterClusterMapping.cc.

References overlap(), and query::result.

                                                                                                                                         {
  int result=-1;
  unsigned nsc=sc.size();
  
  for(unsigned isc=0;isc<nsc;++isc) {
    if(overlap(pfc,*(sc[isc]),minfrac,debug))
      return isc;
  }
  return result;
}

int ClusterClusterMapping::checkOverlap ( const reco::PFCluster &  pfc, const std::vector< reco::SuperClusterRef > &  sc, float  minfrac = 0.01, bool  debug = false  )

static

Definition at line 59 of file ClusterClusterMapping.cc.

References overlap(), and query::result.

                                                                                                                                     {
  int result=-1;
  unsigned nsc=sc.size();
  
  for(unsigned isc=0;isc<nsc;++isc) {
    if(overlap(pfc,*sc[isc],minfrac,debug))
      return isc;
  }
  return result;
}

int ClusterClusterMapping::checkOverlap ( const reco::PFClusterRef &  pfc, const std::vector< reco::SuperClusterRef > &  sc, const edm::ValueMap< reco::CaloClusterPtr > &  pfclusassoc  )

static

Definition at line 70 of file ClusterClusterMapping.cc.

References overlap(), and query::result.

                                                                                                                                                                  {
  int result=-1;
  unsigned nsc=sc.size();
  
  for(unsigned isc=0;isc<nsc;++isc) {
    if(overlap(pfc,*sc[isc],pfclusassoc))
      return isc;
  }
  return result;
}

bool ClusterClusterMapping::overlap ( const reco::CaloCluster &  sc1, const reco::CaloCluster &  sc, float  minfrac = 0.01, bool  debug = false  )

static

Definition at line 4 of file ClusterClusterMapping.cc.

References gather_cfg::cout, first, reco::CaloCluster::hitsAndFractions(), and edm::second().

Referenced by checkOverlap(), ECALClusterImporter::importToBlock(), PFElecTkProducer::isSharingEcalEnergyWithEgSC(), and SCAndECALLinker::testLink().

                                                                                                                       {
  const std::vector< std::pair<DetId, float> > & hits1 = sc1.hitsAndFractions();
  const std::vector< std::pair<DetId, float> > & hits2 = sc2.hitsAndFractions();
  unsigned nhits1=hits1.size();
  unsigned nhits2=hits2.size();
  
  for(unsigned i1=0;i1<nhits1;++i1)
    { 
      // consider only with a minimum fraction of minfrac (default 1%) of the RecHit
      if(hits1[i1].second<minfrac) {
      if(debug) std::cout << " Discarding " << hits1[i1].first << " with " << hits1[i1].second << std::endl;
      continue;
    }
      for(unsigned i2=0;i2<nhits2;++i2)
    {
      // consider only with a minimum fraction of minfract (default 1%) of the RecHit
      if(hits2[i2].second<minfrac ) {
        if(debug) std::cout << " Discarding " << hits2[i2].first << " with " << hits2[i2].second << std::endl;
        continue;
      }
      if(hits1[i1].first==hits2[i2].first)
        {
          if(debug)
        {
          std::cout << " Matching hits " << hits1[i1].first << " with " <<  hits1[i1].second << " and " <<  hits2[i2].first;
          std::cout << " with " << hits2[i2].second << std::endl;         
        }
          return true;
        }
    }
    }
  return false;
}

bool ClusterClusterMapping::overlap ( const reco::PFClusterRef &  pfclustest, const reco::SuperCluster &  sc, const edm::ValueMap< reco::CaloClusterPtr > &  pfclusassoc  )

static

Definition at line 38 of file ClusterClusterMapping.cc.

References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), edm::Ptr< T >::id(), edm::Ref< C, T, F >::id(), edm::Ptr< T >::key(), and edm::Ref< C, T, F >::key().

                                                                                                                                                      {
    
  for (reco::CaloCluster_iterator caloclus = sc.clustersBegin(); caloclus!=sc.clustersEnd(); ++caloclus) {
    const reco::CaloClusterPtr &pfclus = pfclusassoc[*caloclus];
    //printf("pfclus prodidx = %i, key = %i, pfclustest prodidx = %i, key = %i\n",pfclus.id().productIndex(),int(pfclus.key()),pfclustest.id().productIndex(),int(pfclustest.key()));
    if (pfclus.id()==pfclustest.id() && pfclus.key()==pfclustest.key()) return true;    
  }
  return false;
}