#include <TruthD2Decorator.h>

Inheritance diagram for DerivationFramework::TruthD2Decorator:

Collaboration diagram for DerivationFramework::TruthD2Decorator:

Public Member Functions
	TruthD2Decorator (const std::string &t, const std::string &n, const IInterface *p)

StatusCode	initialize ()

virtual StatusCode	addBranches () const

Private Attributes
SG::ReadHandleKey< xAOD::JetContainer >	m_jetContainerKey {this, "JetContainerKey", "AntiKt10TruthTrimmedPtFrac5SmallR20Jets", "Name of jet container key for input"}

SG::WriteDecorHandleKey< xAOD::JetContainer >	m_decorationName {this, "DecorationName", "AntiKt10TruthTrimmedPtFrac5SmallR20Jets.D2", "Decoration Name"}

Detailed Description

Definition at line 22 of file TruthD2Decorator.h.

Constructor & Destructor Documentation

◆ TruthD2Decorator()

DerivationFramework::TruthD2Decorator::TruthD2Decorator	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 18 of file TruthD2Decorator.cxx.

                            : 
     base_class(t,n,p)
   {
  
   }

Member Function Documentation

◆ addBranches()

StatusCode DerivationFramework::TruthD2Decorator::addBranches ( ) const

virtual

Definition at line 38 of file TruthD2Decorator.cxx.

   {
  
       // Event context
       const EventContext& ctx = Gaudi::Hive::currentContext();
  
       // Set up the decorators 
       SG::WriteDecorHandle< xAOD::JetContainer, float > decoratorD2(m_decorationName, ctx); 
  
       // Get the Large-R jet Container
       SG::ReadHandle<xAOD::JetContainer> largeRjets(m_jetContainerKey, ctx);
       
       if(!largeRjets.isValid()) {
         ATH_MSG_ERROR ("Couldn't retrieve JetContainer with key " << m_jetContainerKey.key());
         return StatusCode::FAILURE;
       }
  
       // loop over jet collection
       for( const auto *jet: *largeRjets){
         //get ECF
         float ecf1 = jet->getAttribute<float>("ECF1");
         float ecf2 = jet->getAttribute<float>("ECF2");
         float ecf3 = jet->getAttribute<float>("ECF3");
         
         //calculate D2 and decorate
         float D2=-999;
         if(fabs(ecf2)>1e-8)
           D2=ecf3 * pow(ecf1, 3.0) / pow(ecf2, 3.0);
         decoratorD2(*jet) = D2;
       }
  
       return StatusCode::SUCCESS;
   }

◆ initialize()

StatusCode DerivationFramework::TruthD2Decorator::initialize ( )

Definition at line 26 of file TruthD2Decorator.cxx.

   {
  
       ATH_CHECK(m_jetContainerKey.initialize()); 
       m_decorationName = m_jetContainerKey.key()+".D2";
       ATH_CHECK(m_decorationName.initialize());
  
       return StatusCode::SUCCESS;
  
   }

Member Data Documentation

◆ m_decorationName

SG::WriteDecorHandleKey<xAOD::JetContainer> DerivationFramework::TruthD2Decorator::m_decorationName {this, "DecorationName", "AntiKt10TruthTrimmedPtFrac5SmallR20Jets.D2", "Decoration Name"}

private

Definition at line 31 of file TruthD2Decorator.h.

◆ m_jetContainerKey

SG::ReadHandleKey<xAOD::JetContainer> DerivationFramework::TruthD2Decorator::m_jetContainerKey {this, "JetContainerKey", "AntiKt10TruthTrimmedPtFrac5SmallR20Jets", "Name of jet container key for input"}