Particle Flow Algorithm. More...

#include <PFBlockAlgo.h>

Public Types
typedef reco::PFBlockCollection::const_iterator	IBC

typedef std::list < reco::PFBlockElement * > ::iterator	IE
	define these in *Fwd files in DataFormats/ParticleFlowReco? More...

typedef std::list < reco::PFBlockElement * > ::const_iterator	IEC

typedef std::vector< bool >	Mask

Public Member Functions
const std::auto_ptr < reco::PFBlockCollection > &	blocks () const

void	findBlocks ()
	build blocks More...

	PFBlockAlgo ()

void	setDebug (bool debug)
	sets debug printout flag More...

void	setHOTag (bool ho)

template<template< typename > class T>
void	setInput (const T< reco::PFRecTrackCollection > &trackh, const T< reco::GsfPFRecTrackCollection > &gsftrackh, const T< reco::GsfPFRecTrackCollection > &convbremgsftrackh, const T< reco::MuonCollection > &muonh, const T< reco::PFDisplacedTrackerVertexCollection > &nuclearh, const T< reco::PFRecTrackCollection > &nucleartrackh, const T< reco::PFConversionCollection > &conv, const T< reco::PFV0Collection > &v0, const T< reco::PFClusterCollection > &ecalh, const T< reco::PFClusterCollection > &hcalh, const T< reco::PFClusterCollection > &hoh, const T< reco::PFClusterCollection > &hfemh, const T< reco::PFClusterCollection > &hfhadh, const T< reco::PFClusterCollection > &psh, const T< reco::PhotonCollection > &egphh, const T< reco::SuperClusterCollection > &sceb, const T< reco::SuperClusterCollection > &scee, const T< edm::ValueMap< reco::CaloClusterPtr > > &pfclusterassoc, const Mask &trackMask=dummyMask_, const Mask &gsftrackMask=dummyMask_, const Mask &ecalMask=dummyMask_, const Mask &hcalMask=dummyMask_, const Mask &hoMask=dummyMask_, const Mask &hfemMask=dummyMask_, const Mask &hfhadMask=dummyMask_, const Mask &psMask=dummyMask_, const Mask &phMask=dummyMask_, const Mask &scMask=dummyMask_)
	set input collections of tracks and clusters More...

template<template< typename > class T>
void	setInput (const T< reco::PFRecTrackCollection > &trackh, const T< reco::MuonCollection > &muonh, const T< reco::PFClusterCollection > &ecalh, const T< reco::PFClusterCollection > &hcalh, const T< reco::PFClusterCollection > &hoh, const T< reco::PFClusterCollection > &hfemh, const T< reco::PFClusterCollection > &hfhadh, const T< reco::PFClusterCollection > &psh, const Mask &trackMask=dummyMask_, const Mask &ecalMask=dummyMask_, const Mask &hcalMask=dummyMask_, const Mask &hoMask=dummyMask_, const Mask &psMask=dummyMask_)
	COLIN: I think this is for particle flow at HLT... More...

template<template< typename > class T>
void	setInput (const T< reco::PFRecTrackCollection > &trackh, const T< reco::GsfPFRecTrackCollection > &gsftrackh, const T< reco::PFClusterCollection > &ecalh, const T< reco::PFClusterCollection > &hcalh, const T< reco::PFClusterCollection > &hoh, const T< reco::PFClusterCollection > &psh, const Mask &trackMask=dummyMask_, const Mask &gsftrackMask=dummyMask_, const Mask &ecalMask=dummyMask_, const Mask &hcalMask=dummyMask_, const Mask &hoMask=dummyMask_, const Mask &psMask=dummyMask_)
	COLIN: what is this setinput function for? can it be removed? More...

void	setParameters (std::vector< double > &DPtovPtCut, std::vector< unsigned > &NHitCut, bool useConvBremPFRecTracks, bool useIterTracking, int nuclearInteractionsPurity, bool useEGPhotons, std::vector< double > &photonSelectionCuts, bool useSuperClusters, bool superClusterMatchByRef)

void	setUseOptimization (bool useKDTreeTrackEcalLinker)

std::auto_ptr < reco::PFBlockCollection >	transferBlocks ()

	~PFBlockAlgo ()

Private Member Functions
IE	associate (IE next, IE last, std::vector< PFBlockLink > &links)

void	buildGraph ()

void	checkDisplacedVertexLinks (reco::PFBlock &block) const
	remove extra links between primary track and clusters More...

void	checkMaskSize (const reco::PFRecTrackCollection &tracks, const reco::GsfPFRecTrackCollection &gsftracks, const reco::PFClusterCollection &ecals, const reco::PFClusterCollection &hcals, const reco::PFClusterCollection &hos, const reco::PFClusterCollection &hfems, const reco::PFClusterCollection &hfhads, const reco::PFClusterCollection &pss, const reco::PhotonCollection &egphh, const reco::SuperClusterCollection &sceb, const reco::SuperClusterCollection &scee, const Mask &trackMask, const Mask &gsftrackMask, const Mask &ecalMask, const Mask &hcalMask, const Mask &hoMask, const Mask &hfemMask, const Mask &hfhadMask, const Mask &psMask, const Mask &phMask, const Mask &scMask) const

template<template< typename > class T>
void	fillFromPhoton (const T< reco::PhotonCollection > &, unsigned isc, reco::PFBlockElementSuperCluster *pfbe)

bool	goodPtResolution (const reco::TrackRef &trackref)
	open a resolution map More...

void	link (const reco::PFBlockElement el1, const reco::PFBlockElement el2, PFBlockLink::Type &linktype, reco::PFBlock::LinkTest &linktest, double &dist) const
	check whether 2 elements are linked. Returns distance and linktype More...

bool	linkPrefilter (const reco::PFBlockElement last, const reco::PFBlockElement next) const
	Avoid to check links when not useful. More...

int	muAssocToTrack (const reco::TrackRef &trackref, const edm::Handle< reco::MuonCollection > &muonh) const

int	muAssocToTrack (const reco::TrackRef &trackref, const edm::OrphanHandle< reco::MuonCollection > &muonh) const

void	packLinks (reco::PFBlock &block, const std::vector< PFBlockLink > &links) const

double	testECALAndHCAL (const reco::PFCluster &ecal, const reco::PFCluster &hcal) const

double	testHCALAndHO (const reco::PFCluster &hcal, const reco::PFCluster &ho) const

double	testLinkBySuperCluster (const reco::PFClusterRef &elt1, const reco::PFClusterRef &elt2) const
	test association by Supercluster between two ECAL More...

double	testLinkByVertex (const reco::PFBlockElement elt1, const reco::PFBlockElement elt2) const

double	testPS1AndPS2 (const reco::PFCluster &ps1, const reco::PFCluster &ps2) const

double	testSuperClusterPFCluster (const reco::SuperClusterRef &sct1, const reco::PFClusterRef &elt2) const
	test association between SuperClusters and ECAL More...

double	testTrackAndPS (const reco::PFRecTrack &track, const reco::PFCluster &ps) const

Private Attributes
std::auto_ptr < reco::PFBlockCollection >	blocks_

bool	bNoSuperclus_

bool	debug_
	if true, debug printouts activated More...

std::vector< double >	DPtovPtCut_
	DPt/Pt cut for creating atrack element. More...

std::list< reco::PFBlockElement * >	elements_
	actually, particles will be created by a separate producer More...

std::vector< unsigned >	NHitCut_
	Number of layers crossed cut for creating atrack element. More...

int	nuclearInteractionsPurity_

const edm::ValueMap < reco::CaloClusterPtr > *	pfclusterassoc_

std::vector< int >	pfcSCVec_
	SC corresponding to the PF cluster. More...

const PhotonSelectorAlgo *	photonSelector_
	PhotonSelector. More...

KDTreeLinkerPSEcal	PSELinker_

std::vector< std::vector < reco::PFClusterRef > >	scpfcRefs_
	PF clusters corresponding to a given SC. More...

bool	superClusterMatchByRef_

std::vector < reco::SuperClusterRef >	superClusters_
	list of superclusters More...

KDTreeLinkerTrackEcal	TELinker_

KDTreeLinkerTrackHcal	THLinker_

bool	useConvBremPFRecTracks_
	switch on/off Conversions Brem Recovery with KF Tracks More...

bool	useEGPhotons_
	Flag to turn off the import of EG Photons. More...

bool	useHO_

bool	useIterTracking_
	Flag to turn off quality cuts which require iterative tracking (for heavy-ions) More...

bool	useKDTreeTrackEcalLinker_

bool	useSuperClusters_
	Flag to turn off the import of SuperCluster collections. More...

Static Private Attributes
static const Mask	dummyMask_

Friends
std::ostream &	operator<< (std::ostream &, const PFBlockAlgo &)

Detailed Description

Particle Flow Algorithm.

Author: Colin Bernet

Date: January 2006

Definition at line 71 of file PFBlockAlgo.h.

Member Typedef Documentation

typedef reco::PFBlockCollection::const_iterator PFBlockAlgo::IBC

Definition at line 203 of file PFBlockAlgo.h.

typedef std::list< reco::PFBlockElement* >::iterator PFBlockAlgo::IE

define these in *Fwd files in DataFormats/ParticleFlowReco?

Definition at line 201 of file PFBlockAlgo.h.

typedef std::list< reco::PFBlockElement* >::const_iterator PFBlockAlgo::IEC

Definition at line 202 of file PFBlockAlgo.h.

typedef std::vector<bool> PFBlockAlgo::Mask

Definition at line 95 of file PFBlockAlgo.h.

Constructor & Destructor Documentation

PFBlockAlgo::PFBlockAlgo ( )

Definition at line 21 of file PFBlockAlgo.cc.

                          : 
   blocks_( new reco::PFBlockCollection ),
   DPtovPtCut_(std::vector<double>(4,static_cast<double>(999.))),
   NHitCut_(std::vector<unsigned int>(4,static_cast<unsigned>(0))), 
   useIterTracking_(true),
   photonSelector_(0),
   debug_(false) {}

PFBlockAlgo::~PFBlockAlgo ( )

Definition at line 73 of file PFBlockAlgo.cc.

References gather_cfg::cout, debug_, elements_, and photonSelector_.

                           {
 
 #ifdef PFLOW_DEBUG
   if(debug_)
     cout<<"~PFBlockAlgo - number of remaining elements: "
     <<elements_.size()<<endl;
 #endif
 
   if(photonSelector_) delete photonSelector_;
 
 }

Member Function Documentation

PFBlockAlgo::IE PFBlockAlgo::associate	(	IE	next,
		IE	last,
		std::vector< PFBlockLink > &	links
	)

private

recursive procedure which adds elements from elements_ to the current block, ie blocks_->back(). the resulting links between elements are stored in links, not in the block. afterwards, packLinks( reco::PFBlock& block, const vector<PFBlockLink>& links) has to be called in order to pack the link information in the block.

Definition at line 128 of file PFBlockAlgo.cc.

References blocks_, gather_cfg::cout, debug_, elements_, link(), linkPrefilter(), GetRecoTauVFromDQM_MC_cff::next, and PFBlockLink::NONE.

Referenced by findBlocks().

                                          {
 
     
 #ifdef PFLOW_DEBUG
   if(debug_ ) cout<<"PFBlockAlgo::associate start ----"<<endl;
 #endif
 
   if( last!= elements_.end() ) {
     PFBlockLink::Type linktype = PFBlockLink::NONE;
     double dist = -1;
     PFBlock::LinkTest linktest = PFBlock::LINKTEST_RECHIT;
     link( *last, *next, linktype, linktest, dist ); 
 
     
     if(dist<-0.5) {
 #ifdef PFLOW_DEBUG
       if(debug_ ) cout<<"link failed"<<endl;
 #endif
       return ++next; // association failed
     }
     else {
       // add next element to the current pflowblock
       blocks_->back().addElement( *next );  
 
       // (*next)->setIndex( blocks_->back()->indexToLastElement() );
       
       // this is not necessary? 
       // next->setPFBlock(this);
       
       // create a link between next and last      
       links.push_back( PFBlockLink(linktype, 
                    linktest,
                    dist,
                    (*last)->index(), 
                    (*next)->index() ) );
       // not necessary ?
       //       next->connect( links_.size()-1 );
       //       last->connect( links_.size()-1 );      
     }
   }
   else {
     // add next element to this eflowblock
 #ifdef PFLOW_DEBUG
     if(debug_ ) cout<<"adding to block element "<<(**next)<<endl;
 #endif
     blocks_->back().addElement( *next );
     // (*next)->setIndex( blocks_->back()->indexToLastElement() );   
     // next->setPFBlock(this);
   }
 
   // recursive call: associate next and other unused elements 
   
   //   IE afterNext = next;
   //   ++afterNext;
   //  cout<<"last "<<**last<<" next "<<**next<<endl;
   
   for(IE ie = elements_.begin(); 
       ie != elements_.end();) {
     
     if( ie == last || ie == next ) {
       ++ie;
       continue;
     } 
     bool bTestLink = linkPrefilter(*next, *ie);
 
     if(!bTestLink) {
       ++ie;
       continue;
     }
     
     // *ie already included to a block
     if( (*ie)->locked() ) {
 #ifdef PFLOW_DEBUG
       if(debug_ ) cout<<"element "<<(**ie)<<"already used"<<endl;
 #endif
       ++ie;
       continue;
     }    
     
     
 #ifdef PFLOW_DEBUG
     if(debug_ ) cout<<"calling associate "<<(**next)<<" & "<<(**ie)<<endl;
 #endif
     ie = associate(next, ie, links);
   }       
 
 #ifdef PFLOW_DEBUG
   if(debug_ ) {
     cout<<"**** deleting element "<<endl;
     cout<<**next<<endl;
   }
 #endif
   delete *next;
 
 #ifdef PFLOW_DEBUG
   if(debug_ ) {
     cout<<"**** removing element "<<endl;
   }
 #endif
 
   IE iteratorToNextFreeElement = elements_.erase( next );
 
 #ifdef PFLOW_DEBUG
   if(debug_ ) cout<<"PFBlockAlgo::associate stop ----"<<endl;
 #endif
 
   return iteratorToNextFreeElement;
 }

const std::auto_ptr< reco::PFBlockCollection >& PFBlockAlgo::blocks ( ) const

inline

Returns: collection of blocks

Definition at line 194 of file PFBlockAlgo.h.

References blocks_.

Referenced by operator<<().

195 {return blocks_;}

PFBlockAlgo::blocks_

std::auto_ptr< reco::PFBlockCollection > blocks_

Definition: PFBlockAlgo.h:311

void PFBlockAlgo::buildGraph ( )

private

COLIN: not used. Could be removed. Could also be implemented, to produce a graph of a block, Showing the connections between the elements

Definition at line 316 of file PFBlockAlgo.cc.

                         {
   // loop on all blocks and create a big graph
 }

void PFBlockAlgo::checkDisplacedVertexLinks ( reco::PFBlock & block ) const

private

remove extra links between primary track and clusters

void PFBlockAlgo::checkMaskSize	(	const reco::PFRecTrackCollection &	tracks,
		const reco::GsfPFRecTrackCollection &	gsftracks,
		const reco::PFClusterCollection &	ecals,
		const reco::PFClusterCollection &	hcals,
		const reco::PFClusterCollection &	hos,
		const reco::PFClusterCollection &	hfems,
		const reco::PFClusterCollection &	hfhads,
		const reco::PFClusterCollection &	pss,
		const reco::PhotonCollection &	egphh,
		const reco::SuperClusterCollection &	sceb,
		const reco::SuperClusterCollection &	scee,
		const Mask &	trackMask,
		const Mask &	gsftrackMask,
		const Mask &	ecalMask,
		const Mask &	hcalMask,
		const Mask &	hoMask,
		const Mask &	hfemMask,
		const Mask &	hfhadMask,
		const Mask &	psMask,
		const Mask &	phMask,
		const Mask &	scMask
	)		const

private

checks size of the masks with respect to the vectors they refer to. throws std::length_error if one of the masks has the wrong size

Definition at line 1179 of file PFBlockAlgo.cc.

Referenced by setInput().

                                           {
 
   if( !trackMask.empty() && 
       trackMask.size() != tracks.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the track mask is different ";
     err += "from the size of the track vector.";
     throw std::length_error( err.c_str() );
   }
 
   if( !gsftrackMask.empty() && 
       gsftrackMask.size() != gsftracks.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the gsf track mask is different ";
     err += "from the size of the gsftrack vector.";
     throw std::length_error( err.c_str() );
   }
 
   if( !ecalMask.empty() && 
       ecalMask.size() != ecals.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the ecal mask is different ";
     err += "from the size of the ecal clusters vector.";
     throw std::length_error( err.c_str() );
   }
   
   if( !hcalMask.empty() && 
       hcalMask.size() != hcals.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the hcal mask is different ";
     err += "from the size of the hcal clusters vector.";
     throw std::length_error( err.c_str() );
   }
 
   if( !hoMask.empty() && 
       hoMask.size() != hos.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the ho mask is different ";
     err += "from the size of the ho clusters vector.";
     throw std::length_error( err.c_str() );
   }
 
 
   if( !hfemMask.empty() && 
       hfemMask.size() != hfems.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the hfem mask is different ";
     err += "from the size of the hfem clusters vector.";
     throw std::length_error( err.c_str() );
   }
 
   if( !hfhadMask.empty() && 
       hfhadMask.size() != hfhads.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the hfhad mask is different ";
     err += "from the size of the hfhad clusters vector.";
     throw std::length_error( err.c_str() );
   }
 
   if( !psMask.empty() && 
       psMask.size() != pss.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the ps mask is different ";
     err += "from the size of the ps clusters vector.";
     throw std::length_error( err.c_str() );
   }
   
   if( !phMask.empty() && 
       phMask.size() != egphh.size() ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the photon mask is different ";
     err += "from the size of the photon vector.";
     throw std::length_error( err.c_str() );
   }
   
   if( !scMask.empty() && 
       scMask.size() != (sceb.size() + scee.size()) ) {
     string err = "PFBlockAlgo::setInput: ";
     err += "The size of the SC mask is different ";
     err += "from the size of the SC vectors.";
     throw std::length_error( err.c_str() );
   }  
 
 }

template<template< typename > class T>

void PFBlockAlgo::fillFromPhoton	(	const T< reco::PhotonCollection > &	egh,
		unsigned	isc,
		reco::PFBlockElementSuperCluster *	pfbe
	)

private

Definition at line 1228 of file PFBlockAlgo.h.

References reco::PFBlockElementSuperCluster::setEcalIso(), reco::PFBlockElementSuperCluster::setFromPhoton(), reco::PFBlockElementSuperCluster::setHcalIso(), reco::PFBlockElementSuperCluster::setHoE(), reco::PFBlockElementSuperCluster::setPhotonRef(), and reco::PFBlockElementSuperCluster::setTrackIso().

Referenced by setInput().

                                                                                                                            {
   reco::PhotonRef photonRef(egh,isc);
     pfbe->setTrackIso(photonRef->trkSumPtHollowConeDR04());
     pfbe->setEcalIso(photonRef->ecalRecHitSumEtConeDR04());
     pfbe->setHcalIso(photonRef->hcalTowerSumEtConeDR04());
     pfbe->setHoE(photonRef->hadronicOverEm());
     pfbe->setPhotonRef(photonRef);
     pfbe->setFromPhoton(true);
   }

void PFBlockAlgo::findBlocks ( )

build blocks

Definition at line 86 of file PFBlockAlgo.cc.

References associate(), blocks_, gather_cfg::cout, debug_, elements_, electronStore::links, packLinks(), KDTreeLinkerBase::process(), PSELinker_, TELinker_, THLinker_, and useKDTreeTrackEcalLinker_.

Referenced by PFRootEventManager::particleFlow().

                         {
 
   // Glowinski & Gouzevitch
   if (useKDTreeTrackEcalLinker_) {
     TELinker_.process();
     THLinker_.process();
     PSELinker_.process();
   }
   // !Glowinski & Gouzevitch
 
   // the blocks have not been passed to the event, and need to be cleared
   if(blocks_.get() )blocks_->clear();
   else 
     blocks_.reset( new reco::PFBlockCollection );
 
   blocks_->reserve(elements_.size());
   for(IE ie = elements_.begin(); 
       ie != elements_.end();) {
     
 #ifdef PFLOW_DEBUG
     if(debug_) {
       cout<<" PFBlockAlgo::findBlocks() ----------------------"<<endl;
       cout<<" element "<<**ie<<endl;
       cout<<" creating new block"<<endl;
     }
 #endif
 
     blocks_->push_back( PFBlock() );
     
     vector< PFBlockLink > links;
 
     //    list< IE > used;
     ie = associate( elements_.end() , ie, links );
 
     packLinks( blocks_->back(), links );
   }
 }

bool PFBlockAlgo::goodPtResolution ( const reco::TrackRef & trackref )

private

open a resolution map

check the Pt resolution

Definition at line 1321 of file PFBlockAlgo.cc.

References gather_cfg::cout, debug_, DPtovPtCut_, NHitCut_, P, reco::tau::disc::Pt(), mathSSE::sqrt(), and useIterTracking_.

Referenced by setInput().

                                                            {
 
   double P = trackref->p();
   double Pt = trackref->pt();
   double DPt = trackref->ptError();
   unsigned int NHit = trackref->hitPattern().trackerLayersWithMeasurement();
   unsigned int NLostHit = trackref->hitPattern().trackerLayersWithoutMeasurement();
   unsigned int LostHits = trackref->numberOfLostHits();
   double sigmaHad = sqrt(1.20*1.20/P+0.06*0.06) / (1.+LostHits);
 
   // iteration 1,2,3,4,5 correspond to algo = 1/4,5,6,7,8,9
   unsigned int Algo = 0; 
   switch (trackref->algo()) {
   case TrackBase::ctf:
   case TrackBase::iter0:
   case TrackBase::iter1:
   case TrackBase::iter2:
     Algo = 0;
     break;
   case TrackBase::iter3:
     Algo = 1;
     break;
   case TrackBase::iter4:
     Algo = 2;
     break;
   case TrackBase::iter5:
     Algo = 3;
     break;
   case TrackBase::iter6:
     Algo = 4;
     break;
   default:
     Algo = useIterTracking_ ? 5 : 0;
     break;
   }
 
   // Protection against 0 momentum tracks
   if ( P < 0.05 ) return false;
 
   // Temporary : Reject all tracking iteration beyond 5th step. 
   if ( Algo > 4 ) return false;
  
   if (debug_) cout << " PFBlockAlgo: PFrecTrack->Track Pt= "
            << Pt << " DPt = " << DPt << endl;
   if ( ( DPtovPtCut_[Algo] > 0. && 
      DPt/Pt > DPtovPtCut_[Algo]*sigmaHad ) || 
        NHit < NHitCut_[Algo] ) { 
     // (Algo >= 3 && LostHits != 0) ) {
     if (debug_) cout << " PFBlockAlgo: skip badly measured track"
              << ", P = " << P 
              << ", Pt = " << Pt 
              << " DPt = " << DPt 
              << ", N(hits) = " << NHit << " (Lost : " << LostHits << "/" << NLostHit << ")"
              << ", Algo = " << Algo
              << endl;
     if (debug_) cout << " cut is DPt/Pt < " << DPtovPtCut_[Algo] * sigmaHad << endl;
     if (debug_) cout << " cut is NHit >= " << NHitCut_[Algo] << endl;
     /*
     std::cout << "Track REJECTED : ";
     std::cout << ", P = " << P 
           << ", Pt = " << Pt 
           << " DPt = " << DPt 
           << ", N(hits) = " << NHit << " (Lost : " << LostHits << "/" << NLostHit << ")"
           << ", Algo = " << Algo
           << std::endl;
     */
     return false;
   }
 
   /*
   std::cout << "Track Accepted : ";
   std::cout << ", P = " << P 
        << ", Pt = " << Pt 
        << " DPt = " << DPt 
        << ", N(hits) = " << NHit << " (Lost : " << LostHits << "/" << NLostHit << ")"
        << ", Algo = " << Algo
        << std::endl;
   */
   return true;
 }

void PFBlockAlgo::link	(	const reco::PFBlockElement *	el1,
		const reco::PFBlockElement *	el2,
		PFBlockLink::Type &	linktype,
		reco::PFBlock::LinkTest &	linktest,
		double &	dist
	)		const

private

check whether 2 elements are linked. Returns distance and linktype

Definition at line 323 of file PFBlockAlgo.cc.

Referenced by associate(), and packLinks().

                                {
   
   // ACHTUNG!!!! If you introduce new links check that they are not desabled in linkPrefilter!!!!
 
 
   dist=-1.;
   linktest = PFBlock::LINKTEST_RECHIT; //rechit by default 
 
   PFBlockElement::Type type1 = el1->type();
   PFBlockElement::Type type2 = el2->type();
 
   linktype = static_cast<PFBlockLink::Type>
     ((1<< (type1-1) ) | (1<< (type2-1) ));
 
   if(debug_ ) std::cout << " PFBlockAlgo links type1 " << type1 << " type2 " << type2 << std::endl;
 
   const PFBlockElement* lowEl = el1;
   const PFBlockElement* highEl = el2;
   
   if(type1>type2) {
     lowEl = el2;
     highEl = el1;
   }
   
   switch(linktype) {
   // Track and preshower cluster links are not used for now - disable
   case PFBlockLink::PS1andECAL:
   case PFBlockLink::PS2andECAL:
     {
       // if(debug_ ) cout<< "PSandECAL" <<endl;
       PFClusterRef  psref = lowEl->clusterRef();
       PFClusterRef  ecalref = highEl->clusterRef();
       assert( !psref.isNull() );
       assert( !ecalref.isNull() );
 
       // Check if the linking has been done using the KDTree algo
       // Glowinski & Gouzevitch
       if ( useKDTreeTrackEcalLinker_ && lowEl->isMultilinksValide() ) { // KDTree algo
     const reco::PFMultilinksType& multilinks = lowEl->getMultilinks();
     
     double ecalPhi = ecalref->positionREP().Phi();
     double ecalEta = ecalref->positionREP().Eta();
     
     // Check if the link PS/Ecal exist
     reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
     for (; mlit != multilinks.end(); ++mlit)
       if ((mlit->first == ecalPhi) && (mlit->second == ecalEta))
         break;
     
     // If the link exist, we fill dist and linktest. We use old algorithme method.
     if (mlit != multilinks.end()){
       double xPS = psref->position().X();
       double yPS = psref->position().Y();
       double xECAL  = ecalref->position().X();
       double yECAL  = ecalref->position().Y();
 
       dist = LinkByRecHit::computeDist(xECAL/1000.,yECAL/1000.,xPS/1000.  ,yPS/1000, false);
     }
 
       } else { //Old algorithm
     dist = LinkByRecHit::testECALAndPSByRecHit( *ecalref, *psref ,debug_);
       }
 
       //      linktest = PFBlock::LINKTEST_RECHIT;
       
       break;
     }
   case PFBlockLink::TRACKandECAL:
     {
       if(debug_ ) cout<<"TRACKandECAL"<<endl;
 
       PFRecTrackRef trackref = lowEl->trackRefPF();
       PFClusterRef  clusterref = highEl->clusterRef();
       assert( !trackref.isNull() );
       assert( !clusterref.isNull() );
 
       if(debug_ ) std::cout << " Track pt " << trackref->trackRef()->pt() << std::endl;
 
       // Check if the linking has been done using the KDTree algo
       // Glowinski & Gouzevitch
       if ( useKDTreeTrackEcalLinker_ && lowEl->isMultilinksValide() ) { //KDTree Algo
 
     const reco::PFMultilinksType& multilinks = lowEl->getMultilinks();
     double ecalphi = clusterref->positionREP().Phi();
     double ecaleta = clusterref->positionREP().Eta();
 
     // Check if the link Track/Ecal exist
     reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
     for (; mlit != multilinks.end(); ++mlit)
       if ((mlit->first == ecalphi) && (mlit->second == ecaleta))
         break;
 
     
     // If the link exist, we fill dist and linktest. We use old algorithme method.
     if (mlit != multilinks.end()){
 
 
       //Should be something like this :
       //      const reco::PFRecTrack& track = *trackref;
       //instead of this :
       /*
       reco::PFRecTrack track (*trackref);
       const reco::PFTrajectoryPoint& atECAL_tmp = 
         (*trackref).extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax );
       if(std::abs(atECAL_tmp.positionREP().Eta())<1E-9 &&
          std::abs(atECAL_tmp.positionREP().Phi())<1E-9 &&
          atECAL_tmp.positionREP().R()<1E-9) 
         track.calculatePositionREP();
       */
 
       const reco::PFTrajectoryPoint& atECAL = 
         trackref->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax );
       
       double tracketa = atECAL.positionREP().Eta();
       double trackphi = atECAL.positionREP().Phi();
 
       dist = LinkByRecHit::computeDist(ecaleta, ecalphi, tracketa, trackphi);
     }
 
       } else {// Old algorithm
     if ( trackref->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() )
       dist = LinkByRecHit::testTrackAndClusterByRecHit( *trackref, *clusterref, false, debug_ );
     else
       dist = -1.;
       }
 
       if ( debug_ ) { 
     if( dist > 0. ) { 
       std::cout << " Here a link has been established"
             << " between a track an Ecal with dist  " 
             << dist <<  std::endl;
     } else
       std::cout << " No link found " << std::endl;
       }
       
       //     linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::TRACKandHCAL:
     {
       //      if(debug_ ) cout<<"TRACKandHCAL"<<endl;
 
       PFRecTrackRef trackref = lowEl->trackRefPF();
       PFClusterRef  clusterref = highEl->clusterRef();
       assert( !trackref.isNull() );
       assert( !clusterref.isNull() );
 
       // Check if the linking has been done using the KDTree algo
       // Glowinski & Gouzevitch
       if ( useKDTreeTrackEcalLinker_ && highEl->isMultilinksValide() ) { //KDTree Algo
     
     const reco::PFMultilinksType& multilinks = highEl->getMultilinks();
 
     /*
     reco::PFRecTrack track (*trackref);
     const reco::PFTrajectoryPoint& atHCALEntrance_tmp = 
       (*trackref).extrapolatedPoint( reco::PFTrajectoryPoint::HCALEntrance);
     if (std::abs(atHCALEntrance_tmp.positionREP().Eta())<1E-9 &&
         std::abs(atHCALEntrance_tmp.positionREP().Phi())<1E-9 &&
         atHCALEntrance_tmp.positionREP().R()<1E-9)   
       track.calculatePositionREP();
     */
 
     const reco::PFTrajectoryPoint& atHCAL = 
       trackref->extrapolatedPoint(reco::PFTrajectoryPoint::HCALEntrance);
     
       
     double tracketa = atHCAL.positionREP().Eta();
     double trackphi = atHCAL.positionREP().Phi();
 
     // Check if the link Track/Ecal exist
     reco::PFMultilinksType::const_iterator mlit = multilinks.begin();
     for (; mlit != multilinks.end(); ++mlit)
       if ((mlit->first == trackphi) && (mlit->second == tracketa))
         break;
 
     // If the link exist, we fill dist and linktest. We use old algorithme method.
     if (mlit != multilinks.end()){
 
       const reco::PFTrajectoryPoint& atHCALExit = 
         trackref->extrapolatedPoint(reco::PFTrajectoryPoint::HCALExit);
       double dHEta = 0.0;
       double dHPhi = 0.0;
       if (atHCALExit.position().R()>atHCAL.position().R()) {
         dHEta = atHCALExit.positionREP().Eta()-atHCAL.positionREP().Eta();
         dHPhi = atHCALExit.positionREP().Phi()-atHCAL.positionREP().Phi(); 
         if ( dHPhi > M_PI ) dHPhi = dHPhi - 2.*M_PI;
         else if ( dHPhi < -M_PI ) dHPhi = dHPhi + 2.*M_PI; 
       } else {
         std::cout << "Qu'est ce que c'est que ce gag ? " 
               << atHCALExit.position().R() << " is larger than " 
               << atHCAL.position().R() << " !" << std::endl;
       }
 
       tracketa += 0.1 * dHEta;
       trackphi += 0.1 * dHPhi;
 
       double clusterphi = clusterref->positionREP().Phi();
       double clustereta = clusterref->positionREP().Eta();
       
       dist = LinkByRecHit::computeDist(clustereta, clusterphi, tracketa, trackphi);
     }       
 
       } else {// Old algorithm
     if ( trackref->extrapolatedPoint( reco::PFTrajectoryPoint::HCALEntrance ).isValid() )
       dist = LinkByRecHit::testTrackAndClusterByRecHit( *trackref, *clusterref, false, debug_ );
     else
       dist = -1.;
       }
       
       //      linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::TRACKandHO:
     {
       if(debug_ ) cout<<"TRACKandHO"<<endl;
 
       PFRecTrackRef trackref = lowEl->trackRefPF();
       PFClusterRef  clusterref = highEl->clusterRef();
       
       assert( !trackref.isNull() );
       assert( !clusterref.isNull() );
 
       
       // Old algorithm
       //      cout<<"TRACKandHO1"<<trackref->pt()<<" "<<trackref->eta()<<" "<<trackref->phi()<<endl;
       //Same value is used in PFTrackTransformer::addPoints() for HOLayer, but allow for some rounding precision
       if ( lowEl->trackRef()->pt() > 3.00001 && trackref->extrapolatedPoint( reco::PFTrajectoryPoint::HOLayer ).isValid() ) {
     //  cout<<"TRACKandHO2"<<endl;
     dist = LinkByRecHit::testTrackAndClusterByRecHit( *trackref, *clusterref, false, debug_ );
     
     //  cout <<"dist TRACKandHO "<<dist<<endl;
       } else {
     dist = -1.;
       }
       //      linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::ECALandHCAL:
     {
       //       cout<<"ECALandHCAL"<<endl;
       PFClusterRef  ecalref = lowEl->clusterRef();
       PFClusterRef  hcalref = highEl->clusterRef();
       assert( !ecalref.isNull() );
       assert( !hcalref.isNull() );
       // PJ - 14-May-09 : A link by rechit is needed here !
       dist = testECALAndHCAL( *ecalref, *hcalref );
       // dist = -1.;
       //     linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::HCALandHO:
     {
       //       cout<<"HCALandH0"<<endl;
       PFClusterRef  hcalref = lowEl->clusterRef();
       PFClusterRef  horef = highEl->clusterRef();
       assert( !hcalref.isNull() );
       assert( !horef.isNull() );
       // PJ - 14-May-09 : A link by rechit is needed here !
       dist = testHCALAndHO( *hcalref, *horef );
       // dist = -1.;
       //      cout <<"Dist "<<dist<<endl;
       //     linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::HFEMandHFHAD:
     {
       // cout<<"HFEMandHFHAD"<<endl;
       PFClusterRef  eref = lowEl->clusterRef();
       PFClusterRef  href = highEl->clusterRef();
       assert( !eref.isNull() );
       assert( !href.isNull() );
       dist = LinkByRecHit::testHFEMAndHFHADByRecHit( *eref, *href, debug_ );
       //    linktest = PFBlock::LINKTEST_RECHIT;
       break;      
     }
 
   case PFBlockLink::TRACKandTRACK:
     {
       if(debug_ ) 
     cout<<"TRACKandTRACK"<<endl;
       dist = testLinkByVertex(lowEl, highEl);
       if(debug_ ) 
     std::cout << " PFBlockLink::TRACKandTRACK dist " << dist << std::endl;
       //   linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
 
   case PFBlockLink::ECALandECAL:
       {
     
     PFClusterRef  ecal1ref = lowEl->clusterRef();
     PFClusterRef  ecal2ref = highEl->clusterRef();
     assert( !ecal1ref.isNull() );
     assert( !ecal2ref.isNull() );
     if(debug_)
       cout << " PFBlockLink::ECALandECAL" << endl;
     dist = testLinkBySuperCluster(ecal1ref,ecal2ref);
     break;
       }
 
   case PFBlockLink::ECALandGSF:
     {
       PFClusterRef  clusterref = lowEl->clusterRef();
       assert( !clusterref.isNull() );
       const reco::PFBlockElementGsfTrack *  GsfEl =  dynamic_cast<const reco::PFBlockElementGsfTrack*>(highEl);
       const PFRecTrack * myTrack =  &(GsfEl->GsftrackPF());
       if ( myTrack->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() )
     dist = LinkByRecHit::testTrackAndClusterByRecHit( *myTrack, *clusterref, false, debug_ );
       else
     dist = -1.;
       //   linktest = PFBlock::LINKTEST_RECHIT;
       
       if ( debug_ ) {
     if ( dist > 0. ) {
       std::cout << " Here a link has been established" 
             << " between a GSF track an Ecal with dist  " 
             << dist <<  std::endl;
     } else {
       if(debug_ ) std::cout << " No link found " << std::endl;
     }
       }
       break;
     }
   case PFBlockLink::TRACKandGSF:
     {
       PFRecTrackRef trackref = lowEl->trackRefPF();
       assert( !trackref.isNull() );
       const reco::PFBlockElementGsfTrack *  GsfEl =  
     dynamic_cast<const reco::PFBlockElementGsfTrack*>(highEl);
       GsfPFRecTrackRef gsfref = GsfEl->GsftrackRefPF();
       reco::TrackRef kftrackref= (*trackref).trackRef();
       assert( !gsfref.isNull() );
       PFRecTrackRef refkf = (*gsfref).kfPFRecTrackRef();
       if(refkf.isNonnull()) {
     reco::TrackRef gsftrackref = (*refkf).trackRef();
     if (gsftrackref.isNonnull()&&kftrackref.isNonnull()) {
       if (kftrackref == gsftrackref) { 
         dist = 0.001;
       } else { 
         dist = -1.;
       }
     } else { 
       dist = -1.;
     }
       } else {
     dist = -1.;
       }
       
       
       if(useConvBremPFRecTracks_) {
     if(lowEl->isLinkedToDisplacedVertex()){
       vector<PFRecTrackRef> pfrectrack_vec = GsfEl->GsftrackRefPF()->convBremPFRecTrackRef();
       if(pfrectrack_vec.size() > 0){
         for(unsigned int iconv = 0; iconv <  pfrectrack_vec.size(); ++iconv) {
           if( lowEl->trackType(reco::PFBlockElement::T_FROM_GAMMACONV)) {
         // use track ref
         if(kftrackref == (*pfrectrack_vec[iconv]).trackRef()) {     
           dist = 0.001;
         }
           } 
           else{
         // use the track base ref
         reco::TrackBaseRef newTrackBaseRef((*pfrectrack_vec[iconv]).trackRef());
         reco::TrackBaseRef elemTrackBaseRef(kftrackref);          
         if(newTrackBaseRef == elemTrackBaseRef){
           dist = 0.001;
         } 
           }
         }
       }
     }
       }
  
 
       break;      
     }
      
   case PFBlockLink::GSFandBREM:
     {
       const reco::PFBlockElementGsfTrack * GsfEl  =  dynamic_cast<const reco::PFBlockElementGsfTrack*>(lowEl);
       const reco::PFBlockElementBrem * BremEl =  dynamic_cast<const reco::PFBlockElementBrem*>(highEl);
       GsfPFRecTrackRef gsfref = GsfEl->GsftrackRefPF();
       GsfPFRecTrackRef bremref = BremEl->GsftrackRefPF();
       assert( !gsfref.isNull() );
       assert( !bremref.isNull() );
       if (gsfref == bremref)  { 
     dist = 0.001;
       } else { 
     dist = -1.;
       }
       break;
     }
   case PFBlockLink::GSFandGSF:
     {
       const reco::PFBlockElementGsfTrack * lowGsfEl  =  
     dynamic_cast<const reco::PFBlockElementGsfTrack*>(lowEl);
       const reco::PFBlockElementGsfTrack * highGsfEl  =  
     dynamic_cast<const reco::PFBlockElementGsfTrack*>(highEl);
       
       GsfPFRecTrackRef lowgsfref = lowGsfEl->GsftrackRefPF();
       GsfPFRecTrackRef highgsfref = highGsfEl->GsftrackRefPF();
       assert( !lowgsfref.isNull() );
       assert( !highgsfref.isNull() );
       
       if( (lowGsfEl->trackType(reco::PFBlockElement::T_FROM_GAMMACONV) == false && 
        highGsfEl->trackType(reco::PFBlockElement::T_FROM_GAMMACONV)) ||
       (highGsfEl->trackType(reco::PFBlockElement::T_FROM_GAMMACONV) == false && 
        lowGsfEl->trackType(reco::PFBlockElement::T_FROM_GAMMACONV))) {
     if(lowgsfref->trackId() == highgsfref->trackId()) {
       dist = 0.001;
     }
     else {
       dist = -1.;
     }
       }
       break;
     }
   case PFBlockLink::ECALandBREM:
     {
       PFClusterRef  clusterref = lowEl->clusterRef();
       assert( !clusterref.isNull() );
       const reco::PFBlockElementBrem * BremEl =  
     dynamic_cast<const reco::PFBlockElementBrem*>(highEl);
       const PFRecTrack * myTrack = &(BremEl->trackPF());
       /*
       double DP = (BremEl->DeltaP())*(-1.);
       double SigmaDP = BremEl->SigmaDeltaP();
       double SignBremDp = DP/SigmaDP;
       */
       bool isBrem = true;
       if ( myTrack->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() )
     dist = LinkByRecHit::testTrackAndClusterByRecHit( *myTrack, *clusterref, isBrem, debug_);
       else
     dist = -1.;
       if( debug_ && dist > 0. ) 
     std::cout << "ECALandBREM: dist testTrackAndClusterByRecHit " 
           << dist << std::endl;
       //   linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::HCALandGSF:
     {
       PFClusterRef  clusterref = lowEl->clusterRef();
       assert( !clusterref.isNull() );
       const reco::PFBlockElementGsfTrack *  GsfEl =  dynamic_cast<const reco::PFBlockElementGsfTrack*>(highEl);
       const PFRecTrack * myTrack =  &(GsfEl->GsftrackPF());
       if ( myTrack->extrapolatedPoint( reco::PFTrajectoryPoint::HCALEntrance ).isValid() )
     dist = LinkByRecHit::testTrackAndClusterByRecHit( *myTrack, *clusterref, false, debug_ );
       else
     dist = -1.;
 
       //    linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::HCALandBREM:
     {
       PFClusterRef  clusterref = lowEl->clusterRef();
       assert( !clusterref.isNull() );
       const reco::PFBlockElementBrem * BremEl =  dynamic_cast<const reco::PFBlockElementBrem*>(highEl);
       const PFRecTrack * myTrack = &(BremEl->trackPF());
       bool isBrem = true;
       if ( myTrack->extrapolatedPoint( reco::PFTrajectoryPoint::HCALEntrance ).isValid() )
     dist = LinkByRecHit::testTrackAndClusterByRecHit( *myTrack, *clusterref, isBrem, debug_);
       else
     dist = -1.;
       break;
     }
   case PFBlockLink::SCandECAL:
     {
       PFClusterRef  clusterref = lowEl->clusterRef();
 
       assert( !clusterref.isNull() );
       
       const reco::PFBlockElementSuperCluster * scEl = 
     dynamic_cast<const reco::PFBlockElementSuperCluster*>(highEl);
       assert (!scEl->superClusterRef().isNull());
       dist = testSuperClusterPFCluster(scEl->superClusterRef(),
                        clusterref);
       break;
     }
     /*
   // Links between the two preshower layers are not used for now - disable
   case PFBlockLink::PS1andPS2:
     {
       PFClusterRef  ps1ref = lowEl->clusterRef();
       PFClusterRef  ps2ref = highEl->clusterRef();
       assert( !ps1ref.isNull() );
       assert( !ps2ref.isNull() );
       // PJ - 14-May-09 : A link by rechit is needed here !
       // dist = testPS1AndPS2( *ps1ref, *ps2ref );
       dist = -1.;
       linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::TRACKandPS1:
   case PFBlockLink::TRACKandPS2:
     {
       //cout<<"TRACKandPS"<<endl;
       PFRecTrackRef trackref = lowEl->trackRefPF();
       PFClusterRef  clusterref = highEl->clusterRef();
       assert( !trackref.isNull() );
       assert( !clusterref.isNull() );
       // PJ - 14-May-09 : A link by rechit is needed here !
       dist = testTrackAndPS( *trackref, *clusterref );
       linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
     // GSF Track/Brem Track and preshower cluster links are not used for now - disable
   case PFBlockLink::PS1andGSF:
   case PFBlockLink::PS2andGSF:
     {
       PFClusterRef  psref = lowEl->clusterRef();
       assert( !psref.isNull() );
       const reco::PFBlockElementGsfTrack *  GsfEl =  dynamic_cast<const reco::PFBlockElementGsfTrack*>(highEl);
       const PFRecTrack * myTrack =  &(GsfEl->GsftrackPF());
       // PJ - 14-May-09 : A link by rechit is needed here !
       dist = testTrackAndPS( *myTrack, *psref );
       linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
   case PFBlockLink::PS1andBREM:
   case PFBlockLink::PS2andBREM:
     {
       PFClusterRef  psref = lowEl->clusterRef();
       assert( !psref.isNull() );
       const reco::PFBlockElementBrem * BremEl =  dynamic_cast<const reco::PFBlockElementBrem*>(highEl);
       const PFRecTrack * myTrack = &(BremEl->trackPF());
       // PJ - 14-May-09 : A link by rechit is needed here !
       dist = testTrackAndPS( *myTrack, *psref );
       linktest = PFBlock::LINKTEST_RECHIT;
       break;
     }
     */
 
   default:
     dist = -1.;
     //   linktest = PFBlock::LINKTEST_RECHIT;
     // cerr<<"link type not implemented yet: 0x"<<hex<<linktype<<dec<<endl;
     // assert(0);
     return;
   }
 }

bool PFBlockAlgo::linkPrefilter	(	const reco::PFBlockElement *	last,
		const reco::PFBlockElement *	next
	)		const

inlineprivate

Avoid to check links when not useful.

Definition at line 1432 of file PFBlockAlgo.cc.

References bNoSuperclus_, ECAL, reco::PFBlockElement::getMultilinks(), reco::PFBlockElement::isMultilinksValide(), prof2calltree::last, GetRecoTauVFromDQM_MC_cff::next, and useKDTreeTrackEcalLinker_.

Referenced by associate(), and packLinks().

                                                                                                {
 
   PFBlockElement::Type type1 = (last)->type();
   PFBlockElement::Type type2 = (next)->type();
 
   if( type1==type2 ) {
     // cannot link 2 elements of the same type. 
     // except if the elements are 2 tracks or 2 ECAL
     if( type1!=PFBlockElement::TRACK && type1!=PFBlockElement::GSF &&
     type1!=PFBlockElement::ECAL) { // && type1!=PFBlockElement::HCAL) {
       return false;
     }
 
     if (type1==PFBlockElement::ECAL && bNoSuperclus_) return false;
 
     // cannot link two primary tracks  (except if they come from a V0)
     if( type1 ==PFBlockElement::TRACK) {
       if ( !((last)->isLinkedToDisplacedVertex()) || !((next)->isLinkedToDisplacedVertex())) 
       return false;
     }
   }
 
   if ((type1 == PFBlockElement::PS1 || type1 == PFBlockElement::PS2) && (type2 != PFBlockElement::ECAL)) return false;
   if ((type2 == PFBlockElement::PS1 || type2 == PFBlockElement::PS2) && (type1 != PFBlockElement::ECAL)) return false;
   if ((type1 == PFBlockElement::HFEM && type2 != PFBlockElement::HFHAD) || (type1 == PFBlockElement::HFHAD && type2 != PFBlockElement::HFEM)) return false;
 
   if (useKDTreeTrackEcalLinker_){ 
   
     if ( type1 == PFBlockElement::TRACK && type2 == PFBlockElement::ECAL)
       if ( last->isMultilinksValide()  && last->getMultilinks().size()==0 ) return false;
     if ( type2 == PFBlockElement::TRACK && type1 == PFBlockElement::ECAL)
       if ( next->isMultilinksValide() && next->getMultilinks().size()==0 ) return false;
     if ( type1 == PFBlockElement::PS1 || type1 == PFBlockElement::PS2)
       if ( last->isMultilinksValide()  && last->getMultilinks().size()==0 ) return false;
     if ( type2 == PFBlockElement::PS1 || type2 == PFBlockElement::PS2)
       if ( next->isMultilinksValide() && next->getMultilinks().size()==0 ) return false;
 
   }
 
   return true;
 
 }

int PFBlockAlgo::muAssocToTrack	(	const reco::TrackRef &	trackref,
		const edm::Handle< reco::MuonCollection > &	muonh
	)		const

private

Definition at line 1403 of file PFBlockAlgo.cc.

References edm::Ref< C, T, F >::isNonnull(), edm::HandleBase::isValid(), and j.

Referenced by setInput().

                                                                    {
   if(muonh.isValid() ) {
     for(unsigned j=0;j<muonh->size(); ++j) {
       reco::MuonRef muonref( muonh, j );
       if (muonref->track().isNonnull()) 
     if( muonref->track() == trackref ) return j;
     }
   }
   return -1; // not found
 }

int PFBlockAlgo::muAssocToTrack	(	const reco::TrackRef &	trackref,
		const edm::OrphanHandle< reco::MuonCollection > &	muonh
	)		const

private

Definition at line 1416 of file PFBlockAlgo.cc.

References edm::Ref< C, T, F >::isNonnull(), edm::OrphanHandleBase::isValid(), and j.

                                                                          {
   if(muonh.isValid() ) {
     for(unsigned j=0;j<muonh->size(); ++j) {
       reco::MuonRef muonref( muonh, j );
       if (muonref->track().isNonnull())
     if( muonref->track() == trackref ) return j;
     }
   }
   return -1; // not found
 }

void PFBlockAlgo::packLinks	(	reco::PFBlock &	block,
		const std::vector< PFBlockLink > &	links
	)		const

private

compute missing links in the blocks (the recursive procedure does not build all links)

Definition at line 242 of file PFBlockAlgo.cc.

References reco::PFBlock::bookLinkData(), gather_cfg::cout, debug_, reco::PFBlock::elements(), asciidump::els, link(), reco::PFBlock::linkData(), linkPrefilter(), PFBlockLink::NONE, reco::PFBlock::setLink(), and edm::OwnVector< T, P >::size().

Referenced by findBlocks().

                                                      {
   
   
   const edm::OwnVector< reco::PFBlockElement >& els = block.elements();
   
   block.bookLinkData();
   unsigned elsize = els.size();
   unsigned ilStart = 0;
   //First Loop: update all link data
   for( unsigned i1=0; i1<elsize; ++i1 ) {
     for( unsigned i2=0; i2<i1; ++i2 ) {
       
       // no reflexive link
       //if( i1==i2 ) continue;
       
       double dist = -1;
       
       bool linked = false;
       PFBlock::LinkTest linktest 
     = PFBlock::LINKTEST_RECHIT; 
 
       // are these elements already linked ?
       // this can be optimized
       unsigned linksize = links.size();
       for( unsigned il = ilStart; il<linksize; ++il ) {
     // The following three lines exploits the increasing-element2 ordering of links.
     if ( links[il].element2() < i1 ) ilStart = il;
     if ( links[il].element2() > i1 ) break;
     if( (links[il].element1() == i2 &&
              links[il].element2() == i1) ) {  // yes
       
       dist = links[il].dist();
       linked = true;
 
       //modif-beg     
       //retrieve type of test used to get distance
       linktest = links[il].test();
 #ifdef PFLOW_DEBUG
       if( debug_ )
         cout << "Reading link vector: linktest used=" 
          << linktest 
          << " distance = " << dist 
          << endl; 
 #endif
       //modif-end
       
       break;
     } 
       }
       
       if(!linked) {
     PFBlockLink::Type linktype = PFBlockLink::NONE;
     bool bTestLink = linkPrefilter(&els[i1], &els[i2]);
     if (bTestLink) link( & els[i1], & els[i2], linktype, linktest, dist);
       }
 
       //loading link data according to link test used: RECHIT 
       //block.setLink( i1, i2, chi2, block.linkData() );
 #ifdef PFLOW_DEBUG
       if( debug_ )
     cout << "Setting link between elements " << i1 << " and " << i2
          << " of dist =" << dist << " computed from link test "
          << linktest << endl;
 #endif
       block.setLink( i1, i2, dist, block.linkData(), linktest );
     }
   }
 
 }

void PFBlockAlgo::setDebug ( bool debug )

inline

sets debug printout flag

Definition at line 186 of file PFBlockAlgo.h.

References debug, and debug_.

Referenced by PFRootEventManager::readOptions().

186 {debug_ = debug;}

debug

#define debug

Definition: HDRShower.cc:19

PFBlockAlgo::debug_

bool debug_

if true, debug printouts activated

Definition: PFBlockAlgo.h:377

void PFBlockAlgo::setHOTag ( bool ho )

inline

Definition at line 205 of file PFBlockAlgo.h.

References useHO_.

Referenced by PFRootEventManager::readOptions().

205 { useHO_ = ho;}

PFBlockAlgo::useHO_

bool useHO_

Definition: PFBlockAlgo.h:382

template<template< typename > class T>

void PFBlockAlgo::setInput	(	const T< reco::PFRecTrackCollection > &	trackh,
		const T< reco::GsfPFRecTrackCollection > &	gsftrackh,
		const T< reco::GsfPFRecTrackCollection > &	convbremgsftrackh,
		const T< reco::MuonCollection > &	muonh,
		const T< reco::PFDisplacedTrackerVertexCollection > &	nuclearh,
		const T< reco::PFRecTrackCollection > &	nucleartrackh,
		const T< reco::PFConversionCollection > &	conv,
		const T< reco::PFV0Collection > &	v0,
		const T< reco::PFClusterCollection > &	ecalh,
		const T< reco::PFClusterCollection > &	hcalh,
		const T< reco::PFClusterCollection > &	hoh,
		const T< reco::PFClusterCollection > &	hfemh,
		const T< reco::PFClusterCollection > &	hfhadh,
		const T< reco::PFClusterCollection > &	psh,
		const T< reco::PhotonCollection > &	egphh,
		const T< reco::SuperClusterCollection > &	sceb,
		const T< reco::SuperClusterCollection > &	scee,
		const T< edm::ValueMap< reco::CaloClusterPtr > > &	pfclusterassoc,
		const Mask &	trackMask = `dummyMask_`,
		const Mask &	gsftrackMask = `dummyMask_`,
		const Mask &	ecalMask = `dummyMask_`,
		const Mask &	hcalMask = `dummyMask_`,
		const Mask &	hoMask = `dummyMask_`,
		const Mask &	hfemMask = `dummyMask_`,
		const Mask &	hfhadMask = `dummyMask_`,
		const Mask &	psMask = `dummyMask_`,
		const Mask &	phMask = `dummyMask_`,
		const Mask &	scMask = `dummyMask_`
	)

set input collections of tracks and clusters

–———— GSF Primary tracks and brems ———————

get tracks from converted brems

Loop over the photons

Loop over the SuperClusters

–———— conversions ———————

The tracks from conversions are filled into the elements collection

–———— V0 ———————

The tracks from V0 are filled into the elements collection

One need to cross check if those tracks was not already filled from the conversion collection

This is a new track not yet included into the elements collection

–———— Displaced Vertices ———————

The tracks from Displaced Vertices are filled into the elements collection

One need to cross check if those tracks was not already filled from the conversion or V0 collections

This is a new track not yet included into the elements collection

Fill the displaced vertex ref

--------—— Tracks ---------------——

Mask the tracks in trackh collection already included from Conversions V0 and Displaced Vertices. Take care that all those collections come from the same "generalTracks" collection.

Definition at line 394 of file PFBlockAlgo.h.

Referenced by PFRootEventManager::particleFlow().

                                   {
 
 
   checkMaskSize( *trackh,
          *gsftrackh,
                  *ecalh,
                  *hcalh,
          *hoh,
          *hfemh,
          *hfhadh,
                  *psh,
          *egphh,
          *sceb,
          *scee,      
                  trackMask,
          gsftrackMask,
                  ecalMask,
                  hcalMask,
          hoMask,
          hfemMask,
          hfhadMask,
                  psMask,
          phMask,
          scMask);
 
   /*
   if (nucleartrackh.isValid()){
     for(unsigned i=0;i<nucleartrackh->size(); i++) {
       reco::PFRecTrackRef trackRef(nucleartrackh,i);
       std::cout << *trackRef << std::endl;
     }
   }
   */
 
   if (superClusterMatchByRef_) {
     pfclusterassoc_ = pfclusterassoc.product();
   }
   
   std::vector<reco::PFRecTrackRef> convBremPFRecTracks;
   convBremPFRecTracks.clear();
   // Super cluster mapping
   superClusters_.clear();
   scpfcRefs_.clear();
   pfcSCVec_.clear();
 
   if(gsftrackh.isValid() ) {
     const  reco::GsfPFRecTrackCollection PFGsfProd = *(gsftrackh.product());
     for(unsigned i=0;i<gsftrackh->size(); ++i) {
       if( !gsftrackMask.empty() &&
           !gsftrackMask[i] ) continue;
       reco::GsfPFRecTrackRef refgsf(gsftrackh,i );   
    
       if((refgsf).isNull()) continue;
       reco::GsfTrackRef gsf=refgsf->gsfTrackRef();
 
       // retrieve and save the SC if ECAL-driven - Florian
       if(gsf->extra().isAvailable() && gsf->extra()->seedRef().isAvailable()) {
     reco::ElectronSeedRef seedRef=  gsf->extra()->seedRef().castTo<reco::ElectronSeedRef>();
     // check that the seed is valid
     if(seedRef.isAvailable() && seedRef->isEcalDriven()) {
       reco::SuperClusterRef scRef = seedRef->caloCluster().castTo<reco::SuperClusterRef>();
       if(scRef.isNonnull())   {      
         std::vector<reco::SuperClusterRef>::iterator itcheck=find(superClusters_.begin(),superClusters_.end(),scRef);
         // not present, add it
         if(itcheck==superClusters_.end())
           {
         superClusters_.push_back(scRef);
         reco::PFBlockElementSuperCluster * sce =
           new reco::PFBlockElementSuperCluster(scRef);
         sce->setFromGsfElectron(true);
                 sce->setFromPFSuperCluster(superClusterMatchByRef_);
         elements_.push_back(sce);
           }
         else // it is already present, update the PFBE
           {
         PFBlockElementSCEqual myEqual(scRef);
         std::list<reco::PFBlockElement*>::iterator itcheck=find_if(elements_.begin(),elements_.end(),myEqual);
         if(itcheck!=elements_.end())
           {
             reco::PFBlockElementSuperCluster* thePFBE=dynamic_cast<reco::PFBlockElementSuperCluster*>(*itcheck);
             thePFBE->setFromGsfElectron(true);
             //          std::cout << " Updating element to add electron information" << std::endl;
           }
 //      else
 //        {
 //          std::cout << " Missing element " << std::endl;
 //        }
           }
       }
     }
       }
 
       reco::PFBlockElement* gsfEl;
       
       const  std::vector<reco::PFTrajectoryPoint> 
     PfGsfPoint =  PFGsfProd[i].trajectoryPoints();
   
       unsigned int c_gsf=0;
       bool PassTracker = false;
       bool GetPout = false;
       unsigned int IndexPout = 0;
       
       typedef std::vector<reco::PFTrajectoryPoint>::const_iterator IP;
       for(IP itPfGsfPoint =  PfGsfPoint.begin();  
       itPfGsfPoint!= PfGsfPoint.end();++itPfGsfPoint) {
     
     if (itPfGsfPoint->isValid()){
       int layGsfP = itPfGsfPoint->layer();
       if (layGsfP == -1) PassTracker = true;
       if (PassTracker && layGsfP > 0 && GetPout == false) {
         IndexPout = c_gsf-1;
         GetPout = true;
       }
       //const math::XYZTLorentzVector GsfMoment = itPfGsfPoint->momentum();
       ++c_gsf;
     }
       }
       math::XYZTLorentzVector pin = PfGsfPoint[0].momentum();      
       math::XYZTLorentzVector pout = PfGsfPoint[IndexPout].momentum();
 
       if(useConvBremPFRecTracks_) {
     const std::vector<reco::PFRecTrackRef>& temp_convBremPFRecTracks(refgsf->convBremPFRecTrackRef());
     if(temp_convBremPFRecTracks.size() > 0) {
       for(unsigned int iconv = 0; iconv <temp_convBremPFRecTracks.size(); ++iconv) {
         convBremPFRecTracks.push_back(temp_convBremPFRecTracks[iconv]);
       }
     }
       }
       
       gsfEl = new reco::PFBlockElementGsfTrack(refgsf, pin, pout);
       
       elements_.push_back( gsfEl);
 
       std::vector<reco::PFBrem> pfbrem = refgsf->PFRecBrem();
       
       for (unsigned i2=0;i2<pfbrem.size(); ++i2) {
     const double DP = pfbrem[i2].DeltaP();
     const double SigmaDP =  pfbrem[i2].SigmaDeltaP(); 
     const unsigned int TrajP = pfbrem[i2].indTrajPoint();
     if(TrajP == 99) continue;
 
     reco::PFBlockElement* bremEl;
     bremEl = new reco::PFBlockElementBrem(refgsf,DP,SigmaDP,TrajP);
     elements_.push_back(bremEl);
     
       }
     }
 
   }
 
   if(useEGPhotons_ && egphh.isValid()) {
     unsigned size=egphh->size();
     for(unsigned isc=0; isc<size; ++isc) {
       if(!phMask.empty() && !(phMask)[isc] ) continue;
       if(!photonSelector_->passPhotonSelection((*egphh)[isc])) continue;
       //      std::cout << " Selected a supercluster" << std::endl;
       // Add only the super clusters not already included 
       reco::SuperClusterRef scRef((*egphh)[isc].superCluster());
       std::vector<reco::SuperClusterRef>::iterator itcheck=find(superClusters_.begin(),superClusters_.end(),(*egphh)[isc].superCluster());
       if(itcheck==superClusters_.end())
     {
       superClusters_.push_back(scRef);
       reco::PFBlockElementSuperCluster * sce =
         new reco::PFBlockElementSuperCluster((*egphh)[isc].superCluster());
       fillFromPhoton(egphh,isc,sce);
           sce->setFromPFSuperCluster(superClusterMatchByRef_);
       elements_.push_back(sce);
     }
       else
     {
       PFBlockElementSCEqual myEqual(scRef);
       std::list<reco::PFBlockElement*>::iterator itcheck=find_if(elements_.begin(),elements_.end(),myEqual);
       if(itcheck!=elements_.end())
         {
           reco::PFBlockElementSuperCluster* thePFBE=dynamic_cast<reco::PFBlockElementSuperCluster*>(*itcheck);
           fillFromPhoton(egphh,isc,thePFBE);
           thePFBE->setFromPhoton(true);
           //          std::cout << " Updating element to add Photon information " << photonSelector_->passPhotonSelection((*egphh)[isc]) << std::endl;
 
         }
 //    else
 //      {
 //        std::cout << " Missing element " << std::endl;
 //      }
     }
     }
   }
   
   if(useSuperClusters_ && sceb.isValid() && scee.isValid()) {
     std::vector<reco::SuperClusterRef> screfs;
     screfs.reserve(sceb->size()+scee->size());
     
     for (unsigned int isc=0; isc<sceb->size(); ++isc) {
       reco::SuperClusterRef scRef(sceb,isc);
       screfs.push_back(scRef);
     }
 
     for (unsigned int isc=0; isc<scee->size(); ++isc) {
       reco::SuperClusterRef scRef(scee,isc);
       screfs.push_back(scRef);
     }    
     
     unsigned size=screfs.size();
     for(unsigned isc=0; isc<size; ++isc) {
       if(!scMask.empty() && !(scMask)[isc] ) continue;
       //if(!photonSelector_->passPhotonSelection((*egphh)[isc])) continue;
       
       //      std::cout << " Selected a supercluster" << std::endl;
       // Add only the super clusters not already included 
       reco::SuperClusterRef &scRef = screfs[isc];
       std::vector<reco::SuperClusterRef>::iterator itcheck=find(superClusters_.begin(),superClusters_.end(),scRef);
       if(itcheck==superClusters_.end())
     {
       superClusters_.push_back(scRef);
       reco::PFBlockElementSuperCluster * sce =
         new reco::PFBlockElementSuperCluster(scRef);
           sce->setFromPFSuperCluster(superClusterMatchByRef_);
       //fillFromPhoton(egphh,isc,sce);
       //sce->setFromPhoton(true);
       elements_.push_back(sce);
     }
 //       else
 //  {
 //    PFBlockElementSCEqual myEqual(scRef);
 //    std::list<reco::PFBlockElement*>::iterator itcheck=find_if(elements_.begin(),elements_.end(),myEqual);
 //    if(itcheck!=elements_.end())
 //      {
 //        reco::PFBlockElementSuperCluster* thePFBE=dynamic_cast<reco::PFBlockElementSuperCluster*>(*itcheck);
 //        //fillFromPhoton(egphh,isc,thePFBE);
 //        //thePFBE->setFromPhoton(true);
 //        //          std::cout << " Updating element to add Photon information " << photonSelector_->passPhotonSelection((*egphh)[isc]) << std::endl;
 // 
 //      }
 // //     else
 // //       {
 // //         std::cout << " Missing element " << std::endl;
 // //       }
 //  }
     }
   }  
   
   // set the vector to the right size so to allow random access
   scpfcRefs_.resize(superClusters_.size());
 
 
 
   if(convh.isValid() ) {
     reco::PFBlockElement* trkFromConversionElement;
     for(unsigned i=0;i<convh->size(); ++i) {
       reco::PFConversionRef convRef(convh,i);
 
       unsigned int trackSize=(convRef->pfTracks()).size();
       if ( convRef->pfTracks().size() < 2) continue;
       for(unsigned iTk=0;iTk<trackSize; ++iTk) {
     
     reco::PFRecTrackRef compPFTkRef = convRef->pfTracks()[iTk]; 
     trkFromConversionElement = new reco::PFBlockElementTrack(convRef->pfTracks()[iTk]);
     trkFromConversionElement->setConversionRef( convRef->originalConversion(), reco::PFBlockElement::T_FROM_GAMMACONV);
 
     elements_.push_back( trkFromConversionElement );
 
 
     if (debug_){
       std::cout << "PF Block Element from Conversion electron " << 
         (*trkFromConversionElement).trackRef().key() << std::endl;
       std::cout << *trkFromConversionElement << std::endl;
     }
        
       }     
     }  
   }
   
   
   
   
   if(v0.isValid() ) {
     reco::PFBlockElement* trkFromV0Element = 0;
     for(unsigned i=0;i<v0->size(); ++i) {
       reco::PFV0Ref v0Ref( v0, i );
       unsigned int trackSize=(v0Ref->pfTracks()).size();
       for(unsigned iTk=0;iTk<trackSize; ++iTk) {
 
     reco::PFRecTrackRef newPFRecTrackRef = (v0Ref->pfTracks())[iTk]; 
     reco::TrackBaseRef newTrackBaseRef(newPFRecTrackRef->trackRef());
     bool bNew = true;
     
     for(IE iel = elements_.begin(); iel != elements_.end(); ++iel){
       reco::TrackBaseRef elemTrackBaseRef((*iel)->trackRef());
       if (newTrackBaseRef == elemTrackBaseRef){     
         trkFromV0Element = *iel;
         bNew = false;
         continue;
       }
     } 
 
     if (bNew) {
       trkFromV0Element = new reco::PFBlockElementTrack(v0Ref->pfTracks()[iTk]);
       elements_.push_back( trkFromV0Element );
     }
 
     trkFromV0Element->setV0Ref( v0Ref->originalV0(),
                     reco::PFBlockElement::T_FROM_V0 );
       
     if (debug_){
       std::cout << "PF Block Element from V0 track New = " << bNew 
             << (*trkFromV0Element).trackRef().key() << std::endl;
       std::cout << *trkFromV0Element << std::endl;
     }
 
     
       }
     }
   }
   
 
 
   if(nuclearh.isValid()) {
     reco::PFBlockElement* trkFromDisplacedVertexElement = 0;
     for(unsigned i=0;i<nuclearh->size(); ++i) {
 
       const reco::PFDisplacedTrackerVertexRef dispacedVertexRef( nuclearh, i );
 
       //      std::cout << "Nuclear Interactions Purity " <<  nuclearInteractionsPurity_ << std::endl;
       //     dispacedVertexRef->displacedVertexRef()->Dump();
       //bool bIncludeVertices = true;
       // We add a cut at rho > 2.7 since this corresponds to the lower edge of the beam pipe
       // This cut have to be changer when a new beam pipe would be installed
       
       bool bIncludeVertices = false; 
       bool bNucl = dispacedVertexRef->displacedVertexRef()->isNucl()
     && dispacedVertexRef->displacedVertexRef()->position().rho()>  2.7;
       bool bNucl_Loose = dispacedVertexRef->displacedVertexRef()->isNucl_Loose();
       bool bNucl_Kink = dispacedVertexRef->displacedVertexRef()->isNucl_Kink();
 
       if (nuclearInteractionsPurity_ >= 1) bIncludeVertices = bNucl;
       if (nuclearInteractionsPurity_ >= 2) bIncludeVertices = bIncludeVertices || bNucl_Loose;
       if (nuclearInteractionsPurity_ >= 3) bIncludeVertices = bIncludeVertices || bNucl_Kink;
 
       if (bIncludeVertices){
 
     unsigned int trackSize= dispacedVertexRef->pfRecTracks().size();
     if (debug_){
       std::cout << "" << std::endl;
       std::cout << "Displaced Vertex " << i << std::endl;
       dispacedVertexRef->displacedVertexRef()->Dump();
     }
     for(unsigned iTk=0;iTk < trackSize; ++iTk) {
 
 
       // This peace of code looks weired at first but it seems to be necessary to let 
       // PFRooTEvent work properly. Since the track position called REPPoint is transient 
       // it has to be calculated and the PFRecTrack collection associted to Displaced Vertex is not
       // anymore the original collection. So here we match both collections if possible
       reco::PFRecTrackRef newPFRecTrackRef = dispacedVertexRef->pfRecTracks()[iTk]; 
       reco::TrackBaseRef constTrackBaseRef(newPFRecTrackRef->trackRef());
 
       
       if (nucleartrackh.isValid()){
         for(unsigned i=0;i<nucleartrackh->size(); ++i) {
           reco::PFRecTrackRef transientPFRecTrackRef(nucleartrackh,i);
           reco::TrackBaseRef transientTrackBaseRef(transientPFRecTrackRef->trackRef());
           if (constTrackBaseRef==transientTrackBaseRef){
         newPFRecTrackRef = transientPFRecTrackRef;
         break;
           }
         }
       }
       reco::TrackBaseRef newTrackBaseRef(newPFRecTrackRef->trackRef());
       
 
 
       bool bNew = true;
       reco::PFBlockElement::TrackType blockType;
 
       for(IE iel = elements_.begin(); iel != elements_.end(); ++iel){
         reco::TrackBaseRef elemTrackBaseRef((*iel)->trackRef());
         if (newTrackBaseRef == elemTrackBaseRef){
           trkFromDisplacedVertexElement = *iel;
           bNew = false;
           continue;
         }
       }
 
 
       if (bNew) { 
         
 
 
         trkFromDisplacedVertexElement = new reco::PFBlockElementTrack(newPFRecTrackRef);
         elements_.push_back( trkFromDisplacedVertexElement );
       }
 
       if (dispacedVertexRef->isIncomingTrack(newPFRecTrackRef)) 
         blockType = reco::PFBlockElement::T_TO_DISP;
       else if (dispacedVertexRef->isOutgoingTrack(newPFRecTrackRef)) 
         blockType = reco::PFBlockElement::T_FROM_DISP;
       else 
         blockType = reco::PFBlockElement::DEFAULT;
 
       trkFromDisplacedVertexElement->setDisplacedVertexRef( dispacedVertexRef, blockType );
 
 
       if (debug_){
         std::cout << "PF Block Element from DisplacedTrackingVertex track New = " << bNew
               << (*trkFromDisplacedVertexElement).trackRef().key() << std::endl;
         std::cout << *trkFromDisplacedVertexElement << std::endl;
       }
     
     
     }
       }
     }  
 
     if (debug_) std::cout << "" << std::endl;
 
   }
 
 
 
   if(trackh.isValid() ) {
 
     if (debug_) std::cout << "Tracks already in from Displaced Vertices " << std::endl;
 
     Mask trackMaskVertex;
 
     for(unsigned i=0;i<trackh->size(); ++i) {
       reco::PFRecTrackRef pfRefTrack( trackh,i );
       reco::TrackRef trackRef = pfRefTrack->trackRef();
 
       bool bMask = true;
       for(IE iel = elements_.begin(); iel != elements_.end(); ++iel){
     reco::TrackRef elemTrackRef = (*iel)->trackRef();
     if( trackRef == elemTrackRef ) {
       if (debug_) std::cout << " " << trackRef.key();
       bMask = false; continue;
     }
       }
     
       trackMaskVertex.push_back(bMask);
     }
 
     if (debug_) std::cout << "" << std::endl;
 
     if (debug_) std::cout << "Additionnal tracks from main collection " << std::endl;
 
     for(unsigned i=0;i<trackh->size(); ++i) {
 
 
       // this track has been disabled
       if( !trackMask.empty() && !trackMask[i] ) continue;
       
       reco::PFRecTrackRef ref( trackh,i );
 
       if (debug_) std::cout << " " << ref->trackRef().key();
 
       // Get the eventual muon associated to this track
       int muId_ = muAssocToTrack( ref->trackRef(), muonh );
       bool thisIsAPotentialMuon = false;
       if( muId_ != -1 ) {
     reco::MuonRef muonref( muonh, muId_ );
     thisIsAPotentialMuon = 
       PFMuonAlgo::isLooseMuon(muonref) || 
       PFMuonAlgo::isMuon(muonref);
       }
       // Reject bad tracks (except if identified as muon
       if( !thisIsAPotentialMuon && !goodPtResolution( ref->trackRef() ) ) continue;
 
       if (thisIsAPotentialMuon && debug_) std::cout << "Potential Muon P " <<  ref->trackRef()->p() 
                             << " pt " << ref->trackRef()->p() << std::endl; 
 
 
 
       reco::PFBlockElement* primaryElement = new reco::PFBlockElementTrack( ref );
 
       if( muId_ != -1 ) {
     // if a muon has been found
     reco::MuonRef muonref( muonh, muId_ );
 
     // If this track was already added to the collection, we just need to find the associated element and 
     // attach to it the reference
     if (!trackMaskVertex.empty() && !trackMaskVertex[i]){
       reco::TrackRef primaryTrackRef = ref->trackRef();
       for(IE iel = elements_.begin(); iel != elements_.end(); ++iel){
         reco::TrackRef elemTrackRef = (*iel)->trackRef();
         if( primaryTrackRef == elemTrackRef ) {
           (*iel)->setMuonRef( muonref );
           if (debug_) std::cout << "One of the tracks identified in displaced vertices collections was spotted as muon" <<std:: endl;
         }
       }
     } else primaryElement->setMuonRef( muonref );
 
       }
 
       if (!trackMaskVertex.empty() && !trackMaskVertex[i]) continue;
 
       
       // set track type T_FROM_GAMMA for pfrectracks associated to conv brems
       if(useConvBremPFRecTracks_) {
     if(convBremPFRecTracks.size() > 0.) {
       for(unsigned int iconv = 0; iconv < convBremPFRecTracks.size(); ++iconv) {
         if((*ref).trackRef() == (*convBremPFRecTracks[iconv]).trackRef()) {
           bool value = true;
           primaryElement->setTrackType(reco::PFBlockElement::T_FROM_GAMMACONV, value);
         }
       }
     }
       }
       elements_.push_back( primaryElement );
     }
 
     if (debug_) std::cout << " " << std::endl;
 
   }
 
  
   // -------------- GSF tracks and brems for Conversion Recovery ----------
    
   if(convbremgsftrackh.isValid() ) {
     
  
     const  reco::GsfPFRecTrackCollection ConvPFGsfProd = *(convbremgsftrackh.product());
     for(unsigned i=0;i<convbremgsftrackh->size(); ++i) {
 
       reco::GsfPFRecTrackRef refgsf(convbremgsftrackh,i );   
       
       if((refgsf).isNull()) continue;
       
       reco::PFBlockElement* gsfEl;
       
       const  std::vector<reco::PFTrajectoryPoint> 
     PfGsfPoint =  ConvPFGsfProd[i].trajectoryPoints();
       
       unsigned int c_gsf=0;
       bool PassTracker = false;
       bool GetPout = false;
       unsigned int IndexPout = -1;
       
       typedef std::vector<reco::PFTrajectoryPoint>::const_iterator IP;
       for(IP itPfGsfPoint =  PfGsfPoint.begin();  
       itPfGsfPoint!= PfGsfPoint.end();++itPfGsfPoint) {
     
     if (itPfGsfPoint->isValid()){
       int layGsfP = itPfGsfPoint->layer();
       if (layGsfP == -1) PassTracker = true;
       if (PassTracker && layGsfP > 0 && GetPout == false) {
         IndexPout = c_gsf-1;
         GetPout = true;
       }
       //const math::XYZTLorentzVector GsfMoment = itPfGsfPoint->momentum();
       ++c_gsf;
     }
       }
       math::XYZTLorentzVector pin = PfGsfPoint[0].momentum();      
       math::XYZTLorentzVector pout = PfGsfPoint[IndexPout].momentum();
       
       
     
       gsfEl = new reco::PFBlockElementGsfTrack(refgsf, pin, pout);
       
       bool valuegsf = true;
       // IMPORTANT SET T_FROM_GAMMACONV trackType() FOR CONVERSIONS
       gsfEl->setTrackType(reco::PFBlockElement::T_FROM_GAMMACONV, valuegsf);
 
       
 
       elements_.push_back( gsfEl);
       std::vector<reco::PFBrem> pfbrem = refgsf->PFRecBrem();
       
       for (unsigned i2=0;i2<pfbrem.size(); ++i2) {
     const double DP = pfbrem[i2].DeltaP();
     const double SigmaDP =  pfbrem[i2].SigmaDeltaP(); 
     const unsigned int TrajP = pfbrem[i2].indTrajPoint();
     if(TrajP == 99) continue;
 
     reco::PFBlockElement* bremEl;
     bremEl = new reco::PFBlockElementBrem(refgsf,DP,SigmaDP,TrajP);
     elements_.push_back(bremEl);
     
       }
     }
   }
 
   
   // -------------- ECAL clusters ---------------------
 
 
   if(ecalh.isValid() ) {
     //  pfcSCVec_.resize(ecalh->size(),-1);
  
     bNoSuperclus_  = (superClusters_.size() == 0);
     if (!bNoSuperclus_) pfcSCVec_.resize(ecalh->size(),-1);
 
 
     for(unsigned i=0;i<ecalh->size(); ++i)  {
 
       // this ecal cluster has been disabled
       if( !ecalMask.empty() &&
           !ecalMask[i] ) continue;
 
       reco::PFClusterRef ref( ecalh,i );
       reco::PFBlockElement* te
         = new reco::PFBlockElementCluster( ref,
                        reco::PFBlockElement::ECAL);
       elements_.push_back( te );
 
       if (!bNoSuperclus_) {
 
     // Now mapping with Superclusters
     int scindex = -1;
     if (superClusterMatchByRef_) {
           scindex = ClusterClusterMapping::checkOverlap(ref,superClusters_,*pfclusterassoc_);
         }
     else {
           scindex= ClusterClusterMapping::checkOverlap(*ref,superClusters_);
         }
 
     if(scindex>=0)  {
       pfcSCVec_[ref.key()]=scindex;
       scpfcRefs_[scindex].push_back(ref);
     }
       }
 
     }
 
     bNoSuperclus_ = (scpfcRefs_.size() == 0);
 
   }
 
   // -------------- HCAL clusters ---------------------
 
   if(hcalh.isValid() ) {
     
     for(unsigned i=0;i<hcalh->size(); ++i)  {
       
       // this hcal cluster has been disabled
       if( !hcalMask.empty() &&
           !hcalMask[i] ) continue;
       
       reco::PFClusterRef ref( hcalh,i );
       reco::PFBlockElement* th
         = new reco::PFBlockElementCluster( ref,
                        reco::PFBlockElement::HCAL );
       elements_.push_back( th );
     }
   }
 
   // -------------- HO clusters ---------------------
 
   if(useHO_ && hoh.isValid() ) {
     
     for(unsigned i=0;i<hoh->size(); ++i)  {
       
       // this hcal cluster has been disabled
       if( !hoMask.empty() &&
           !hoMask[i] ) continue;
       
       reco::PFClusterRef ref( hoh,i );
       reco::PFBlockElement* th
         = new reco::PFBlockElementCluster( ref,
                        reco::PFBlockElement::HO );
       elements_.push_back( th );
     }
   }
 
   // -------------- HFEM clusters ---------------------
 
   if(hfemh.isValid() ) {
     
     for(unsigned i=0;i<hfemh->size(); ++i)  {
       
       // this hfem cluster has been disabled
       if( !hfemMask.empty() &&
           !hfemMask[i] ) continue;
       
       reco::PFClusterRef ref( hfemh,i );
       reco::PFBlockElement* th
         = new reco::PFBlockElementCluster( ref,
                        reco::PFBlockElement::HFEM );
       elements_.push_back( th );
     }
   }
 
 
   // -------------- HFHAD clusters ---------------------
 
   if(hfhadh.isValid() ) {
     
     for(unsigned i=0;i<hfhadh->size(); ++i)  {
       
       // this hfhad cluster has been disabled
       if( !hfhadMask.empty() &&
           !hfhadMask[i] ) continue;
       
       reco::PFClusterRef ref( hfhadh,i );
       reco::PFBlockElement* th
         = new reco::PFBlockElementCluster( ref,
                        reco::PFBlockElement::HFHAD );
       elements_.push_back( th );
     }
   }
 
 
 
 
   // -------------- PS clusters ---------------------
 
   if(psh.isValid() ) {
     for(unsigned i=0;i<psh->size(); ++i)  {
 
       // this ps cluster has been disabled
       if( !psMask.empty() &&
           !psMask[i] ) continue;
       reco::PFBlockElement::Type type = reco::PFBlockElement::NONE;
       reco::PFClusterRef ref( psh,i );
       // two types of elements:  PS1 (V) and PS2 (H) 
       // depending on layer:  PS1 or PS2
       switch(ref->layer()){
       case PFLayer::PS1:
         type = reco::PFBlockElement::PS1;
         break;
       case PFLayer::PS2:
         type = reco::PFBlockElement::PS2;
         break;
       default:
         break;
       }
       reco::PFBlockElement* tp
         = new reco::PFBlockElementCluster( ref,
                        type );
       elements_.push_back( tp );
     }
   }
 
 
   // -------------- Loop over block elements ---------------------
 
   // Here we provide to all KDTree linkers the collections to link.
   // Glowinski & Gouzevitch
   
   for (std::list< reco::PFBlockElement* >::iterator it = elements_.begin();
        it != elements_.end(); ++it) {
     switch ((*it)->type()){
     
     case reco::PFBlockElement::TRACK:
       if (useKDTreeTrackEcalLinker_) {
     if ( (*it)->trackRefPF()->extrapolatedPoint( reco::PFTrajectoryPoint::ECALShowerMax ).isValid() )
       TELinker_.insertTargetElt(*it);
     if ( (*it)->trackRefPF()->extrapolatedPoint( reco::PFTrajectoryPoint::HCALEntrance ).isValid() )
       THLinker_.insertTargetElt(*it);
       }
       
       break;
 
     case reco::PFBlockElement::PS1:
     case reco::PFBlockElement::PS2:
       if (useKDTreeTrackEcalLinker_)
     PSELinker_.insertTargetElt(*it);
       break;
 
     case reco::PFBlockElement::HCAL:
       if (useKDTreeTrackEcalLinker_)
     THLinker_.insertFieldClusterElt(*it);
       break;
 
     case reco::PFBlockElement::HO: 
       if (useHO_ && useKDTreeTrackEcalLinker_) {
     // THLinker_.insertFieldClusterElt(*it);
       }
       break;
 
     
     case reco::PFBlockElement::ECAL:
       if (useKDTreeTrackEcalLinker_) {
     TELinker_.insertFieldClusterElt(*it);
     PSELinker_.insertFieldClusterElt(*it);
       }
       break;
 
     default:
       break;
     }
   }
 }

template<template< typename > class T>

void PFBlockAlgo::setInput	(	const T< reco::PFRecTrackCollection > &	trackh,
		const T< reco::MuonCollection > &	muonh,
		const T< reco::PFClusterCollection > &	ecalh,
		const T< reco::PFClusterCollection > &	hcalh,
		const T< reco::PFClusterCollection > &	hoh,
		const T< reco::PFClusterCollection > &	hfemh,
		const T< reco::PFClusterCollection > &	hfhadh,
		const T< reco::PFClusterCollection > &	psh,
		const Mask &	trackMask = `dummyMask_`,
		const Mask &	ecalMask = `dummyMask_`,
		const Mask &	hcalMask = `dummyMask_`,
		const Mask &	hoMask = `dummyMask_`,
		const Mask &	psMask = `dummyMask_`
	)

inline

COLIN: I think this is for particle flow at HLT...

Definition at line 130 of file PFBlockAlgo.h.

                                            {
     T<reco::GsfPFRecTrackCollection> gsftrackh;
     T<reco::GsfPFRecTrackCollection> convbremgsftrackh;
     //T<reco::MuonCollection> muonh;
     T<reco::PFDisplacedTrackerVertexCollection> nuclearh;
     T<reco::PFRecTrackCollection>    nucleartrackh;
     T<reco::PFConversionCollection> convh;
     T<reco::PFV0Collection> v0;
     T<reco::PhotonCollection> phh;
     T<reco::SuperClusterCollection> scebh;
     T<reco::SuperClusterCollection> sceeh;    
     T<edm::ValueMap<reco::CaloClusterPtr> > pfclusterassoc;
     setInput<T>( trackh, gsftrackh, convbremgsftrackh, muonh, nuclearh, nucleartrackh, convh, v0, 
          ecalh, hcalh, hoh, hfemh, hfhadh, psh, phh, scebh, sceeh, pfclusterassoc,
          trackMask, ecalMask, hcalMask, hoMask, psMask); 
   }

template<template< typename > class T>

void PFBlockAlgo::setInput	(	const T< reco::PFRecTrackCollection > &	trackh,
		const T< reco::GsfPFRecTrackCollection > &	gsftrackh,
		const T< reco::PFClusterCollection > &	ecalh,
		const T< reco::PFClusterCollection > &	hcalh,
		const T< reco::PFClusterCollection > &	hoh,
		const T< reco::PFClusterCollection > &	psh,
		const Mask &	trackMask = `dummyMask_`,
		const Mask &	gsftrackMask = `dummyMask_`,
		const Mask &	ecalMask = `dummyMask_`,
		const Mask &	hcalMask = `dummyMask_`,
		const Mask &	hoMask = `dummyMask_`,
		const Mask &	psMask = `dummyMask_`
	)

inline

COLIN: what is this setinput function for? can it be removed?

Definition at line 161 of file PFBlockAlgo.h.

                                            {
     T<reco::GsfPFRecTrackCollection> convbremgsftrackh;
     T<reco::MuonCollection> muonh;
     T<reco::PFDisplacedTrackerVertexCollection>  nuclearh;
     T<reco::PFRecTrackCollection>    nucleartrackh;
     T<reco::PFConversionCollection> convh;
     T<reco::PFV0Collection> v0;
     T<reco::PhotonCollection> egphh;
     setInput<T>( trackh, gsftrackh, convbremgsftrackh, muonh, nuclearh, nucleartrackh, convh, v0, ecalh, hcalh, hoh, psh, egphh,
          trackMask, gsftrackMask,ecalMask, hcalMask, hoMask, psMask); 
   }

void PFBlockAlgo::setParameters	(	std::vector< double > &	DPtovPtCut,
		std::vector< unsigned > &	NHitCut,
		bool	useConvBremPFRecTracks,
		bool	useIterTracking,
		int	nuclearInteractionsPurity,
		bool	useEGPhotons,
		std::vector< double > &	photonSelectionCuts,
		bool	useSuperClusters,
		bool	superClusterMatchByRef
	)

Definition at line 31 of file PFBlockAlgo.cc.

References DPtovPtCut_, NHitCut_, nuclearInteractionsPurity_, photonSelector_, superClusterMatchByRef_, useConvBremPFRecTracks_, useEGPhotons_, useIterTracking_, and useSuperClusters_.

Referenced by PFRootEventManager::readOptions().

                                  {
   
   DPtovPtCut_    = DPtovPtCut;
   NHitCut_       = NHitCut;
   useIterTracking_ = useIterTracking;
   useConvBremPFRecTracks_ = useConvBremPFRecTracks;
   nuclearInteractionsPurity_ = nuclearInteractionsPurity;
   useEGPhotons_ = useEGPhotons;
   // Pt cut; Track iso (constant + slope), Ecal iso (constant + slope), HCAL iso (constant+slope), H/E
   if(useEGPhotons_)
     photonSelector_ = new PhotonSelectorAlgo(photonSelectionCuts[0],   
                          photonSelectionCuts[1], photonSelectionCuts[2],   
                          photonSelectionCuts[3], photonSelectionCuts[4],    
                          photonSelectionCuts[5], photonSelectionCuts[6],    
                          photonSelectionCuts[7],
                          photonSelectionCuts[8],
                          photonSelectionCuts[9],
                          photonSelectionCuts[10]
                          );
 
 
   useSuperClusters_ = useSuperClusters;
   superClusterMatchByRef_ = superClusterMatchByRef;
 }

void PFBlockAlgo::setUseOptimization ( bool useKDTreeTrackEcalLinker )

Definition at line 66 of file PFBlockAlgo.cc.

References useKDTreeTrackEcalLinker_.

Referenced by PFRootEventManager::readOptions().

 {
   useKDTreeTrackEcalLinker_ = useKDTreeTrackEcalLinker;
 }

double PFBlockAlgo::testECALAndHCAL	(	const reco::PFCluster &	ecal,
		const reco::PFCluster &	hcal
	)		const

private

tests association between an ECAL and an HCAL cluster

Returns: distance

Definition at line 936 of file PFBlockAlgo.cc.

References LinkByRecHit::computeDist(), gather_cfg::cout, debug_, and reco::PFCluster::positionREP().

Referenced by link().

                                                {
   
   //   cout<<"entering testECALAndHCAL"<<endl;
   
   double dist = fabs(ecal.positionREP().Eta()) > 2.5 ?
     LinkByRecHit::computeDist( ecal.positionREP().Eta(),
                    ecal.positionREP().Phi(), 
                    hcal.positionREP().Eta(), 
                    hcal.positionREP().Phi() )
     : 
     -1.;
 
 #ifdef PFLOW_DEBUG
   if(debug_) cout<<"testECALAndHCAL "<< dist <<" "<<endl;
   if(debug_){
     cout<<" ecaleta " << ecal.positionREP().Eta()
     <<" ecalphi " << ecal.positionREP().Phi()
     <<" hcaleta " << hcal.positionREP().Eta()
     <<" hcalphi " << hcal.positionREP().Phi()
   }
 #endif
 
   if ( dist < 0.2 ) return dist; 
  
   // Need to implement a link by RecHit
   return -1.;
 }

double PFBlockAlgo::testHCALAndHO	(	const reco::PFCluster &	hcal,
		const reco::PFCluster &	ho
	)		const

private

tests association between an HCAL and an HO cluster

Returns: distance

Definition at line 966 of file PFBlockAlgo.cc.

References LinkByRecHit::computeDist(), gather_cfg::cout, debug_, and reco::PFCluster::positionREP().

Referenced by link().

                                              {
   
   double dist = fabs(hcal.positionREP().Eta()) < 1.5 ?
     LinkByRecHit::computeDist( hcal.positionREP().Eta(),
                    hcal.positionREP().Phi(), 
                    ho.positionREP().Eta(), 
                    ho.positionREP().Phi() )
     : 
     -1.;
 
 #ifdef PFLOW_DEBUG
   if(debug_) cout<<"testHCALAndHO "<< dist <<" "<<endl;
   if(debug_){
     cout<<" hcaleta " << hcal.positionREP().Eta()
     <<" hcalphi " << hcal.positionREP().Phi()
     <<" hoeta " << ho.positionREP().Eta()
     <<" hophi " << ho.positionREP().Phi()
     <<" dist " << dist<<endl;
   }
 #endif
 
   if ( dist < 0.20 ) return dist; 
  
   // Need to implement a link by RecHit
   return -1.;
 }

double PFBlockAlgo::testLinkBySuperCluster	(	const reco::PFClusterRef &	elt1,
		const reco::PFClusterRef &	elt2
	)		const

private

test association by Supercluster between two ECAL

Definition at line 997 of file PFBlockAlgo.cc.

References LinkByRecHit::computeDist(), i, edm::Ref< C, T, F >::key(), pfcSCVec_, and scpfcRefs_.

Referenced by link().

                                                       {
   
   //  cout<<"entering testECALAndECAL "<< pfcRefSCMap_.size() << endl;
   
   double dist = -1;
   
   // the first one is not in any super cluster
   int testindex=pfcSCVec_[ecal1.key()];
   if(testindex == -1.) return dist;
   //  if(itcheck==pfcRefSCMap_.end()) return dist;
   // now retrieve the of PFclusters in this super cluster  
 
   const std::vector<reco::PFClusterRef> & thePFClusters(scpfcRefs_[testindex]);
   
   unsigned npf=thePFClusters.size();
   for(unsigned i=0;i<npf;++i)
     {
       if(thePFClusters[i]==ecal2) // yes they are in the same SC 
     {
       dist=LinkByRecHit::computeDist( ecal1->positionREP().Eta(),
                       ecal1->positionREP().Phi(), 
                       ecal2->positionREP().Eta(), 
                       ecal2->positionREP().Phi() );
 //    std::cout << " DETA " << fabs(ecal1->positionREP().Eta()-ecal2->positionREP().Eta()) << std::endl;
 //    if(fabs(ecal1->positionREP().Eta()-ecal2->positionREP().Eta())>0.2)
 //      {
 //        std::cout <<  " Super Cluster " <<  *(superClusters_[testindex]) << std::endl;
 //        std::cout <<  " Cluster1 " <<  *ecal1 << std::endl;
 //        std::cout <<  " Cluster2 " <<  *ecal2 << std::endl;
 //        ClusterClusterMapping::checkOverlap(*ecal1,superClusters_,0.01,true);
 //        ClusterClusterMapping::checkOverlap(*ecal2,superClusters_,0.01,true);
 //      }
       return dist;
     }
     }
   return dist;
 }

double PFBlockAlgo::testLinkByVertex	(	const reco::PFBlockElement *	elt1,
		const reco::PFBlockElement *	elt2
	)		const

private

Definition at line 1120 of file PFBlockAlgo.cc.

References reco::PFBlockElement::convRef(), gather_cfg::cout, debug_, reco::PFBlockElement::displacedVertexRef(), edm::Ref< C, T, F >::isNonnull(), query::result, reco::PFBlockElement::T_FROM_DISP, reco::PFBlockElement::T_FROM_GAMMACONV, reco::PFBlockElement::T_FROM_V0, reco::PFBlockElement::T_TO_DISP, reco::PFBlockElement::trackType(), and reco::PFBlockElement::V0Ref().

Referenced by link().

                                                          {
 
   //  cout << "Test link by vertex between" << endl << *elt1 << endl << " and " << endl << *elt2 << endl;
 
   double result=-1.;
 
   reco::PFBlockElement::TrackType T_TO_DISP = reco::PFBlockElement::T_TO_DISP;
   reco::PFBlockElement::TrackType T_FROM_DISP = reco::PFBlockElement::T_FROM_DISP;
   PFDisplacedTrackerVertexRef ni1_TO_DISP = elt1->displacedVertexRef(T_TO_DISP);
   PFDisplacedTrackerVertexRef ni2_TO_DISP = elt2->displacedVertexRef(T_TO_DISP);
   PFDisplacedTrackerVertexRef ni1_FROM_DISP = elt1->displacedVertexRef(T_FROM_DISP);
   PFDisplacedTrackerVertexRef ni2_FROM_DISP = elt2->displacedVertexRef(T_FROM_DISP);
   
   if( ni1_TO_DISP.isNonnull() && ni2_FROM_DISP.isNonnull())
     if( ni1_TO_DISP == ni2_FROM_DISP ) { result = 1.0; return result; }
 
   if( ni1_FROM_DISP.isNonnull() && ni2_TO_DISP.isNonnull())
     if( ni1_FROM_DISP == ni2_TO_DISP ) { result = 1.0; return result; }
 
   if( ni1_FROM_DISP.isNonnull() && ni2_FROM_DISP.isNonnull())
     if( ni1_FROM_DISP == ni2_FROM_DISP ) { result = 1.0; return result; }
     
   
   if (  elt1->trackType(reco::PFBlockElement::T_FROM_GAMMACONV)  &&
          elt2->trackType(reco::PFBlockElement::T_FROM_GAMMACONV)  ) {
     
     if(debug_ ) std::cout << " testLinkByVertex On Conversions " << std::endl;
     
     if ( elt1->convRef().isNonnull() && elt2->convRef().isNonnull() ) {
       if(debug_ ) std::cout << " PFBlockAlgo.cc testLinkByVertex  Cconversion Refs are non null  " << std::endl;      
       if ( elt1->convRef() ==  elt2->convRef() ) {
     result=1.0;
     if(debug_ ) std::cout << " testLinkByVertex  Cconversion Refs are equal  " << std::endl;      
     return result;
       }
     } 
     
   }
   
   if (  elt1->trackType(reco::PFBlockElement::T_FROM_V0)  &&
              elt2->trackType(reco::PFBlockElement::T_FROM_V0)  ) {
     if(debug_ ) std::cout << " testLinkByVertex On V0 " << std::endl;
     if ( elt1->V0Ref().isNonnull() && elt2->V0Ref().isNonnull() ) {
       if(debug_ ) std::cout << " PFBlockAlgo.cc testLinkByVertex  V0 Refs are non null  " << std::endl;
       if ( elt1->V0Ref() ==  elt2->V0Ref() ) {
     result=1.0;
     if(debug_ ) std::cout << " testLinkByVertex  V0 Refs are equal  " << std::endl;
     return result;
       }
     }
   }
 
   return result;
 }

double PFBlockAlgo::testPS1AndPS2	(	const reco::PFCluster &	ps1,
		const reco::PFCluster &	ps2
	)		const

private

tests association between a PS1 v cluster and a PS2 h cluster returns distance

Definition at line 1068 of file PFBlockAlgo.cc.

References gather_cfg::cout, debug_, reco::CaloCluster::position(), pileupReCalc_HLTpaths::scale, and mathSSE::sqrt().

                                             {
   
 #ifdef PFLOW_DEBUG
   //   cout<<"entering testPS1AndPS2"<<endl;
   
   // compute chi2 in y, z using swimming formulae
   // y2 = y1 * z2/z1   and x2 = x1 *z2/z1
   
   // ps position1  x, y, z
   double x1 = ps1.position().X();
   double y1 = ps1.position().Y();
   double z1 = ps1.position().Z();
   double x2 = ps2.position().X();
   double y2 = ps2.position().Y();
   double z2 = ps2.position().Z();
   // MDN Bug correction Jan 09: check that z1 and z2 have the same sign!
   if (z1*z2<0.) -1.;
   // swim to PS2
   double scale = z2/z1;
   double x1atPS2 = x1*scale;
   double y1atPS2 = y1*scale;
   // resolution of PS cluster dxdx and dydy from strip pitch and length
   // vertical strips in PS1, measure x with pitch precision
   double dx1dx1 = resPSpitch_*resPSpitch_*scale*scale;
   double dy1dy1 = resPSlength_*resPSlength_*scale*scale;
   // horizontal strips in PS2 , measure y with pitch precision
   double dy2dy2 = resPSpitch_*resPSpitch_;
   double dx2dx2 = resPSlength_*resPSlength_;
   
   // double chi2 = (x2-x1atPS2)*(x2-x1atPS2)/(dx1dx1 + dx2dx2) 
   //  + (y2-y1atPS2)*(y2-y1atPS2)/(dy1dy1 + dy2dy2);
   
   double dist = std::sqrt( (x2-x1atPS2)*(x2-x1atPS2)
              + (y2-y1atPS2)*(y2-y1atPS2));
     
   if(debug_) cout<<"testPS1AndPS2 "<<dist<<" "<<endl;
   if(debug_){
     cout<<" x1atPS2 "<< x1atPS2 << " dx1 "<<resPSpitch_*scale
     <<" y1atPS2 "<< y1atPS2 << " dy1 "<<resPSlength_*scale<< endl
     <<" x2 " <<x2  << " dx2 "<<resPSlength_
     <<" y2 " << y2 << " dy2 "<<resPSpitch_<< endl;
   }
 #endif
 
   // Need a link by rechit here
   return -1.; 
 }

double PFBlockAlgo::testSuperClusterPFCluster	(	const reco::SuperClusterRef &	sct1,
		const reco::PFClusterRef &	elt2
	)		const

private

test association between SuperClusters and ECAL

Definition at line 1038 of file PFBlockAlgo.cc.

References LinkByRecHit::computeDist(), ClusterClusterMapping::overlap(), muon::overlap(), pfclusterassoc_, and superClusterMatchByRef_.

Referenced by link().

                                                          {
   
   //  cout<<"entering testECALAndECAL "<< pfcRefSCMap_.size() << endl;
   
   double dist = -1;
   
   if(superClusterMatchByRef_) {
     //loop over supercluster CaloClusters, look up PFCluster ptrs in value map, and match by ref
     for (reco::CaloCluster_iterator caloclus = ecal1->clustersBegin(); caloclus!=ecal1->clustersEnd(); ++caloclus) {
       bool overlap = ClusterClusterMapping::overlap(ecal2, *ecal1,*pfclusterassoc_);
       if (overlap) dist = 0.001;
     }
   }
   else {
     bool overlap=ClusterClusterMapping::overlap(*ecal1,*ecal2);
     
     if(overlap)     {
       dist=LinkByRecHit::computeDist( ecal1->position().eta(),
                                       ecal1->position().phi(), 
                                       ecal2->positionREP().Eta(), 
                                       ecal2->positionREP().Phi() );
     }
   }
   return dist;
 }

double PFBlockAlgo::testTrackAndPS	(	const reco::PFRecTrack &	track,
		const reco::PFCluster &	ps
	)		const

private

tests association between a track and a PS cluster returns distance

Definition at line 872 of file PFBlockAlgo.cc.

References gather_cfg::cout, debug_, reco::PFTrack::extrapolatedPoint(), reco::PFTrajectoryPoint::isValid(), reco::PFCluster::layer(), reco::PFTrajectoryPoint::position(), reco::CaloCluster::position(), PFLayer::PS1, reco::PFTrajectoryPoint::PS1, PFLayer::PS2, reco::PFTrajectoryPoint::PS2, and mathSSE::sqrt().

                                             {
 
 #ifdef PFLOW_DEBUG
   //   cout<<"entering testTrackAndPS"<<endl;
   // resolution of PS cluster dxdx and dydy from strip pitch and length
   double dx=0.;
   double dy=0.;
   
   unsigned layerid =0;
   // PS1: vertical strips  PS2: horizontal strips
   switch (ps.layer()) {
   case PFLayer::PS1:
     layerid = reco::PFTrajectoryPoint::PS1;
     
     // vertical strips in PS1, measure x with pitch precision
     dx = resPSpitch_;
     dy = resPSlength_; 
     break;
   case PFLayer::PS2:
     layerid = reco::PFTrajectoryPoint::PS2;
     // horizontal strips in PS2, measure y with pitch precision
     dy = resPSpitch_;
     dx = resPSlength_;
     break;
   default:
     break;
   }
   const reco::PFTrajectoryPoint& atPS
     = track.extrapolatedPoint( layerid );  
   // did not reach PS, cannot be associated with a cluster.
   if( ! atPS.isValid() ) return -1.;   
   
   double trackx = atPS.position().X();
   double tracky = atPS.position().Y();
   double trackz = atPS.position().Z(); // MDN jan 09
   
   // ps position  x, y
   double psx = ps.position().X();
   double psy = ps.position().Y();
   // MDN Jan 09: check that trackz and psz have the same sign
   double psz = ps.position().Z();
   if( trackz*psz < 0.) return -1.; 
   
   // double chi2 = (psx-trackx)*(psx-trackx)/(dx*dx + trackresolx*trackresolx)
   //  + (psy-tracky)*(psy-tracky)/(dy*dy + trackresoly*trackresoly);
 
   double dist = std::sqrt( (psx-trackx)*(psx-trackx)
              + (psy-tracky)*(psy-tracky));  
   if(debug_) cout<<"testTrackAndPS "<< dist <<" "<<endl;
   if(debug_){
     cout<<" trackx " << trackx 
     <<" tracky " << tracky 
     <<" psx "    <<  psx   
     <<" psy "    << psy    
     << endl;
   }
 #endif
   
   // Return -1. as long as no link by rechit is available
   return -1.;
 }

std::auto_ptr< reco::PFBlockCollection > PFBlockAlgo::transferBlocks ( )

inline

Returns: auto_ptr to collection of blocks

Definition at line 198 of file PFBlockAlgo.h.

References blocks_.

Referenced by PFRootEventManager::particleFlow().

198 {return blocks_;}

PFBlockAlgo::blocks_

std::auto_ptr< reco::PFBlockCollection > blocks_

Definition: PFBlockAlgo.h:311

Friends And Related Function Documentation

std::ostream& operator<<	(	std::ostream &	out,
		const PFBlockAlgo &	a
	)

friend

Definition at line 1285 of file PFBlockAlgo.cc.

                                                               {
   if(! out) return out;
   
   out<<"====== Particle Flow Block Algorithm ======= ";
   out<<endl;
   out<<"number of unassociated elements : "<<a.elements_.size()<<endl;
   out<<endl;
   
   for(PFBlockAlgo::IEC ie = a.elements_.begin(); 
       ie != a.elements_.end(); ++ie) {
     out<<"\t"<<**ie <<endl;
   }
 
   
   //   const PFBlockCollection& blocks = a.blocks();
 
   const std::auto_ptr< reco::PFBlockCollection >& blocks
     = a.blocks(); 
     
   if(!blocks.get() ) {
     out<<"blocks already transfered"<<endl;
   }
   else {
     out<<"number of blocks : "<<blocks->size()<<endl;
     out<<endl;
     
     for(PFBlockAlgo::IBC ib=blocks->begin(); 
     ib != blocks->end(); ++ib) {
       out<<(*ib)<<endl;
     }
   }
 
   return out;
 }