DerivationFramework::SUSYGenFilterTool Class Reference

#include <SUSYGenFilterTool.h>

Inheritance diagram for DerivationFramework::SUSYGenFilterTool:

Collaboration diagram for DerivationFramework::SUSYGenFilterTool:

Public Member Functions
	SUSYGenFilterTool (const std::string &t, const std::string &n, const IInterface *p)
	~SUSYGenFilterTool ()
virtual StatusCode	initialize () override
virtual StatusCode	addBranches () const override
StatusCode	getGenFiltVars (const xAOD::TruthParticleContainer *tpc, float &genFiltHT, float &genFiltMET) const
bool	isPrompt (const xAOD::TruthParticle *tp) const
MCTruthPartClassifier::ParticleOrigin	getPartOrigin (const xAOD::TruthParticle *tp) const

Private Attributes
std::string	m_eventInfoName
std::string	m_mcName
std::string	m_truthJetsName
float	m_MinJetPt
	Min pT for the truth jets. More...
float	m_MaxJetEta
	Max eta for the truth jets. More...
float	m_MinLepPt
	Min pT for the truth leptons. More...
float	m_MaxLepEta
	Max eta for the truth leptons. More...
ToolHandle< IMCTruthClassifier >	m_classif

Detailed Description

Definition at line 28 of file SUSYGenFilterTool.h.

Constructor & Destructor Documentation

SUSYGenFilterTool()

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

Definition at line 19 of file SUSYGenFilterTool.cxx.

                                                                                                 :
    base_class(t,n,p),
    m_classif("MCTruthClassifier/SUSYGenFilt_MCTruthClassifier")
  {
    
    
    declareProperty("EventInfoName",m_eventInfoName="EventInfo");
    declareProperty("MCCollectionName",m_mcName="TruthParticles");
    declareProperty("TruthJetCollectionName",m_truthJetsName="AntiKt4TruthWZJets");
    declareProperty("MinJetPt",m_MinJetPt = 35e3);  
    declareProperty("MaxJetEta",m_MaxJetEta = 2.5);
    declareProperty("MinLeptonPt",m_MinLepPt = 25e3);
    declareProperty("MaxLeptonEta",m_MaxLepEta = 2.5);
  }

~SUSYGenFilterTool()

DerivationFramework::SUSYGenFilterTool::~SUSYGenFilterTool

(

)

Definition at line 34 of file SUSYGenFilterTool.cxx.

{}

Member Function Documentation

addBranches()

StatusCode DerivationFramework::SUSYGenFilterTool::addBranches ( ) const

overridevirtual

Definition at line 72 of file SUSYGenFilterTool.cxx.

                                                 {
    ATH_MSG_VERBOSE("SUSYGenFilterTool::addBranches()");
    
    const xAOD::EventInfo* eventInfo;
    if (evtStore()->retrieve(eventInfo,m_eventInfoName).isFailure()) {
      ATH_MSG_ERROR("could not retrieve event info " <<m_eventInfoName);
      return StatusCode::FAILURE;
    }
    
    const xAOD::TruthParticleContainer* truthPC = 0;
    if (evtStore()->retrieve(truthPC,m_mcName).isFailure()) {
      ATH_MSG_ERROR("WARNING could not retrieve TruthParticleContainer " <<m_mcName);
      return StatusCode::FAILURE;
    }
 
    float genFiltHT(0.), genFiltMET(0.);
    ATH_CHECK( getGenFiltVars(truthPC, genFiltHT, genFiltMET) );
 
    ATH_MSG_DEBUG("Computed generator filter quantities: HT " << genFiltHT/1e3 << ", MET " << genFiltMET/1e3 );
 
    dec_genFiltHT(*eventInfo) = genFiltHT;
    dec_genFiltMET(*eventInfo) = genFiltMET;
 
    return StatusCode::SUCCESS;
  }

getGenFiltVars()

StatusCode DerivationFramework::SUSYGenFilterTool::getGenFiltVars	(	const xAOD::TruthParticleContainer *	tpc,
		float &	genFiltHT,
		float &	genFiltMET
	)		const

Definition at line 98 of file SUSYGenFilterTool.cxx.

                                                                                                                               {
    // Get jet container out
    const xAOD::JetContainer* truthjets = 0;
    if ( evtStore()->retrieve( truthjets, m_truthJetsName).isFailure() || !truthjets ){
      ATH_MSG_ERROR( "No xAOD::JetContainer found in StoreGate with key " << m_truthJetsName ); 
      return StatusCode::FAILURE;
    }
 
    // Get HT
    genFiltHT = -1;
    for (const auto tj : *truthjets) {
      if ( tj->pt()>m_MinJetPt && fabs(tj->eta())<m_MaxJetEta ) {
    ATH_MSG_VERBOSE("Adding truth jet with pt " << tj->pt()
            << ", eta " << tj->eta()
            << ", phi " << tj->phi()
            << ", nconst = " << tj->numConstituents());
    genFiltHT += tj->pt();
      }
    }
 
    float MEx(0.), MEy(0.);
    for (const auto tp : *tpc){
      int pdgid = tp->pdgId();
      if (HepMC::is_simulation_particle(tp)) continue; // Particle is from G4
      if (MC::isZeroEnergyPhoton(tp)) continue; // Work around for an old generator bug
      if ( !MC::isStable(tp) ) continue; // Stable!
 
      if ((MC::isElectron(pdgid) || MC::isMuon(pdgid)) && tp->pt()>m_MinLepPt && std::fabs(tp->eta())<m_MaxLepEta) {
    if( isPrompt(tp) ) {
      ATH_MSG_VERBOSE("Adding prompt lepton " << tp);
      genFiltHT += tp->pt();
    }
      }
 
      if (MC::isSpecialNonInteracting(tp) && isPrompt(tp) ) {
    ATH_MSG_VERBOSE("Found prompt nonInteracting particle " << tp);
    MEx += tp->px();
    MEy += tp->py();
      }
    }
    genFiltMET = std::sqrt(MEx*MEx+MEy*MEy);
 
    return StatusCode::SUCCESS;
  }

getPartOrigin()

MCTruthPartClassifier::ParticleOrigin DerivationFramework::SUSYGenFilterTool::getPartOrigin

(

const xAOD::TruthParticle *

tp

)

const

initialize()

StatusCode DerivationFramework::SUSYGenFilterTool::initialize ( )

overridevirtual

Definition at line 36 of file SUSYGenFilterTool.cxx.

                                          {
    ATH_MSG_INFO("Initialize " );
    return StatusCode::SUCCESS;
  }

isPrompt()

bool DerivationFramework::SUSYGenFilterTool::isPrompt

(

const xAOD::TruthParticle *

tp

)

const

Definition at line 41 of file SUSYGenFilterTool.cxx.

  {
    ParticleOrigin orig = m_classif->particleTruthClassifier( tp ).second;
    ATH_MSG_VERBOSE("Particle has origin " << orig);
 
    switch(orig) {
    case PhotonConv:
    case DalitzDec:
    case ElMagProc:
    case Mu:
    case TauLep:
    case LightMeson:
    case StrangeMeson:
    case CharmedMeson:
    case BottomMeson:
    case CCbarMeson:
    case JPsi:
    case BBbarMeson: 
    case LightBaryon:
    case StrangeBaryon:
    case CharmedBaryon: 
    case BottomBaryon:
    case PionDecay:
    case KaonDecay:
      return false;
    default:
      break;
    }
    return true;
  }

Member Data Documentation

m_classif

ToolHandle<IMCTruthClassifier> DerivationFramework::SUSYGenFilterTool::m_classif

private

Definition at line 52 of file SUSYGenFilterTool.h.

m_eventInfoName

std::string DerivationFramework::SUSYGenFilterTool::m_eventInfoName

private

Definition at line 43 of file SUSYGenFilterTool.h.

m_MaxJetEta

float DerivationFramework::SUSYGenFilterTool::m_MaxJetEta

private

Max eta for the truth jets.

Definition at line 48 of file SUSYGenFilterTool.h.

m_MaxLepEta

float DerivationFramework::SUSYGenFilterTool::m_MaxLepEta

private

Max eta for the truth leptons.

Definition at line 50 of file SUSYGenFilterTool.h.

m_mcName

std::string DerivationFramework::SUSYGenFilterTool::m_mcName

private

Definition at line 44 of file SUSYGenFilterTool.h.

m_MinJetPt

float DerivationFramework::SUSYGenFilterTool::m_MinJetPt

private

Min pT for the truth jets.

Definition at line 47 of file SUSYGenFilterTool.h.

m_MinLepPt

float DerivationFramework::SUSYGenFilterTool::m_MinLepPt

private

Min pT for the truth leptons.

Definition at line 49 of file SUSYGenFilterTool.h.

m_truthJetsName

std::string DerivationFramework::SUSYGenFilterTool::m_truthJetsName

private

Definition at line 45 of file SUSYGenFilterTool.h.

---

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

• SUSYGenFilterTool.h
• SUSYGenFilterTool.cxx