DerivationFramework::SUSYSignalTagger Class Reference

#include <SUSYSignalTagger.h>

Inheritance diagram for DerivationFramework::SUSYSignalTagger:

Collaboration diagram for DerivationFramework::SUSYSignalTagger:

Public Member Functions
	SUSYSignalTagger (const std::string &t, const std::string &n, const IInterface *p)
	~SUSYSignalTagger ()=default
virtual StatusCode	initialize () override final
virtual StatusCode	addBranches (const EventContext &ctx) const override

Private Member Functions
bool	FindSusyHardProc (const xAOD::TruthParticleContainer &truthP, int &pdgid1, int &pdgid2) const

Private Attributes
SG::ReadHandleKey< xAOD::TruthParticleContainer >	m_mcName { this,"MCCollectionName", "TruthParticles", "MC Collection Key"}
SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfoName { this, "EventInfoName", "EventInfo", "Event Info Key"}
SG::WriteDecorHandleKey< xAOD::EventInfo >	m_dec_procIDKey {this, "SUSY_procIDKey", m_eventInfoName, "SUSY_procID"}
SG::WriteDecorHandleKey< xAOD::EventInfo >	m_dec_pdgId1Key {this, "SUSY_pid1Key", m_eventInfoName, "SUSY_pid1"}
SG::WriteDecorHandleKey< xAOD::EventInfo >	m_dec_pdgId2Key {this, "SUSY_pid2Key", m_eventInfoName, "SUSY_pid2"}

Detailed Description

Definition at line 28 of file SUSYSignalTagger.h.

Constructor & Destructor Documentation

SUSYSignalTagger()

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

Definition at line 18 of file SUSYSignalTagger.cxx.

                                                                                               :
    base_class(t,n,p){
  }

~SUSYSignalTagger()

DerivationFramework::SUSYSignalTagger::~SUSYSignalTagger ( )

default

Member Function Documentation

addBranches()

StatusCode DerivationFramework::SUSYSignalTagger::addBranches ( const EventContext &  ctx ) const

overridevirtual

Definition at line 31 of file SUSYSignalTagger.cxx.

                                                                       {
    SG::ReadHandle<xAOD::EventInfo> eventInfo(m_eventInfoName, ctx);
    if (!eventInfo.isValid()) {
      ATH_MSG_ERROR("could not retrieve event info " <<m_eventInfoName);
      return StatusCode::FAILURE;
    }
 
    SG::ReadHandle<xAOD::TruthParticleContainer> truthPC(m_mcName, ctx);
    if (!truthPC.isValid()) {
      ATH_MSG_ERROR("WARNING could not retrieve TruthParticleContainer " <<m_mcName);
      return StatusCode::FAILURE;
    }
 
    SG::WriteDecorHandle<xAOD::EventInfo, int> dec_procID(m_dec_procIDKey, ctx);
    SG::WriteDecorHandle<xAOD::EventInfo, int> dec_pdgId1(m_dec_pdgId1Key, ctx);
    SG::WriteDecorHandle<xAOD::EventInfo, int> dec_pdgId2(m_dec_pdgId2Key, ctx);
 
    //Identify SUSY hard proc
    int pdgId1(0);
    int pdgId2(0);
    bool found = FindSusyHardProc( *truthPC, pdgId1, pdgId2);
    if (!found) {
      ATH_MSG_DEBUG("could not identify SUSY process! ");
      dec_procID(*eventInfo) = 0;
      dec_pdgId1(*eventInfo) = -99;
      dec_pdgId2(*eventInfo) = -99;
      return StatusCode::SUCCESS;
    }
    //Get SUSY proc ID
    unsigned int procID = finalStateID(pdgId1, pdgId2);
    if (procID == 0)  ATH_MSG_WARNING("could not identify SUSY procID! ");
    dec_procID(*eventInfo) = procID;
    dec_pdgId1(*eventInfo) = pdgId1;
    dec_pdgId2(*eventInfo) = pdgId2;
    return StatusCode::SUCCESS;
  }

FindSusyHardProc()

bool DerivationFramework::SUSYSignalTagger::FindSusyHardProc ( const xAOD::TruthParticleContainer &  truthP, int &  pdgid1, int &  pdgid2  ) const

private

Definition at line 70 of file SUSYSignalTagger.cxx.

                                   {
    pdgid1 = 0;
    pdgid2 = 0;
    const xAOD::TruthParticle* firstsp(0);
    const xAOD::TruthParticle* secondsp(0);
    if (truthP.empty()) {
      return false;
    }
    for (const auto tp : truthP) {
      if (MC::isSquark(tp) || MC::isSlepton(tp) || MC::isGaugino(tp)) {
        if (tp->nParents() != 0) {
          if ( !MC::isSUSY(tp->parent(0))) {
            if (!firstsp) {
              firstsp = tp;
            } else if (!secondsp) {
              secondsp = tp;
            } else {
              if (firstsp->nChildren() != 0 && HepMC::is_same_particle(tp,firstsp->child(0))) {
                firstsp = tp;
              }
              else if (secondsp->nChildren() != 0 && HepMC::is_same_particle(tp,secondsp->child(0))) {
                secondsp = tp;
              }
              else if (firstsp->nChildren() != 0 && HepMC::is_same_particle(firstsp->child(0),secondsp)) {
                firstsp = secondsp;
                secondsp = tp;
              }
              else if (secondsp->nChildren() != 0 && HepMC::is_same_particle(secondsp->child(0),firstsp)) {
                secondsp = firstsp;
                firstsp = tp;
              }
            }
          }
        }
      }
    }
    // quit if no sparticles found
    if (!firstsp && !secondsp) return false; // should find none or two
 
    if (firstsp && firstsp->nChildren() == 1) {
      for (const auto tp : truthP) {
        if (HepMC::is_same_particle(firstsp->child(0),tp) && tp->pdgId() != firstsp->pdgId()) {
          firstsp = tp;
          break;
        }
      }
    }
    if (secondsp && secondsp->nChildren() == 1) {
      for (const auto tp : truthP) {
        if (HepMC::is_same_particle(secondsp->child(0),tp) && tp->pdgId() != secondsp->pdgId()) {
          secondsp = tp;
          break;
        }
      }
    }
    if (firstsp && abs(firstsp->pdgId()) > 1000000) pdgid1 = firstsp->pdgId();  // Replace with (firstsp && MC::isSUSY(firstsp)) ?
    if (secondsp && abs(secondsp->pdgId()) > 1000000) pdgid2 = secondsp->pdgId(); // Replace with (secondsp && MC::isSUSY(secondsp)) ?
    // Return gracefully:
    return true;
  }

initialize()

StatusCode DerivationFramework::SUSYSignalTagger::initialize ( )

finaloverridevirtual

Definition at line 22 of file SUSYSignalTagger.cxx.

                                          {
    ATH_CHECK(m_eventInfoName.initialize());
    ATH_CHECK(m_mcName.initialize());
    ATH_CHECK(m_dec_procIDKey.initialize());
    ATH_CHECK(m_dec_pdgId1Key.initialize());
    ATH_CHECK(m_dec_pdgId2Key.initialize());
    return StatusCode::SUCCESS;
  }

Member Data Documentation

m_dec_pdgId1Key

SG::WriteDecorHandleKey<xAOD::EventInfo> DerivationFramework::SUSYSignalTagger::m_dec_pdgId1Key {this, "SUSY_pid1Key", m_eventInfoName, "SUSY_pid1"}

private

Definition at line 42 of file SUSYSignalTagger.h.

m_dec_pdgId2Key

SG::WriteDecorHandleKey<xAOD::EventInfo> DerivationFramework::SUSYSignalTagger::m_dec_pdgId2Key {this, "SUSY_pid2Key", m_eventInfoName, "SUSY_pid2"}

private

Definition at line 43 of file SUSYSignalTagger.h.

m_dec_procIDKey

SG::WriteDecorHandleKey<xAOD::EventInfo> DerivationFramework::SUSYSignalTagger::m_dec_procIDKey {this, "SUSY_procIDKey", m_eventInfoName, "SUSY_procID"}

private

Definition at line 41 of file SUSYSignalTagger.h.

m_eventInfoName

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::SUSYSignalTagger::m_eventInfoName { this, "EventInfoName", "EventInfo", "Event Info Key"}

private

Definition at line 40 of file SUSYSignalTagger.h.

m_mcName

SG::ReadHandleKey<xAOD::TruthParticleContainer> DerivationFramework::SUSYSignalTagger::m_mcName { this,"MCCollectionName", "TruthParticles", "MC Collection Key"}

private

Definition at line 39 of file SUSYSignalTagger.h.

---

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

• SUSYSignalTagger.h
• SUSYSignalTagger.cxx