TrigHitDVHypoTool.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2022 CERN for the benefit of the ATLAS collaboration
*/
#ifndef TRIGLONGLIVEDPARTICLESHYPO_TRIGHITDVHYPOTOOL_H
#define TRIGLONGLIVEDPARTICLESHYPO_TRIGHITDVHYPOTOOL_H
 
#include <string>
 
#include "Gaudi/Property.h"
#include "AthenaBaseComps/AthAlgTool.h"
#include "AthenaMonitoringKernel/GenericMonitoringTool.h"
#include "TrigCompositeUtils/TrigCompositeUtils.h"
#include "xAODTracking/TrackParticlexAODHelpers.h"
#include "xAODTracking/TrackParticleContainer.h"
#include "xAODJet/JetContainer.h"
 
 
 
class TrigHitDVHypoTool : virtual public ::AthAlgTool
{
public:
   TrigHitDVHypoTool( const std::string& type,
              const std::string& name,
              const IInterface* parent );
 
   virtual ~TrigHitDVHypoTool();
   virtual StatusCode initialize() override;
 
   struct HitDVHypoInfo {
      TrigCompositeUtils::Decision*       decision{nullptr};
      bool                                isSPOverflow{false};
      float                               averageMu{0};
      const xAOD::TrigComposite*          hitDV{nullptr};
      const TrigCompositeUtils::DecisionIDContainer previousDecisionsIDs;
   };
 
   // seed type enum
   enum SeedType {
      HLTJet = 2,
      SP = 4
   };
 
   StatusCode decide( std::vector<HitDVHypoInfo>& )  const;
 
private:
 
   HLT::Identifier m_decisionId;
   Gaudi::Property< std::vector<float> >  m_cutJetPtGeV { this, "cutJetPtGeV", { 20.0 }, "Jet pT requirement in GeV" };
   Gaudi::Property< std::vector<float> >  m_cutJetEta   { this, "cutJetEta",   { 2.5 },  "Jet Eta requirement" };
   Gaudi::Property< std::vector<float> >  m_effBDT      { this, "effBDT",      { 0.7 },  "Efficiency for BDT cut." };
   Gaudi::Property< std::vector<bool>  >  m_doSPseed    { this, "doSPseed",    { true },  "Switch to do SP seeding" };
 
   //
   StatusCode inclusiveSelection(std::vector<HitDVHypoInfo>& ) const;
   StatusCode multiplicitySelection(std::vector<HitDVHypoInfo>& ) const;
   bool decideOnSingleObject( HitDVHypoInfo&, size_t, bool ) const;
 
   //
   float      deltaR(float, float, float, float) const;
   int        getSPLayer(int, float) const;
   StatusCode findSPSeeds( const xAOD::TrigCompositeContainer*, std::vector<float>&, std::vector<float>& ) const;
   StatusCode findJetSeeds(const xAOD::JetContainer*, const float, const float, std::vector<float>&, std::vector<float>&) const;
   StatusCode calculateBDT(const xAOD::TrigCompositeContainer*, const xAOD::TrigCompositeContainer*,
               const std::vector<float>&, const std::vector<float>&, const float, const int, xAOD::TrigCompositeContainer*, int&) const;
   float      getBDTthreshold_0eta1(float,float) const;
   float      getBDTthreshold_1eta2(float,float) const;
   float      getBDTthreshold(float) const;
};
 
#endif //> !TRIGLONGLIVEDPARTICLESHYPO_TRIGHITDVHYPOTOOL_H