TrigMuonEFIdtpInvMassHypoTool Class Reference

#include <TrigMuonEFIdtpInvMassHypoTool.h>

Inheritance diagram for TrigMuonEFIdtpInvMassHypoTool:

Collaboration diagram for TrigMuonEFIdtpInvMassHypoTool:

Public Member Functions
	TrigMuonEFIdtpInvMassHypoTool (const std::string &type, const std::string &name, const IInterface *parent)
	~TrigMuonEFIdtpInvMassHypoTool ()=default
virtual StatusCode	initialize () override
virtual StatusCode	decide (Combo::LegDecisionsMap &passingLegs, const EventContext &) const override
	retrieves the decisions associated to this decId, make their combinations and apply the algorithm More...
virtual HLT::Identifier	decisionId () const
	retrieves this ComboHypoTool's chain's decision ID More...
StatusCode	setLegMultiplicity (const Combo::MultiplicityReqMap &multiplicityRequiredMap)
	Sets the number of legs and the multiplicity required on each leg. More...
const std::vector< int > &	legMultiplicity () const
	Gets the number of legs and the multiplicity required on each leg. More...
HLT::Identifier	legDecisionId (size_t i) const
	Retrieves this ComboHypoTool's chain's decision ID for a given leg. More...
const std::vector< HLT::Identifier > &	legDecisionIds () const
	Retrieves this ComboHypoTool's chain's decision IDs for all legs. More...
virtual StatusCode	decideOnSingleObject (TrigCompositeUtils::Decision *, const std::vector< const TrigCompositeUtils::DecisionIDContainer * > &) const
	Alternate method called by BPhysics ComboHypoAlgs instead of the base method decide(...). More...

Protected Member Functions
StatusCode	selectLegs (const Combo::LegDecisionsMap &IDCombMap, std::vector< std::vector< Combo::LegDecision >> &leg_decisions) const
	Creates the per-leg vectors of Decision objects starting from the initial LegDecision map, storing only those concerning this HypoTool's chain Pack the Decision objects in std::pair<DecisionID, ElementLink<Decision>> so the derived class' executeAlg function knows which leg each object is on. More...
void	updateLegDecisionsMap (const std::vector< std::vector< Combo::LegDecision >> &passing_comb, Combo::LegDecisionsMap &passingLegs) const
	For when the tool accepts some/all combinations. More...
void	eraseFromLegDecisionsMap (Combo::LegDecisionsMap &passingLegs) const
	For when the tool rejects all combinations. More...
StatusCode	printDebugInformation (const Combo::LegDecisionsMap &passingLegs) const
	Print the output of the tool, after having removed failed Decision Objects. More...

Protected Attributes
Gaudi::Property< size_t >	m_combinationsThresholdWarn
Gaudi::Property< size_t >	m_combinationsThresholdBreak
Gaudi::Property< bool >	m_modeOR
Gaudi::Property< bool >	m_enableOverride

Private Member Functions
bool	executeAlg (const std::vector< Combo::LegDecision > &combinaiton) const override
	Only a dummy implementation exists in ComboHypoToolBase. More...
StatusCode	doTPIdperf (const xAOD::TrackParticle *metrack, const std::vector< const xAOD::TrackParticle * > &tracks_pt, const std::vector< const xAOD::TrackParticle * > &tracks_ftf) const
float	qOverPMatching (const xAOD::TrackParticle *metrack, const xAOD::TrackParticle *idtrack) const
bool	passedCBQualityCuts (const xAOD::Muon *muon) const
bool	passedSAQualityCuts (const xAOD::Muon *muon) const

Private Attributes
Gaudi::Property< double >	m_invMassLow
Gaudi::Property< double >	m_invMassHigh
Gaudi::Property< bool >	m_acceptAll
Gaudi::Property< bool >	m_selOS
Gaudi::Property< bool >	m_muonqualityCut
ToolHandle< GenericMonitoringTool >	m_monTool {this, "MonTool", "", "Monitoring tool" }
ToolHandle< CP::IMuonSelectionTool >	m_muonSelTool {this, "MuonSelectionTool", "CP::MuonSelectionTool/MuonSelectionTool", "Tool for muon quality selection"}
HLT::Identifier	m_decisionId
	The DecisionID of the chain, obtained from the Tool's name. More...
std::vector< HLT::Identifier >	m_legDecisionIds
	The DecisionIDs of the individual legs, derived from both m_decisionId and m_legMultiplicities. More...
std::vector< int >	m_legMultiplicities
	The number of legs, and the required multiplicity on each leg. More...

Detailed Description

Definition at line 19 of file TrigMuonEFIdtpInvMassHypoTool.h.

Constructor & Destructor Documentation

TrigMuonEFIdtpInvMassHypoTool()

TrigMuonEFIdtpInvMassHypoTool::TrigMuonEFIdtpInvMassHypoTool	(	const std::string &	type,
		const std::string &	name,
		const IInterface *	parent
	)

Definition at line 16 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

                                                                                                                                    :
   ComboHypoToolBase(type, name, parent) {
}

~TrigMuonEFIdtpInvMassHypoTool()

TrigMuonEFIdtpInvMassHypoTool::~TrigMuonEFIdtpInvMassHypoTool ( )

default

Member Function Documentation

decide()

StatusCode ComboHypoToolBase::decide ( Combo::LegDecisionsMap &  passingLegs, const EventContext &    ) const

overridevirtualinherited

retrieves the decisions associated to this decId, make their combinations and apply the algorithm

Parameters

[in]	LegDecisionsMap	that lists all the passing decisions, to be updated by the tool depending on the outcome of executeAlg
[in]	Event	Context, currently unused

Reimplemented in DisplacedJetRankComboHypoTool.

Definition at line 48 of file ComboHypoToolBase.cxx.

                                                                                                   {
  if (m_legMultiplicities.size() == 0) {
    ATH_MSG_ERROR("ComboHypoTool for " << m_decisionId << " has not been properly configured. setLegMultiplicity should be called by the parent alg in initalize");
    return StatusCode::FAILURE;
  }
 
  // if no combinations passed, then exit 
  if (passingLegs.size() == 0) {
    return StatusCode::SUCCESS;
  }
 
  ATH_MSG_DEBUG("Looking for legs from " << decisionId() << " in the map. Map contains features for " << passingLegs.size() << " legs, which may be data for many chains.");
 
  // select the leg decisions from the map with this ID:
  std::vector<std::vector<Combo::LegDecision>> legDecisions;
  ATH_CHECK(selectLegs(passingLegs, legDecisions));
 
  for (size_t legIndex = 0; legIndex < m_legMultiplicities.size(); ++legIndex) {
    // The static cast is because this gets read as a Gaudi property which only seems to work with vector<int> rather than vector<size_t> 
    if (legDecisions.at(legIndex).size() < static_cast<size_t>(m_legMultiplicities.at(legIndex))) {
      ATH_MSG_DEBUG("Too few features are found for " << m_decisionId << " on leg " << legIndex <<", require:" << m_legMultiplicities.at(legIndex) << " have:" << legDecisions.at(legIndex).size());
      ATH_MSG_DEBUG("This ComboHypoTool cannot run in this event, this chain **REJECTS** this event.");
      eraseFromLegDecisionsMap(passingLegs);
      ATH_CHECK(printDebugInformation(passingLegs));
      return StatusCode::SUCCESS;
    }
  }
 
  // Create and initialise the combinations generator for the requirements of this chain, given the objects available in this event.
  HLT::NestedUniqueCombinationGenerator nucg;
  for (size_t legIndex = 0; legIndex < m_legMultiplicities.size(); ++legIndex) {
    const size_t choose_any = m_legMultiplicities.at(legIndex);
    const size_t out_of = legDecisions.at(legIndex).size();
    nucg.add({out_of, choose_any});
    ATH_MSG_DEBUG("For leg " << legIndex << " we will be choosing any " << choose_any << " Decision Objects out of " << out_of);
  }
 
  std::vector<std::vector<Combo::LegDecision>> passingCombinations;
 
  size_t warnings = 0, iterations = 0;
  do {
 
    const std::vector<size_t> combination = nucg();
    ++nucg;
    std::vector<Combo::LegDecision> combinationToCheck;
 
    size_t location_in_combination = 0;
    for (size_t legIndex = 0; legIndex < m_legMultiplicities.size(); ++legIndex) {
      // We loop over however many objects are required on the leg,
      // but we take their index from the 'combination'. Hence 'object' is not used directly
      for (size_t object = 0; object < static_cast<size_t>(m_legMultiplicities.at(legIndex)); ++object) {
        const size_t objectIndex = combination.at(location_in_combination++);
        combinationToCheck.push_back( legDecisions.at(legIndex).at(objectIndex) );
      }
    }
 
    ++iterations;
 
    try {
      if (executeAlg(combinationToCheck)) {
        ATH_MSG_DEBUG("Combination " << (iterations - 1) << " decided to be passing");
        passingCombinations.push_back(combinationToCheck);
        if (m_modeOR == true and m_enableOverride) {
          break;
        }
      } else { // the combination failed
        if (m_modeOR == false and m_enableOverride) {
          break;
        }
      }
    } catch (std::exception& e) {
      ATH_MSG_ERROR(e.what());
      return StatusCode::FAILURE;
    }
 
    if ((iterations >= m_combinationsThresholdWarn && warnings == 0) or (iterations >= m_combinationsThresholdBreak)) {
      ATH_MSG_WARNING("Have so far processed " << iterations << " combinations for " << m_decisionId << " in this event, " << passingCombinations.size() << " passing.");
      ++warnings;
      if (iterations >= m_combinationsThresholdBreak) {
        ATH_MSG_WARNING("Too many combinations! Breaking the loop at this point.");
        break;
      }
    }
 
  } while (nucg);
 
 
  if (m_modeOR) {
 
    ATH_MSG_DEBUG("Passing " << passingCombinations.size() << " combinations out of " << iterations << ", " 
      << m_decisionId << (passingCombinations.size() ? " **ACCEPTS**" : " **REJECTS**") << " this event based on OR logic.");
 
    if (m_enableOverride) {
      ATH_MSG_DEBUG("Note: stopped after the first successful combination due to the EnableOverride flag.");  
    }
 
  } else {  // modeAND
 
    const bool passAll = (passingCombinations.size() == iterations);
 
    ATH_MSG_DEBUG("Passing " << passingCombinations.size() << " combinations out of " << iterations << ", " 
      << m_decisionId << (passAll ? " **ACCEPTS**" : " **REJECTS**") << " this event based on AND logic.");
 
    if (m_enableOverride) {
      ATH_MSG_DEBUG("Note: stopped after the first failed combination due to the EnableOverride flag.");  
    }
 
    if (not passAll) {
      passingCombinations.clear();
    }
 
  }
 
  if (not passingCombinations.empty()) { // need partial erasure of the decsions (only those not present in any combination)
    updateLegDecisionsMap(passingCombinations, passingLegs);
  } else { // need complete erasure of input decisions
    eraseFromLegDecisionsMap(passingLegs);
  }
 
  ATH_CHECK(printDebugInformation(passingLegs));
  return StatusCode::SUCCESS;
}

decideOnSingleObject()

StatusCode ComboHypoToolBase::decideOnSingleObject ( TrigCompositeUtils::Decision *  , const std::vector< const TrigCompositeUtils::DecisionIDContainer * > &    ) const

virtualinherited

Alternate method called by BPhysics ComboHypoAlgs instead of the base method decide(...).

This function should be considered a specialist use-case only. It must be over-ridden to do anything useful.

Reimplemented in TrigMultiTrkComboHypoTool, and TrigBmumuxComboHypoTool.

Definition at line 284 of file ComboHypoToolBase.cxx.

                                                                                                                                  {
  ATH_MSG_ERROR("Do not use ComboHypoToolBase on its own, inherit this class and override decideOnSingleObject.");
  ATH_MSG_ERROR("NOTE: Only if you are also supplying your own decide(...) implimentation, or similar.");
  ATH_MSG_ERROR("NOTE: Most uses cases should only need to override executeAlg(...).");
  return StatusCode::FAILURE;
}

decisionId()

virtual HLT::Identifier ComboHypoToolBase::decisionId ( ) const

inlinevirtualinherited

retrieves this ComboHypoTool's chain's decision ID

Definition at line 41 of file ComboHypoToolBase.h.

{ return m_decisionId; }

doTPIdperf()

StatusCode TrigMuonEFIdtpInvMassHypoTool::doTPIdperf ( const xAOD::TrackParticle *  metrack, const std::vector< const xAOD::TrackParticle * > &  tracks_pt, const std::vector< const xAOD::TrackParticle * > &  tracks_ftf  ) const

private

Definition at line 163 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

{
   ATH_MSG_VERBOSE("----- doTPIdperf -----");
 
   // monitored variables
   std::vector<float>  mnt_PT_dr,  mnt_PT_qovp;
   std::vector<float>  mnt_FTF_dr, mnt_FTF_qovp;
   auto mon_PT_dr     = Monitored::Collection("PT_dr",    mnt_PT_dr);
   auto mon_PT_qovp   = Monitored::Collection("PT_qovp",  mnt_PT_qovp);
   auto mon_FTF_dr    = Monitored::Collection("FTF_dr",   mnt_FTF_dr);
   auto mon_FTF_qovp  = Monitored::Collection("FTF_qovp", mnt_FTF_qovp);
   auto mon_phi_effi  = Monitored::Scalar<float>("probePhiEfficiency", -999.0);
   auto mon_eta_effi  = Monitored::Scalar<float>("probeEtaEfficiency", -999.0);
   auto mon_pt_found  = Monitored::Scalar<float>("PTfound",  -1);
   auto mon_ftf_found = Monitored::Scalar<float>("FTFfound", -1);
   auto mon_pt_phi    = Monitored::Scalar<float>("PTphi", -999.0);
   auto mon_pt_pix    = Monitored::Scalar<float>("PTpixelFound", -1);
   auto mon_pt_pixnext= Monitored::Scalar<float>("PTpixelNextToFound", -1);
   auto monitorIt     = Monitored::Group(m_monTool, mon_PT_dr, mon_PT_qovp, mon_FTF_dr, mon_FTF_qovp, mon_phi_effi, mon_eta_effi, mon_pt_found, mon_ftf_found, mon_pt_phi, mon_pt_pix, mon_pt_pixnext);
 
   // probe values
   auto mnt_probe_pt  = Monitored::Scalar<float>("probe_pt",  0);
   auto mnt_probe_eta = Monitored::Scalar<float>("probe_eta", 0);
   auto monProbe      = Monitored::Group(m_monTool, mnt_probe_pt, mnt_probe_eta);
   mnt_probe_pt  = metrack->pt()/Gaudi::Units::GeV;
   mnt_probe_eta = metrack->eta();
 
   // PT 
   const float PT_DR_CUT   = 0.1;
   const float PT_QOVP_CUT = 2.0;
   ATH_MSG_VERBOSE("PT: n size="<<tracks_pt.size());
   int n_pt_matched = 0;
   float min_dr_pt_matched = 999;
   uint8_t pt_expectInnermost=0;
   uint8_t pt_numberInnermost=0;
   uint8_t pt_expectNextToInnermost=0;
   uint8_t pt_numberNextToInnermost=0;
   for(auto idtrack : tracks_pt) {
      float dr_pt   = xAOD::P4Helpers::deltaR(metrack,idtrack);
      float qovp_pt = TrigMuonEFIdtpCommon::qOverPMatching(metrack,idtrack);
      mnt_PT_dr.push_back(dr_pt);
      mnt_PT_qovp.push_back(qovp_pt);
      ATH_MSG_VERBOSE("... dr (ME-PT): "<<dr_pt);
      ATH_MSG_VERBOSE("... Q/p match (ME-PT): "<<qovp_pt);
      if( dr_pt > PT_DR_CUT ) continue;
      if( qovp_pt > PT_QOVP_CUT ) continue;
      ++n_pt_matched;
      if( dr_pt < min_dr_pt_matched ) {
     min_dr_pt_matched = dr_pt;
     mon_pt_phi = idtrack->phi();
     idtrack->summaryValue(pt_expectInnermost,xAOD::expectInnermostPixelLayerHit);
     idtrack->summaryValue(pt_numberInnermost,xAOD::numberOfInnermostPixelLayerHits);
     idtrack->summaryValue(pt_expectNextToInnermost,xAOD::expectNextToInnermostPixelLayerHit);
     idtrack->summaryValue(pt_numberNextToInnermost,xAOD::numberOfNextToInnermostPixelLayerHits);    
      }
   }
   if( n_pt_matched > 1 ) n_pt_matched=1;
 
   // FTF
   const float FTF_DR_CUT   = 0.2;
   const float FTF_QOVP_CUT = 3.0;
   ATH_MSG_VERBOSE("FTF: n size="<<tracks_ftf.size());
   int n_ftf_matched = 0;
   for(auto idtrack : tracks_ftf) {
      float dr_pt   = xAOD::P4Helpers::deltaR(metrack,idtrack);
      float qovp_pt = TrigMuonEFIdtpCommon::qOverPMatching(metrack,idtrack);
      mnt_FTF_dr.push_back(dr_pt);
      mnt_FTF_qovp.push_back(qovp_pt);
      ATH_MSG_VERBOSE("... dr (ME-FTF): "<<dr_pt);
      ATH_MSG_VERBOSE("... Q/p match (ME-FTF): "<<qovp_pt);
      if( dr_pt > FTF_DR_CUT ) continue;
      if( qovp_pt > FTF_QOVP_CUT ) continue;
      ++n_ftf_matched;
   }
   if( n_ftf_matched > 1 ) n_ftf_matched=1;
   
   // efficiency
   float me_eta    = metrack->eta();
   float me_phi    = metrack->phi();
   float me_pt_gev = metrack->pt()/Gaudi::Units::GeV;
 
   int eta_nr = 0;
   if( std::abs(me_eta) > 2.0 ) { eta_nr = 2; }
   else if( std::abs(me_eta) > 1.0 ) { eta_nr = 1; }
   int   pt_nr = 0;
   if( me_pt_gev > 20.0 ) { pt_nr = 1; }
 
   std::stringstream ss;
   ss << "PT_effi_pt" << pt_nr << "_eta" << eta_nr;
   auto mnt_PT_effi = Monitored::Scalar<int>(ss.str(), 0);
   ss.clear();
   ss << "PT_effi_pt" << pt_nr << "_eta" << eta_nr;
   auto mnt_FTF_effi = Monitored::Scalar<int>(ss.str(), 0);
   auto monEffi = Monitored::Group(m_monTool, mnt_PT_effi, mnt_FTF_effi);
   mnt_PT_effi  = n_pt_matched;
   mnt_FTF_effi = n_ftf_matched;
 
   // TnP monitoring
   mon_phi_effi  = me_phi;
   mon_eta_effi  = me_eta;
   mon_pt_found  = n_pt_matched;
   mon_ftf_found = n_ftf_matched;
   if( n_pt_matched != 0 ) {
      if( pt_expectInnermost == 1 ) {
     if( pt_numberInnermost>0 ) {
        mon_pt_pix=1;
     } else {
        mon_pt_pix=0;
     }
      }
      if( pt_expectNextToInnermost == 1 ) {
     if( pt_numberNextToInnermost>0 ) {
        mon_pt_pixnext=1;
     } else {
        mon_pt_pixnext=0;
     }
      }
   }
 
   //
   return StatusCode::SUCCESS;
}

eraseFromLegDecisionsMap()

void ComboHypoToolBase::eraseFromLegDecisionsMap ( Combo::LegDecisionsMap &  passingLegs ) const

protectedinherited

For when the tool rejects all combinations.

Remove all Decision Objects from all the legs of this HypoTool's chain.

Definition at line 253 of file ComboHypoToolBase.cxx.

                                                                                        {
  for (auto& it : passingLegs) {
    DecisionID id = it.first;
    if (id == m_decisionId or (isLegId(id) and getIDFromLeg(id) == m_decisionId)) {
      const size_t nDecisionObjects = it.second.size();
      it.second.clear();
      ATH_MSG_VERBOSE("-- Removed " << nDecisionObjects << " from " << id);
    }
  }
}

executeAlg()

bool TrigMuonEFIdtpInvMassHypoTool::executeAlg ( const std::vector< Combo::LegDecision > &  combination ) const

overrideprivatevirtual

Only a dummy implementation exists in ComboHypoToolBase.

This should be over-ridden by a derived class. The derived class should return a boolean pass/fail for each possible combination in the event. param[in] combination A single combination of objects to be discriminated against. Vector contains the required number of objects over all legs. Use the pair.first to tell which leg a given pair.second decision object belongs to in the current combination.

Reimplemented from ComboHypoToolBase.

Definition at line 53 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

{
   ATH_MSG_VERBOSE("in executeAlg");
   bool result = false;
 
   if( m_acceptAll ) {
      ATH_MSG_DEBUG("Accept property is set: taking all the events");
      return true;
   }
 
   // monitored variables
   std::vector<float>  mnt_invMass;
   auto mon_invMass   = Monitored::Collection("Mass", mnt_invMass);
   auto monitorIt     = Monitored::Group(m_monTool, mon_invMass);
 
   // retrieve muon and tracks
   std::vector<const xAOD::TrackParticle*> tracks_pt;
   std::vector<const xAOD::TrackParticle*> tracks_ftf;
   std::vector<const xAOD::Muon*> muons_cb;
   std::vector<const xAOD::Muon*> muons_sa;
   int i_leg_idperf = -1;
   for(auto leg: combination) {
      auto decision= (*(leg.second));
      auto i_leg = TrigCompositeUtils::getIndexFromLeg(leg.first);
      ATH_MSG_VERBOSE("i_leg="<<i_leg);
      auto muonLinks = TrigCompositeUtils::findLinks<xAOD::MuonContainer>( decision, TrigCompositeUtils::featureString(), TrigDefs::lastFeatureOfType);
      if( muonLinks.size() != 1 )        continue;
      if( ! muonLinks.at(0).isValid() )  continue;
      const xAOD::Muon *mu = *(muonLinks.at(0).link);
      if( ! mu->primaryTrackParticle() ) continue;
 
      bool is_idperf_muon = false;
      ATH_MSG_VERBOSE("... selected: muonType="<<mu->muonType()<<", pT="<<mu->pt()/Gaudi::Units::GeV<<", eta="<<mu->eta()<<", phi="<<mu->phi()<<", author="<<mu->author());
      if(mu->author()==xAOD::Muon::Author::MuidCo && mu->muonType()==xAOD::Muon::MuonType::Combined) {
     muons_cb.push_back(mu);
      }
      else if(mu->author()==xAOD::Muon::Author::MuidSA && mu->muonType()==xAOD::Muon::MuonType::MuonStandAlone) {
     muons_sa.push_back(mu);
     is_idperf_muon = true;
     if( i_leg_idperf == -1 ) { i_leg_idperf = i_leg; }
     else if( i_leg_idperf != (int)i_leg ) {
        ATH_MSG_ERROR("i_leg for idperf looks inconsistent: i_leg_idperf / i_leg="<<i_leg_idperf<<" / "<<i_leg);
        return result;
     }
      }
      
      // if this leg is idperf leg (muon is SA) get ID tracks
      if( is_idperf_muon ) {
     const std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrackParticleContainer> > ptLinks = TrigCompositeUtils::findLinks< xAOD::TrackParticleContainer >( decision, "HLT_IDTrack_Muon_IDTrig", TrigDefs::lastFeatureOfType);
     ATH_MSG_VERBOSE("PT TrackParticleContainer Links size = "<<ptLinks.size());
     if( ptLinks.size() == 1 && ptLinks.at(0).isValid() ) {
        const xAOD::TrackParticle* track = *(ptLinks.at(0).link);
        tracks_pt.push_back(track);
        float pt  = track->pt();
        float eta = track->eta();
        float phi = track->phi();
        ATH_MSG_VERBOSE("... pt / eta / phi = "<<pt/Gaudi::Units::GeV << " / " << eta << " / " << phi<<", fitter="<<track->trackFitter());
     }
     const std::vector< TrigCompositeUtils::LinkInfo<xAOD::TrackParticleContainer> > ftfLinks = TrigCompositeUtils::findLinks< xAOD::TrackParticleContainer >( decision, "HLT_IDTrack_Muon_FTF", TrigDefs::lastFeatureOfType);
     ATH_MSG_VERBOSE("FTF TrackParticleContainer Links size = "<< ftfLinks.size());
     if( ftfLinks.size() == 1 && ftfLinks.at(0).isValid() ) {
        const xAOD::TrackParticle* track = *(ftfLinks.at(0).link);
        tracks_ftf.push_back(track);
        float pt  = track->pt();
        float eta = track->eta();
        float phi = track->phi();
        ATH_MSG_VERBOSE("... pt / eta / phi = "<<pt/Gaudi::Units::GeV << " / " << eta << " / " << phi<<", fitter="<<track->trackFitter());
     }
      }
   } // end of combination loop
 
   // mass between CB and SA
   for(auto muon_cb : muons_cb) {
      if( m_muonqualityCut && ! passedCBQualityCuts(muon_cb) ) continue;
 
      for(auto muon_sa : muons_sa) {
     if( m_muonqualityCut && ! passedSAQualityCuts(muon_sa) ) continue;
 
     const xAOD::TrackParticle* tr_cb = muon_cb->trackParticle(xAOD::Muon::TrackParticleType::CombinedTrackParticle);
     const xAOD::TrackParticle* tr_sa = muon_sa->trackParticle(xAOD::Muon::TrackParticleType::ExtrapolatedMuonSpectrometerTrackParticle);
     if (!tr_cb || !tr_sa) {
        ATH_MSG_ERROR("Either CB or SA TrackParticle not found.");
        continue;
     }
     
     if( m_selOS && (muon_cb->charge()*muon_sa->charge() > 0) ) continue;
 
     float diMuMass = (tr_cb->p4()+tr_sa->p4()).M()/Gaudi::Units::GeV;
     mnt_invMass.push_back(diMuMass);
     ATH_MSG_VERBOSE("pt CB / pt SA = " << (*tr_cb).pt()/Gaudi::Units::GeV << " / " << (*tr_sa).pt()/Gaudi::Units::GeV << ", mass =" << diMuMass);
 
     if( ((m_invMassLow >=0 && diMuMass>m_invMassLow ) || m_invMassLow <0) &&
         ((m_invMassHigh>=0 && diMuMass<m_invMassHigh) || m_invMassHigh<0) ) {
        result = true;
        StatusCode sc = doTPIdperf(tr_sa, tracks_pt, tracks_ftf);
        if( sc != StatusCode::SUCCESS ) {
           ATH_MSG_ERROR("doTPIdperf failed with StatuCode="<<sc);
        }
     }
      }
   }
   
   //
   ATH_MSG_VERBOSE("idperf TP overall result is: "<<(result?"true":"false"));
   return result;   
}

initialize()

StatusCode TrigMuonEFIdtpInvMassHypoTool::initialize ( )

overridevirtual

Definition at line 23 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

{
   if( m_muonqualityCut ) {
      if( m_muonSelTool.retrieve().isFailure() ) {
     ATH_MSG_ERROR("Unable to retrieve " << m_muonSelTool);
     return StatusCode::FAILURE;
      }
   } else m_muonSelTool.disable();
 
   if( m_acceptAll ) {
      ATH_MSG_DEBUG("Accepting all the events!");
   } else {
      if( m_invMassLow<0 && m_invMassHigh<0 ) { 
     ATH_MSG_ERROR("Both mass cuts are <0. This is probably a configuration mistake.");
     return StatusCode::FAILURE;
      }
   }
 
   if ( not m_monTool.name().empty() ) {
      ATH_CHECK( m_monTool.retrieve() );
      ATH_MSG_DEBUG("MonTool name: " << m_monTool);
   }
 
   ATH_MSG_DEBUG("Initialization completed successfully");
   return StatusCode::SUCCESS;
}

legDecisionId()

HLT::Identifier ComboHypoToolBase::legDecisionId ( size_t  i ) const

inlineinherited

Retrieves this ComboHypoTool's chain's decision ID for a given leg.

Only populated for chains with more than one leg. For chains with one leg, use decisionId()

Definition at line 60 of file ComboHypoToolBase.h.

{ return m_legDecisionIds.at(i); }

legDecisionIds()

const std::vector<HLT::Identifier>& ComboHypoToolBase::legDecisionIds ( ) const

inlineinherited

Retrieves this ComboHypoTool's chain's decision IDs for all legs.

Only populated for chains with more than one leg. For chains with one leg, use decisionId()

Definition at line 66 of file ComboHypoToolBase.h.

{ return m_legDecisionIds; }

legMultiplicity()

const std::vector<int>& ComboHypoToolBase::legMultiplicity ( ) const

inlineinherited

Gets the number of legs and the multiplicity required on each leg.

Definition at line 54 of file ComboHypoToolBase.h.

{ return m_legMultiplicities; }

passedCBQualityCuts()

bool TrigMuonEFIdtpInvMassHypoTool::passedCBQualityCuts ( const xAOD::Muon *  muon ) const

private

Definition at line 289 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

{
   bool passCut = false;
 
   const xAOD::TrackParticle* idtrack = muon->trackParticle( xAOD::Muon::InnerDetectorTrackParticle );
   const xAOD::TrackParticle* metrack = muon->trackParticle( xAOD::Muon::ExtrapolatedMuonSpectrometerTrackParticle );
 
   const float CHI2_CUT = 8.0;
   const float QOVP_CUT = 7.0;
    
   if( idtrack && metrack ) {
      float qOverPsignif = TrigMuonEFIdtpCommon::qOverPMatching(metrack,idtrack);
      float reducedChi2  = muon->primaryTrackParticle()->chiSquared()/muon->primaryTrackParticle()->numberDoF(); 
      // Selection criteria based on the requirements that are part of the muon quality working points (offline)
      if(std::abs(reducedChi2) < CHI2_CUT && !m_muonSelTool->isBadMuon(*muon) && qOverPsignif<QOVP_CUT && muon->author()==xAOD::Muon::MuidCo) passCut = true;
   }
   
   return passCut;
}

passedSAQualityCuts()

bool TrigMuonEFIdtpInvMassHypoTool::passedSAQualityCuts ( const xAOD::Muon *  muon ) const

private

Definition at line 310 of file TrigMuonEFIdtpInvMassHypoTool.cxx.

{
   bool passCut = false;
 
   const xAOD::TrackParticle* metrack = muon->trackParticle( xAOD::Muon::ExtrapolatedMuonSpectrometerTrackParticle );
 
   const float CHI2_CUT = 8.0;
 
   if( metrack ) {
      float reducedChi2 = muon->primaryTrackParticle()->chiSquared()/muon->primaryTrackParticle()->numberDoF(); 
      if( std::abs(reducedChi2) < CHI2_CUT && !m_muonSelTool->isBadMuon(*muon) ) passCut = true;
   }
   
   return passCut;
}

printDebugInformation()

StatusCode ComboHypoToolBase::printDebugInformation ( const Combo::LegDecisionsMap &  passingLegs ) const

protectedinherited

Print the output of the tool, after having removed failed Decision Objects.

Restricted to the ComboHypoTool's chain's legs.

Definition at line 264 of file ComboHypoToolBase.cxx.

                                                                                                 {
  ATH_MSG_DEBUG("ComboHypoToolBase: End of " << m_decisionId << ", passing elements are: ");
  for (const auto& [id, ELV] : passingLegs) {
    // Only print for this chain
    if (id == m_decisionId or (isLegId(id) and m_decisionId == getIDFromLeg(id))) { 
      ATH_MSG_DEBUG("-- " << HLT::Identifier(id) << " with " << ELV.size() << " elements");
      for (const auto& EL : ELV) {
        ATH_MSG_DEBUG("-- -- container:" << EL.dataID() << ", index:" << EL.index());
      }
    }
  }
  return StatusCode::SUCCESS;
}

qOverPMatching()

float TrigMuonEFIdtpInvMassHypoTool::qOverPMatching ( const xAOD::TrackParticle *  metrack, const xAOD::TrackParticle *  idtrack  ) const

private

selectLegs()

StatusCode ComboHypoToolBase::selectLegs ( const Combo::LegDecisionsMap &  IDCombMap, std::vector< std::vector< Combo::LegDecision >> &  leg_decisions  ) const

protectedinherited

Creates the per-leg vectors of Decision objects starting from the initial LegDecision map, storing only those concerning this HypoTool's chain Pack the Decision objects in std::pair<DecisionID, ElementLink<Decision>> so the derived class' executeAlg function knows which leg each object is on.

Definition at line 172 of file ComboHypoToolBase.cxx.

{
  /*
    legDecisions: nested vector like:
    [(leg0, el0), (leg0, el1), (leg0, el2), (leg0, el3)], <-- All leg0 objects are in legDecisions[0][X]
    [(leg1, mu0), (leg1, mu1)] <-- All leg1 objects are in legDecisions[1][X]
    We keep the legID in the std::pair inside the inner vector as these pairs will be flattened into a single vector
    when individual combinations of objects are passed to executeAlg, pair.first then contains the leg, with pair.second containing the Decision Object.
  */
 
  // Extract from IDCombMap the features for this chain. Wrap the features up in a pair along with their leg ID
  for (size_t legIndex = 0; legIndex < m_legMultiplicities.size(); ++ legIndex) {
 
    // If the chain has more than one leg, then we have per-leg IDs. Otherwise we just use the chain's ID
    const HLT::Identifier& legIdentifier = (m_legMultiplicities.size() > 1 ? m_legDecisionIds.at(legIndex) : m_decisionId);
 
    // Combo::LegDecision is a pair of <DecisionID, ElementLink<Decision>>
    std::vector<Combo::LegDecision> decisionObjectsOnLeg;
 
    // Find physics objects on this leg. May be zero.
    const Combo::LegDecisionsMap::const_iterator it = IDCombMap.find(legIdentifier.numeric());
 
    if (it != IDCombMap.end()) {
      for (const ElementLink<DecisionContainer>& el : it->second) {
        decisionObjectsOnLeg.emplace_back(legIdentifier, el);
      }
    }
 
    legDecisions.push_back(std::move(decisionObjectsOnLeg));
  }
 
  if (msgLvl(MSG::DEBUG)) {
    ATH_MSG_DEBUG("Getting " << legDecisions.size() << " legs to combine, for ID: " << decisionId());
    size_t count = 0;
    for (const auto& leg : legDecisions) {
      ATH_MSG_DEBUG("Leg " << count++ << " --");
      for (const auto& dEL : leg) {
        ATH_MSG_DEBUG("-- " << HLT::Identifier(dEL.first) << " container:" << dEL.second.dataID() << ", index:" << dEL.second.index());
      }
    }
  }
  return StatusCode::SUCCESS;
}

setLegMultiplicity()

StatusCode ComboHypoToolBase::setLegMultiplicity ( const Combo::MultiplicityReqMap &  multiplicityRequiredMap )

inherited

Sets the number of legs and the multiplicity required on each leg.

This should be called when the Tool is retrieved by its parent ComboHypo alg. This also sets the leg Decision IDs at the same time param[in] multiplicityRequiredMap: Mapping of chains to required multiplicity per leg.

Definition at line 16 of file ComboHypoToolBase.cxx.

                                                                                                       {
  const std::string nameOfToolsChain = m_decisionId.name();
  const Combo::MultiplicityReqMap::const_iterator it = multiplicityRequiredMap.find(nameOfToolsChain);
 
  if (it == multiplicityRequiredMap.end()) {
    ATH_MSG_ERROR("ComboHypoTool for " << m_decisionId << " could not find its required multiplcity data in the map supplied by its parent alg.");
    return StatusCode::FAILURE;
  }
 
  m_legMultiplicities = it->second;
  if (m_legMultiplicities.size() == 0) {
    ATH_MSG_ERROR("ComboHypoTool for " << m_decisionId << " was listed in the supplied multiplicityRequiredMap, but data was not supplied for any legs.");
    return StatusCode::FAILURE;
  }
 
  m_legDecisionIds.clear();
  for (size_t i = 0; i < m_legMultiplicities.size(); ++i) {
    if(m_legMultiplicities.at(i) <= 0) {
      ATH_MSG_ERROR("ComboHypoTool for " << m_decisionId << " has been configured with an impossible multiplicity requirement of " << m_legMultiplicities.at(i) << " on leg " << i);
      return StatusCode::FAILURE;
    }
    ATH_MSG_DEBUG("ComboHypoTool for " << m_decisionId << " will require multiplicity " << m_legMultiplicities.at(i) << " on leg " << i);
    if (m_legMultiplicities.size() > 1) { // We only have per-leg IDs when there is more than one leg
      m_legDecisionIds.push_back(createLegName(m_decisionId, i));
    } else {  // For one leg, just repete the chain's ID
      m_legDecisionIds.push_back(m_decisionId);
    }
  }
 
  return StatusCode::SUCCESS;
}

updateLegDecisionsMap()

void ComboHypoToolBase::updateLegDecisionsMap ( const std::vector< std::vector< Combo::LegDecision >> &  passing_comb, Combo::LegDecisionsMap &  passingLegs  ) const

protectedinherited

For when the tool accepts some/all combinations.

Remove Decision Objects from legs of this HypoTool's chain which participated in NONE of combinations which were flagged as accepting the event.

Definition at line 216 of file ComboHypoToolBase.cxx.

                                                                                                                                                {
  if (msgLvl(MSG::DEBUG)) {
    size_t count = 0;
    for (const std::vector<Combo::LegDecision>& comb : passingComb) {
      ATH_MSG_DEBUG("-- Passing combination " << count++ << " of " << passingComb.size());
      for (const auto& [id, EL] : comb) {
        ATH_MSG_DEBUG("-- -- " << HLT::Identifier(id) << " container:" << EL.dataID() << ", index:" << EL.index());
      }
    }
  }
 
  // remove combinations that didn't pass from the final map passingLegs
  for (auto& it : passingLegs) {
    DecisionID legId = it.first;
    if (not(legId == m_decisionId or (isLegId(legId) and getIDFromLeg(legId) == m_decisionId))) {
      continue; // Some other chain, ignore it to get faster execution.
    } 
    std::vector<ElementLink<DecisionContainer>> updatedDecisionObjectsOnLeg;
    bool update = false;
    // Loop over all passing combinations, and all Decision Objects in each passing combination. Find Decision Objects on this leg.
    for (const std::vector<Combo::LegDecision>& comb : passingComb) {
      for (const auto& [id, EL] : comb) {
        // Check that this Decision Object is on the correct leg, and that we haven't collated it already from another combination.
        if (id == legId and std::find(updatedDecisionObjectsOnLeg.begin(), updatedDecisionObjectsOnLeg.end(), EL) == updatedDecisionObjectsOnLeg.end()) {
          ATH_MSG_VERBOSE("Keeping on leg " << HLT::Identifier(id) << " the Decision Object container:" << EL.dataID() << ", index:" << EL.index());
          updatedDecisionObjectsOnLeg.push_back(EL);
          update = true;
        }
      }
    }
    // only update those concerning this tool ID
    if (update){ 
      it.second = updatedDecisionObjectsOnLeg;
    }
  }
}

Member Data Documentation

m_acceptAll

Gaudi::Property< bool > TrigMuonEFIdtpInvMassHypoTool::m_acceptAll

private

Initial value:

{
      this, "AcceptAll", false, "Ignore selection" }

Definition at line 42 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_combinationsThresholdBreak

Gaudi::Property<size_t> ComboHypoToolBase::m_combinationsThresholdBreak

protectedinherited

Initial value:

{this, "CombinationsThresholdBreak", 10000,
    "Events processing this many combinations will generate a second WARNING message, and the loop over combinations will be terminated at this point."}

Definition at line 109 of file ComboHypoToolBase.h.

m_combinationsThresholdWarn

Gaudi::Property<size_t> ComboHypoToolBase::m_combinationsThresholdWarn

protectedinherited

Initial value:

{this, "CombinationsThresholdWarn", 1000,
    "Events processing this many combinations will generate a WARNING message."}

Definition at line 106 of file ComboHypoToolBase.h.

m_decisionId

HLT::Identifier ComboHypoToolBase::m_decisionId

privateinherited

The DecisionID of the chain, obtained from the Tool's name.

Definition at line 122 of file ComboHypoToolBase.h.

m_enableOverride

Gaudi::Property<bool> ComboHypoToolBase::m_enableOverride

protectedinherited

Initial value:

{this, "EnableOverride", false,
    "Stops processing combinations as soon as a valid combination is found in OR mode, or as soon as an invalid combination is found in AND mode. This is to save CPU."}

Definition at line 115 of file ComboHypoToolBase.h.

m_invMassHigh

Gaudi::Property< double> TrigMuonEFIdtpInvMassHypoTool::m_invMassHigh

private

Initial value:

{
      this, "InvMassHigh", -1.0, "High threshold for invariant mass cut" }

Definition at line 40 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_invMassLow

Gaudi::Property< double> TrigMuonEFIdtpInvMassHypoTool::m_invMassLow

private

Initial value:

{
      this, "InvMassLow", -1.0, "Low threshold for invariant mass cut" }

Definition at line 38 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_legDecisionIds

std::vector<HLT::Identifier> ComboHypoToolBase::m_legDecisionIds

privateinherited

The DecisionIDs of the individual legs, derived from both m_decisionId and m_legMultiplicities.

Definition at line 123 of file ComboHypoToolBase.h.

m_legMultiplicities

std::vector<int> ComboHypoToolBase::m_legMultiplicities

privateinherited

The number of legs, and the required multiplicity on each leg.

Definition at line 124 of file ComboHypoToolBase.h.

m_modeOR

Gaudi::Property<bool> ComboHypoToolBase::m_modeOR

protectedinherited

Initial value:

{this, "ModeOR", true,
    "Accepts based on the logical OR over all calls to executeAlg. If this flag is set to false then the logical AND is required instead."}

Definition at line 112 of file ComboHypoToolBase.h.

m_monTool

ToolHandle< GenericMonitoringTool > TrigMuonEFIdtpInvMassHypoTool::m_monTool {this, "MonTool", "", "Monitoring tool" }

private

Definition at line 50 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_muonqualityCut

Gaudi::Property< bool > TrigMuonEFIdtpInvMassHypoTool::m_muonqualityCut

private

Initial value:

{
      this, "MuonQualityCut", false, "Ignore selection" }

Definition at line 46 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_muonSelTool

ToolHandle<CP::IMuonSelectionTool> TrigMuonEFIdtpInvMassHypoTool::m_muonSelTool {this, "MuonSelectionTool", "CP::MuonSelectionTool/MuonSelectionTool", "Tool for muon quality selection"}

private

Definition at line 51 of file TrigMuonEFIdtpInvMassHypoTool.h.

m_selOS

Gaudi::Property< bool > TrigMuonEFIdtpInvMassHypoTool::m_selOS

private

Initial value:

{
      this, "SelectOppositeSign", false, "Select only oppositly charged pairs" }

Definition at line 44 of file TrigMuonEFIdtpInvMassHypoTool.h.

---

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

• TrigMuonEFIdtpInvMassHypoTool.h
• TrigMuonEFIdtpInvMassHypoTool.cxx