Inheritance diagram for JetPartonMatcher:

Classes
struct	WorkingVariables

Public Member Functions
	JetPartonMatcher (const edm::ParameterSet &)

	~JetPartonMatcher ()

Public Member Functions inherited from edm::global::EDProducer<>
	EDProducer ()=default

Public Member Functions inherited from edm::global::EDProducerBase
	EDProducerBase ()

ModuleDescription const &	moduleDescription () const

virtual	~EDProducerBase ()

Public Member Functions inherited from edm::ProducerBase
void	callWhenNewProductsRegistered (std::function< void(BranchDescription const &)> const &func)

std::vector< edm::ProductResolverIndex > const &	indiciesForPutProducts (BranchType iBranchType) const

	ProducerBase ()

void	registerProducts (ProducerBase , ProductRegistry , ModuleDescription const &)

std::function< void(BranchDescription const &)>	registrationCallback () const
	used by the fwk to register list of products More...

void	resolvePutIndicies (BranchType iBranchType, std::unordered_multimap< std::string, edm::ProductResolverIndex > const &iIndicies, std::string const &moduleLabel)

virtual	~ProducerBase () noexcept(false)

Public Member Functions inherited from edm::EDConsumerBase
std::vector< ConsumesInfo >	consumesInfo () const

	EDConsumerBase ()

	EDConsumerBase (EDConsumerBase const &)=delete

	EDConsumerBase (EDConsumerBase &&)=default

ProductResolverIndexAndSkipBit	indexFrom (EDGetToken, BranchType, TypeID const &) const

void	itemsMayGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

void	itemsToGet (BranchType, std::vector< ProductResolverIndexAndSkipBit > &) const

std::vector< ProductResolverIndexAndSkipBit > const &	itemsToGetFrom (BranchType iType) const

void	labelsForToken (EDGetToken iToken, Labels &oLabels) const

void	modulesWhoseProductsAreConsumed (std::vector< ModuleDescription const * > &modules, ProductRegistry const &preg, std::map< std::string, ModuleDescription const * > const &labelsToDesc, std::string const &processName) const

EDConsumerBase const &	operator= (EDConsumerBase const &)=delete

EDConsumerBase &	operator= (EDConsumerBase &&)=default

bool	registeredToConsume (ProductResolverIndex, bool, BranchType) const

bool	registeredToConsumeMany (TypeID const &, BranchType) const

void	updateLookup (BranchType iBranchType, ProductResolverIndexHelper const &, bool iPrefetchMayGet)

virtual	~EDConsumerBase () noexcept(false)

Private Member Functions
int	fillAlgoritDefinition (const Jet &, WorkingVariables &) const

int	fillPhysicsDefinition (const Jet &, WorkingVariables &) const

virtual void	produce (edm::StreamID, edm::Event &, const edm::EventSetup &) const override

Private Attributes
double	coneSizeToAssociate

bool	doPriority

edm::EDGetTokenT< edm::View< reco::Jet > >	m_jetsSrcToken

edm::EDGetTokenT< GenParticleRefVector >	m_ParticleSrcToken

bool	physDefinition

vector< int >	priorityList

Additional Inherited Members
Public Types inherited from edm::global::EDProducerBase
typedef EDProducerBase	ModuleType

Public Types inherited from edm::ProducerBase
typedef ProductRegistryHelper::TypeLabelList	TypeLabelList

Public Types inherited from edm::EDConsumerBase
typedef ProductLabels	Labels

Static Public Member Functions inherited from edm::global::EDProducerBase
static const std::string &	baseType ()

static void	fillDescriptions (ConfigurationDescriptions &descriptions)

static void	prevalidate (ConfigurationDescriptions &descriptions)

Protected Member Functions inherited from edm::EDConsumerBase
template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	consumes (edm::InputTag const &tag)

EDGetToken	consumes (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	consumes (TypeToGet const &id, edm::InputTag const &tag)

ConsumesCollector	consumesCollector ()
	Use a ConsumesCollector to gather consumes information from helper functions. More...

template<typename ProductType , BranchType B = InEvent>
void	consumesMany ()

void	consumesMany (const TypeToGet &id)

template<BranchType B>
void	consumesMany (const TypeToGet &id)

template<typename ProductType , BranchType B = InEvent>
EDGetTokenT< ProductType >	mayConsume (edm::InputTag const &tag)

EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

template<BranchType B>
EDGetToken	mayConsume (const TypeToGet &id, edm::InputTag const &tag)

Detailed Description

Definition at line 94 of file JetPartonMatcher.cc.

Constructor & Destructor Documentation

JetPartonMatcher::JetPartonMatcher ( const edm::ParameterSet & iConfig )

Definition at line 125 of file JetPartonMatcher.cc.

References AK4CaloJetsMCFlavour_cff::coneSizeToAssociate, edm::ParameterSet::exists(), and edm::ParameterSet::getParameter().

 {
     produces<JetMatchedPartonsCollection>();
     m_jetsSrcToken           = consumes<edm::View <reco::Jet> >(iConfig.getParameter<edm::InputTag>("jets"));
     m_ParticleSrcToken       = consumes<GenParticleRefVector>(iConfig.getParameter<edm::InputTag>("partons"));
     coneSizeToAssociate = iConfig.getParameter<double>("coneSizeToAssociate");
     if ( iConfig.exists("doPriority") ) {
       doPriority = iConfig.getParameter<bool>("doPriority");
       priorityList = iConfig.getParameter<vector<int> >("priorityList");
     } else {
       doPriority = false;
       priorityList.clear();
     }
 }

JetPartonMatcher::~JetPartonMatcher ( )

Definition at line 142 of file JetPartonMatcher.cc.

143 {

144 }

Member Function Documentation

int JetPartonMatcher::fillAlgoritDefinition	(	const Jet &	theJet,
		WorkingVariables &	wv
	)		const

private

Definition at line 215 of file JetPartonMatcher.cc.

References funct::abs(), AK4CaloJetsMCFlavour_cff::coneSizeToAssociate, HiRegitMuonDetachedTripletStep_cff::DeltaR, spr::find(), funct::m, MuonErrorMatrixAnalyzer_cfi::maxPt, reco::Candidate::p4(), reco::LeafCandidate::p4(), JetPartonMatcher::WorkingVariables::particles, reco::Candidate::pdgId(), reco::Candidate::pt(), reco::Candidate::status(), JetPartonMatcher::WorkingVariables::theHardest, JetPartonMatcher::WorkingVariables::theHeaviest, and JetPartonMatcher::WorkingVariables::theNearest2.

                                                                                            {
 
   int tempParticle = -1;
   int tempPartonHighestPt = -1;
   int tempNearest = -1;
   float maxPt = 0;
   float minDr = 1000;
   bool foundPriority = false;
 
   // Loop over the particles in question until we find a priority
   // "hit", or if we find none in the priority list (or we don't want
   // to consider priority), then loop through all particles and fill
   // standard definition.
 
   //
   // Matching:
   //
   // 1) First try to match by hand. The "priority list" is given
   //    by the user. The algorithm finds any particles in that
   //    "priority list" that are within the specified cone size.
   //    If it finds one, it counts the object as associating to that
   //    particle.
   //    NOTE! The objects in the priority list are given in order
   //    of priority. So if the user specifies:
   //    6, 24, 21,
   //    then it will first look for top quarks, then W bosons, then gluons.
   // 2) If no priority items are found, do the default "standard"
   //    matching.
   for( size_t m = 0; m != wv.particles->size() && !foundPriority; ++ m ) {
     const Candidate & aParton = *(wv.particles->at(m).get());
 
     // "Priority" behavoir:
     // Associate to the first particle found in the priority list, regardless
     // of delta R.
     if ( doPriority ) {
       int ipdgId = aParton.pdgId();
       vector<int>::const_iterator ipriority = find( priorityList.begin(), priorityList.end(), abs(ipdgId) );
       // Check to see if this particle is in our priority list
       if ( ipriority != priorityList.end() ) {
     // This particle is on our priority list. If it matches,
     // we break, since we take in order of priority, not deltaR
     double dist = DeltaR( theJet.p4(), aParton.p4() );
     if( dist <= coneSizeToAssociate ) {
       tempParticle = m;
       foundPriority = true;
     }
       }
     }
     // Here we do not want to do priority matching. Ensure the "foundPriority" swtich
     // is turned off.
     else {
       foundPriority = false;
     }
 
     // Default behavior:
     // Look for partons before the color string to associate.
     // Do this if we don't want to do priority matching, or if
     // we didn't find a priority object.
 
     if( !foundPriority && aParton.status() != 3 ) {
       double dist = DeltaR( theJet.p4(), aParton.p4() );
       if( dist <= coneSizeToAssociate ) {
     if( dist < minDr ) {
       minDr = dist;
       tempNearest = m;
     }
     if( tempParticle == -1 && ( abs( aParton.pdgId() ) == 4 )  ) tempParticle = m;
     if(                         abs( aParton.pdgId() ) == 5    ) tempParticle = m;
     if( aParton.pt() > maxPt ) {
       maxPt = aParton.pt();
       tempPartonHighestPt = m;
     }
       }
     }
   }
 
   if ( foundPriority ) {
     wv.theHeaviest = tempParticle; // The priority-matched particle
     wv.theHardest  = -1;  //  set the rest to -1
     wv.theNearest2 = -1;  // "                  "
   } else {
     wv.theHeaviest = tempParticle;
     wv.theHardest  = tempPartonHighestPt;
     wv.theNearest2 = tempNearest;
     if ( tempParticle == -1 ) tempParticle = tempPartonHighestPt;
   }
   return tempParticle;
 }

int JetPartonMatcher::fillPhysicsDefinition	(	const Jet &	theJet,
		WorkingVariables &	wv
	)		const

private

Definition at line 326 of file JetPartonMatcher.cc.

References funct::abs(), AK4CaloJetsMCFlavour_cff::coneSizeToAssociate, DEFINE_FWK_MODULE, HiRegitMuonDetachedTripletStep_cff::DeltaR, spr::find(), GenHFHadronMatcher_cff::flavour, funct::m, reco::Candidate::numberOfDaughters(), reco::Candidate::p4(), reco::LeafCandidate::p4(), JetPartonMatcher::WorkingVariables::particles, cosmictrackingParticleSelector_cfi::pdgId, reco::Candidate::pdgId(), reco::Candidate::status(), JetPartonMatcher::WorkingVariables::theHardest, JetPartonMatcher::WorkingVariables::theHeaviest, JetPartonMatcher::WorkingVariables::theNearest2, and JetPartonMatcher::WorkingVariables::theNearest3.

                                                                                            {
 
   float TheBiggerConeSize = 0.7; // In HepMC it's 0.3 --> it's a mistake: value has to be 0.7
   int tempParticle = -1;
   int nInTheCone = 0;
   int tempNearest = -1;
   float minDr = 1000;
   bool foundPriority = false;
 
   vector<const reco::Candidate *> theContaminations;
   theContaminations.clear();
 
   // Loop over the particles in question until we find a priority
   // "hit", or if we find none in the priority list (or we don't want
   // to consider priority), then loop through all particles and fill
   // standard definition.
 
   //
   // Matching:
   //
   // 1) First try to match by hand. The "priority list" is given
   //    by the user. The algorithm finds any particles in that
   //    "priority list" that are within the specified cone size.
   //    If it finds one, it counts the object as associating to that
   //    particle.
   //    NOTE! The objects in the priority list are given in order
   //    of priority. So if the user specifies:
   //    6, 24, 21,
   //    then it will first look for top quarks, then W bosons, then gluons.
   // 2) If no priority items are found, do the default "standard"
   //    matching.
   for( size_t m = 0; m != wv.particles->size() && !foundPriority; ++ m ) {
 
     const Candidate & aParticle = *(wv.particles->at(m).get());
 
     // "Priority" behavoir:
     // Associate to the first particle found in the priority list, regardless
     // of delta R.
     if ( doPriority ) {
       int ipdgId = aParticle.pdgId();
       vector<int>::const_iterator ipriority = find( priorityList.begin(), priorityList.end(), abs(ipdgId) );
       // Check to see if this particle is in our priority list
       if ( ipriority != priorityList.end() ) {
     // This particle is on our priority list. If it matches,
     // we break, since we take in order of priority, not deltaR
     double dist = DeltaR( theJet.p4(), aParticle.p4() );
     if( dist <= coneSizeToAssociate ) {
       tempParticle = m;
       foundPriority = true;
     }
       }
     }
     // Here we do not want to do priority matching. Ensure the "foundPriority" swtich
     // is turned off.
     else {
       foundPriority = false;
     }
 
     // Default behavior:
     if ( !foundPriority ) {
 
       // skipping all particle but udscbg (is this correct/enough?!?!)
       bool isAParton = false;
       int flavour = abs(aParticle.pdgId());
       if(flavour == 1 ||
      flavour == 2 ||
      flavour == 3 ||
      flavour == 4 ||
      flavour == 5 ||
      flavour == 21 ) isAParton = true;
       if(!isAParton) continue;
       double dist = DeltaR( theJet.p4(), aParticle.p4() );
       if( aParticle.status() == 3 && dist < minDr ) {
     minDr = dist;
     tempNearest = m;
       }
       if( aParticle.status() == 3 && dist <= coneSizeToAssociate ) {
     //cout << "particle in small cone=" << aParticle.pdgId() << endl;
     tempParticle = m;
     nInTheCone++;
       }
       // Look for heavy partons in TheBiggerConeSize now
       if( aParticle.numberOfDaughters() > 0 && aParticle.status() != 3 ) {
     if( flavour ==  1 ||
         flavour ==  2 ||
         flavour ==  3 ||
         flavour == 21 ) continue;
     if( dist < TheBiggerConeSize ) theContaminations.push_back( &aParticle );
       }
     }
   }
 
 
   // Here's the default behavior for assignment if there is no priority.
   if ( !foundPriority ) {
     wv.theNearest3 = tempNearest;
 
     if(nInTheCone != 1) return -1; // rejected --> only one initialParton requested
     if(theContaminations.size() == 0 ) return tempParticle; //no contamination
     int initialPartonFlavour = abs( (wv.particles->at(tempParticle).get()) ->pdgId() );
 
     vector<const Candidate *>::const_iterator itCont = theContaminations.begin();
     for( ; itCont != theContaminations.end(); itCont++ ) {
       int contaminatingFlavour = abs( (*itCont)->pdgId() );
       if( (*itCont)->numberOfMothers()>0 && (*itCont)->mother(0) == wv.particles->at(tempParticle).get() ) continue; // mother is the initialParton --> OK
       if( initialPartonFlavour == 4 ) {
     if( contaminatingFlavour == 4 ) continue; // keep association --> the initialParton is a c --> the contaminated parton is a c
     tempParticle = -1; // all the other cases reject!
     return tempParticle;
       }
     }
   }
   // If there is priority, then just set the heaviest to priority, the rest are -1.
   else {
     wv.theHeaviest = tempParticle; // Set the heaviest to tempParticle
     wv.theNearest2 = -1; //  Set the rest to -1
     wv.theNearest3 = -1; // "                  "
     wv.theHardest = -1;  // "                  "
   }
 
   return tempParticle;
 }

void JetPartonMatcher::produce	(	edm::StreamID	,
		edm::Event &	iEvent,
		const edm::EventSetup &	iEs
	)		const

overrideprivatevirtual

Definition at line 148 of file JetPartonMatcher.cc.

References reco::LeafCandidate::eta(), edm::Event::getByToken(), funct::m, eostools::move(), JetPartonMatcher::WorkingVariables::particles, reco::LeafCandidate::pdgId(), reco::LeafCandidate::phi(), reco::LeafCandidate::pt(), edm::Event::put(), reco::LeafCandidate::status(), JetPartonMatcher::WorkingVariables::theHeaviest, JetPartonMatcher::WorkingVariables::theNearest2, and JetPartonMatcher::WorkingVariables::theNearest3.

 {
   WorkingVariables wv;
   edm::Handle <edm::View <reco::Jet> > jets_h;
   iEvent.getByToken(m_jetsSrcToken,     jets_h    );
   iEvent.getByToken(m_ParticleSrcToken, wv.particles );
 
   edm::LogVerbatim("JetPartonMatcher") << "=== Partons size:" << wv.particles->size();
 
   for( size_t m = 0; m != wv.particles->size(); ++ m ) {
     const GenParticle & aParton = *(wv.particles->at(m).get());
     edm::LogVerbatim("JetPartonMatcher") <<  aParton.status() << " " <<
                                              aParton.pdgId()  << " " <<
                                              aParton.pt()     << " " <<
                                              aParton.eta()    << " " <<
                                              aParton.phi()    << endl;
   }
 
   auto jetMatchedPartons = std::make_unique<JetMatchedPartonsCollection>(reco::JetRefBaseProd(jets_h));
 
   for (size_t j = 0; j < jets_h->size(); j++) {
 
     const int theMappedPartonAlg = fillAlgoritDefinition( (*jets_h)[j], wv );
     const int theMappedPartonPhy = fillPhysicsDefinition( (*jets_h)[j], wv );
 
     GenParticleRef pHV;
     GenParticleRef pN2;
     GenParticleRef pN3;
     GenParticleRef pPH;
     GenParticleRef pAL;
 
     if(wv.theHeaviest>=0)        pHV = wv.particles->at( wv.theHeaviest        );
     if(wv.theNearest2>=0)        pN2 = wv.particles->at( wv.theNearest2        );
     if(wv.theNearest3>=0)        pN3 = wv.particles->at( wv.theNearest3        );
     if(theMappedPartonPhy>=0) pPH = wv.particles->at( theMappedPartonPhy );
     if(theMappedPartonAlg>=0) pAL = wv.particles->at( theMappedPartonAlg );
 
     (*jetMatchedPartons)[jets_h->refAt(j)]=MatchedPartons(pHV,pN2,pN3,pPH,pAL);
   }
 
   iEvent.put(std::move(jetMatchedPartons));
 
 }

Member Data Documentation

double JetPartonMatcher::coneSizeToAssociate

private

Definition at line 116 of file JetPartonMatcher.cc.

bool JetPartonMatcher::doPriority

private

Definition at line 118 of file JetPartonMatcher.cc.

edm::EDGetTokenT<edm::View <reco::Jet> > JetPartonMatcher::m_jetsSrcToken

private

Definition at line 114 of file JetPartonMatcher.cc.

edm::EDGetTokenT<GenParticleRefVector> JetPartonMatcher::m_ParticleSrcToken

private

Definition at line 115 of file JetPartonMatcher.cc.

bool JetPartonMatcher::physDefinition

private

Definition at line 117 of file JetPartonMatcher.cc.

vector<int> JetPartonMatcher::priorityList

private

Definition at line 119 of file JetPartonMatcher.cc.

Classes

Public Member Functions

Private Member Functions

Private Attributes

Additional Inherited Members

Detailed Description

Constructor & Destructor Documentation

Member Function Documentation

Member Data Documentation