GenJetClosestMatchSelectorDefinition Struct Reference

#include <GenJetClosestMatchSelectorDefinition.h>

Public Types
typedef reco::GenJetCollection	collection
typedef container::const_iterator	const_iterator
typedef std::vector< reco::GenJet * >	container
typedef edm::Handle< collection >	HandleToCollection

Public Member Functions
const_iterator	begin () const
const_iterator	end () const
	GenJetClosestMatchSelectorDefinition (const edm::ParameterSet &cfg, edm::ConsumesCollector &&iC)
void	select (const HandleToCollection &hc, const edm::Event &e, const edm::EventSetup &s)
size_t	size () const

Private Attributes
edm::EDGetTokenT< edm::View< reco::Candidate > >	matchTo_
container	selected_

Detailed Description

Definition at line 12 of file GenJetClosestMatchSelectorDefinition.h.

Member Typedef Documentation

typedef reco::GenJetCollection GenJetClosestMatchSelectorDefinition::collection

Definition at line 15 of file GenJetClosestMatchSelectorDefinition.h.

typedef container::const_iterator GenJetClosestMatchSelectorDefinition::const_iterator

Definition at line 18 of file GenJetClosestMatchSelectorDefinition.h.

typedef std::vector< reco::GenJet *> GenJetClosestMatchSelectorDefinition::container

Definition at line 17 of file GenJetClosestMatchSelectorDefinition.h.

typedef edm::Handle< collection > GenJetClosestMatchSelectorDefinition::HandleToCollection

Definition at line 16 of file GenJetClosestMatchSelectorDefinition.h.

Constructor & Destructor Documentation

GenJetClosestMatchSelectorDefinition::GenJetClosestMatchSelectorDefinition ( const edm::ParameterSet &  cfg, edm::ConsumesCollector &&  iC  )

inline

Definition at line 20 of file GenJetClosestMatchSelectorDefinition.h.

References edm::ParameterSet::getParameter(), and matchTo_.

                                                                                                   {

    matchTo_ = iC.consumes< edm::View<reco::Candidate> >(cfg.getParameter< edm::InputTag >( "MatchTo" ));
  }

Member Function Documentation

const_iterator GenJetClosestMatchSelectorDefinition::begin ( void  ) const

inline

Definition at line 25 of file GenJetClosestMatchSelectorDefinition.h.

References selected_.

{ return selected_.begin(); }

const_iterator GenJetClosestMatchSelectorDefinition::end ( void  ) const

inline

Definition at line 27 of file GenJetClosestMatchSelectorDefinition.h.

References selected_.

{ return selected_.end(); }

void GenJetClosestMatchSelectorDefinition::select ( const HandleToCollection &  hc, const edm::Event &  e, const edm::EventSetup &  s  )

inline

Definition at line 29 of file GenJetClosestMatchSelectorDefinition.h.

References reco::deltaR2(), edm::Event::getByToken(), crabWrapper::key, matchTo_, and selected_.

  {

    selected_.clear();

    edm::Handle< edm::View<reco::Candidate> > matchCandidates;
    e.getByToken( matchTo_, matchCandidates);


    unsigned key=0;

    //    std::cout<<"number of candidates "<<matchCandidates->size()<<std::endl;

    typedef edm::View<reco::Candidate>::const_iterator IC;
    for( IC ic = matchCandidates->begin();
     ic!= matchCandidates->end(); ++ic ) {

      double eta2 = ic->eta();
      double phi2 = ic->phi();

      //      std::cout<<"cand "<<eta2<<" "<<phi2<<std::endl;


      // look for the closest gen jet
      double deltaR2Min = 9999;
      collection::const_iterator closest = hc->end();
      for( collection::const_iterator genjet = hc->begin();
       genjet != hc->end();
       ++genjet, ++key) {

    reco::GenJetRef genJetRef(hc, key);

    // is it matched?

    double eta1 = genjet->eta();
    double phi1 = genjet->phi();


    double deltaR2 = reco::deltaR2(eta1, phi1, eta2, phi2);

    // std::cout<<"  genjet "<<eta1<<" "<<phi1<<" "<<deltaR2<<std::endl;

    // cut should be a parameter
    if( deltaR2<deltaR2Min ) {
      deltaR2Min = deltaR2;
      closest = genjet;
    }
      }

      if(deltaR2Min<0.01 ) {
    // std::cout<<deltaR2Min<<std::endl;
    selected_.push_back( new reco::GenJet(*closest) );
      }
    } // end collection iteration

    // std::cout<<selected_.size()<<std::endl;
  } // end select()

size_t GenJetClosestMatchSelectorDefinition::size ( void  ) const

inline

Definition at line 89 of file GenJetClosestMatchSelectorDefinition.h.

References selected_.

Referenced by ntupleDataFormat._Collection::__iter__(), and ntupleDataFormat._Collection::__len__().

{ return selected_.size(); }

Member Data Documentation

edm::EDGetTokenT<edm::View<reco::Candidate> > GenJetClosestMatchSelectorDefinition::matchTo_

private

Definition at line 93 of file GenJetClosestMatchSelectorDefinition.h.

Referenced by GenJetClosestMatchSelectorDefinition(), and select().

container GenJetClosestMatchSelectorDefinition::selected_

private

Definition at line 92 of file GenJetClosestMatchSelectorDefinition.h.

Referenced by begin(), end(), select(), and size().