TrigBmumuxComboHypo.h

Go to the documentation of this file.

/*
  Copyright (C) 2002-2024 CERN for the benefit of the ATLAS collaboration
*/
 
#ifndef TRIG_TrigBmumuxComboHypo_H
#define TRIG_TrigBmumuxComboHypo_H
 
#include <string>
#include <vector>
#include <utility>
 
#include "Gaudi/Property.h"
#include "xAODTracking/TrackParticleContainer.h"
#include "xAODTracking/VertexContainer.h"
#include "xAODTracking/VertexAuxContainer.h"
#include "xAODMuon/MuonContainer.h"
#include "xAODTrigger/TrigComposite.h"
#include "xAODTrigBphys/TrigBphys.h"
#include "xAODTrigBphys/TrigBphysContainer.h"
 
#include "StoreGate/ReadHandleKey.h"
#include "StoreGate/WriteHandleKey.h"
 
#include "TrigCompositeUtils/TrigCompositeUtils.h"
#include "DecisionHandling/ComboHypo.h"
 
#include "TrkVKalVrtFitter/TrkVKalVrtFitter.h"
#include "InDetConversionFinderTools/VertexPointEstimator.h"
#include "ITrackToVertex/ITrackToVertex.h"
 
#include "AthenaMonitoringKernel/Monitored.h"
#include "AthenaMonitoringKernel/GenericMonitoringTool.h"
 
#include "ITrigBphysState.h"
#include "TrigBmumuxComboHypoTool.h"
 
#include "Constants.h"
typedef struct PDG20 PDG;
 

class TrigBmumuxState: public ::ITrigBphysState {
 public:
  TrigBmumuxState() = delete;
  TrigBmumuxState(const EventContext& context,
                  const TrigCompositeUtils::DecisionContainer& previousDecisions,
                  TrigCompositeUtils::DecisionContainer& decisions,
                  xAOD::TrigBphysContainer* trigBphysCollection = nullptr,
                  const InDet::BeamSpotData* beamSpotData = nullptr)
      : ITrigBphysState(context, previousDecisions, decisions, trigBphysCollection, beamSpotData) {
    dimuons.setStore(&dimuonsStore);
  }
  virtual ~TrigBmumuxState() = default;
 
  // EFCB muon candidates from mergeMuonsFromDecisions()
  struct Muon {
    ElementLink<xAOD::MuonContainer> link;
    std::vector<ElementLink<TrigCompositeUtils::DecisionContainer>> decisionLinks;
    TrigCompositeUtils::DecisionIDContainer decisionIDs;
  };
  std::vector<Muon> muons;
 
  // tracks from mergeTracksFromViews()
  std::vector<ElementLink<xAOD::TrackParticleContainer>> tracks;
  bool isCompositeRoI = false;
 
  // dimuon candidates from findDimuonCandidates()
  xAOD::VertexContainer dimuons;
  xAOD::VertexAuxContainer dimuonsStore;
  std::vector<std::array<size_t, 2>> trigBphysMuonIndices;  // {mu1, mu2} indices in state.muons for dimuon in state.dimuons
 
  // tracks to be fitted with given dimuon candidate from findBmumuxCandidates_selectTracks()
  std::vector<ElementLink<xAOD::TrackParticleContainer>> selectedTracks;
  std::map<const xAOD::TrackParticle*, double> selectedTrackZ0;  // track z0 impact parameters wrt dimuon vertex
 
  // track combinations from findBmumuxCandidates_fit(makeCombinations = true)
  std::map<size_t, size_t> trackCombinations;  // key = track index for dimuon + track case OR state.getCombinationKey(trk1, trk2) for track + track case : number of found combinations
 
  // sorted vector of track combinations which failed vertex fit from findBmumuxCandidates_fastFit()
  std::vector<size_t> badTrackCombinations;
 
  size_t getCombinationKey(size_t i1, size_t i2) const { size_t n = selectedTracks.size(); return (i1 < i2 ? i1 + n * i2 : i2 + n * i1); }
  void addTrackCombination(size_t i1) { trackCombinations[i1]++; }
  void addTrackCombination(size_t i1, size_t i2) { trackCombinations[getCombinationKey(i1, i2)]++; }
  void addTrackCombination(size_t i1, size_t i2, size_t i3) { addTrackCombination(i1, i2); addTrackCombination(i1, i3); addTrackCombination(i2, i3); }
  bool isBadCombination(size_t i1) const { return std::binary_search(badTrackCombinations.begin(), badTrackCombinations.end(), i1); }
  bool isBadCombination(size_t i1, size_t i2, size_t i3) const { return (isBadCombination(getCombinationKey(i1, i2)) || isBadCombination(getCombinationKey(i1, i3)) || isBadCombination(getCombinationKey(i2, i3))); }
 
  StatusCode addTriggerObject(xAOD::TrigBphys* triggerObject) {
    if (!triggerObject) {
      return StatusCode::FAILURE;
    }
    trigBphysCollection().push_back(triggerObject);
    return StatusCode::SUCCESS;
  }
 
};
 

class TrigBmumuxComboHypo: public ::ComboHypo {
 public:
  TrigBmumuxComboHypo(const std::string& name, ISvcLocator* pSvcLocator);
  TrigBmumuxComboHypo() = delete;
 
  virtual StatusCode initialize() override;
  virtual StatusCode execute(const EventContext& context) const override;
 
  enum Decay : size_t {
    kPsi_2mu,         // psi -> mu+ mu-
    kB_2mu1trk,       // B -> mu+ mu- trk1
    kB_2mu2trk,       // B -> mu+ mu- trk1 trk2
    kDs,              // D_s+ -> K+ K- pi+
    kDplus,           // D+ -> K- pi+ pi+
    kD0,              // D0 -> K- pi+
    kB_PsiPi,         // psi + pion from D*+
    kFastFit_2trk,    // trk1 + trk2
    kFastFit_2mu1trk  // dimuon + trk1
  };
 
 private:

  StatusCode mergeMuonsFromDecisions(TrigBmumuxState&) const;

  StatusCode mergeTracksFromViews(TrigBmumuxState&) const;

  StatusCode findDimuonCandidates(TrigBmumuxState&) const;

  StatusCode findBmumuxCandidates(TrigBmumuxState&) const;

  StatusCode findBmumuxCandidates_selectTracks(TrigBmumuxState&, size_t dimuonIndex) const;

  StatusCode findBmumuxCandidates_fastFit(TrigBmumuxState&, size_t dimuonIndex) const;

  StatusCode findBmumuxCandidates_fit(TrigBmumuxState&, size_t dimuonIndex, bool makeCombinations = false) const;

  StatusCode createDecisionObjects(TrigBmumuxState&) const;

  std::unique_ptr<xAOD::Vertex> fit(
      const EventContext& context,
      const std::vector<ElementLink<xAOD::TrackParticleContainer>>& trackParticleLinks,
      Decay decay = kPsi_2mu,
      const xAOD::Vertex* dimuon = nullptr) const;

  xAOD::TrigBphys* makeTriggerObject(
      TrigBmumuxState& state,
      const xAOD::Vertex& vertex,
      xAOD::TrigBphys::pType type = xAOD::TrigBphys::MULTIMU,
      const std::vector<double>& trkMass = {PDG::mMuon, PDG::mMuon},
      const ElementLink<xAOD::TrigBphysContainer>& dimuonLink = ElementLink<xAOD::TrigBphysContainer>()) const;

  bool isIdenticalTracks(const xAOD::TrackParticle* lhs, const xAOD::TrackParticle* rhs) const;
  bool isIdenticalTracks(const xAOD::Muon* lhs, const xAOD::Muon* rhs) const;
  bool isInSameRoI(const xAOD::Muon*, const xAOD::TrackParticle*) const;
  bool passDimuonTrigger(const std::vector<const TrigCompositeUtils::DecisionIDContainer*>& previousDecisionIDs) const;

  bool isInMassRange(double mass, const std::pair<double, double>& range) const { return (mass > range.first && mass < range.second); }
  double Lxy(const Amg::Vector3D& productionVertex, const xAOD::Vertex& decayVertex) const;

  xAOD::TrackParticle::GenVecFourMom_t momentum(const xAOD::Vertex& vertex, const std::vector<double>& trkMass) const;
 
  SG::ReadHandleKey<xAOD::TrackParticleContainer> m_trackParticleContainerKey {this,
    "TrackCollectionKey", "InDetTrackParticles", "input TrackParticle container name"};
  SG::ReadHandleKey<xAOD::MuonContainer> m_muonContainerKey {this,
    "MuonCollectionKey", "Muons", "input EF Muon container name"};
  SG::WriteHandleKey<xAOD::TrigBphysContainer> m_trigBphysContainerKey {this,
    "TrigBphysCollectionKey", "TrigBphysContainer", "output TrigBphysContainer name"};
  SG::ReadCondHandleKey<InDet::BeamSpotData>
    m_beamSpotKey {this, "BeamSpotKey", "BeamSpotData", "SG key for beam spot"};
 
  // general properties
  Gaudi::Property<double> m_deltaR {this,
    "DeltaR", 0.01, "minimum deltaR between same-sign tracks (overlap removal)"};
  Gaudi::Property<double> m_trkZ0 {this,
    "TrkZ0", 50., "maximum z0 impact parameter of the track wrt the fitted dimuon vertex; no preselection if negative"};
  Gaudi::Property<double> m_roiEtaWidth {this,
    "RoiEtaWidth", 0.75, "extent of the RoI in eta from initial muon, to be check for SuperRoI"};
  Gaudi::Property<double> m_roiPhiWidth {this,
    "RoiPhiWidth", 0.75, "extent of the RoI in phi from initial muon, to be check for SuperRoI"};
  Gaudi::Property<size_t> m_fitAttemptsWarningThreshold {this,
    "FitAttemptsWarningThreshold", 200, "Events processing this many calls of the vertex fitter will generate a WARNING message (time-out protect)"};
  Gaudi::Property<size_t> m_fitAttemptsBreakThreshold {this,
    "FitAttemptsBreakThreshold", 1000, "Events processing this many calls of the vertex fitter will generate a second WARNING message and the loop over combinations will be terminated at this point (time-out protect)"};
  Gaudi::Property<double> m_fastFit_2trk_chi2 {this,
    "FastFit_2trk_chi2", 50., "maximum chi2 for fast fit of trk1 + trk2"};
   Gaudi::Property<double> m_fastFit_2mu1trk_chi2 {this,
    "FastFit_2mu1trk_chi2", 60., "maximum chi2 for fast fit of dimuon + trk1"};
 
  // dimuon properties
  Gaudi::Property<bool> m_dimuon_rejectSameChargeTracks {this,
    "Dimuon_rejectSameChargeTracks", true, "if true, the only (mu+, mu-) pairs will be kept (no wrong-charge combinations)"};
  Gaudi::Property<std::pair<double, double>> m_dimuon_massRange {this,
    "Dimuon_massRange", {100., 5500.}, "dimuon mass range"};
  Gaudi::Property<double> m_dimuon_chi2 {this,
    "Dimuon_chi2", 20., "maximum chi2 of the dimuon vertex"};
 
  // B+ -> mu+ mu- K+
  Gaudi::Property<bool> m_BplusToMuMuKaon {this,
    "BplusToMuMuKaon", true, "switch on/off B+ -> mu+ mu- K+ decay"};
  Gaudi::Property<double> m_BplusToMuMuKaon_minKaonPt {this,
    "BplusToMuMuKaon_minKaonPt", 100., "minimum pT of kaon track"};
  Gaudi::Property<std::pair<double, double>> m_BplusToMuMuKaon_massRange {this,
    "BplusToMuMuKaon_massRange", {4500., 5900.}, "B+ mass range"};
  Gaudi::Property<float> m_BplusToMuMuKaon_chi2 {this,
    "BplusToMuMuKaon_chi2", 50., "maximum chi2 of the fitted B+ vertex"};
  Gaudi::Property<bool> m_BplusToMuMuKaon_useFastFit {this,
    "BplusToMuMuKaon_useFastFit", false, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B_c+ -> J/psi(-> mu+ mu-) pi+
  Gaudi::Property<bool> m_BcToMuMuPion {this,
    "BcToMuMuPion", true, "switch on/off B_c+ -> J/psi(-> mu+ mu-) pi+ decay"};
  Gaudi::Property<double> m_BcToMuMuPion_minPionPt {this,
    "BcToMuMuPion_minPionPt", 2000., "minimum pT of pion track"};
  Gaudi::Property<std::pair<double, double>> m_BcToMuMuPion_dimuonMassRange {this,
    "BcToMuMuPion_dimuonMassRange", {2500., 4300.}, "dimuon mass range for B_c+ decay"};
  Gaudi::Property<std::pair<double, double>> m_BcToMuMuPion_massRange {this,
    "BcToMuMuPion_massRange", {5500., 7300.}, "B_c+ mass range"};
  Gaudi::Property<float> m_BcToMuMuPion_chi2 {this,
    "BcToMuMuPion_chi2", 50., "maximum chi2 of the fitted B_c+ vertex"};
  Gaudi::Property<bool> m_BcToMuMuPion_useFastFit {this,
    "BcToMuMuPion_useFastFit", false, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B_s0 -> mu+ mu- phi(-> K+ K-)
  Gaudi::Property<bool> m_BsToMuMuPhi1020 {this,
    "BsToMuMuPhi1020", true, "switch on/off B_s0 -> mu+ mu- phi(-> K+ K-) decay"};
  Gaudi::Property<bool> m_BsToMuMuPhi1020_rejectSameChargeTracks {this,
    "BsToMuMuPhi1020_rejectSameChargeTracks", true, "if true, the only (K+, K-) pairs will be kept (no wrong-charge combinations)"};
  Gaudi::Property<double> m_BsToMuMuPhi1020_minKaonPt {this,
    "BsToMuMuPhi1020_minKaonPt", 100., "minimum pT of kaon tracks"};
  Gaudi::Property<std::pair<double, double>> m_BsToMuMuPhi1020_massRange {this,
    "BsToMuMuPhi1020_massRange", {4800., 5800.}, "B_s0 mass range"};
  Gaudi::Property<std::pair<double, double>> m_BsToMuMuPhi1020_phiMassRange {this,
    "BsToMuMuPhi1020_phiMassRange", {940., 1100.}, "phi1020 mass range"};
  Gaudi::Property<float> m_BsToMuMuPhi1020_chi2 {this,
    "BsToMuMuPhi1020_chi2", 60., "maximum chi2 of the fitted B+ vertex"};
  Gaudi::Property<bool> m_BsToMuMuPhi1020_useFastFit {this,
    "BsToMuMuPhi1020_useFastFit", false, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B0 -> mu+ mu- K*0(-> K+ pi-)
  Gaudi::Property<bool> m_BdToMuMuKstar0 {this,
    "BdToMuMuKstar0", true, "switch on/off B0 -> mu+ mu- K*0(-> K+ pi-) decay"};
  Gaudi::Property<bool> m_BdToMuMuKstar0_rejectSameChargeTracks {this,
    "BdToMuMuKstar0_rejectSameChargeTracks", true, "if true, the only (K+, pi-) and (K-, pi+) pairs will be kept (no wrong-charge combinations)"};
  Gaudi::Property<double> m_BdToMuMuKstar0_minKaonPt {this,
    "BdToMuMuKstar0_minKaonPt", 100., "minimum pT of kaon track"};
  Gaudi::Property<double> m_BdToMuMuKstar0_minPionPt {this,
    "BdToMuMuKstar0_minPionPt", 100., "minimum pT of pion track"};
  Gaudi::Property<std::pair<double, double>> m_BdToMuMuKstar0_massRange {this,
    "BdToMuMuKstar0_massRange", {4600., 5900.}, "B0 mass range"};
  Gaudi::Property<std::pair<double, double>> m_BdToMuMuKstar0_KstarMassRange {this,
    "BdToMuMuKstar0_KstarMassRange", {700., 1100.}, "K*0 mass range"};
  Gaudi::Property<float> m_BdToMuMuKstar0_chi2 {this,
    "BdToMuMuKstar0_chi2", 60., "maximum chi2 of the fitted B0 vertex"};
  Gaudi::Property<bool> m_BdToMuMuKstar0_useFastFit {this,
    "BdToMuMuKstar0_useFastFit", true, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // Lambda_b0 -> J/psi(-> mu+ mu-) p K-
  Gaudi::Property<bool> m_LambdaBToMuMuProtonKaon {this,
    "LambdaBToMuMuProtonKaon", true, "switch on/off Lambda_b0 -> J/psi(-> mu+ mu-) p K- decay"};
  Gaudi::Property<bool> m_LambdaBToMuMuProtonKaon_rejectSameChargeTracks {this,
    "LambdaBToMuMuProtonKaon_rejectSameChargeTracks", false, "if true, the only (p, K-) and (anti-p, K+) pairs will be kept (no wrong-charge combinations)"};
  Gaudi::Property<double> m_LambdaBToMuMuProtonKaon_minProtonPt {this,
    "LambdaBToMuMuProtonKaon_minProtonPt", 1000., "minimum pT of proton track"};
  Gaudi::Property<double> m_LambdaBToMuMuProtonKaon_minKaonPt {this,
    "LambdaBToMuMuProtonKaon_minKaonPt", 1000., "minimum pT of kaon track"};
  Gaudi::Property<double> m_LambdaBToMuMuProtonKaon_minKstarMass {this,
    "LambdaBToMuMuProtonKaon_minKstarMass", 1300., "min value for both mass(trk1=kaon, trk2=pion) and mass(trk1=pion, trk2=kaon)"};
  Gaudi::Property<std::pair<double, double>> m_LambdaBToMuMuProtonKaon_dimuonMassRange {this,
    "LambdaBToMuMuProtonKaon_dimuonMassRange", {2500., 4300.}, "dimuon mass range for Lambda_b0 decay"};
  Gaudi::Property<std::pair<double, double>> m_LambdaBToMuMuProtonKaon_massRange {this,
    "LambdaBToMuMuProtonKaon_massRange", {4800., 6400.}, "Lambda_b0 mass range"};
  Gaudi::Property<float> m_LambdaBToMuMuProtonKaon_chi2 {this,
    "LambdaBToMuMuProtonKaon_chi2", 60., "maximum chi2 of the fitted Lambda_b0 vertex"};
  Gaudi::Property<bool> m_LambdaBToMuMuProtonKaon_useFastFit {this,
    "LambdaBToMuMuProtonKaon_useFastFit", true, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B_c+ -> J/psi(-> mu+ mu-) D_s+(->phi(-> K+ K-) pi+)
  Gaudi::Property<bool> m_BcToDsMuMu {this,
    "BcToDsMuMu", true, "switch on/off B_c+ -> J/psi(-> mu+ mu-) D_s+(->phi(-> K+ K-) pi+) decay"};
  Gaudi::Property<double> m_BcToDsMuMu_minKaonPt {this,
    "BcToDsMuMu_minKaonPt", 1000., "minimum pT of kaon track from phi(1020)"};
  Gaudi::Property<double> m_BcToDsMuMu_minPionPt {this,
    "BcToDsMuMu_minPionPt", 1000., "minimum pT of pion track from D_s+"};
  Gaudi::Property<std::pair<double, double>> m_BcToDsMuMu_massRange {this,
    "BcToDsMuMu_massRange", {5500., 7300.}, "B_c+ mass range"};
  Gaudi::Property<std::pair<double, double>> m_BcToDsMuMu_dimuonMassRange {this,
    "BcToDsMuMu_dimuonMassRange", {2500., 4300.}, "dimuon mass range for B_c+ -> J/psi D_s+ decay"};
  Gaudi::Property<std::pair<double, double>> m_BcToDsMuMu_phiMassRange {this,
    "BcToDsMuMu_phiMassRange", {940., 1100.}, "phi(1020) mass range"};
  Gaudi::Property<std::pair<double, double>> m_BcToDsMuMu_DsMassRange {this,
    "BcToDsMuMu_DsMassRange", {1750., 2100.}, "D_s+ mass range"};
  Gaudi::Property<float> m_BcToDsMuMu_chi2 {this,
    "BcToDsMuMu_chi2", 60., "maximum chi2 of the fitted B_c+ vertex"};
  Gaudi::Property<bool> m_BcToDsMuMu_useFastFit {this,
    "BcToDsMuMu_useFastFit", true, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B_c+ -> J/psi(-> mu+ mu-) D+(-> K- pi+ pi+)
  Gaudi::Property<bool> m_BcToDplusMuMu {this,
    "BcToDplusMuMu", true, "switch on/off B_c+ -> J/psi(-> mu+ mu-) D+(-> K- pi+ pi+) decay"};
  Gaudi::Property<double> m_BcToDplusMuMu_minKaonPt {this,
    "BcToDplusMuMu_minKaonPt", 1000., "minimum pT of kaon track from D+"};
  Gaudi::Property<double> m_BcToDplusMuMu_minPionPt {this,
    "BcToDplusMuMu_minPionPt", 1000., "minimum pT of pion track from D+"};
  Gaudi::Property<std::pair<double, double>> m_BcToDplusMuMu_massRange {this,
    "BcToDplusMuMu_massRange", {5500., 7300.}, "B_c+ mass range"};
  Gaudi::Property<std::pair<double, double>> m_BcToDplusMuMu_dimuonMassRange {this,
    "BcToDplusMuMu_dimuonMassRange", {2500., 4300.}, "dimuon mass range for B_c+ -> J/psi D+ decay"};
  Gaudi::Property<std::pair<double, double>> m_BcToDplusMuMu_DplusMassRange {this,
    "BcToDplusMuMu_DplusMassRange", {1750., 2000.}, "D+ mass range"};
  Gaudi::Property<float> m_BcToDplusMuMu_chi2 {this,
    "BcToDplusMuMu_chi2", 60., "maximum chi2 of the fitted B_c+ vertex"};
  Gaudi::Property<bool> m_BcToDplusMuMu_useFastFit {this,
    "BcToDplusMuMu_useFastFit", true, "true: perform vertex fit depending on only if fast fit; false: always perform vertex fit"};
 
  // B_c+ -> J/psi(-> mu+ mu-) D*+(-> D0(-> K- pi+) pi+)
  Gaudi::Property<bool> m_BcToDstarMuMu {this,
    "BcToDstarMuMu", true, "switch on/off partial reconstruction of B_c+ -> J/psi(-> mu+ mu-) D0(-> K- pi+) X decay"};
  Gaudi::Property<bool> m_BcToDstarMuMu_makeDstar {this,
    "BcToDstarMuMu_makeDstar", true, "switch on/off full reconstruction of B_c+ -> J/psi(-> mu+ mu-) D*+(-> D0(-> K- pi+) pi+) decay"};
  Gaudi::Property<double> m_BcToDstarMuMu_minD0KaonPt {this,
    "BcToDstarMuMu_minD0KaonPt", 1000., "minimum pT of kaon track from D0"};
  Gaudi::Property<double> m_BcToDstarMuMu_minD0PionPt {this,
    "BcToDstarMuMu_minD0PionPt", 1000., "minimum pT of pion track from D0"};
  Gaudi::Property<double> m_BcToDstarMuMu_minDstarPionPt {this,
    "BcToDstarMuMu_minDstarPionPt", 500., "minimum pT of pion track from D*+"};
  Gaudi::Property<double> m_BcToDstarMuMu_maxDstarPionZ0 {this,
    "BcToDstarMuMu_maxDstarPionZ0", 5., "maximum z0 impact parameter of the pion track from D*+ wrt the fitted dimuon vertex; no preselection if negative"};
  Gaudi::Property<std::pair<double, double>> m_BcToDstarMuMu_massRange {this,
    "BcToDstarMuMu_massRange", {5500., 7300.}, "B_c+ mass range"};
  Gaudi::Property<std::pair<double, double>> m_BcToDstarMuMu_dimuonMassRange {this,
    "BcToDstarMuMu_dimuonMassRange", {2500., 4300.}, "dimuon mass range for B_c+ -> J/psi D+ decay"};
  Gaudi::Property<std::pair<double, double>> m_BcToDstarMuMu_D0MassRange {this,
    "BcToDstarMuMu_D0MassRange", {1750., 2000.}, "D0 mass range"};
  Gaudi::Property<std::pair<double, double>> m_BcToDstarMuMu_DstarMassRange {this,
    "BcToDstarMuMu_DstarMassRange", {-1., 2110.}, "D*+ mass range"};
  Gaudi::Property<float> m_BcToDstarMuMu_chi2 {this,
    "BcToDstarMuMu_chi2", 60., "maximum chi2 of the fitted B_c+ vertex"};
 
  // external tools
  ToolHandle<InDet::VertexPointEstimator> m_vertexPointEstimator {this,
    "VertexPointEstimator", "", "tool to find starting point for the vertex fitter"};
  ToolHandle<Trk::TrkVKalVrtFitter> m_vertexFitter {this,
    "VertexFitter", "", "VKalVrtFitter tool to fit tracks into the common vertex"};
  ToolHandle<Reco::ITrackToVertex> m_trackToVertexTool {this,
    "TrackToVertexTool", "", "tool to extrapolate track to vertex or beamspot"};
  ToolHandle<GenericMonitoringTool> m_monTool {this,
    "MonTool", "", "monitoring tool"};
 
  TrigCompositeUtils::DecisionIDContainer m_allowedIDs;
 
  const static std::vector<std::vector<double>> s_trkMass;
 
};
 
#endif  // TRIG_TrigBmumuxComboHypo_H