---

PFTopProjector Class Reference

Author:

Colin Bernet

More...

#include <PhysicsTools/PFCandProducer/interface/PFTopProjector.h>

Inheritance diagram for PFTopProjector:

List of all members.

Public Member Functions
virtual void	beginJob (const edm::EventSetup &c)
	PFTopProjector (const edm::ParameterSet &)
virtual void	produce (edm::Event &, const edm::EventSetup &)
	~PFTopProjector ()
Private Member Functions
void	maskAncestors (const reco::CandidatePtrVector &ancestors, std::vector< bool > &masked) const
	ancestors is a RefToBase vector.
template<class T>
void	printAncestors (const reco::CandidatePtrVector &ancestors, const edm::Handle< std::vector< T > > &allPFCandidates) const
template<class T, class U>
void	processCollection (const edm::Handle< std::vector< T > > &handle, const edm::Handle< std::vector< U > > &allPFCandidates, std::vector< bool > &masked, const char *objectName) const
void	ptrToAncestor (reco::CandidatePtr candRef, reco::CandidatePtrVector &ancestors, const edm::ProductID &ancestorsID) const
	fills ancestors with ptrs to the PFCandidates that in one way or another contribute to the candidate pointed to by candPtr
Private Attributes
edm::InputTag	inputTagIsolatedElectrons_
	optional collection of electrons
edm::InputTag	inputTagIsolatedMuons_
	optional collection of muons
edm::InputTag	inputTagPFCandidates_
	ancestor PFCandidates
edm::InputTag	inputTagPFJets_
	optional collection of jets
edm::InputTag	inputTagPFTaus_
	optional collection of taus
edm::InputTag	inputTagPileUpPFCandidates_
	optional collection of PileUpPFCandidates
bool	verbose_
	verbose ?
Static Private Attributes
static const char *	pfCandidatesOutLabel_ = "PFCandidates"
	label for output PFCandidates collection
static const char *	pfJetsOutLabel_ = "PFJets"
	label for output PFJet collection

---

Detailed Description

Author:

Colin Bernet

Date:

february 2008

Definition at line 37 of file PFTopProjector.h.

---

Constructor & Destructor Documentation

PFTopProjector::PFTopProjector

(

const edm::ParameterSet &

iConfig

)

[explicit]

Definition at line 31 of file PFTopProjector.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), edm::ParameterSet::getParameter(), edm::ParameterSet::getUntrackedParameter(), inputTagIsolatedElectrons_, inputTagIsolatedMuons_, inputTagPFCandidates_, inputTagPFJets_, inputTagPFTaus_, inputTagPileUpPFCandidates_, alivecheck_mergeAndRegister::msg, pfCandidatesOutLabel_, pfJetsOutLabel_, and verbose_.

                                                             {
  
  inputTagPFCandidates_ 
    = iConfig.getParameter<InputTag>("PFCandidates");

  inputTagPileUpPFCandidates_ 
    = iConfig.getParameter<InputTag>("PileUpPFCandidates");

  inputTagIsolatedElectrons_ 
    = iConfig.getParameter<InputTag>("IsolatedElectrons");

  inputTagIsolatedMuons_ 
    = iConfig.getParameter<InputTag>("IsolatedMuons");

  inputTagPFJets_ 
    = iConfig.getParameter<InputTag>("PFJets");

  inputTagPFTaus_ 
    = iConfig.getParameter<InputTag>("PFTaus");

  verbose_ = 
    iConfig.getUntrackedParameter<bool>("verbose",false);


  // produces<reco::PFCandidateCollection>("PFCandidates");
  //  produces<reco::PFJetCollection>("PFJets");

  produces<reco::PFCandidateCollection>(pfCandidatesOutLabel_); 
  produces<reco::PFJetCollection>(pfJetsOutLabel_);
  
  if(verbose_) {    
    ostringstream  msg;
    msg<<"input PFCandidateCollection         : "
       <<inputTagPFCandidates_<<endl
       <<"input PileUpPFCandidateCollection   : "
       <<inputTagPileUpPFCandidates_<<endl 
       <<"input IsolatedPFCandidates of type muon : "
       <<inputTagIsolatedMuons_<<endl  
       <<"input IsolatedPFCandidates of type electron : "
       <<inputTagIsolatedElectrons_<<endl  
       <<"input PFJetCollection : "
       <<inputTagPFJets_<<endl
       <<"input PFTauCollection : "
       <<inputTagPFTaus_<<endl;  
    cout<<msg.str()<<endl;
    //     LogInfo("PFTopProjector")<<msg.str();
  }
}

PFTopProjector::~PFTopProjector

(

)

Definition at line 82 of file PFTopProjector.cc.

{ }

---

Member Function Documentation

void PFTopProjector::beginJob

(

const edm::EventSetup &

c

)

[virtual]

Reimplemented from edm::EDProducer.

Definition at line 86 of file PFTopProjector.cc.

{ }

void PFTopProjector::maskAncestors	(	const reco::CandidatePtrVector &	ancestors,
		std::vector< bool > &	masked
	)			const [private]

ancestors is a RefToBase vector.

For each object in this vector get the index and set the corresponding slot to true in the masked vector

Definition at line 280 of file PFTopProjector.cc.

References i, index, and edm::PtrVectorBase::size().

Referenced by processCollection().

                                                                    {
  
  for(unsigned i=0; i<ancestors.size(); i++) {
    unsigned index = ancestors[i].key();
    assert( index<masked.size() );
    
    //     if(verbose_) {
    //       ProductID id = ancestors[i].id();
    //       cout<<"\tmasking "<<index<<", ancestor "<<id<<"/"<<index<<endl;
    //     }
    masked[index] = true;
  }
}

template<class T>

void PFTopProjector::printAncestors	(	const reco::CandidatePtrVector &	ancestors,
		const edm::Handle< std::vector< T > > &	allPFCandidates
	)			const [inline, private]

Definition at line 150 of file PFTopProjector.h.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), i, edm::PtrVectorBase::id(), index, and edm::PtrVectorBase::size().

Referenced by processCollection().

                                                                                               {
  
  std::vector<T> pfs = *allPFCandidates;

  for(unsigned i=0; i<ancestors.size(); i++) {

    edm::ProductID id = ancestors[i].id();
    assert( id == allPFCandidates.id() );
 
    unsigned index = ancestors[i].key();
    assert( index < pfs.size() );
    
    std::cout<<"\t\t"<<pfs[index]<<std::endl;
  }
}

template<class T, class U>

void PFTopProjector::processCollection	(	const edm::Handle< std::vector< T > > &	handle,
		const edm::Handle< std::vector< U > > &	allPFCandidates,
		std::vector< bool > &	masked,
		const char *	objectName
	)			const [inline, private]

Definition at line 107 of file PFTopProjector.h.

References collection, GenMuonPlsPt100GeV_cfg::cout, lat::endl(), i, edm::PtrVectorBase::id(), maskAncestors(), printAncestors(), ptr, ptrToAncestor(), and verbose_.

Referenced by produce().

                                                                      {

  if( handle.isValid() && allPFCandidates.isValid() ) {
    const std::vector<T>& collection = *handle;
    
    if(verbose_) 
      std::cout<<" Collection: "<<objectName
               <<" size = "<<collection.size()<<std::endl;
    
    for(unsigned i=0; i<collection.size(); i++) {
      
      
      edm::Ptr<T>   ptr( handle, i);
      reco::CandidatePtr basePtr( ptr );
 
      
      reco::CandidatePtrVector ancestors;
      ptrToAncestor( basePtr,
                     ancestors,
                     allPFCandidates.id() );
      
      if(verbose_) {
/*      std::cout<<"\t"<<objectName<<" "<<i */
/*               <<" pt,eta,phi = " */
/*               <<basePtr->pt()<<"," */
/*               <<basePtr->eta()<<"," */
/*               <<basePtr->phi()<<std::endl; */
        
        std::cout<<"\t"<<collection[i]<<std::endl;
        printAncestors( ancestors, allPFCandidates );
      }
  
      maskAncestors( ancestors, masked );
    }
  }

}

void PFTopProjector::produce	(	edm::Event &	iEvent,
		const edm::EventSetup &	iSetup
	)			[virtual]

Implements edm::EDProducer.

Definition at line 89 of file PFTopProjector.cc.

References GenMuonPlsPt100GeV_cfg::cout, lat::endl(), err, edm::EventID::event(), Exception, pfpat::fetchCollection(), i, edm::Event::id(), edm::Handle< T >::id(), inputTagIsolatedElectrons_, inputTagIsolatedMuons_, inputTagPFCandidates_, inputTagPFJets_, inputTagPFTaus_, inputTagPileUpPFCandidates_, edm::Handle< T >::isValid(), pfCandidatesOutLabel_, pfJetBenchmark_cfi::pfJets, pfJetsOutLabel_, processCollection(), edm::Event::put(), and verbose_.

                                                       {
  
  
  if( verbose_)
    cout<<"Event -------------------- "<<iEvent.id().event()<<endl;
  
  // get the various collections

  // this collection is the collection of PFCandidates to 
  // be masked by the other ones
  Handle<PFCandidateCollection> pfCandidates;
  pfpat::fetchCollection(pfCandidates, 
                         inputTagPFCandidates_, 
                         iEvent );

  if( !pfCandidates.isValid() ) {
    std::ostringstream  err;
    err<<"The collection of input PFCandidates must be supplied."<<endl
       <<"It is now set to : "<<inputTagPFCandidates_<<endl;
    edm::LogError("PFPAT")<<err.str();
    throw cms::Exception( "MissingProduct", err.str());
  }

  edm::ProductID pfCandidatesID = pfCandidates.id();

  // for each object in the following collections, the 
  // PFTopProjector will find and mask all PFCandidates in pfCandidates
  // involved in the construction of this object
  
  // empty input tags can be specified, and will be recognized 
  // by the fetchCollection function
  
  Handle<PileUpPFCandidateCollection> pfPileUpCandidates;
  pfpat::fetchCollection(pfPileUpCandidates, 
                         inputTagPileUpPFCandidates_, 
                         iEvent );

  Handle<PFCandidateCollection> pfIsolatedElectrons;
  pfpat::fetchCollection(pfIsolatedElectrons, 
                         inputTagIsolatedElectrons_, 
                         iEvent );

  Handle<PFCandidateCollection> pfIsolatedMuons;
  pfpat::fetchCollection(pfIsolatedMuons, 
                         inputTagIsolatedMuons_, 
                         iEvent );

  Handle<PFJetCollection> pfJets;
  pfpat::fetchCollection(pfJets, 
                         inputTagPFJets_, 
                         iEvent );

  Handle<PFTauCollection> pfTaus;
  pfpat::fetchCollection(pfTaus, 
                         inputTagPFTaus_, 
                         iEvent );

  if(verbose_) {
    cout<<"Top projector: event "<<iEvent.id().event()<<endl;
    cout<<"product Ids --------------------"<<endl;
    cout<<"PF       :  "<<pfCandidates.id()<<endl
        <<"PFPU     :  "<<pfPileUpCandidates.id()<<endl
        <<"PFIso e  :  "<<pfIsolatedElectrons.id()<<endl
        <<"PFIso mu :  "<<pfIsolatedMuons.id()<<endl
        <<"PFJets   :  "<<pfJets.id()<<endl
        <<"PFTaus   :  "<<pfTaus.id()<<endl;
  }

  // output PFCandidate collection
  // will contain a copy of each PFCandidate in pfCandidates
  // that remains unmasked. 
  auto_ptr< reco::PFCandidateCollection > 
    pPFCandidateOutput( new reco::PFCandidateCollection ); 
  
  auto_ptr< reco::PFJetCollection > 
    pPFJetOutput( new reco::PFJetCollection ); 
  
  // mask for each PFCandidate.
  // at the beginning, all PFCandidates are unmasked.
  vector<bool> masked( pfCandidates->size(), false);
    
  assert( pfCandidates.isValid() );


  processCollection( pfTaus, pfCandidates, masked, 
                     "PFTau");
  processCollection( pfJets, pfCandidates, masked, 
                     "PFJet");
  processCollection( pfPileUpCandidates, pfCandidates, masked, 
                     "PileUpParticle");
  processCollection( pfIsolatedElectrons, pfCandidates, masked, 
                     "IsoElectron");
  processCollection( pfIsolatedMuons, pfCandidates, masked, 
                     "IsoMuon");


  const PFCandidateCollection& inCands = *pfCandidates;

  if(verbose_) 
    cout<<" Remaining PFCandidates ------ "<<endl;
  
  for(unsigned i=0; i<inCands.size(); i++) {
    
    if(masked[i]) {
      if(verbose_)
        cout<<"X "<<i<<" "<<inCands[i]<<endl;
      continue;
    }
    else {
      if(verbose_)
        cout<<"O "<<i<<" "<<inCands[i]<<endl;
      PFCandidatePtr motherPtr( pfCandidates, i );
      pPFCandidateOutput->push_back( inCands[i] );
      pPFCandidateOutput->back().setSourcePtr(motherPtr);
    }
  }

  iEvent.put( pPFCandidateOutput, pfCandidatesOutLabel_ );


  // now mask the jets with the taus (if the jet collection has been provided)
  
  if( pfJets.isValid() ) {
    vector<bool> maskedJets( pfJets->size(), false);

    
    processCollection( pfTaus, pfJets, maskedJets, 
                       "PFTau masking PFJets");


    const PFJetCollection& inJets = *pfJets;
    
    if(verbose_) 
      cout<<" Remaining PFJets ------ "<<endl;
    
    for(unsigned i=0; i<inJets.size(); i++) {
      
      if(maskedJets[i]) {
        if(verbose_)
          cout<<"X "<<i<<" "<<inJets[i]<<endl;
        continue;
      }
      else {
        if(verbose_)
          cout<<"O "<<i<<" "<<inJets[i]<<endl;
        pPFJetOutput->push_back( inJets[i] );
      }
    }
  }



  iEvent.put( pPFJetOutput, pfJetsOutLabel_);
  
  //   LogDebug("PFTopProjector")<<"STOP event: "<<iEvent.id().event()
  //                     <<" in run "<<iEvent.id().run()<<endl;
}

void PFTopProjector::ptrToAncestor	(	reco::CandidatePtr	candRef,
		reco::CandidatePtrVector &	ancestors,
		const edm::ProductID &	ancestorsID
	)			const [private]

fills ancestors with ptrs to the PFCandidates that in one way or another contribute to the candidate pointed to by candPtr

Definition at line 250 of file PFTopProjector.cc.

References i, edm::Ptr< T >::id(), and edm::PtrVector< T >::push_back().

Referenced by processCollection().

                                                                         {

  
 
  unsigned nSources = candPtr->numberOfSourceCandidatePtrs();

//   cout<<"going down from "<<candPtr.id()
//       <<"/"<<candPtr.key()<<" #mothers "<<nSources
//       <<" ancestor id "<<ancestorsID<<endl;
  
  for(unsigned i=0; i<nSources; i++) {
    
    CandidatePtr mother = candPtr->sourceCandidatePtr(i);
//     cout<<"  mother id "<<mother.id()<<endl;
    
    if(  mother.id() != ancestorsID ) {
      // the mother is not yet at lowest level
      ptrToAncestor( mother, ancestors, ancestorsID);
    }
    else {
      // adding mother to the list of ancestors
      ancestors.push_back( mother ); 
    }
  }
}

---

Member Data Documentation

edm::InputTag PFTopProjector::inputTagIsolatedElectrons_ [private]

optional collection of electrons

Definition at line 83 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

edm::InputTag PFTopProjector::inputTagIsolatedMuons_ [private]

optional collection of muons

Definition at line 86 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

edm::InputTag PFTopProjector::inputTagPFCandidates_ [private]

ancestor PFCandidates

Definition at line 77 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

edm::InputTag PFTopProjector::inputTagPFJets_ [private]

optional collection of jets

Definition at line 89 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

edm::InputTag PFTopProjector::inputTagPFTaus_ [private]

optional collection of taus

Definition at line 92 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

edm::InputTag PFTopProjector::inputTagPileUpPFCandidates_ [private]

optional collection of PileUpPFCandidates

Definition at line 80 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

const char * PFTopProjector::pfCandidatesOutLabel_ = "PFCandidates" [static, private]

label for output PFCandidates collection

Definition at line 101 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

const char * PFTopProjector::pfJetsOutLabel_ = "PFJets" [static, private]

label for output PFJet collection

Definition at line 98 of file PFTopProjector.h.

Referenced by PFTopProjector(), and produce().

bool PFTopProjector::verbose_ [private]

verbose ?

Definition at line 95 of file PFTopProjector.h.

Referenced by PFTopProjector(), processCollection(), and produce().

---

The documentation for this class was generated from the following files:

• PhysicsTools/PFCandProducer/interface/PFTopProjector.h
• PhysicsTools/PFCandProducer/src/PFTopProjector.cc

---