reco::PFBlock Class Reference

Block of elements. More...

#include <PFBlock.h>

Classes
struct	Link

Public Types
typedef edm::OwnVector < reco::PFBlockElement > ::const_iterator	IE
typedef std::map< unsigned int, Link >	LinkData
enum	LinkTest { LINKTEST_RECHIT, LINKTEST_NLINKTEST, LINKTEST_ALL }

Public Member Functions
void	addElement (reco::PFBlockElement *element)
void	associatedElements (unsigned i, const LinkData &linkData, std::multimap< double, unsigned > &sortedAssociates, reco::PFBlockElement::Type type=PFBlockElement::NONE, LinkTest test=LINKTEST_RECHIT) const
void	bookLinkData ()
double	dist (unsigned ie1, unsigned ie2, const LinkData &linkData, LinkTest test) const
double	dist (unsigned ie1, unsigned ie2, const LinkData &linkData) const
const edm::OwnVector < reco::PFBlockElement > &	elements () const
const LinkData &	linkData () const
LinkData &	linkData ()
bool	matrix2vector (unsigned i, unsigned j, unsigned &index) const
	PFBlock ()
void	setLink (unsigned i1, unsigned i2, double dist, LinkData &linkData, LinkTest test=LINKTEST_RECHIT) const

Private Member Functions
unsigned	linkDataSize () const

Private Attributes
edm::OwnVector < reco::PFBlockElement >	elements_
	all elements More...
LinkData	linkData_
	link data (permanent) More...

Friends
std::ostream &	operator<< (std::ostream &out, const PFBlock &co)

Detailed Description

Block of elements.

Author

Colin Bernet

Date

January 2006

A PFBlock is:

• a set of topologically connected elements.
• a set of links between these elements

Definition at line 30 of file PFBlock.h.

Member Typedef Documentation

typedef edm::OwnVector< reco::PFBlockElement >::const_iterator reco::PFBlock::IE

Definition at line 42 of file PFBlock.h.

typedef std::map< unsigned int, Link > reco::PFBlock::LinkData

Definition at line 46 of file PFBlock.h.

Member Enumeration Documentation

enum reco::PFBlock::LinkTest

Enumerator
LINKTEST_RECHIT
LINKTEST_NLINKTEST
LINKTEST_ALL

Definition at line 48 of file PFBlock.h.

                  {
      LINKTEST_RECHIT,
      LINKTEST_NLINKTEST,
      LINKTEST_ALL
    };

Constructor & Destructor Documentation

reco::PFBlock::PFBlock ( )

inline

Definition at line 54 of file PFBlock.h.

{}

Member Function Documentation

void PFBlock::addElement

(

reco::PFBlockElement *

element

)

add an element to the current PFBlock the block will keep a copy.

Definition at line 12 of file PFBlock.cc.

References reco::PFBlockElement::clone(), reco::PFBlockElement::lock(), and reco::PFBlockElement::setIndex().

                                                   {
  element->setIndex( elements_.size() );
  element->lock();
  elements_.push_back( element->clone() ); 

}

void PFBlock::associatedElements	(	unsigned	i,
		const LinkData &	linkData,
		std::multimap< double, unsigned > &	sortedAssociates,
		reco::PFBlockElement::Type	type = PFBlockElement::NONE,
		LinkTest	test = LINKTEST_RECHIT
	)		const

lock an element ( unlink it from the others ) Colin: this function is misleading void lock(unsigned i, LinkData& linkData ) const; fills a map with the elements associated to element i. elements are sorted by increasing distance. if specified, only the elements of type "type" will be considered if specified, only the link calculated from a certain "test" will be considered: distance test, etc..

Definition at line 75 of file PFBlock.cc.

References cmsHarvester::index, NONE, and test::SmallWORMDict::test().

Referenced by PFAlgo::associatePSClusters(), PFPhotonAlgo::EvaluateSingleLegMVA(), PFEGammaAlgo::EvaluateSingleLegMVA(), PFEGammaFilters::isElectronSafeForJetMET(), PFAlgoTestBenchElectrons::processBlock(), PFElectronAlgo::SetActive(), PFElectronAlgo::SetCandidates(), PFElectronAlgo::SetIDOutputs(), and PFElectronAlgo::SetLinks().

                                        {

  sortedAssociates.clear();
  
  // i is too large
  if( i > elements_.size() ) return;
  // assert(i>=0); i >= 0, since i is unsigned
  
  for(unsigned ie=0; ie<elements_.size(); ie++) {
    
    // considered element itself
    if( ie == i ) {
      continue;
    }
    // not the right type
    if(type !=  PFBlockElement::NONE && 
       elements_[ie].type() != type ) {
      continue;
    }

    // Order the elements by increasing distance !

    unsigned index = 0;
    if( !matrix2vector(i, ie, index) ) continue;

    double c2=-1;
    LinkData::const_iterator it =  linkData.find(index);
    if ( it!=linkData.end() && 
     ( ( (1 << test ) & it->second.test) !=0 || (test == LINKTEST_ALL) ) ) 
      c2= it->second.distance;

    // not associated
    if( c2 < 0 ) { 
      continue;
    }

    sortedAssociates.insert( pair<double,unsigned>(c2, ie) );
  }
} 

void PFBlock::bookLinkData

(

)

Definition at line 20 of file PFBlock.cc.

                           {

}

double reco::PFBlock::dist ( unsigned  ie1, unsigned  ie2, const LinkData &  linkData, LinkTest  test  ) const

inline

Returns

distance of link

Definition at line 94 of file PFBlock.h.

Referenced by PFElectronAlgo::SetLinks().

                               {
      return dist(ie1,ie2,linkData);
    }

double PFBlock::dist	(	unsigned	ie1,
		unsigned	ie2,
		const LinkData &	linkData
	)		const

Returns

distance of link

Definition at line 146 of file PFBlock.cc.

References cmsHarvester::index.

                                                      {

  double Dist = -1;

  unsigned index = 0;
  if( !matrix2vector(ie1, ie2, index) ) return Dist;
  LinkData::const_iterator it =  linkData.find(index);
  if( it!=linkData.end() ) Dist= it->second.distance;

  return Dist;

}

const edm::OwnVector< reco::PFBlockElement >& reco::PFBlock::elements ( ) const

inline

Returns

elements

Definition at line 107 of file PFBlock.h.

References elements_.

Referenced by PFRecoTauAlgorithm::buildPFTau(), FWPFBlockProxyBuilder::buildViewType(), PFPhotonAlgo::EvaluateSingleLegMVA(), PFEGammaAlgo::EvaluateSingleLegMVA(), CalibratableTest::extractCandidate(), PFEGammaFilters::isElectronSafeForJetMET(), PFEGammaFilters::isPhotonSafeForJetMET(), reco::tau::RecoTauElectronRejectionPlugin::operator()(), PFAlgoTestBenchElectrons::processBlock(), PFAlgo::reconstructParticles(), PFElectronAlgo::SetActive(), PFElectronAlgo::SetCandidates(), PFElectronAlgo::SetIDOutputs(), and PFElectronAlgo::SetLinks().

                                                               {
      return elements_;
    }

const LinkData& reco::PFBlock::linkData ( ) const

inline

Returns

link data

Definition at line 112 of file PFBlock.h.

References linkData_.

Referenced by PFPhotonAlgo::EvaluateSingleLegMVA(), PFEGammaAlgo::EvaluateSingleLegMVA(), PFEGammaFilters::isElectronSafeForJetMET(), PFAlgoTestBenchElectrons::processBlock(), PFElectronAlgo::SetActive(), PFElectronAlgo::SetCandidates(), PFElectronAlgo::SetIDOutputs(), and PFElectronAlgo::SetLinks().

                                     {
      return linkData_;
    }

LinkData& reco::PFBlock::linkData ( )

inline

Returns

link data

Definition at line 117 of file PFBlock.h.

References linkData_.

                         {
      return linkData_;
    }

unsigned PFBlock::linkDataSize ( ) const

private

Returns

size of linkData_, calculated from the number of elements

Definition at line 349 of file PFBlock.cc.

References gen::n.

                                     {
  unsigned n = elements_.size();
  
  // number of possible undirected links between n elements.
  // reflective links impossible.
 
  return n*(n-1)/2; 
}

bool PFBlock::matrix2vector	(	unsigned	i,
		unsigned	j,
		unsigned &	index
	)		const

makes the correspondance between a 2d element matrix and the 1D vector which is the most compact way to store the matrix

Definition at line 119 of file PFBlock.cc.

References combine::missing, findQualityFiles::size, and std::swap().

                                                     {

  unsigned size = elements_.size();
  if( iindex == jindex || 
      iindex >=  size ||
      jindex >=  size ) {
    return false;
  }
  
  if( iindex > jindex ) 
    std::swap( iindex, jindex);

  
  index = jindex-iindex-1;

  if(iindex>0) {
    index += iindex*size;
    unsigned missing = iindex*(iindex+1)/2;
    index -= missing;
  }
  
  return true;
}

void PFBlock::setLink	(	unsigned	i1,
		unsigned	i2,
		double	dist,
		LinkData &	linkData,
		LinkTest	test = LINKTEST_RECHIT
	)		const

set a link between elements of indices i1 and i2, of "distance" dist the link is set in the linkData vector provided as an argument. As indicated by the 'const' statement, 'this' is not modified.

Definition at line 26 of file PFBlock.cc.

References assert(), reco::PFBlock::Link::distance, cmsHarvester::index, cmsLHEtoEOSManager::l, convertSQLiteXML::ok, test::SmallWORMDict::test(), and reco::PFBlock::Link::test.

                                   {
  
  assert( test<LINKTEST_ALL );
  
  unsigned index = 0;
  bool ok =  matrix2vector(i1,i2, index);

  if(ok) {
    //ignore the  -1, -1 pair
    if ( Dist > -0.5 ) {
      Link & l = linkData[index];
      l.distance=Dist;
      l.test |= (1 << test);
    }     else  //delete if existing
      {
    LinkData::iterator it = linkData.find(index);
    if(it!=linkData.end()) linkData.erase(it);
      }

  } else {
    assert(0);
  }
  
}

Friends And Related Function Documentation

std::ostream& operator<< ( std::ostream &  out, const PFBlock &  co  )

friend

Member Data Documentation

edm::OwnVector< reco::PFBlockElement > reco::PFBlock::elements_

private

all elements

Definition at line 132 of file PFBlock.h.

Referenced by elements().

LinkData reco::PFBlock::linkData_

private

link data (permanent)

Definition at line 129 of file PFBlock.h.

Referenced by linkData().