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 54 of file ClusterClusterMapping.cc.

References overlap(), and mps_fire::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 68 of file ClusterClusterMapping.cc.

References overlap(), and mps_fire::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 82 of file ClusterClusterMapping.cc.

References overlap(), and mps_fire::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(), SpecialClusterImporter< T >::importToBlock(), PFElecTkProducer::isSharingEcalEnergyWithEgSC(), SCAndECALLinker::testLink(), and SCAndHGCalLinker::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 40 of file ClusterClusterMapping.cc.

References reco::SuperCluster::clustersBegin(), reco::SuperCluster::clustersEnd(), edm::ValueMap< T >::contains(), 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) {
    if (pfclusassoc.contains(caloclus->id())) {
      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;
}