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	addBranches () const override

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

Private Attributes
Gaudi::Property< std::string >	m_eventInfoName { this, "EventInfoName", "EventInfo", "Event Info Key"}
Gaudi::Property< std::string >	m_mcName { this,"MCCollectionName", "TruthParticles", "MC Collection Key"}

Detailed Description

Definition at line 26 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 23 of file SUSYSignalTagger.cxx.

                                                                                               :
    base_class(t,n,p){
  }  

~SUSYSignalTagger()

DerivationFramework::SUSYSignalTagger::~SUSYSignalTagger ( )

default

Member Function Documentation

addBranches()

StatusCode DerivationFramework::SUSYSignalTagger::addBranches ( ) const

overridevirtual

Definition at line 27 of file SUSYSignalTagger.cxx.

                                                {
    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_DEBUG("WARNING could not retrieve TruthParticleContainer " <<m_mcName);
      return StatusCode::FAILURE;
    }
    //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 61 of file SUSYSignalTagger.cxx.

                                   {
    pdgid1 = 0;
    pdgid2 = 0;
    const xAOD::TruthParticle* firstsp(0);
    const xAOD::TruthParticle* secondsp(0);
    if (!truthP || 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;
  }

Member Data Documentation

m_eventInfoName

Gaudi::Property<std::string> DerivationFramework::SUSYSignalTagger::m_eventInfoName { this, "EventInfoName", "EventInfo", "Event Info Key"}

private

Definition at line 34 of file SUSYSignalTagger.h.

m_mcName

Gaudi::Property<std::string> DerivationFramework::SUSYSignalTagger::m_mcName { this,"MCCollectionName", "TruthParticles", "MC Collection Key"}

private

Definition at line 35 of file SUSYSignalTagger.h.

---

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

• SUSYSignalTagger.h
• SUSYSignalTagger.cxx