TrigmuCombHypoTool Class Reference

#include <TrigmuCombHypoTool.h>

Inheritance diagram for TrigmuCombHypoTool:

Collaboration diagram for TrigmuCombHypoTool:

Classes
struct	CombinedMuonInfo

Public Member Functions
	TrigmuCombHypoTool (const std::string &type, const std::string &name, const IInterface *parent)
virtual	~TrigmuCombHypoTool ()
virtual StatusCode	initialize () override
virtual StatusCode	decide (std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &toolInput) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	evtStore () const
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc. More...
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc. More...
virtual StatusCode	sysInitialize () override
	Perform system initialization for an algorithm. More...
virtual StatusCode	sysStart () override
	Handle START transition. More...
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles. More...
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles. More...
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T > &t)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKey &hndl, const std::string &doc, const SG::VarHandleKeyType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleBase &hndl, const std::string &doc, const SG::VarHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, SG::VarHandleKeyArray &hndArr, const std::string &doc, const SG::VarHandleKeyArrayType &)
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc, const SG::NotHandleType &)
	Declare a new Gaudi property. More...
Gaudi::Details::PropertyBase *	declareProperty (const std::string &name, T &property, const std::string &doc="none")
	Declare a new Gaudi property. More...
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
MsgStream &	msg (const MSG::Level lvl) const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution More...
std::enable_if_t< std::is_void_v< std::result_of_t< decltype(&T::renounce)(T)> > &&!std::is_base_of_v< SG::VarHandleKeyArray, T > &&std::is_base_of_v< Gaudi::DataHandle, T >, void >	renounce (T &h)
void	extraDeps_update_handler (Gaudi::Details::PropertyBase &ExtraDeps)
	Add StoreName to extra input/output deps as needed. More...

Private Types
enum	{ MaxNumberTools = 20 }
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
bool	decideOnSingleObject (TrigmuCombHypoTool::CombinedMuonInfo &input, size_t cutIndex) const
StatusCode	inclusiveSelection (std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &input) const
StatusCode	multiplicitySelection (std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &input) const
StatusCode	applyOverlapRemoval (std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &toolInput) const
StatusCode	checkOverlap (std::vector< TrigmuCombHypoTool::CombinedMuonInfo * > &input) const
bool	isOverlap (const xAOD::L2CombinedMuon *mf1, const xAOD::L2CombinedMuon *mf2) const
double	dR (double eta1, double phi1, double eta2, double phi2) const
double	invMass (double m1, double pt1, double eta1, double phi1, double m2, double pt2, double eta2, double phi2) const
StatusCode	chooseBestMuon (std::vector< TrigmuCombHypoTool::CombinedMuonInfo * > &input, const std::vector< unsigned int > &mucombResult) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleKeyArrayType &)
	specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &hndl, const SG::VarHandleType &)
	specialization for handling Gaudi::Property<SG::VarHandleBase> More...
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T > &t, const SG::NotHandleType &)
	specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...

Private Attributes
HLT::Identifier	m_decisionId
ToolHandle< GenericMonitoringTool >	m_monTool { this, "MonTool", "", "Monitoring tool" }
Gaudi::Property< std::vector< std::vector< double > > >	m_ptBins
Gaudi::Property< std::vector< std::vector< double > > >	m_ptThresholds
Gaudi::Property< bool >	m_acceptAll
Gaudi::Property< bool >	m_pikCuts
Gaudi::Property< double >	m_maxPtToApplyPik
Gaudi::Property< double >	m_chi2MaxID
Gaudi::Property< size_t >	m_numMuonThreshold
Gaudi::Property< bool >	m_strategydependent
Gaudi::Property< std::vector< double > >	m_strategyDependentPtCuts
Gaudi::Property< bool >	m_decisionPerCluster
Gaudi::Property< bool >	m_applyOR
Gaudi::Property< bool >	m_requireDR
Gaudi::Property< bool >	m_requireMufastDR
Gaudi::Property< bool >	m_requireMass
Gaudi::Property< bool >	m_requireSameSign
Gaudi::Property< std::vector< float > >	m_etaBins
Gaudi::Property< std::vector< float > >	m_dRThres
Gaudi::Property< std::vector< float > >	m_mufastDRThres
Gaudi::Property< std::vector< float > >	m_massThres
Gaudi::Property< float >	m_d0min
std::vector< size_t >	m_bins = {0}
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default) More...
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default) More...
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 31 of file TrigmuCombHypoTool.h.

Member Typedef Documentation

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< AlgTool > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Member Enumeration Documentation

anonymous enum

anonymous enum

private

Enumerator
MaxNumberTools

Definition at line 32 of file TrigmuCombHypoTool.h.

{ MaxNumberTools = 20 }; 

Constructor & Destructor Documentation

TrigmuCombHypoTool()

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

Definition at line 18 of file TrigmuCombHypoTool.cxx.

   : AthAlgTool( type, name, parent ),
     m_decisionId( HLT::Identifier::fromToolName( name ) ) 
{
}

~TrigmuCombHypoTool()

TrigmuCombHypoTool::~TrigmuCombHypoTool ( )

virtual

Definition at line 26 of file TrigmuCombHypoTool.cxx.

                                       {
}

Member Function Documentation

applyOverlapRemoval()

StatusCode TrigmuCombHypoTool::applyOverlapRemoval ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &  toolInput ) const

private

Definition at line 359 of file TrigmuCombHypoTool.cxx.

                                                                                                                 {
 
  ATH_MSG_DEBUG("Running Overlap Removal for muComb");
 
  std::vector<TrigmuCombHypoTool::CombinedMuonInfo*> input;
 
  // set pT threshold for events where so many muons
  // in such events, muons are removed if pT < pTthreshold
  // pT threshold is set to the pT value of m_numMuonThreshold-th leading muon
  float pTthreshold = 0;
  std::vector<float> pTvec;
  for ( auto& i: toolInput ) {
    if ( TrigCompositeUtils::passed( m_decisionId.numeric(), i.previousDecisionIDs) &&
     i.muComb!=nullptr ){
      pTvec.emplace_back(i.muComb->pt());
    }
  }
  if(pTvec.size() > m_numMuonThreshold) {
    std::sort(pTvec.begin(),pTvec.end(), std::greater<float>{});
    pTthreshold = pTvec.at(m_numMuonThreshold);
  }
 
  for ( auto& i: toolInput ) {
    if ( TrigCompositeUtils::passed( m_decisionId.numeric(), i.previousDecisionIDs) &&
     i.muComb!=nullptr ){
      if(i.muComb->pt() > pTthreshold)
    input.emplace_back(&i);
      else // set isOR for removed muons
    i.isOR.insert(m_decisionId.numeric());
    }
  }
 
  size_t numMuon = input.size();
 
  auto mucombNrAllEVs     = Monitored::Scalar("NrAllEVs", -9999.);
  auto mucombNrActiveEVs  = Monitored::Scalar("NrActiveEVs", -9999.);
  auto monitorIt          = Monitored::Group(m_monTool, mucombNrAllEVs, mucombNrActiveEVs);
  if ( numMuon == 0) {
    ATH_MSG_DEBUG( "No positive previous hypo decision. Not need overlap removal." );
    mucombNrActiveEVs = numMuon;
    mucombNrAllEVs = numMuon;
    return StatusCode::SUCCESS;
  }
  else if ( numMuon == 1 ) {
    ATH_MSG_DEBUG("Number of muon event = " << numMuon );
    ATH_MSG_DEBUG("no overlap Removal necessary. exitting with all EventViews active." );
    mucombNrActiveEVs = numMuon;
    mucombNrAllEVs = numMuon;
    return StatusCode::SUCCESS;
  } else {
    ATH_MSG_DEBUG("Number of muon event = " << numMuon );
    mucombNrAllEVs = numMuon;
    ATH_CHECK(checkOverlap(input));
    return StatusCode::SUCCESS;
  }
 
   return StatusCode::SUCCESS;
}

checkOverlap()

StatusCode TrigmuCombHypoTool::checkOverlap ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo * > &  input ) const

private

Definition at line 421 of file TrigmuCombHypoTool.cxx.

                                                                                                       {
 
  size_t numMuon = input.size();
  unsigned int i,j;
  std::vector<unsigned int> mucombResult;
 
  bool errorWhenIdentifyingOverlap = false;
 
  for(i=0; i<numMuon; i++) {mucombResult.emplace_back(i); }
  for(i=0; i<numMuon-1; i++){
    for(j=i+1; j<numMuon; j++){
      ATH_MSG_DEBUG("++ i=" << i << " vs j=" << j);
      bool overlapped = isOverlap((*input[i]).muComb, (*input[j]).muComb);
      if( ! overlapped ){ // judged as different
    ATH_MSG_DEBUG("   judged as: differenr objects");
    if( mucombResult[i] == mucombResult[j] ) { // but marked as same by someone
      ATH_MSG_DEBUG( "inconsistentency in muComb overlap removal for more than two objects" );
      ATH_MSG_DEBUG( "two objects are judged as different but both were already marked as identical by someone else as: " );
      ATH_MSG_DEBUG( "i/j/result[i]/result[j]=" << i << " / " << j << " / " << mucombResult[i] << " / "  << mucombResult[j] );
      auto mucombError  = Monitored::Scalar("MucombError", -9999.);
      auto monitorIt    = Monitored::Group(m_monTool, mucombError);
      mucombError = TrigmuCombHypoToolConsts::errorCode_inconsistent_overlap1;
      errorWhenIdentifyingOverlap = true;
    }
      }
      else{ // judged as overlap
    if( (mucombResult[j] != j && mucombResult[i] != mucombResult[j]) || (mucombResult[j] == j && mucombResult[i] != i) ){
      ATH_MSG_DEBUG( "inconsistentency in muComb based overlap removal for more than two objects" );
      ATH_MSG_DEBUG( "two objects are judged as overlap but only either was already marked as overlap to someone else: " );
      ATH_MSG_DEBUG( "i/j/result[i]/result[j]=" << i << " / " << j << " / " << mucombResult[i] << " / "  << mucombResult[j] );
      auto mucombError  = Monitored::Scalar("MucombError", -9999.);
      auto monitorIt    = Monitored::Group(m_monTool, mucombError);
      mucombError = TrigmuCombHypoToolConsts::errorCode_inconsistent_overlap2;
      errorWhenIdentifyingOverlap = true;
    }
    ATH_MSG_DEBUG("   judged as: overlapped objects");
    if( mucombResult[i] == i ) {
      ATH_MSG_DEBUG( "   i is not yet marked as overlap. so, it is a newly found overlap" );
      ATH_MSG_DEBUG( "   -> marking mucombResult[j] as i..." );
      mucombResult[j] = i;
    } else {
      ATH_MSG_DEBUG( "   both i/j already marked as overlap by: mucombResult[i]=" << mucombResult[i] );
      ATH_MSG_DEBUG( "   -> do nothing..." );
    }
      }
    }
  }
 
  if( errorWhenIdentifyingOverlap ) {
    ATH_MSG_WARNING( "error when resolving overlap. exitting with all EVs active..." );
    auto mucombNrActiveEVs  = Monitored::Scalar("NrActiveEVs", -9999.);
    auto monitorIt          = Monitored::Group(m_monTool, mucombNrActiveEVs);
    mucombNrActiveEVs = numMuon;
    // for(i=0; i<numMuon; i++) TrigCompositeUtils::addDecisionID( m_decisionId, toolInput[i].decision );
    return StatusCode::SUCCESS;
  }
 
  unsigned int n_uniqueMuon = 0;
  for(i=0; i<numMuon; i++) {
    ATH_MSG_DEBUG( "muComb based results: i=" << i << ": ");
    if( mucombResult[i] != i ) { ATH_MSG_DEBUG( "      overlap to j=" << mucombResult[i] ); }
    else {
      n_uniqueMuon++;
      ATH_MSG_DEBUG( "      unique" );
    }
  }
 
  ATH_MSG_DEBUG( "nr of unique Muons after muComb-based removal=" << n_uniqueMuon );
 
  if( numMuon != n_uniqueMuon ){
    ATH_CHECK(chooseBestMuon(input, mucombResult));
  } else {
    ATH_MSG_DEBUG( "no overlap identified. exitting with all EventViews active" );
    auto mucombNrActiveEVs  = Monitored::Scalar("NrActiveEVs", -9999.);
    auto monitorIt          = Monitored::Group(m_monTool, mucombNrActiveEVs);
    mucombNrActiveEVs = n_uniqueMuon;
  }
 
  return StatusCode::SUCCESS;
}

chooseBestMuon()

StatusCode TrigmuCombHypoTool::chooseBestMuon ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo * > &  input, const std::vector< unsigned int > &  mucombResult  ) const

private

Definition at line 674 of file TrigmuCombHypoTool.cxx.

{
  const double ZERO_LIMIT = 1e-4;
  size_t numMuon = input.size();
  unsigned int i,j,k;
 
  auto mucombNrActiveEVs    = Monitored::Scalar("NrActiveEVs", -9999.);
  auto mucombNrOverlapped   = Monitored::Scalar("NrOverlapped", 0);
  auto mucombOverlappedEta  = Monitored::Scalar("OverlappedEta", -9999.);
  auto mucombOverlappedPhi  = Monitored::Scalar("OverlappedPhi", -9999.);
  auto mucombOverlappedPt   = Monitored::Scalar("OverlappedPt", -9999.);
 
  auto monitorIt       = Monitored::Group(m_monTool, mucombNrActiveEVs, mucombNrOverlapped,
                      mucombOverlappedPt, mucombOverlappedEta, mucombOverlappedPhi);
 
  ATH_MSG_DEBUG( "--- choose best among overlaps & disable EVs (muComb based) ---" );
  for(i=0; i<numMuon; i++) {
    ATH_MSG_DEBUG( "++ i=" << i << ": result=" << mucombResult[i] );
    if( mucombResult[i] != i ) {
      ATH_MSG_DEBUG( "   overlap to some one. already the best one was chosen. skip." );
      continue;
    }
    std::vector<unsigned int> others;
    for(j=0; j<numMuon; j++) {
      if( mucombResult[j] == mucombResult[i] ) others.emplace_back(j);
    }
    if( others.size() == 1 ) {
      ATH_MSG_DEBUG( "   unique object. keep it active." );
      continue;
    }
    else {// must choose one best
      ATH_MSG_DEBUG( "   overlapped objects among: " << others );
      unsigned int best_ev = 0;
      float maxPtCombMf  = 0;
      float mindRRoadRoI = 999;
      for(k=0; k<others.size(); k++) {
    j=others[k];
 
    float ptCombMf  = 0.;
    const xAOD::L2CombinedMuon* combMf = (*input[j]).muComb;
    ptCombMf  = fabs(combMf->pt()/Gaudi::Units::GeV);
 
    const xAOD::L2StandAloneMuon* mf = (*input[j]).muComb->muSATrack();
    const float roadPhiP = atan2(mf->dirPhiMS(),1.);
    const float roadPhiM = atan2(-1*mf->dirPhiMS(),-1.);
    const float roadPhi = (std::abs(CxxUtils::deltaPhi(roadPhiP, mf->roiPhi())) < std::abs(CxxUtils::deltaPhi(roadPhiM, mf->roiPhi())))? roadPhiP : roadPhiM;
    float roadAw = 0;
    if(std::abs(mf->roiEta()) < 1.05) { // barrel
      if( std::abs(mf->roadAw(1,0)) > ZERO_LIMIT )      roadAw = mf->roadAw(1,0);
      else if( std::abs(mf->roadAw(2,0)) > ZERO_LIMIT ) roadAw = mf->roadAw(2,0);
      else if( std::abs(mf->roadAw(0,0)) > ZERO_LIMIT ) roadAw = mf->roadAw(0,0);
    }
    else { // endcap
      if( std::abs(mf->roadAw(4,0)) > ZERO_LIMIT )      roadAw = mf->roadAw(4,0);
      else if( std::abs(mf->roadAw(5,0)) > ZERO_LIMIT ) roadAw = mf->roadAw(5,0);
      else if( std::abs(mf->roadAw(3,0)) > ZERO_LIMIT ) roadAw = mf->roadAw(3,0);
    }
    float roadEta = 999;
    if(std::abs(roadAw) > ZERO_LIMIT)
      roadEta = -std::log(std::tan(0.5*std::atan(std::abs(roadAw))));
    if(roadAw < 0) roadEta *= -1.;
    const double dRRoadRoI = dR(roadEta, roadPhi, mf->roiEta(), mf->roiPhi());
    ATH_MSG_DEBUG("     j="<< j << " , ptCombMf=" << ptCombMf << ", dRRoadRoI=" << dRRoadRoI);
 
    if( (ptCombMf > maxPtCombMf) ||
        (std::abs(ptCombMf - maxPtCombMf) < ZERO_LIMIT &&
         dRRoadRoI < mindRRoadRoI) ) {
      maxPtCombMf  = ptCombMf;
      mindRRoadRoI = dRRoadRoI;
      best_ev  = j;
    }
      }
      ATH_MSG_DEBUG( "      best is: best_ev/maxPtCombMf=" << best_ev << " / " <<  maxPtCombMf );
 
      for(k=0; k<others.size(); k++) {
    j=others[k];
    if( j != best_ev ) {
      ATH_MSG_DEBUG( "      EventView( j=" << j << " ) is not active" );
 
      (*input[j]).isOR.insert(m_decisionId.numeric());
 
      // monitoring
      const xAOD::L2CombinedMuon* CombMf = (*input[j]).muComb;
      ++mucombNrOverlapped;
      mucombOverlappedPt = CombMf->pt()* CombMf->charge() /Gaudi::Units::GeV;
      mucombOverlappedEta = CombMf->eta();
      mucombOverlappedPhi = CombMf->phi();
    }
    if( j == best_ev ){
      ATH_MSG_DEBUG( "      EventView( j=" << j << " ) is best one" );
    }
      }
    }
  }
  mucombNrActiveEVs = numMuon - mucombNrOverlapped;
 
  return StatusCode::SUCCESS;
}

decide()

StatusCode TrigmuCombHypoTool::decide ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &  toolInput ) const

virtual

Definition at line 246 of file TrigmuCombHypoTool.cxx.

{
   size_t numTrigger = m_ptBins.size();
   size_t numMuon = toolInput.size();
   ATH_MSG_DEBUG("Retrieved from TrigmuCombHypoAlg and Running TrigmuCombHypoTool for selections.");
  
   if ( numTrigger == 1 ) { // in case of HLT_mu4, HLT_mu6 and so on.
      ATH_MSG_DEBUG("Number of muon event = " << numMuon );
      ATH_MSG_DEBUG("Applying selection of single << " << m_decisionId );
      return inclusiveSelection(toolInput); 
   } else {         // in case of HLT_2mu6 and so on.
     ATH_MSG_DEBUG("Number of muon event = " << numMuon );
     ATH_MSG_DEBUG("Applying selection of multiplicity << " << m_decisionId );
 
     if(m_applyOR)
       ATH_CHECK(applyOverlapRemoval(toolInput));
 
     return multiplicitySelection(toolInput);
   }
 
   return StatusCode::SUCCESS;
}

decideOnSingleObject()

bool TrigmuCombHypoTool::decideOnSingleObject ( TrigmuCombHypoTool::CombinedMuonInfo &  input, size_t  cutIndex  ) const

private

Definition at line 105 of file TrigmuCombHypoTool.cxx.

{
 
   // defined Monitoring variables
   auto fexPt   = Monitored::Scalar("Pt", -9999.);
   auto ptFL    = Monitored::Scalar("PtFL", -9999.);
   auto Strategy    = Monitored::Scalar("StrategyFlag", 0);
   auto idEta   = Monitored::Scalar("Eta", -9999.);
   auto idPhi   = Monitored::Scalar("Phi", -9999.);
   auto idZ0    = Monitored::Scalar("Z0", -9999.);
   auto idA0    = Monitored::Scalar("A0", -9999.);
 
   auto monitorIt = Monitored::Group(m_monTool, fexPt, ptFL, Strategy, 
                       idEta, idPhi, idZ0, idA0);
 
   ATH_MSG_VERBOSE( "Cut index " << cutIndex );
 
   // flag to pass or not
   bool result = false;
 
   if(m_acceptAll) {
      result = true;
      ATH_MSG_DEBUG("Accept property is set: taking all the events");
      return result;
   } else {
      result = false;
      ATH_MSG_DEBUG("Accept property not set: applying selection!");
   }
 
   //Retrieve combined muon
   //xAOD::L2CombinedMuon
   auto pMuon = input.muComb;
   if (!pMuon) {
      result = false;
      ATH_MSG_ERROR("Retrieval of xAOD::L2CombinedMuon from vector failed");
      return result;
   }
 
   auto ptValue = pMuon->pt() * pMuon->charge() / Gaudi::Units::GeV;
 
 
   if (pMuon->pt() == 0) {
      ATH_MSG_DEBUG("L2CombinedMuon pt == 0, empty container -> rejected");
      result = false;
      return result;
   }
 
   fexPt    = ptValue;
   ptFL      = ptValue;
   idEta    = pMuon->eta();
   idPhi    = pMuon->phi();
   int usealgo      = pMuon->strategy();
   float ptresComb = pMuon->sigmaPt() / Gaudi::Units::GeV;
   Strategy  = usealgo;
   ATH_MSG_DEBUG("combined muon pt (GeV)/ sigma_pt (GeV)/ eta / phi / usedalgo: " 
               << fexPt << " (GeV) / " << ptresComb << " (GeV) / " << idEta << " / " << idPhi 
               << " / " << usealgo);
 
   // check the pointers to the L2StandAloneMuon
   if (!pMuon->muSATrack()) {
      ATH_MSG_DEBUG("L2CombinedMuon has no valid xaOD::L2StandaloneMuon -> rejected");
      result = false;
      return result;
   }
 
   // check the pointer to the ID track
   if (!pMuon->idTrack()) {
      ATH_MSG_DEBUG("L2CombinedMuon has no valid xAOD:TrackParticle IDtrack -> rejected");
      result = false;
      return result;
   }
 
   idA0     = pMuon->idTrack()->d0();
   idZ0     = pMuon->idTrack()->z0();
 
   //Get the Pt cut for that eta bin
   float threshold = (idEta != -9999) ? 0 : 999999;
   float absEta = fabs(idEta);
   for (std::vector<float>::size_type i = 0; i < m_bins[cutIndex]; ++i) {
      if (absEta > m_ptBins[cutIndex][i] && absEta < m_ptBins[cutIndex][i + 1]) threshold = m_ptThresholds[cutIndex][i];
   }
 
   // Check pt threshold for hypothesis and pi/k rejection cuts,
   // convert units since Muonfeature is in GeV
 
   //Kpi rejection
   bool pikCut = true;
   if (m_pikCuts && (std::abs(fexPt) < m_maxPtToApplyPik)) {
      if (pMuon->idTrack()->chiSquared() > m_chi2MaxID) pikCut = false;
   }
 
   //Std Pt cut
   bool stdCut = true;
   if (std::abs(fexPt) <= (threshold / Gaudi::Units::GeV)) stdCut = false;
   ATH_MSG_DEBUG("REGTEST muon pt is " << fexPt
              << " GeV and threshold cut is " << threshold / Gaudi::Units::GeV
              << " GeV and pik_cut is " << (pikCut ? "true" : "false"));
 
   //Strategy dependent Pt cuts
   bool sdpCut = true;
   if (m_strategydependent && usealgo > 0) {
      if (usealgo >= 1 && usealgo <= 4) {
         double tmpcut = m_strategyDependentPtCuts.value()[usealgo - 1];
         if (std::abs(fexPt) <= std::abs(tmpcut)) sdpCut = false;
         if (tmpcut < 0) stdCut = true; //Do not apply std Pt cut
         ATH_MSG_DEBUG("REGTEST muon pt is " << fexPt << " GeV"
                    << " and threshold for strategy dependent cut is " << tmpcut
                    << " GeV and strategy dependent / std cuts are " << (sdpCut ? "true" : "false") << " / " << (stdCut ? "true" : "false"));
      } else {
         ATH_MSG_DEBUG("usealgo out of range, is: " << usealgo << " while should be in [1, 4]");
      }
   }
 
 
   //d0 cut
   bool d0Cut = true;
   if (m_d0min>0. && std::abs(idA0)<m_d0min) d0Cut = false;
 
   result = stdCut && pikCut && sdpCut && d0Cut;
 
   if (result) ptFL = -9999.;
 
   if (m_d0min>0.) {
     ATH_MSG_DEBUG("REGTEST: Muon passed pt threshold: " << (stdCut ? "true" : "false")
              << " and pik_cut is " << (pikCut ? "true" : "false")
              << " and strategy dependent cuts is " << (sdpCut ? "true" : "false")
          << " and result of d0min cut is "<< (d0Cut ? "true" : "false")
              << " so hypothesis is " << (result ? "true" : "false"));
   } else {
     ATH_MSG_DEBUG("REGTEST: Muon passed pt threshold: " << (stdCut ? "true" : "false")
              << " and pik_cut is " << (pikCut ? "true" : "false")
              << " and strategy dependent cuts is " << (sdpCut ? "true" : "false")
              << " so hypothesis is " << (result ? "true" : "false"));
   }
   return result;
}

declareGaudiProperty() [1/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyArrayType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKeyArray>

Definition at line 170 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [2/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleKeyType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleKey>

Definition at line 156 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
 
  }

declareGaudiProperty() [3/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  hndl, const SG::VarHandleType &    )

inlineprivateinherited

specialization for handling Gaudi::Property<SG::VarHandleBase>

Definition at line 184 of file AthCommonDataStore.h.

  {
    return *AthCommonDataStore<PBASE>::declareProperty(hndl.name(),
                                                       hndl.value(), 
                                                       hndl.documentation());
  }

declareGaudiProperty() [4/4]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareGaudiProperty ( Gaudi::Property< T > &  t, const SG::NotHandleType &    )

inlineprivateinherited

specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray>

Definition at line 199 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(t);
  }

declareProperty() [1/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleBase &  hndl, const std::string &  doc, const SG::VarHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleBase. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 245 of file AthCommonDataStore.h.

  {
    this->declare(hndl.vhKey());
    hndl.vhKey().setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [2/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleKey &  hndl, const std::string &  doc, const SG::VarHandleKeyType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
hndl	Object holding the property value.
doc	Documentation string for the property.

This is the version for types that derive from SG::VarHandleKey. The property value object is put on the input and output lists as appropriate; then we forward to the base class.

Definition at line 221 of file AthCommonDataStore.h.

  {
    this->declare(hndl);
    hndl.setOwner(this);
 
    return PBASE::declareProperty(name,hndl,doc);
  }

declareProperty() [3/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, SG::VarHandleKeyArray &  hndArr, const std::string &  doc, const SG::VarHandleKeyArrayType &    )

inlineinherited

Definition at line 259 of file AthCommonDataStore.h.

  {
 
    // std::ostringstream ost;
    // ost << Algorithm::name() << " VHKA declareProp: " << name 
    //     << " size: " << hndArr.keys().size() 
    //     << " mode: " << hndArr.mode() 
    //     << "  vhka size: " << m_vhka.size()
    //     << "\n";
    // debug() << ost.str() << endmsg;
 
    hndArr.setOwner(this);
    m_vhka.push_back(&hndArr);
 
    Gaudi::Details::PropertyBase* p =  PBASE::declareProperty(name, hndArr, doc);
    if (p != 0) {
      p->declareUpdateHandler(&AthCommonDataStore<PBASE>::updateVHKA, this);
    } else {
      ATH_MSG_ERROR("unable to call declareProperty on VarHandleKeyArray " 
                    << name);
    }
 
    return p;
 
  }

declareProperty() [4/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc, const SG::NotHandleType &    )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This is the generic version, for types that do not derive from SG::VarHandleKey. It just forwards to the base class version of declareProperty.

Definition at line 333 of file AthCommonDataStore.h.

  {
    return PBASE::declareProperty(name, property, doc);
  }

declareProperty() [5/6]

Gaudi::Details::PropertyBase* AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( const std::string &  name, T &  property, const std::string &  doc = "none"  )

inlineinherited

Declare a new Gaudi property.

Parameters

name	Name of the property.
property	Object holding the property value.
doc	Documentation string for the property.

This dispatches to either the generic declareProperty or the one for VarHandle/Key/KeyArray.

Definition at line 352 of file AthCommonDataStore.h.

  {
    typedef typename SG::HandleClassifier<T>::type htype;
    return declareProperty (name, property, doc, htype());
  }

declareProperty() [6/6]

Gaudi::Details::PropertyBase& AthCommonDataStore< AthCommonMsg< AlgTool > >::declareProperty ( Gaudi::Property< T > &  t )

inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

                                                                 {
    typedef typename SG::HandleClassifier<T>::type htype;
    return AthCommonDataStore<PBASE>::declareGaudiProperty(t, htype());
  }

detStore()

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::detStore ( ) const

inlineinherited

The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

dR()

double TrigmuCombHypoTool::dR ( double  eta1, double  phi1, double  eta2, double  phi2  ) const

private

Definition at line 628 of file TrigmuCombHypoTool.cxx.

{
  const double deta = eta1 - eta2;
  const double dphi = CxxUtils::deltaPhi(phi1, phi2);
  return std::sqrt(deta*deta + dphi*dphi);
}

evtStore() [1/2]

ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( )

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

evtStore() [2/2]

const ServiceHandle<StoreGateSvc>& AthCommonDataStore< AthCommonMsg< AlgTool > >::evtStore ( ) const

inlineinherited

The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.

Definition at line 90 of file AthCommonDataStore.h.

{ return m_evtStore; }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::extraDeps_update_handler ( Gaudi::Details::PropertyBase &  ExtraDeps )

protectedinherited

Add StoreName to extra input/output deps as needed.

use the logic of the VarHandleKey to parse the DataObjID keys supplied via the ExtraInputs and ExtraOuputs Properties to add the StoreName if it's not explicitly given

inclusiveSelection()

StatusCode TrigmuCombHypoTool::inclusiveSelection ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &  input ) const

private

Definition at line 270 of file TrigmuCombHypoTool.cxx.

{
   for ( auto& i: input) {
     // If muon event has difference DecisionID, it shouldn't apply.
     if (TrigCompositeUtils::passed(m_decisionId.numeric(), i.previousDecisionIDs)) {
       if ( decideOnSingleObject(i, 0)==true ) {
         ATH_MSG_DEBUG("Pass through selection " << m_decisionId );
         TrigCompositeUtils::addDecisionID(m_decisionId, i.decision);
       } else {
         ATH_MSG_DEBUG("Not pass through selection " << m_decisionId );
       } 
     } else {
       ATH_MSG_DEBUG("Not match DecisionID:" << m_decisionId );
     } 
   }  
 
   return StatusCode::SUCCESS;
}

initialize()

StatusCode TrigmuCombHypoTool::initialize ( )

overridevirtual

Definition at line 32 of file TrigmuCombHypoTool.cxx.

{
   ATH_MSG_DEBUG( "Tool configured for chain/id: " << m_decisionId );
 
   if (m_acceptAll) {
      ATH_MSG_DEBUG("Accepting all the events!");
   } else {
      ATH_MSG_DEBUG("AcceptAll = False");
      m_bins.resize (m_ptBins.size());
      for ( size_t j=0; j<m_ptBins.size(); ++j) {
         m_bins[j] = m_ptBins[j].size() - 1;
         if (m_bins[j] != m_ptThresholds[j].size()) {
            ATH_MSG_ERROR("bad thresholds setup .... exiting!");
            return StatusCode::FAILURE;
         }
 
         for (std::vector<float>::size_type i = 0; i < m_bins[j]; ++i) {
 
            ATH_MSG_DEBUG("bin[" << j << "] " << m_ptBins[j][i] << " - " <<  m_ptBins[j][i + 1]
                  << " with Pt Threshold of " << (m_ptThresholds[j][i]) / Gaudi::Units::GeV << " GeV");
         }
      }  
   }
 
   if ( not m_monTool.name().empty() ) {
      ATH_CHECK( m_monTool.retrieve() );
      ATH_MSG_DEBUG("MonTool name: " << m_monTool);
   }
 
   // Overlap Removal
   if( m_applyOR ) {
     ATH_MSG_DEBUG( "--- overlap removal as: ---"     );
     if( m_requireDR ) {
       ATH_MSG_DEBUG( "+ dR cut:" );
       if( (m_etaBins.size()-1) != m_dRThres.size() ) {
     ATH_MSG_ERROR( "bad thresholds setup .... exiting!" );
     return StatusCode::FAILURE;
       }
       for(unsigned int i=0; i<m_dRThres.size(); i++) {
     ATH_MSG_DEBUG( "     EtaBin " << m_etaBins[i] << " - " << m_etaBins[i+1]
            << " : dR < " << m_dRThres[i] );
       }
     }
     if( m_requireMufastDR ) {
       ATH_MSG_DEBUG( "+ dr(by mF) cut:" );
       if( (m_etaBins.size()-1) != m_mufastDRThres.size() ) {
     ATH_MSG_ERROR( "bad thresholds setup .... exiting!" );
     return StatusCode::FAILURE;
       }
       for(unsigned int i=0; i<m_mufastDRThres.size(); i++) {
     ATH_MSG_DEBUG( "     EtaBin " << m_etaBins[i] << " - " << m_etaBins[i+1]
            << " : dR(mF) < " << m_mufastDRThres[i] );
       }
     }
     if( m_requireMass ) {
       ATH_MSG_DEBUG( "+ Mass cut:" );
       if( (m_etaBins.size()-1) != m_massThres.size() ) {
     ATH_MSG_ERROR( "bad thresholds setup .... exiting!" );
     return StatusCode::FAILURE;
       }
       for(unsigned int i=0; i<m_massThres.size(); i++) {
     ATH_MSG_DEBUG( "     EtaBin " << m_etaBins[i] << " - " << m_etaBins[i+1]
            << " : Mass < " << m_massThres[i] );
       }
     }
     if( m_requireSameSign ) ATH_MSG_DEBUG( "+ Same charge sign" );
   }
   // minimum d0 cut for displaced muon triggers
   if (m_d0min>0.) ATH_MSG_DEBUG( " Rejecting muons with abs(d0) < "<<m_d0min<<" mm");
 
   return StatusCode::SUCCESS;
}

inputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::inputHandles ( ) const

overridevirtualinherited

Return this algorithm's input handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

invMass()

double TrigmuCombHypoTool::invMass ( double  m1, double  pt1, double  eta1, double  phi1, double  m2, double  pt2, double  eta2, double  phi2  ) const

private

Definition at line 638 of file TrigmuCombHypoTool.cxx.

{
  const double ZERO_LIMIT = 1e-12;
 
  double theta1 = 2*atan2((double)exp(-eta1),1.);
  double theta2 = 2*atan2((double)exp(-eta2),1.);
 
  double fpt1   = fabs(pt1);
  double fpt2   = fabs(pt2);
 
  double px1    = fpt1*cos(phi1);
  double py1    = fpt1*sin(phi1);
  double pz1    = fpt1/tan(theta1);
  double  e1    = sqrt(px1*px1+py1*py1+pz1*pz1+m1*m1);
 
  double px2    = fpt2*cos(phi2);
  double py2    = fpt2*sin(phi2);
  double pz2    = fpt2/tan(theta2);
  double  e2    = sqrt(px2*px2+py2*py2+pz2*pz2+m2*m2);
 
  double pxsum  = px1 + px2;
  double pysum  = py1 + py2;
  double pzsum  = pz1 + pz2;
  double esum   =  e1 +  e2;
 
  double mass  = 0;
  double mass2 = esum*esum - pxsum*pxsum - pysum*pysum - pzsum*pzsum;
  if( mass2 > ZERO_LIMIT ) mass = sqrt(mass2);
 
  return mass;
}

isOverlap()

bool TrigmuCombHypoTool::isOverlap ( const xAOD::L2CombinedMuon *  mf1, const xAOD::L2CombinedMuon *  mf2  ) const

private

Definition at line 505 of file TrigmuCombHypoTool.cxx.

{
 
  auto mucombDR             = Monitored::Scalar("DR", -9999.);
  auto mucombMass           = Monitored::Scalar("Mass", -9999.);
  auto mucombDRLog10        = Monitored::Scalar("DRLog10", -9999.);
  auto mucombMassLog10      = Monitored::Scalar("MassLog10", -9999.);
 
  auto monitorIt       = Monitored::Group(m_monTool, mucombDR, mucombMass, mucombDRLog10, mucombMassLog10);
 
 
  ATH_MSG_DEBUG( "   ...mF1: pt/eta/phi=" << combMf1->pt()/Gaudi::Units::GeV << " / " << combMf1->eta() << " / " << combMf1->phi() );
  ATH_MSG_DEBUG( "   ...mF2: pt/eta/phi=" << combMf2->pt()/Gaudi::Units::GeV << " / " << combMf2->eta() << " / " << combMf2->phi() );
 
  // if dR or invMass is necessary but (eta,phi) info is not avaiable
  // (i.e. eta,phi=0,0; rec failed)
  const double ZERO_LIMIT_FOR_ETAPHI = 1e-4;
  if( (fabs(combMf1->eta()) <ZERO_LIMIT_FOR_ETAPHI && fabs(combMf1->phi()) < ZERO_LIMIT_FOR_ETAPHI) ||
      (fabs(combMf2->eta()) <ZERO_LIMIT_FOR_ETAPHI && fabs(combMf2->phi()) < ZERO_LIMIT_FOR_ETAPHI) ) {
    ATH_MSG_DEBUG( "   ...-> (eta,phi) info not available (rec at (eta,phi)=(0,0))" );
    if( m_requireDR || m_requireMass ) {
      ATH_MSG_DEBUG( "   ...-> but dR of invMass check is required. cannot judge overlap -> return with false" );
      return false;
    }
  }
 
  // if charge or invMass is necessary but charge(=pT) info is not avaiable
  const double ZERO_LIMIT_FOR_PT = 1e-4;
  if( (fabs(combMf1->pt()) <ZERO_LIMIT_FOR_PT) || (fabs(combMf2->pt()) < ZERO_LIMIT_FOR_PT) ) {
    ATH_MSG_DEBUG( "   ...-> pT info not available (rec at pT=0)" );
    if( m_requireSameSign || m_requireMass ) {
      ATH_MSG_DEBUG( "   ...-> but same sign or invMass check is required. cannot judge overlap -> return with false" );
      return false;
    }
  }
 
  // determine etabin and thresholds
  double absEta = (fabs(combMf1->pt()) > fabs(combMf2->pt())) ? fabs(combMf1->eta()) : fabs(combMf2->eta());
  unsigned int iThres = 0;
  for(unsigned int i=0; i<(m_etaBins.size()-1); i++) {
    if ( m_etaBins[i] <= absEta && absEta < m_etaBins[i+1] ) iThres = i;
  }
  double dRThres     = m_requireDR ? m_dRThres[iThres] : 0.;
  double dRbyMFThres = m_requireMufastDR ? m_mufastDRThres[iThres] : 0.;
  double massThres   = m_requireMass ? m_massThres[iThres] : 0.;
  ATH_MSG_DEBUG( "   ...iThres=" << iThres );
  if(m_requireDR)        ATH_MSG_DEBUG( "   ...dR       threshold=" << dRThres     );
  if(m_requireMufastDR)  ATH_MSG_DEBUG( "   ...dR(byMF) threshold=" << dRbyMFThres );
  if(m_requireMass)      ATH_MSG_DEBUG( "   ...mass     threshold=" << massThres   );
 
  // same sign cut
  bool sameSign = false;
  if( m_requireSameSign ) {
    sameSign = ((combMf1->pt()*combMf2->pt()) > 0) ? true : false;
    ATH_MSG_DEBUG( "   ...-> sameSign=" << sameSign );
  }
 
  // dR cut
  bool dRisClose = false;
  double dr = dR(combMf1->eta(),combMf1->phi(),combMf2->eta(),combMf2->phi());
 
  mucombDR = dr;
  const double monitor_limit = 1e-4;
  double dr_mon = (dr>=monitor_limit) ? dr : monitor_limit;
  mucombDRLog10 = log10(dr_mon);
 
  if( m_requireDR ) {
    if( dr < dRThres ) dRisClose = true;
    ATH_MSG_DEBUG( "   ...-> dR=" << dr << " : dRisClose=" << dRisClose );
  }
 
  // dR(by MF) cut
  bool dRbyMFisClose = false;
  if( m_requireMufastDR ) {
    const xAOD::L2StandAloneMuon* mf1 = combMf1->muSATrack();
    const xAOD::L2StandAloneMuon* mf2 = combMf2->muSATrack();
    if( mf1 == 0 || mf2 == 0 ) {
      ATH_MSG_DEBUG( "mF link from combinedMF broken" );
      ATH_MSG_DEBUG( "   ...-> mF dR is required but mF link broken. cannot judge overlap -> return with false" );
      return false;
    }
    else {
      // here, we do not check (eta,phi) of mF is not (0,0)
      // (i.e. we apply muComb based cut even if muFast rec is failed)
      double dRByMF = dR(mf1->etaMS(),mf1->phiMS(),mf2->etaMS(),mf2->phiMS());
      if( dRByMF < dRbyMFThres ) dRbyMFisClose = true;
      ATH_MSG_DEBUG( "   ...-> dR(by MF)=" << dRByMF << " : dRbyMFisClose=" << dRbyMFisClose );
    }
  }
 
  // mass cut
  const double TRACK_MASS = 0;  // just assume zero mass
  bool massIsClose = false;
  double mass = invMass(TRACK_MASS,combMf1->pt()/Gaudi::Units::GeV,combMf1->eta(),combMf1->phi(),TRACK_MASS,combMf2->pt()/Gaudi::Units::GeV,combMf2->eta(),combMf2->phi());
 
  mucombMass = mass;
  double mass_mon = (mass>=monitor_limit) ? mass : monitor_limit;
  mucombMassLog10 = log10(mass_mon);
 
  if( m_requireMass ) {
    if( mass < massThres ) massIsClose = true;
    ATH_MSG_DEBUG( "   ...-> mass=" << mass << " : massIsClose=" << massIsClose );
  }
 
  // total judge
  bool overlap = false;
  if( ((m_requireSameSign &&   sameSign)     || (! m_requireSameSign)) &&
      ((m_requireDR       &&  dRisClose)     || (! m_requireDR))       &&
      ((m_requireMufastDR &&  dRbyMFisClose) || (! m_requireMufastDR)) &&
      ((m_requireMass     &&  massIsClose)   || (! m_requireMass)) ) {
    overlap = true;
  }
 
  ATH_MSG_DEBUG( "   ...=> isOverlap=" << overlap );
 
  return overlap;
 
}

msg() [1/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msg() [2/2]

MsgStream& AthCommonMsg< AlgTool >::msg ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 27 of file AthCommonMsg.h.

                                                  {
    return this->msgStream(lvl);
  }

msgLvl()

bool AthCommonMsg< AlgTool >::msgLvl ( const MSG::Level  lvl ) const

inlineinherited

Definition at line 30 of file AthCommonMsg.h.

                                               {
    return this->msgLevel(lvl);
  }

multiplicitySelection()

StatusCode TrigmuCombHypoTool::multiplicitySelection ( std::vector< TrigmuCombHypoTool::CombinedMuonInfo > &  input ) const

private

Definition at line 290 of file TrigmuCombHypoTool.cxx.

{
   HLT::Index2DVec passingSelection( m_ptBins.size() );
 
   for ( size_t cutIndex=0; cutIndex < m_ptBins.size(); ++cutIndex ) {
      size_t elementIndex{ 0 };
      for ( auto& i: input ) {
 
    if(!m_acceptAll && m_applyOR && (i.isOR.find(m_decisionId.numeric()) != i.isOR.end())) {
      ATH_MSG_DEBUG("skip due to overap, DecisionID " << m_decisionId );
      elementIndex++;
      continue;
    }
 
         // If muon event has difference DecisionID, it shouldn't apply.
         if ( TrigCompositeUtils::passed( m_decisionId.numeric(), i.previousDecisionIDs ) ) {
            if ( decideOnSingleObject( i, cutIndex ) == true ) {
               ATH_MSG_DEBUG("Pass through selection " << m_decisionId << " : Event[" << elementIndex << "]" );
               passingSelection[cutIndex].push_back( elementIndex );
            } else {
               ATH_MSG_DEBUG("Not pass through selection " << m_decisionId << " : Event[" << elementIndex << "]" );
            } 
         } else {
            ATH_MSG_DEBUG("Not match DecisionID " << m_decisionId );
         }
         elementIndex++;
      }
 
      // If no object passes the selection, multipul selection should stop.
      if ( passingSelection[cutIndex].empty() ) {
         ATH_MSG_DEBUG( "No object passed selection " << cutIndex << " rejecting" );
         return StatusCode::SUCCESS;
      }
   }
 
   std::set<size_t> passingIndices;
   if ( m_decisionPerCluster==true ) {
      auto notFromSameRoI = [&]( const HLT::Index1DVec& comb ) {
         std::set<const xAOD::L2CombinedMuon*> setOfClusters;
         for ( auto index: comb ) {
            setOfClusters.insert( input[index].muComb );
         }
         return setOfClusters.size() == comb.size();
      };
 
      HLT::elementsInUniqueCombinations( passingSelection, passingIndices, std::move(notFromSameRoI) );
 
   } else {
      HLT::elementsInUniqueCombinations( passingSelection, passingIndices );
   }
 
   if ( passingIndices.empty() ) {
      ATH_MSG_DEBUG("No muon event passed through selection " << m_decisionId );
      return StatusCode::SUCCESS;
   }
 
   for ( auto idx: passingIndices ) {
      ATH_MSG_DEBUG("Muon event[" << idx << "] passes through Chain/ID " << m_decisionId 
                    << " with pT = " << input[idx].muComb->pt() << "GeV" );
      TrigCompositeUtils::addDecisionID( m_decisionId.numeric(), input[idx].decision );
   }
 
   return StatusCode::SUCCESS;
}

outputHandles()

virtual std::vector<Gaudi::DataHandle*> AthCommonDataStore< AthCommonMsg< AlgTool > >::outputHandles ( ) const

overridevirtualinherited

Return this algorithm's output handles.

We override this to include handle instances from key arrays if they have not yet been declared. See comments on updateVHKA.

renounce()

std::enable_if_t<std::is_void_v<std::result_of_t<decltype(&T::renounce)(T)> > && !std::is_base_of_v<SG::VarHandleKeyArray, T> && std::is_base_of_v<Gaudi::DataHandle, T>, void> AthCommonDataStore< AthCommonMsg< AlgTool > >::renounce ( T &  h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

  {
    h.renounce();
    PBASE::renounce (h);
  }

renounceArray()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::renounceArray ( SG::VarHandleKeyArray &  handlesArray )

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

sysInitialize()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysInitialize ( )

overridevirtualinherited

Perform system initialization for an algorithm.

We override this to declare all the elements of handle key arrays at the end of initialization. See comments on updateVHKA.

Reimplemented in DerivationFramework::CfAthAlgTool, AthCheckedComponent< AthAlgTool >, AthCheckedComponent<::AthAlgTool >, and asg::AsgMetadataTool.

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< AlgTool > >::sysStart ( )

overridevirtualinherited

Handle START transition.

We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.

updateVHKA()

void AthCommonDataStore< AthCommonMsg< AlgTool > >::updateVHKA ( Gaudi::Details::PropertyBase &  )

inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

                                               {
    // debug() << "updateVHKA for property " << p.name() << " " << p.toString() 
    //         << "  size: " << m_vhka.size() << endmsg;
    for (auto &a : m_vhka) {
      std::vector<SG::VarHandleKey*> keys = a->keys();
      for (auto k : keys) {
        k->setOwner(this);
      }
    }
  }

Member Data Documentation

m_acceptAll

Gaudi::Property< bool > TrigmuCombHypoTool::m_acceptAll

private

Initial value:

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

Definition at line 90 of file TrigmuCombHypoTool.h.

m_applyOR

Gaudi::Property< bool > TrigmuCombHypoTool::m_applyOR

private

Initial value:

{
      this, "ApplyOR", false, "apply overlap removal for mucomb" }

Definition at line 128 of file TrigmuCombHypoTool.h.

m_bins

std::vector<size_t> TrigmuCombHypoTool::m_bins = {0}

private

Definition at line 159 of file TrigmuCombHypoTool.h.

m_chi2MaxID

Gaudi::Property< double > TrigmuCombHypoTool::m_chi2MaxID

private

Initial value:

{
    this, "MaxChi2IDPik", 3.5, "" }

Definition at line 100 of file TrigmuCombHypoTool.h.

m_d0min

Gaudi::Property< float > TrigmuCombHypoTool::m_d0min

private

Initial value:

{
        this, "MinimumD0", 0., "lower d0 cut (mm)"}

Definition at line 155 of file TrigmuCombHypoTool.h.

m_decisionId

HLT::Identifier TrigmuCombHypoTool::m_decisionId

private

Definition at line 80 of file TrigmuCombHypoTool.h.

m_decisionPerCluster

Gaudi::Property<bool> TrigmuCombHypoTool::m_decisionPerCluster

private

Initial value:

{ 
         this, "DecisionPerCluster", true, "Is multiplicity requirement refering to muons ( false ) or RoIs/clusters with muons ( false ), relevant only in when multiplicity > 1" }

Definition at line 123 of file TrigmuCombHypoTool.h.

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_dRThres

Gaudi::Property< std::vector<float> > TrigmuCombHypoTool::m_dRThres

private

Initial value:

{
        this, "DRThres", {0.06, 0.05, 0.05}, "DR threshold"}

Definition at line 146 of file TrigmuCombHypoTool.h.

m_etaBins

Gaudi::Property< std::vector<float> > TrigmuCombHypoTool::m_etaBins

private

Initial value:

{
        this, "EtaBins", {0, 1.9, 2.1, 9.9}, "eta bins of DR and mass thresholds"}

Definition at line 143 of file TrigmuCombHypoTool.h.

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< AlgTool > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_massThres

Gaudi::Property< std::vector<float> > TrigmuCombHypoTool::m_massThres

private

Initial value:

{
        this, "MassThres", {0.20, 0.15, 0.10}, "mass threshold"}

Definition at line 152 of file TrigmuCombHypoTool.h.

m_maxPtToApplyPik

Gaudi::Property< double > TrigmuCombHypoTool::m_maxPtToApplyPik

private

Initial value:

{
    this, "MaxPtToApplyPik", 25., "" }

Definition at line 97 of file TrigmuCombHypoTool.h.

m_monTool

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

private

Definition at line 82 of file TrigmuCombHypoTool.h.

m_mufastDRThres

Gaudi::Property< std::vector<float> > TrigmuCombHypoTool::m_mufastDRThres

private

Initial value:

{
        this, "MufastDRThres", {0.06, 0.05, 0.05}, "mufast DR threshold"}

Definition at line 149 of file TrigmuCombHypoTool.h.

m_numMuonThreshold

Gaudi::Property< size_t > TrigmuCombHypoTool::m_numMuonThreshold

private

Initial value:

{
    this, "numMuonThreshold", 50, "Threshold of the number of muons" }

Definition at line 103 of file TrigmuCombHypoTool.h.

m_pikCuts

Gaudi::Property< bool > TrigmuCombHypoTool::m_pikCuts

private

Initial value:

{
    this, "ApplyPikCuts", false, "Enable pi/k decay rejection" }

Definition at line 94 of file TrigmuCombHypoTool.h.

m_ptBins

Gaudi::Property< std::vector<std::vector<double> > > TrigmuCombHypoTool::m_ptBins

private

Initial value:

{
    this, "PtBins", { {0, 2.5} }, "Bins range of each pT threshold" }

Definition at line 84 of file TrigmuCombHypoTool.h.

m_ptThresholds

Gaudi::Property< std::vector<std::vector<double> > > TrigmuCombHypoTool::m_ptThresholds

private

Initial value:

{
    this, "PtThresholds", { {5.83 * Gaudi::Units::GeV} }, "Track pT requirement ( separate threshold for each muon )" }

Definition at line 87 of file TrigmuCombHypoTool.h.

m_requireDR

Gaudi::Property<bool> TrigmuCombHypoTool::m_requireDR

private

Initial value:

{
        this, "RequireDR", true, "require or not DR cut for overlap removal"}

Definition at line 131 of file TrigmuCombHypoTool.h.

m_requireMass

Gaudi::Property<bool> TrigmuCombHypoTool::m_requireMass

private

Initial value:

{
        this, "RequireMass", true, "require or not mass cut for overlap removal"}

Definition at line 137 of file TrigmuCombHypoTool.h.

m_requireMufastDR

Gaudi::Property<bool> TrigmuCombHypoTool::m_requireMufastDR

private

Initial value:

{
        this, "RequireMufastDR", true, "require or not mufast DR cut for overlap removal"}

Definition at line 134 of file TrigmuCombHypoTool.h.

m_requireSameSign

Gaudi::Property<bool> TrigmuCombHypoTool::m_requireSameSign

private

Initial value:

{
        this, "RequireSameSign", true, "require or not charge cut for overlap removal"}

Definition at line 140 of file TrigmuCombHypoTool.h.

m_strategydependent

Gaudi::Property< bool > TrigmuCombHypoTool::m_strategydependent

private

Initial value:

{
    this, "ApplyStrategyDependentCuts", false, "" }

Definition at line 113 of file TrigmuCombHypoTool.h.

m_strategyDependentPtCuts

Gaudi::Property< std::vector<double> > TrigmuCombHypoTool::m_strategyDependentPtCuts

private

Initial value:

{
    this, "StrategyDependentPtCuts",{ 0.0, 0.0, 0.0, -3.0 }, "" }

Definition at line 116 of file TrigmuCombHypoTool.h.

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< AlgTool > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

---

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

• TrigmuCombHypoTool.h
• TrigmuCombHypoTool.cxx