#include <PseudoJetTranslator.h>

Collaboration diagram for PseudoJetTranslator:

Public Member Functions
	PseudoJetTranslator (bool saveArea, bool saveArea4Vec)

xAOD::Jet &	translate (const fastjet::PseudoJet &pj, const PseudoJetContainer &pjCont, xAOD::JetContainer &jetCont, const xAOD::Vertex *originVertex=nullptr) const

xAOD::Jet &	translate (const fastjet::PseudoJet &pj, const PseudoJetContainer &pjCont, xAOD::JetContainer &jetCont, const xAOD::Jet &parent, const xAOD::Vertex *originVertex=nullptr) const

Protected Attributes
bool	m_saveArea = true

bool	m_saveArea4Vec = true

Detailed Description

Definition at line 12 of file PseudoJetTranslator.h.

Constructor & Destructor Documentation

◆ PseudoJetTranslator()

PseudoJetTranslator::PseudoJetTranslator	(	bool	saveArea,
		bool	saveArea4Vec
	)

inline

Definition at line 16 of file PseudoJetTranslator.h.

16 : m_saveArea(saveArea), m_saveArea4Vec(saveArea4Vec) {}

Member Function Documentation

◆ translate() [1/2]

xAOD::Jet & PseudoJetTranslator::translate	(	const fastjet::PseudoJet &	pj,
		const PseudoJetContainer &	pjCont,
		xAOD::JetContainer &	jetCont,
		const xAOD::Jet &	parent,
		const xAOD::Vertex *	originVertex = `nullptr`
	)		const

Definition at line 43 of file PseudoJetTranslator.cxx.

                                                                                    {
   xAOD::Jet& jet = translate(pj, pjCont, jetCont, originVertex);
  
   const xAOD::JetContainer* parentCont = dynamic_cast<const xAOD::JetContainer*>(parent.container());
   if ( parentCont == nullptr ) { return jet ;}  // can this happen? if so THIS IS an ERROR ! should do something
  
   ElementLink<xAOD::JetContainer> el(*parentCont, parent.index());
   static const SG::AuxElement::Accessor<ElementLink<xAOD::JetContainer> > parentELacc("Parent"); 
   parentELacc(jet) =el;
  
   jet.setInputType(parent.getInputType());
   jet.setAlgorithmType(parent.getAlgorithmType());
   jet.setSizeParameter(parent.getSizeParameter());
   jet.setConstituentsSignalState(parent.getConstituentsSignalState());
  
   float area=0;
   if(parent.getAttribute(xAOD::JetAttribute::JetGhostArea, area)) jet.setAttribute(xAOD::JetAttribute::JetGhostArea, area);
     
   return jet;
 }

◆ translate() [2/2]

xAOD::Jet & PseudoJetTranslator::translate	(	const fastjet::PseudoJet &	pj,
		const PseudoJetContainer &	pjCont,
		xAOD::JetContainer &	jetCont,
		const xAOD::Vertex *	originVertex = `nullptr`
	)		const

Definition at line 7 of file PseudoJetTranslator.cxx.

                                                                                        {
  
   // Create a new jet in place at the end of the container
   jetCont.emplace_back(new xAOD::Jet());
   xAOD::Jet& jet = *jetCont.back();
   jet.setJetP4( xAOD::JetFourMom_t( pj.pt(), pj.eta(), pj.phi(), pj.m() ) );
  
   const static SG::AuxElement::Accessor<const fastjet::PseudoJet*> pjAccessor("PseudoJet");
   pjAccessor(jet) = &pj;
  
   // Record the jet-finding momentum, i.e. the one used to find/groom the jet.
   jet.setJetP4(xAOD::JetConstitScaleMomentum, jet.jetP4());
   
   // save area if needed ---------
   if( pj.has_area() ){
  
     if(m_saveArea) jet.setAttribute(xAOD::JetAttribute::ActiveArea,pj.area());
     if(m_saveArea4Vec){
         fastjet::PseudoJet pja = pj.area_4vector();
         xAOD::JetFourMom_t fvarea(pja.pt(), pja.eta(), pja.phi(), pja.m());
         jet.setAttribute(xAOD::JetAttribute::ActiveArea4vec, fvarea);      
     }    
   }// area -------------
   
   if (originVertex == nullptr){
     pjCont.extractConstituents(jet, pj);
   }
   else{
     pjCont.extractByVertexConstituents(jet, pj, originVertex);
   }
  
   return jet;
 }

Member Data Documentation

◆ m_saveArea

bool PseudoJetTranslator::m_saveArea = true

protected

Definition at line 28 of file PseudoJetTranslator.h.

◆ m_saveArea4Vec

bool PseudoJetTranslator::m_saveArea4Vec = true