#include <MuPlusDsCascade.h>

Inheritance diagram for DerivationFramework::MuPlusDsCascade:

Collaboration diagram for DerivationFramework::MuPlusDsCascade:

Public Member Functions
	MuPlusDsCascade (const std::string &t, const std::string &n, const IInterface *p)

	~MuPlusDsCascade ()

virtual StatusCode	initialize () override

StatusCode	performSearch (std::vector< Trk::VxCascadeInfo * > *cascadeinfoContainer) const

virtual StatusCode	addBranches () const override
	Pass the thinning service More...

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

Static Public Member Functions
static const InterfaceID &	interfaceID ()

static xAOD::Vertex *	FindVertexTrack (const xAOD::MuonContainer c, xAOD::Vertex v)

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
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
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
std::string	m_vertexContainerKey

std::string	m_vertexDxContainerKey

std::vector< std::string >	m_cascadeOutputsKeys

std::string	m_VxPrimaryCandidateName
	Name of primary vertex container. More...

double	m_DxMassLower

double	m_DxMassUpper

double	m_MassLower

double	m_MassUpper

double	m_vtx0MassHypo

double	m_vtx1MassHypo

double	m_vtx0Daug1MassHypo

double	m_vtx1Daug1MassHypo

double	m_vtx1Daug2MassHypo

double	m_vtx1Daug3MassHypo

int	m_Dx_pid

bool	m_constrDx

double	m_chi2cut

SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfo_key {this, "EventInfo", "EventInfo", "Input event information"}

ToolHandle< Trk::TrkVKalVrtFitter >	m_iVertexFitter

ToolHandle< Analysis::PrimaryVertexRefitter >	m_pvRefitter

ToolHandle< Trk::V0Tools >	m_V0Tools

ToolHandle< DerivationFramework::CascadeTools >	m_CascadeTools

std::unique_ptr< InDet::InDetTrackSelectionTool >	m_trackSelectionTools

ServiceHandle< IPartPropSvc >	m_partPropSvc {this, "PartPropSvc", "PartPropSvc"}

bool	m_refitPV

std::string	m_refPVContainerName

std::string	m_hypoName
	name of the mass hypothesis. More...

int	m_PV_max

int	m_DoVertexType

size_t	m_PV_minNTracks

std::string	m_muonCollectionKey

ToolHandle< Trk::ITrackSelectorTool >	m_trkSelector

double	m_thresholdPt

bool	m_mcpCuts

bool	m_combOnly

bool	m_useCombMeasurement

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 41 of file MuPlusDsCascade.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

◆ MuPlusDsCascade()

DerivationFramework::MuPlusDsCascade::MuPlusDsCascade	(	const std::string &	t,
		const std::string &	n,
		const IInterface *	p
	)

Definition at line 518 of file MuPlusDsCascade.cxx.

                                                                                                  : AthAlgTool(t,n,p),
     m_vertexDxContainerKey(""),
     m_cascadeOutputsKeys{ "MuPlusDsCascadeVtx1", "MuPlusDsCascadeVtx2" },
     m_VxPrimaryCandidateName("PrimaryVertices"),
     m_DxMassLower(0.0),
     m_DxMassUpper(10000.0),
     m_MassLower(0.0),
     m_MassUpper(20000.0),
     m_vtx0MassHypo(-1),
     m_vtx1MassHypo(-1),
     m_vtx0Daug1MassHypo(-1),
     m_vtx1Daug1MassHypo(-1),
     m_vtx1Daug2MassHypo(-1),
     m_vtx1Daug3MassHypo(-1),
     m_Dx_pid(431),
     m_constrDx(true),
     m_chi2cut(-1.0),
     m_iVertexFitter("Trk::TrkVKalVrtFitter"),
     m_pvRefitter("Analysis::PrimaryVertexRefitter",this),
     m_V0Tools("Trk::V0Tools"),
     m_CascadeTools("DerivationFramework::CascadeTools"),
     m_muonCollectionKey("StacoMuonCollection"),
     m_trkSelector("InDet::TrackSelectorTool"),
     m_thresholdPt(0.0),
     m_mcpCuts(true),
     m_combOnly(false),
     m_useCombMeasurement(false)
     {
        declareProperty("DxVertices",                m_vertexDxContainerKey);
        declareProperty("VxPrimaryCandidateName",    m_VxPrimaryCandidateName);
        declareProperty("RefPVContainerName",        m_refPVContainerName     = "RefittedPrimaryVertices");
        declareProperty("DxMassLowerCut",            m_DxMassLower);
        declareProperty("DxMassUpperCut",            m_DxMassUpper);
        declareProperty("MassLowerCut",              m_MassLower);
        declareProperty("MassUpperCut",              m_MassUpper);
        declareProperty("HypothesisName",            m_hypoName               = "B");
        declareProperty("Vtx0MassHypo",              m_vtx0MassHypo);
        declareProperty("Vtx1MassHypo",              m_vtx1MassHypo);
        declareProperty("Vtx0Daug1MassHypo",         m_vtx0Daug1MassHypo);
        declareProperty("Vtx1Daug1MassHypo",         m_vtx1Daug1MassHypo);
        declareProperty("Vtx1Daug2MassHypo",         m_vtx1Daug2MassHypo);
        declareProperty("Vtx1Daug3MassHypo",         m_vtx1Daug3MassHypo);
        declareProperty("DxHypothesis",              m_Dx_pid);
        declareProperty("ApplyDxMassConstraint",     m_constrDx);
        declareProperty("Chi2Cut",                   m_chi2cut);
        declareProperty("RefitPV",                   m_refitPV                = true);
        declareProperty("MaxnPV",                    m_PV_max                 = 999);
        declareProperty("MinNTracksInPV",            m_PV_minNTracks          = 0);
        declareProperty("DoVertexType",              m_DoVertexType           = 7);
        declareProperty("TrkVertexFitterTool",       m_iVertexFitter);
        declareProperty("PVRefitter",                m_pvRefitter);
        declareProperty("V0Tools",                   m_V0Tools);
        declareProperty("CascadeTools",              m_CascadeTools);
        declareProperty("CascadeVertexCollections",  m_cascadeOutputsKeys);
        declareProperty("muonCollectionKey",         m_muonCollectionKey);
        declareProperty("TrackSelectorTool",         m_trkSelector);
        declareProperty("muonThresholdPt",           m_thresholdPt);
        declareProperty("useMCPCuts",                m_mcpCuts);
        declareProperty("combOnly",                  m_combOnly);
        declareProperty("useCombinedMeasurement",    m_useCombMeasurement);
     }

◆ ~MuPlusDsCascade()

DerivationFramework::MuPlusDsCascade::~MuPlusDsCascade ( )

Definition at line 580 of file MuPlusDsCascade.cxx.

580 { }

Member Function Documentation

◆ addBranches()

StatusCode DerivationFramework::MuPlusDsCascade::addBranches ( ) const

overridevirtual

Pass the thinning service

Implements DerivationFramework::IAugmentationTool.

Definition at line 95 of file MuPlusDsCascade.cxx.

     {
       std::vector<Trk::VxCascadeInfo*> cascadeinfoContainer;
       constexpr int topoN = 2;
       std::array<xAOD::VertexContainer*, topoN> Vtxwritehandles;
       std::array<xAOD::VertexAuxContainer*, topoN> Vtxwritehandlesaux;
       if(m_cascadeOutputsKeys.size() !=topoN)  { ATH_MSG_FATAL("Incorrect number of VtxContainers"); return StatusCode::FAILURE; }
  
       for(int i =0; i<topoN;i++){
          Vtxwritehandles[i] = new xAOD::VertexContainer();
          Vtxwritehandlesaux[i] = new xAOD::VertexAuxContainer();
          Vtxwritehandles[i]->setStore(Vtxwritehandlesaux[i]);
          ATH_CHECK(evtStore()->record(Vtxwritehandles[i]   , m_cascadeOutputsKeys[i]       ));
          ATH_CHECK(evtStore()->record(Vtxwritehandlesaux[i], m_cascadeOutputsKeys[i] + "Aux."));
       }
  
       //----------------------------------------------------
       // retrieve primary vertices
       //----------------------------------------------------
       const xAOD::Vertex * primaryVertex(nullptr);
       const xAOD::VertexContainer *pvContainer(nullptr);
       ATH_CHECK(evtStore()->retrieve(pvContainer, m_VxPrimaryCandidateName));
       ATH_MSG_DEBUG("Found " << m_VxPrimaryCandidateName << " in StoreGate!");
  
       if (pvContainer->size()==0){
         ATH_MSG_WARNING("You have no primary vertices: " << pvContainer->size());
         return StatusCode::RECOVERABLE;
       } else {
         primaryVertex = (*pvContainer)[0];
       }
  
       //----------------------------------------------------
       // Try to retrieve refitted primary vertices
       //----------------------------------------------------
       xAOD::VertexContainer*    refPvContainer    = nullptr;
       xAOD::VertexAuxContainer* refPvAuxContainer = nullptr;
       if (m_refitPV) {
         if (evtStore()->contains<xAOD::VertexContainer>(m_refPVContainerName)) {
           // refitted PV container exists. Get it from the store gate
           ATH_CHECK(evtStore()->retrieve(refPvContainer   , m_refPVContainerName       ));
           ATH_CHECK(evtStore()->retrieve(refPvAuxContainer, m_refPVContainerName + "Aux."));
         } else {
           // refitted PV container does not exist. Create a new one.
           refPvContainer = new xAOD::VertexContainer;
           refPvAuxContainer = new xAOD::VertexAuxContainer;
           refPvContainer->setStore(refPvAuxContainer);
           ATH_CHECK(evtStore()->record(refPvContainer   , m_refPVContainerName));
           ATH_CHECK(evtStore()->record(refPvAuxContainer, m_refPVContainerName+"Aux."));
         }
       }
  
       ATH_CHECK(performSearch(&cascadeinfoContainer));
       SG::ReadHandle<xAOD::EventInfo> evt(m_eventInfo_key);
       if(!evt.isValid()) {
           ATH_MSG_ERROR("Cannot Retrieve " << m_eventInfo_key.key() );
           return StatusCode::FAILURE;
       }
  
       BPhysPVCascadeTools helper(&(*m_CascadeTools), evt.cptr());
       helper.SetMinNTracksInPV(m_PV_minNTracks);
  
       // Decorators for the main vertex: chi2, ndf, pt and pt error, plus the V0 vertex variables
       SG::AuxElement::Decorator<VertexLinkVector> CascadeLinksDecor("CascadeVertexLinks"); 
       SG::AuxElement::Decorator<MuonsLinkVector> MuonsLinksDecor("MuonsLinks");
       SG::AuxElement::Decorator<VertexLinkVector> DxLinksDecor("DxVertexLinks"); 
       SG::AuxElement::Decorator<float> chi2_decor("ChiSquared");
       SG::AuxElement::Decorator<float> ndof_decor("NumberDoF");
       SG::AuxElement::Decorator<float> Pt_decor("Pt");
       SG::AuxElement::Decorator<float> PtErr_decor("PtErr");
       SG::AuxElement::Decorator<float> Mass_svdecor("Dx_mass");
       SG::AuxElement::Decorator<float> MassErr_svdecor("Dx_massErr");
       SG::AuxElement::Decorator<float> Pt_svdecor("Dx_Pt");
       SG::AuxElement::Decorator<float> PtErr_svdecor("Dx_PtErr");
       SG::AuxElement::Decorator<float> Lxy_svdecor("Dx_Lxy");
       SG::AuxElement::Decorator<float> LxyErr_svdecor("Dx_LxyErr");
       SG::AuxElement::Decorator<float> Tau_svdecor("Dx_Tau");
       SG::AuxElement::Decorator<float> TauErr_svdecor("Dx_TauErr");
       SG::AuxElement::Decorator<float> Dx_chi2_svdecor("Dx_chi2");
       SG::AuxElement::Decorator<float> Dx_ndof_svdecor("Dx_ndof");
  
       //                K X1 X2
       // Ds+/-        : K K  π
       // D+/-         : K π  π
       // Lambda_c+/-  : K π  p
  
       SG::AuxElement::Decorator<float> massKX1_svdecor("KX_mass");
       SG::AuxElement::Decorator<float> massKX1X2_svdecor("KXpi_mass");
       SG::AuxElement::Decorator<float> RapidityKX1X2_svdecor("KXpi_Rapidity");
  
       SG::AuxElement::Decorator<float> MuMass_decor("Mu_mass");
       SG::AuxElement::Decorator<float> MuPt_decor("Mu_pt");
       SG::AuxElement::Decorator<float> MuEta_decor("Mu_eta");
       SG::AuxElement::Decorator<float> MuChi2_decor("Mu_chi2");
       SG::AuxElement::Decorator<float> MunDoF_decor("Mu_nDoF");
       //muon contribution to chi2 of the cascade fit
       SG::AuxElement::Decorator<float> MuChi2B_decor("Mu_chi2_B");
       SG::AuxElement::Decorator<float> MunDoFB_decor("Mu_nDoF_B");
  
       SG::AuxElement::Decorator<float> ChargeK_decor("K_charge");
       SG::AuxElement::Decorator<float> ChargeX1_decor("X_charge"); // K (Ds+) or pi (D+, Lambda_c+)
       SG::AuxElement::Decorator<float> ChargeX2_decor("Pi_charge"); // pi (Ds+, D+) or p (Lambda_c+)
       SG::AuxElement::Decorator<float> ChargeMu_decor("Mu_charge");
  
       ATH_MSG_DEBUG("cascadeinfoContainer size " << cascadeinfoContainer.size());
  
       // Get Muons container
       const xAOD::MuonContainer* muonContainer(nullptr);
       ATH_CHECK( evtStore()->retrieve(muonContainer, m_muonCollectionKey) );
       ATH_MSG_DEBUG("Muon container size "<<muonContainer->size());
         
       // Get D_(s)+/Lambda_c+ container and identify the input D_(s)+/Lambda_c+
       const xAOD::VertexContainer  *dxContainer(nullptr);
       ATH_CHECK(evtStore()->retrieve(dxContainer   , m_vertexDxContainerKey       ));
  
       for (Trk::VxCascadeInfo* x : cascadeinfoContainer) {
         if(x==nullptr){
           ATH_MSG_ERROR("cascadeinfoContainer is null");
           continue;
         }
  
         // the cascade fitter returns:
         // std::vector<xAOD::Vertex*>, each xAOD::Vertex contains the refitted track parameters (perigee at the vertex position)
         //   vertices[iv]              the links to the original TPs and a covariance of size 3+5*NTRK; the chi2 of the total fit
         //                             is split between the cascade vertices as per track contribution
         // std::vector< std::vector<TLorentzVector> >, each std::vector<TLorentzVector> contains the refitted momenta (TLorentzVector)
         //   momenta[iv][...]          of all tracks in the corresponding vertex, including any pseudotracks (from cascade vertices)
         //                             originating in this vertex; the masses are as assigned in the cascade fit
         // std::vector<Amg::MatrixX>,  the corresponding covariance matrices in momentum space
         //   covariance[iv]
         // int nDoF, double Chi2
         //
         // the invariant mass, pt, lifetime etc. errors should be calculated using the covariance matrices in momentum space as these
         // take into account the full track-track and track-vertex correlations
         //
         // in the case of Jpsi+V0: vertices[0] is the V0 vertex, vertices[1] is the B/Lambda_b(bar) vertex, containing the 2 Jpsi tracks.
         // The covariance terms between the two vertices are not stored. In momentum space momenta[0] contains the 2 V0 tracks,
         // their momenta add up to the momentum of the 3rd track in momenta[1], the first two being the Jpsi tracks
  
         const std::vector<xAOD::Vertex*> &cascadeVertices = x->vertices();
  
         if(cascadeVertices.size()!=topoN) ATH_MSG_ERROR("Incorrect number of vertices");
         if(cascadeVertices[0] == nullptr || cascadeVertices[1] == nullptr) ATH_MSG_ERROR("Error null vertex");
         // Keep vertices (bear in mind that they come in reverse order!)
         for(int i =0;i<topoN;i++) Vtxwritehandles[i]->push_back(cascadeVertices[i]);
         
         x->setSVOwnership(false); // Prevent Container from deleting vertices
         const auto mainVertex = cascadeVertices[1];   // this is the B_c+/- vertex
                                                       //and cascadeVertices[0] is Dx vertex
         const std::vector< std::vector<TLorentzVector> > &moms = x->getParticleMoms();
           
  
         // Set links to cascade vertices
         std::vector<const xAOD::Vertex*> verticestoLink;
         verticestoLink.push_back(cascadeVertices[0]);
         if(Vtxwritehandles[1] == nullptr) ATH_MSG_ERROR("Vtxwritehandles[1] is null");
         if(!BPhysPVCascadeTools::LinkVertices(CascadeLinksDecor, verticestoLink, Vtxwritehandles[0], cascadeVertices[1]))
             ATH_MSG_ERROR("Error decorating with cascade vertices");
  
         //----------------------------------- mu
         // Identify muons
         typedef std::vector<const xAOD::Muon*> MuonBag;
         MuonBag selectedMuons; selectedMuons.clear();
  
         for(const xAOD::Muon * mu : *muonContainer){
             const xAOD::TrackParticle* muonTrk = mu->trackParticle(xAOD::Muon::InnerDetectorTrackParticle );
             if (muonTrk == cascadeVertices[1]->trackParticle(0)) selectedMuons.push_back(mu); //there is always only one muon
         }
         ATH_MSG_DEBUG("selectedMuon size "<<selectedMuons.size()); //always only one muon
  
         //-----------------------------------
         // Create link for Muon
         // create tmp vector of preceding vertex links
         MuonsLinkVector preMuLinks;
  
         // loop over input precedingMuons
         auto muItr = selectedMuons.begin();
         for(; muItr != selectedMuons.end(); muItr++) {
             // create element link
             MuonsLink muLink;
             muLink.setElement(*muItr);
             muLink.setStorableObject(*muonContainer);
             
             // sanity check : is the link valid?
             if( !muLink.isValid() ) continue;
                
             // link is OK, store it in the tmp vector
             preMuLinks.push_back( muLink );
           
         } // end of loop over preceding vertices
            
         // all OK: store preceding vertex links in the aux store
         MuonsLinksDecor(*cascadeVertices[1]) = preMuLinks;
         //-----------------------------------
         //----------------------------------- mu
  
         //----------------------------------- Dx / Lambda_c
         // Identify the input D_(s)+ or Lambda_c+
         const xAOD::Vertex* dxVertex = BPhysPVCascadeTools::FindVertex<3>(dxContainer, cascadeVertices[0]);
         ATH_MSG_DEBUG("1 pt D_(s)+/Lambda_c+ tracks " << cascadeVertices[1]->trackParticle(0)->pt()<<", "
                       << cascadeVertices[0]->trackParticle(0)->pt()<< ", "<< cascadeVertices[0]->trackParticle(1)->pt()<< ", "
                       << cascadeVertices[0]->trackParticle(2)->pt());
         if (dxVertex) ATH_MSG_DEBUG("2 pt D_(s)+/Lambda_c+ tracks " << dxVertex->trackParticle(0)->pt() << ", "
                       << dxVertex->trackParticle(1)->pt() << ", " << dxVertex->trackParticle(2)->pt());
         // Set links to input vertices
         std::vector<const xAOD::Vertex*> dxVerticestoLink;
         if (dxVertex) dxVerticestoLink.push_back(dxVertex);
         else ATH_MSG_WARNING("Could not find linking D_(s)+/Lambda_c+");
         if(!BPhysPVCascadeTools::LinkVertices(DxLinksDecor, dxVerticestoLink, dxContainer, cascadeVertices[1]))
             ATH_MSG_ERROR("Error decorating with D_(s)+/Lambda_c+ vertices");
  
         // tag for D- and Lambda_c-
         bool tagDp(true);
         if (dxVertex) {
           if((std::abs(m_Dx_pid)==411 || std::abs(m_Dx_pid)==4122) && (dxVertex->trackParticle(2)->charge()==-1)) tagDp = false;
         }
         //----------------------------------- Dx / Lambda_c
           
         double mass_b = m_vtx0MassHypo;
         double mass_d = m_vtx1MassHypo;
         std::vector<double> massesMu;
         massesMu.push_back(m_vtx0Daug1MassHypo);        // µ
         std::vector<double> massesDx;
         if(tagDp){
           massesDx.push_back(m_vtx1Daug1MassHypo);      // K (Ds+/-) or π (D+, Lambda_c+)
           massesDx.push_back(m_vtx1Daug2MassHypo);      // K
         }else{ // Change the order for D- or Lambda_c-
           massesDx.push_back(m_vtx1Daug2MassHypo);      // K
           massesDx.push_back(m_vtx1Daug1MassHypo);      // π
         }
         massesDx.push_back(m_vtx1Daug3MassHypo);        // π (Ds+/-, D+-) or p (Lambda_c+)
         std::vector<double> Masses;
         Masses.push_back(m_vtx0Daug1MassHypo);          // µ
         Masses.push_back(m_vtx1MassHypo);               // Dx / Lambda_c
  
         // loop over candidates -- Don't apply PV_minNTracks requirement here
         // because it may result in exclusion of the high-pt PV.
         // get good PVs
  
  
         xAOD::BPhysHypoHelper vtx(m_hypoName, mainVertex);
  
         // Get refitted track momenta from all vertices, charged tracks only
         BPhysPVCascadeTools::SetVectorInfo(vtx, x);
  
         // Decorate main vertex
         //
         // 1.a) mass, mass error
         BPHYS_CHECK( vtx.setMass(m_CascadeTools->invariantMass(moms[1])) );
         BPHYS_CHECK( vtx.setMassErr(m_CascadeTools->invariantMassError(moms[1],x->getCovariance()[1])) );
         // 1.b) pt and pT error (the default pt of mainVertex is != the pt of the full cascade fit!)
         Pt_decor(*mainVertex) = m_CascadeTools->pT(moms[1]);
         PtErr_decor(*mainVertex) = m_CascadeTools->pTError(moms[1],x->getCovariance()[1]);
         // 1.c) chi2 and ndof (the default chi2 of mainVertex is != the chi2 of the full cascade fit!)
         chi2_decor(*mainVertex) = x->fitChi2();
         ndof_decor(*mainVertex) = x->nDoF();
         Dx_chi2_svdecor(*mainVertex) = cascadeVertices[0]->chiSquared();
         Dx_ndof_svdecor(*mainVertex) = cascadeVertices[0]->numberDoF();
         
         //--------------------- mu
         float muM = 0., mupt = 0., mueta = 0.;
         if (selectedMuons.size()==1) {
           TLorentzVector  p4_mu1;
           p4_mu1.SetPtEtaPhiM(selectedMuons.at(0)->pt(),
                               selectedMuons.at(0)->eta(),
                               selectedMuons.at(0)->phi(), m_vtx0Daug1MassHypo);
           muM = p4_mu1.M();
           mupt = p4_mu1.Pt();
           mueta = p4_mu1.Eta();
         }
         MuMass_decor(*mainVertex) = muM;
         MuPt_decor(*mainVertex) = mupt;
         MuEta_decor(*mainVertex) = mueta;
         ChargeMu_decor(*mainVertex) = selectedMuons.at(0)->charge();
         MuChi2_decor(*mainVertex) = cascadeVertices[1]->trackParticle(0)->chiSquared();
         MunDoF_decor(*mainVertex) = cascadeVertices[1]->trackParticle(0)->numberDoF();
         
         // track contribution to chi2 of cascasde fit
         std::vector< Trk::VxTrackAtVertex > trkAtB = cascadeVertices[1]->vxTrackAtVertex();
         MuChi2B_decor(*mainVertex) = trkAtB.at(0).trackQuality().chiSquared();
         MunDoFB_decor(*mainVertex) = trkAtB.at(0).trackQuality().numberDoF();
  
         //--------------------- Dx / Lambda_c
         //tagDp = true by default, for Ds+/-, D+ and Lambda_c+
         float massKX1 = 0.;
         float massKX1X2 = 0.;
         float RapidityKX1X2 = 0.;
         if (dxVertex) {
           TLorentzVector  p4_h1, p4_h2, p4_h3;
           if(tagDp){
             p4_h1.SetPtEtaPhiM(dxVertex->trackParticle(0)->pt(), 
                                dxVertex->trackParticle(0)->eta(),
                                dxVertex->trackParticle(0)->phi(), m_vtx1Daug1MassHypo);
             ChargeX1_decor(*mainVertex) = dxVertex->trackParticle(0)->charge();
             p4_h2.SetPtEtaPhiM(dxVertex->trackParticle(1)->pt(), 
                                dxVertex->trackParticle(1)->eta(),
                                dxVertex->trackParticle(1)->phi(), m_vtx1Daug2MassHypo);
             ChargeK_decor(*mainVertex) = dxVertex->trackParticle(1)->charge();
           }else{ // Change the order for D- and Lambda_c-
             p4_h1.SetPtEtaPhiM(dxVertex->trackParticle(0)->pt(), 
                                dxVertex->trackParticle(0)->eta(),
                                dxVertex->trackParticle(0)->phi(), m_vtx1Daug2MassHypo);
             ChargeK_decor(*mainVertex) = dxVertex->trackParticle(0)->charge();
  
             p4_h2.SetPtEtaPhiM(dxVertex->trackParticle(1)->pt(), 
                                dxVertex->trackParticle(1)->eta(),
                                dxVertex->trackParticle(1)->phi(), m_vtx1Daug1MassHypo);
             ChargeX1_decor(*mainVertex) = dxVertex->trackParticle(1)->charge(); // K (Ds+) or π (D+)
           }
             p4_h3.SetPtEtaPhiM(dxVertex->trackParticle(2)->pt(), 
                                dxVertex->trackParticle(2)->eta(),
                                dxVertex->trackParticle(2)->phi(), m_vtx1Daug3MassHypo);
             ChargeX2_decor(*mainVertex) = dxVertex->trackParticle(2)->charge(); // π (Ds+, D+) or p (Lambda_c+)
             
             massKX1 = (p4_h1 + p4_h2).M();
             massKX1X2 = (p4_h1 + p4_h2 + p4_h3).M();
             RapidityKX1X2 = (p4_h1 + p4_h2 + p4_h3).Rapidity();
         }
         massKX1_svdecor(*mainVertex) = massKX1;
         massKX1X2_svdecor(*mainVertex) = massKX1X2;
         RapidityKX1X2_svdecor(*mainVertex) = RapidityKX1X2;
           
  
         ATH_CHECK(helper.FillCandwithRefittedVertices(m_refitPV, pvContainer, 
                                     refPvContainer, &(*m_pvRefitter), m_PV_max, m_DoVertexType, x, 1, mass_b, vtx));
  
  
         // 4) decorate the main vertex with V0 vertex mass, pt, lifetime and lxy values (plus errors) 
         // V0 points to the main vertex, so lifetime and lxy are w.r.t the main vertex
         Mass_svdecor(*mainVertex) = m_CascadeTools->invariantMass(moms[0]);
         MassErr_svdecor(*mainVertex) = m_CascadeTools->invariantMassError(moms[0],x->getCovariance()[0]);
         Pt_svdecor(*mainVertex) = m_CascadeTools->pT(moms[0]);
         PtErr_svdecor(*mainVertex) = m_CascadeTools->pTError(moms[0],x->getCovariance()[0]);
         Lxy_svdecor(*mainVertex) = m_CascadeTools->lxy(moms[0],cascadeVertices[0],cascadeVertices[1]);
         LxyErr_svdecor(*mainVertex) = m_CascadeTools->lxyError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1]);
         Tau_svdecor(*mainVertex) = m_CascadeTools->tau(moms[0],cascadeVertices[0],cascadeVertices[1]);
         TauErr_svdecor(*mainVertex) = m_CascadeTools->tauError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1]);
  
         // Some checks in DEBUG mode
         ATH_MSG_DEBUG("chi2 " << x->fitChi2()
                   << " chi2_1 " << m_V0Tools->chisq(cascadeVertices[0])
                   << " chi2_2 " << m_V0Tools->chisq(cascadeVertices[1])
                   << " vprob " << m_CascadeTools->vertexProbability(x->nDoF(),x->fitChi2()));
         ATH_MSG_DEBUG("ndf " << x->nDoF() << " ndf_1 " << m_V0Tools->ndof(cascadeVertices[0]) << " ndf_2 " << m_V0Tools->ndof(cascadeVertices[1]));
         ATH_MSG_DEBUG("V0Tools mass_d " << m_V0Tools->invariantMass(cascadeVertices[0],massesDx)
                    << " error " << m_V0Tools->invariantMassError(cascadeVertices[0],massesDx)
                    << " mass_J " << m_V0Tools->invariantMass(cascadeVertices[1],massesMu)
                    << " error " << m_V0Tools->invariantMassError(cascadeVertices[1],massesMu));
         // masses and errors, using track masses assigned in the fit
         double Mass_B = m_CascadeTools->invariantMass(moms[1]);
         double Mass_D = m_CascadeTools->invariantMass(moms[0]);
         double Mass_B_err = m_CascadeTools->invariantMassError(moms[1],x->getCovariance()[1]);
         double Mass_D_err = m_CascadeTools->invariantMassError(moms[0],x->getCovariance()[0]);
         ATH_MSG_DEBUG("Mass_B " << Mass_B << " Mass_D " << Mass_D);
         ATH_MSG_DEBUG("Mass_B_err " << Mass_B_err << " Mass_D_err " << Mass_D_err);
         double mprob_B = m_CascadeTools->massProbability(mass_b,Mass_B,Mass_B_err);
         double mprob_D = m_CascadeTools->massProbability(mass_d,Mass_D,Mass_D_err);
         ATH_MSG_DEBUG("mprob_B " << mprob_B << " mprob_D " << mprob_D);
         // masses and errors, assigning user defined track masses
         ATH_MSG_DEBUG("Mass_b " << m_CascadeTools->invariantMass(moms[1],Masses)
                   << " Mass_d " << m_CascadeTools->invariantMass(moms[0],massesDx));
         ATH_MSG_DEBUG("Mass_b_err " << m_CascadeTools->invariantMassError(moms[1],x->getCovariance()[1],Masses)
                   << " Mass_d_err " << m_CascadeTools->invariantMassError(moms[0],x->getCovariance()[0],massesDx));
         ATH_MSG_DEBUG("pt_b " << m_CascadeTools->pT(moms[1])
                   << " pt_d " << m_CascadeTools->pT(moms[0])
                   << " pt_dp " << m_V0Tools->pT(cascadeVertices[0]));
         ATH_MSG_DEBUG("ptErr_b " << m_CascadeTools->pTError(moms[1],x->getCovariance()[1])
                   << " ptErr_d " << m_CascadeTools->pTError(moms[0],x->getCovariance()[0])
                   << " ptErr_dp " << m_V0Tools->pTError(cascadeVertices[0]));
         ATH_MSG_DEBUG("lxy_B " << m_V0Tools->lxy(cascadeVertices[1],primaryVertex) << " lxy_D " << m_V0Tools->lxy(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("lxy_b " << m_CascadeTools->lxy(moms[1],cascadeVertices[1],primaryVertex) << " lxy_d " << m_CascadeTools->lxy(moms[0],cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("lxyErr_b " << m_CascadeTools->lxyError(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex)
                   << " lxyErr_d " << m_CascadeTools->lxyError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1])
                   << " lxyErr_dp " << m_V0Tools->lxyError(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("tau_B " << m_CascadeTools->tau(moms[1],cascadeVertices[1],primaryVertex,mass_b)
                    << " tau_dp " << m_V0Tools->tau(cascadeVertices[0],cascadeVertices[1],massesDx));
         ATH_MSG_DEBUG("tau_b " << m_CascadeTools->tau(moms[1],cascadeVertices[1],primaryVertex)
                    << " tau_d " << m_CascadeTools->tau(moms[0],cascadeVertices[0],cascadeVertices[1])
                    << " tau_D " << m_CascadeTools->tau(moms[0],cascadeVertices[0],cascadeVertices[1],mass_d));
         ATH_MSG_DEBUG("tauErr_b " << m_CascadeTools->tauError(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex)
                   << " tauErr_d " << m_CascadeTools->tauError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1])
                   << " tauErr_dp " << m_V0Tools->tauError(cascadeVertices[0],cascadeVertices[1],massesDx));
         ATH_MSG_DEBUG("TauErr_b " << m_CascadeTools->tauError(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex,mass_b)
                   << " TauErr_d " << m_CascadeTools->tauError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1],mass_d)
                   << " TauErr_dp " << m_V0Tools->tauError(cascadeVertices[0],cascadeVertices[1],massesDx,mass_d));
  
         ATH_MSG_DEBUG("CascadeTools main vert wrt PV " << " CascadeTools SV " << " V0Tools SV");
         ATH_MSG_DEBUG("a0z " << m_CascadeTools->a0z(moms[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0z(moms[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0z(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("a0zErr " << m_CascadeTools->a0zError(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0zError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0zError(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("a0xy " << m_CascadeTools->a0xy(moms[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0xy(moms[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0xy(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("a0xyErr " << m_CascadeTools->a0xyError(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0xyError(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0xyError(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("a0 " << m_CascadeTools->a0(moms[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0(moms[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("a0Err " << m_CascadeTools->a0Error(moms[1],x->getCovariance()[1],cascadeVertices[1],primaryVertex)
                    << ", " << m_CascadeTools->a0Error(moms[0],x->getCovariance()[0],cascadeVertices[0],cascadeVertices[1])
                    << ", " << m_V0Tools->a0Error(cascadeVertices[0],cascadeVertices[1]));
         ATH_MSG_DEBUG("x0 " << m_V0Tools->vtx(cascadeVertices[0]).x() << " y0 " << m_V0Tools->vtx(cascadeVertices[0]).y() << " z0 " << m_V0Tools->vtx(cascadeVertices[0]).z());
         ATH_MSG_DEBUG("x1 " << m_V0Tools->vtx(cascadeVertices[1]).x() << " y1 " << m_V0Tools->vtx(cascadeVertices[1]).y() << " z1 " << m_V0Tools->vtx(cascadeVertices[1]).z());
         ATH_MSG_DEBUG("X0 " << primaryVertex->x() << " Y0 " << primaryVertex->y() << " Z0 " << primaryVertex->z());
         ATH_MSG_DEBUG("rxy0 " << m_V0Tools->rxy(cascadeVertices[0]) << " rxyErr0 " << m_V0Tools->rxyError(cascadeVertices[0]));
         ATH_MSG_DEBUG("rxy1 " << m_V0Tools->rxy(cascadeVertices[1]) << " rxyErr1 " << m_V0Tools->rxyError(cascadeVertices[1]));
         ATH_MSG_DEBUG("Rxy0 wrt PV " << m_V0Tools->rxy(cascadeVertices[0],primaryVertex) << " RxyErr0 wrt PV " << m_V0Tools->rxyError(cascadeVertices[0],primaryVertex));
         ATH_MSG_DEBUG("Rxy1 wrt PV " << m_V0Tools->rxy(cascadeVertices[1],primaryVertex) << " RxyErr1 wrt PV " << m_V0Tools->rxyError(cascadeVertices[1],primaryVertex));
         ATH_MSG_DEBUG("number of covariance matrices " << (x->getCovariance()).size());
           
       } // loop over cascadeinfoContainer
  
       // Deleting cascadeinfo since this won't be stored.
       // Vertices have been kept in m_cascadeOutputs and should be owned by their container
       for (auto x : cascadeinfoContainer)  delete x;
  
       return StatusCode::SUCCESS;
     }

◆ 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.

95 { return m_detStore; }

◆ 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.

85 { 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.

90 { 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

◆ FindVertexTrack()

static xAOD::Vertex* DerivationFramework::MuPlusDsCascade::FindVertexTrack	(	const xAOD::MuonContainer *	c,
		xAOD::Vertex *	v
	)

static

◆ initialize()

StatusCode DerivationFramework::MuPlusDsCascade::initialize ( )

overridevirtual

Definition at line 38 of file MuPlusDsCascade.cxx.

                                            {
  
         // retrieving vertex Fitter
         ATH_CHECK( m_iVertexFitter.retrieve());
         
         // retrieving the V0 tools //from JpsiFinder
         ATH_CHECK( m_V0Tools.retrieve());
         
         // Get the track selector tool from ToolSvc
         if ( m_trkSelector.retrieve().isFailure() ) {
             ATH_MSG_FATAL("Failed to retrieve tool " << m_trkSelector);
             return StatusCode::FAILURE;
         } else {
             ATH_MSG_INFO("Retrieved tool " << m_trkSelector);
         }
  
         //======================== inDetTrack selection tool ==================
         m_trackSelectionTools = std::make_unique<InDet::InDetTrackSelectionTool>("TrackSelector");
         ANA_CHECK(m_trackSelectionTools->setProperty("CutLevel", "LoosePrimary"));
         ANA_CHECK(m_trackSelectionTools->initialize() );
         //=====================================================================
  
         // retrieving the Cascade tools
         ATH_CHECK( m_CascadeTools.retrieve());
  
         ATH_CHECK( m_eventInfo_key.initialize() );
  
         ATH_CHECK( m_partPropSvc.retrieve() );
         auto pdt = m_partPropSvc->PDT();
         
         // Ds+/-        : K K π
         // D+/-         : K π π
         // Lambda_c+/-  : K π p
  
         // retrieve particle masses
         if(m_vtx0MassHypo < 0.)
           m_vtx0MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::BCPLUS);
         if(m_vtx1MassHypo < 0.) {
           if(std::abs(m_Dx_pid) == 411) m_vtx1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::DPLUS);
           if(std::abs(m_Dx_pid) == 431) m_vtx1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::DSPLUS);
           if(std::abs(m_Dx_pid) == 4122) m_vtx1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::LAMBDACPLUS);
         }
  
         if(m_vtx0Daug1MassHypo < 0.) m_vtx0Daug1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::MUON);
         if(m_vtx1Daug1MassHypo < 0.) {
            if(std::abs(m_Dx_pid) == 431) m_vtx1Daug1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::KPLUS); //Ds+
            else m_vtx1Daug1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS); //D+, Lambda_c+
         }
         if(m_vtx1Daug2MassHypo < 0.) m_vtx1Daug2MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::KPLUS);
         if(m_vtx1Daug3MassHypo < 0.) {
             if(std::abs(m_Dx_pid) == 4122) m_vtx1Daug3MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PROTON); //Lambda_c+
             else m_vtx1Daug3MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS); //Ds+, D+
         }
         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.

◆ interfaceID()

static const InterfaceID& DerivationFramework::MuPlusDsCascade::interfaceID ( )

inlinestatic

Definition at line 44 of file MuPlusDsCascade.h.

44 { return IID_MuPlusDsCascade;}

◆ 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);
   }

◆ 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.

◆ performSearch()

StatusCode DerivationFramework::MuPlusDsCascade::performSearch ( std::vector< Trk::VxCascadeInfo * > * cascadeinfoContainer ) const

Definition at line 582 of file MuPlusDsCascade.cxx.

     {
         ATH_MSG_DEBUG( "MuPlusDsCascade::performSearch" );
         assert(cascadeinfoContainer!=nullptr);
  
         // Get TrackParticle container (for setting links to the original tracks)
         const xAOD::TrackParticleContainer  *trackContainer(nullptr);
         ATH_CHECK(evtStore()->retrieve(trackContainer   , "InDetTrackParticles"      ));
  
         // Get V0 container
         const xAOD::VertexContainer  *dxContainer(nullptr);
         ATH_CHECK(evtStore()->retrieve(dxContainer   , m_vertexDxContainerKey       ));
  
         double mass_d = m_vtx1MassHypo; 
         std::vector<const xAOD::TrackParticle*> tracksMu;
         std::vector<const xAOD::TrackParticle*> tracksDx;
         std::vector<double> massesMu;
         massesMu.push_back(m_vtx0Daug1MassHypo); // µ
         std::vector<double> massesDx;
         massesDx.push_back(m_vtx1Daug1MassHypo); // K (Ds+/-) or π (D+, Lambda_c+)
         massesDx.push_back(m_vtx1Daug2MassHypo); // K
         massesDx.push_back(m_vtx1Daug3MassHypo); // π (Ds+/-, D+) ot p (Lambda_c+)
         std::vector<double> massesDm; // Alter the order of masses for D- and Lambda_c-
         massesDm.push_back(m_vtx1Daug2MassHypo); // K
         massesDm.push_back(m_vtx1Daug1MassHypo); // π
         massesDm.push_back(m_vtx1Daug3MassHypo); // π (D-) or p (Lambda_c-)
         std::vector<double> Masses;
         Masses.push_back(m_vtx0Daug1MassHypo);   // µ
         Masses.push_back(m_vtx1MassHypo);        // Dx / Lambda_c
         
         //-------------------------------------------------------------------------------
         // Retrieving and selecting muons (taken from JpsiFinder.cxx)
         // Get the muons from StoreGate
         const xAOD::MuonContainer* importedMuonCollection;
         StatusCode sc = evtStore()->retrieve(importedMuonCollection,m_muonCollectionKey);
  
         if(sc.isFailure()){
             ATH_MSG_WARNING("No muon collection with key " << m_muonCollectionKey << " found in StoreGate");
             return StatusCode::SUCCESS;;
         }else{
             ATH_MSG_DEBUG("Found muon collections with key "<<m_muonCollectionKey);
         }
         ATH_MSG_DEBUG("Muon container size "<<importedMuonCollection->size());
         
         // Typedef for vectors of muons
         typedef std::vector<const xAOD::Muon*> MuonBag;
  
         // Select the muons
         const xAOD::Vertex* vx = 0;
         MuonBag theMuonsAfterSelection;
         for (auto mu : *importedMuonCollection) {
             if ( mu == nullptr ) continue;
             if (!mu->inDetTrackParticleLink().isValid()) continue; // No muons without ID tracks
             const xAOD::TrackParticle* muonTrk = *(mu->inDetTrackParticleLink());
             if ( muonTrk==nullptr) continue;
             if ( !m_trkSelector->decision(*muonTrk, vx) ) continue; // all ID tracks must pass basic tracking cuts
             if ( std::fabs(muonTrk->pt())<m_thresholdPt ) continue; // higher pt cut if needed
             if ( m_mcpCuts && !mu->passesIDCuts()) continue; // cuts of the MCP group recommendation
             if ( m_combOnly && mu->muonType() != xAOD::Muon::Combined ) continue; // require combined muons
             if ( mu->muonType() == xAOD::Muon::SiliconAssociatedForwardMuon && !m_useCombMeasurement) continue;
                 
             theMuonsAfterSelection.push_back(mu);
         }
         if (theMuonsAfterSelection.size() == 0) return StatusCode::SUCCESS;;
         ATH_MSG_DEBUG("Number of muons after selection: " << theMuonsAfterSelection.size());
         //-------------------------------------------------------------------------------
  
         //-------------------------------------------------------------------------------
         // Select the D_s+/D+/Lambda_c+ candidates before calling cascade fit
         std::vector<const xAOD::Vertex*> selectedDxCandidates;
         
  
         for(auto vxcItr : *dxContainer){
  
            // Check the passed flag first
           const xAOD::Vertex* vtx = vxcItr;
            if(std::abs(m_Dx_pid)==431) { // D_s+/-
                SG::AuxElement::Accessor<Char_t> flagAcc1("passed_Ds");
                if(flagAcc1.isAvailable(*vtx)){
                   if(!flagAcc1(*vtx)) continue;
                }
            }
  
            if(std::abs(m_Dx_pid)==411) { // D+/-
                SG::AuxElement::Accessor<Char_t> flagAcc1("passed_Dp");
                SG::AuxElement::Accessor<Char_t> flagAcc2("passed_Dm");
                bool isDp(true);
                bool isDm(true);
                if(flagAcc1.isAvailable(*vtx)){
                   if(!flagAcc1(*vtx)) isDp = false;
                }
                if(flagAcc2.isAvailable(*vtx)){
                   if(!flagAcc2(*vtx)) isDm = false;
                }
                if(!(isDp||isDm)) continue;
            } 
  
            // Track selection - Loose
            if ( !m_trackSelectionTools->accept(vxcItr->trackParticle(0)) ){
                 ATH_MSG_DEBUG(" Original Dx/Lambda_c candidate rejected by the track's cut level - loose ");
                 continue;
            }
            if ( !m_trackSelectionTools->accept(vxcItr->trackParticle(1)) ){
                 ATH_MSG_DEBUG(" Original Dx/Lambda_c candidate rejected by the track's cut level - loose ");
                 continue;
            }
            if ( !m_trackSelectionTools->accept(vxcItr->trackParticle(2)) ){
                 ATH_MSG_DEBUG(" Original Dx/Lambda_c candidate rejected by the track's cut level - loose ");
                 continue;
            }
             
            // Ensure the total charge is correct
            if( std::abs( vxcItr->trackParticle(0)->charge()+vxcItr->trackParticle(1)->charge()+vxcItr->trackParticle(2)->charge() ) != 1 ){
                ATH_MSG_DEBUG(" Original Dx/Lambda_c candidate rejected by the charge requirement: "
                               << vxcItr->trackParticle(0)->charge() << ", " << vxcItr->trackParticle(1)->charge() << ", " << vxcItr->trackParticle(2)->charge() );
                continue;
            }
  
            // Check D_(s)/Lambda_c +/- candidate invariant mass and skip if need be
            double mass_D;
            if( (std::abs(m_Dx_pid) == 411 || std::abs(m_Dx_pid) == 4122) && vxcItr->trackParticle(2)->charge()<0) // D- & Lambda_c-
                mass_D = m_V0Tools->invariantMass(vxcItr,massesDm);
            else // D+, D_s+/-, Lambda_c+
                mass_D = m_V0Tools->invariantMass(vxcItr,massesDx);
            ATH_MSG_DEBUG("D_(s)/Lambda_c mass " << mass_D);
            if(mass_D < m_DxMassLower || mass_D > m_DxMassUpper) {
                ATH_MSG_DEBUG(" Original D_(s)/Lambda_c candidate rejected by the mass cut: mass = "
                               << mass_D << " != (" << m_DxMassLower << ", " << m_DxMassUpper << ")" );
                continue;
            }
  
            // Add loose cut on phi(K+K-) mass for D_s->phi π
            if(std::abs(m_Dx_pid)==431){
               TLorentzVector p4Kp_in, p4Km_in;
               p4Kp_in.SetPtEtaPhiM( vxcItr->trackParticle(0)->pt(),
                                     vxcItr->trackParticle(0)->eta(),
                                     vxcItr->trackParticle(0)->phi(), m_vtx1Daug1MassHypo);
               p4Km_in.SetPtEtaPhiM( vxcItr->trackParticle(1)->pt(),
                                     vxcItr->trackParticle(1)->eta(),
                                     vxcItr->trackParticle(1)->phi(), m_vtx1Daug2MassHypo);
               double mass_phi = (p4Kp_in + p4Km_in).M();
               ATH_MSG_DEBUG("phi mass " << mass_phi);
               if(mass_phi > 1600) { //1st cut on phi(K+K-) mass, loose one
                   ATH_MSG_DEBUG(" Original phi candidate rejected by the mass cut: mass = " << mass_phi );
                   continue;
                   
               }
            }
             
            // Cut on pT(p) for Lambda_c
            if(std::abs(m_Dx_pid)==4122){
                if (vxcItr->trackParticle(2)->pt() < 2400) continue;
            }
             
            // L.G., 26.03.21: Check there are no previous D+/- combinations with inverse pion positions
            if(std::abs(m_Dx_pid)==411 && selectedDxCandidates.size()>0) {
                bool Dpmcopy(false);
                for(auto dxItr : selectedDxCandidates){
                    if(vxcItr->trackParticle(2)/*π*/->charge()>0) { // D+
                        if ( vxcItr->trackParticle(2) == dxItr->trackParticle(0) && vxcItr->trackParticle(0) == dxItr->trackParticle(2)
                            && vxcItr->trackParticle(1) == dxItr->trackParticle(1) ) Dpmcopy = true;
                    } else { // D-
                        if ( vxcItr->trackParticle(2) == dxItr->trackParticle(1) && vxcItr->trackParticle(1) == dxItr->trackParticle(2)
                            && vxcItr->trackParticle(0) == dxItr->trackParticle(0) ) Dpmcopy = true;
                    }
                }
                if (Dpmcopy) continue;
            }
             
            selectedDxCandidates.push_back(vxcItr);
         } // end for(auto vxcItr : dxContainer)
         if(selectedDxCandidates.size()<1) return StatusCode::SUCCESS;
         //-------------------------------------------------------------------------------
         
         //-------------------------------------------------------------------------------
         // Select mu D_(s)+/Lambda_c candidates
         // Iterate over muons
         for(auto muItr : theMuonsAfterSelection){
             tracksMu.clear();
             //Convert to trackParticle base
             auto TrkMuon = muItr->trackParticle( xAOD::Muon::InnerDetectorTrackParticle );
             tracksMu.push_back(TrkMuon);
             if (tracksMu.size() != 1 || massesMu.size() != 1 ) {
               ATH_MSG_WARNING("Problems with muon input");
             }
  
             // Iterate over D_(s)/Lambda_c +/-  vertices
             for(auto dxItr : selectedDxCandidates){
               // Check identical tracks in input
               if(std::find(tracksMu.cbegin(), tracksMu.cend(), dxItr->trackParticle(0)) != tracksMu.cend()) continue;
               if(std::find(tracksMu.cbegin(), tracksMu.cend(), dxItr->trackParticle(1)) != tracksMu.cend()) continue;
               if(std::find(tracksMu.cbegin(), tracksMu.cend(), dxItr->trackParticle(2)) != tracksMu.cend()) continue;
                
               size_t dxTrkNum = dxItr->nTrackParticles();
               tracksDx.clear();
               for( unsigned int it=0; it<dxTrkNum; it++) tracksDx.push_back(dxItr->trackParticle(it));
               if (tracksDx.size() != 3 || massesDx.size() != 3 ) {
                 ATH_MSG_WARNING("Problems with D_(s)/Lambda_c +/- input");
               }
  
               if(std::find(tracksDx.cbegin(), tracksDx.cend(), TrkMuon) != tracksDx.cend()) continue; //looking for the muons in the D+ tracks
                
               ATH_MSG_DEBUG("Using tracks" << tracksMu[0] << ", " << tracksDx[0] << ", " << tracksDx[1] << ", " << tracksDx[2]);
                // Apply the user's settings to the fitter
               // Reset
               std::unique_ptr<Trk::IVKalState> state (m_iVertexFitter->makeState());
               // Robustness
               int robustness = 0;
               m_iVertexFitter->setRobustness(robustness, *state);
               // Build up the topology
               // Vertex list
               std::vector<Trk::VertexID> vrtList;
               // D_(s)/Lambda_c +/- vertex
               Trk::VertexID vID;
               if (m_constrDx) {
                 if((std::abs(m_Dx_pid)==411 || std::abs(m_Dx_pid)==4122 ) && dxItr->trackParticle(2)->charge()<0) // D-, Lambda_c-
                   vID = m_iVertexFitter->startVertex(tracksDx,massesDm, *state,mass_d);
                 else // D+, D_s+/-, Lambda_c+
                   vID = m_iVertexFitter->startVertex(tracksDx,massesDx, *state,mass_d);
               } else {
                 if((std::abs(m_Dx_pid)==411 || std::abs(m_Dx_pid)==4122 ) && dxItr->trackParticle(2)->charge()<0) // D-, Lambda_c-
                   vID = m_iVertexFitter->startVertex(tracksDx,massesDm, *state);
                 else // D+, D_s+/-, Lambda_c+
                   vID = m_iVertexFitter->startVertex(tracksDx,massesDx, *state);
               }
               vrtList.push_back(vID);
                
               // B vertex including muon
               //For mass constraint, not used yet
               Trk::VertexID vID2 = m_iVertexFitter->nextVertex(tracksMu,massesMu,vrtList, *state);
               ATH_MSG_DEBUG(vID2);
                 
               // Do the work
               std::unique_ptr<Trk::VxCascadeInfo> result(m_iVertexFitter->fitCascade(*state));
  
               if (result != nullptr) {
                 // reset links to original tracks
                 BPhysPVCascadeTools::PrepareVertexLinks(result.get(), trackContainer);
                 ATH_MSG_DEBUG("storing tracks " << ((result->vertices())[0] /*D_(s)+*/)->trackParticle(0) << ", "
                                                 << ((result->vertices())[0])->trackParticle(1) << ", "
                                                 << ((result->vertices())[0])->trackParticle(2) << ", "
                                                 << ((result->vertices())[1]/*mu*/)->trackParticle(0));
                                               
                 // necessary to prevent memory leak
                 result->setSVOwnership(true);
                   
                 // Chi2/DOF cut
                 double bChi2DOF = result->fitChi2()/result->nDoF();
                 ATH_MSG_DEBUG("Candidate chi2/DOF is " << bChi2DOF);
                 bool chi2CutPassed = (m_chi2cut <= 0.0 || bChi2DOF < m_chi2cut);
  
                 const std::vector<xAOD::Vertex*> &cascadeVertices = result->vertices();
                 const std::vector< std::vector<TLorentzVector> > &moms = result->getParticleMoms();
                   
                 //----------------------------------------------------
                 // retrieve primary vertices
                 //----------------------------------------------------
                 const xAOD::Vertex * primaryVertex(nullptr);
                 const xAOD::VertexContainer *pvContainer(nullptr);
                 ATH_CHECK(evtStore()->retrieve(pvContainer, m_VxPrimaryCandidateName));
                 ATH_MSG_DEBUG("Found " << m_VxPrimaryCandidateName << " in StoreGate!");
  
                 if (pvContainer->size()==0){
                     ATH_MSG_WARNING("You have no primary vertices: " << pvContainer->size());
                     return StatusCode::RECOVERABLE;
                 } else {
                     primaryVertex = (*pvContainer)[0];
                 }
                   
                 // Add stronger cut on phi(K+K-) mass for D_s->phi π after cascade fit
                 if(std::abs(m_Dx_pid)==431){
                     TLorentzVector p4Kp_in, p4Km_in;
                     p4Kp_in.SetPtEtaPhiM( cascadeVertices[0]->trackParticle(0)->pt(),
                                           cascadeVertices[0]->trackParticle(0)->eta(),
                                           cascadeVertices[0]->trackParticle(0)->phi(), m_vtx1Daug1MassHypo);
                     p4Km_in.SetPtEtaPhiM( cascadeVertices[0]->trackParticle(1)->pt(),
                                           cascadeVertices[0]->trackParticle(1)->eta(),
                                           cascadeVertices[0]->trackParticle(1)->phi(), m_vtx1Daug2MassHypo);
                     double mass_phi = (p4Kp_in + p4Km_in).M();
                     ATH_MSG_DEBUG("phi mass " << mass_phi);
                     if(mass_phi > 1100) {
                         ATH_MSG_DEBUG(" Original phi candidate rejected by the mass cut: mass = " << mass_phi );
                         continue;
                     }
                 }
  
                 double mass = m_CascadeTools->invariantMass(moms[1]);
                 if(chi2CutPassed) {
                   if (mass >= m_MassLower && mass <= m_MassUpper) {
                       if (m_CascadeTools->lxy(moms[1],cascadeVertices[1],primaryVertex) > 0.09){ //B_Lxy
                           if (m_CascadeTools->pT(moms[1]) > 9500){ //B_pT
                               if (m_CascadeTools->lxy(moms[0],cascadeVertices[0],cascadeVertices[1]) > 0){ //Dx_lxy
                                   if (std::abs(m_Dx_pid)!=411) cascadeinfoContainer->push_back(result.release());
                                   else if (m_CascadeTools->lxy(moms[0],cascadeVertices[0],cascadeVertices[1]) > 0.09){ //D+_Lxy
                                       cascadeinfoContainer->push_back(result.release());
                                   }//D+_Lxy
                               }//Dx_lxy
                           } //B_pT
                       }//B_Lxy
                   } else {
                     ATH_MSG_DEBUG("Candidate rejected by the mass cut: mass = "
                                   << mass << " != (" << m_MassLower << ", " << m_MassUpper << ")" );
                   }
                 } // chi2CutPassed
                   
               } // result != nullptr
  
            } //Iterate over D_(s)/Lambda_c +/- vertices
  
         } //Iterate over muons
  
         ATH_MSG_DEBUG("cascadeinfoContainer size " << cascadeinfoContainer->size());
  
         return StatusCode::SUCCESS;
     }

◆ 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_cascadeOutputsKeys

std::vector<std::string> DerivationFramework::MuPlusDsCascade::m_cascadeOutputsKeys

private

Definition at line 58 of file MuPlusDsCascade.h.

◆ m_CascadeTools

ToolHandle< DerivationFramework::CascadeTools > DerivationFramework::MuPlusDsCascade::m_CascadeTools

private

Definition at line 81 of file MuPlusDsCascade.h.

◆ m_chi2cut

double DerivationFramework::MuPlusDsCascade::m_chi2cut

private

Definition at line 75 of file MuPlusDsCascade.h.

◆ m_combOnly

bool DerivationFramework::MuPlusDsCascade::m_combOnly

private

Definition at line 98 of file MuPlusDsCascade.h.

◆ m_constrDx

bool DerivationFramework::MuPlusDsCascade::m_constrDx

private

Definition at line 74 of file MuPlusDsCascade.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_DoVertexType

int DerivationFramework::MuPlusDsCascade::m_DoVertexType

private

Definition at line 91 of file MuPlusDsCascade.h.

◆ m_Dx_pid

int DerivationFramework::MuPlusDsCascade::m_Dx_pid

private

Definition at line 73 of file MuPlusDsCascade.h.

◆ m_DxMassLower

double DerivationFramework::MuPlusDsCascade::m_DxMassLower

private

Definition at line 62 of file MuPlusDsCascade.h.

◆ m_DxMassUpper

double DerivationFramework::MuPlusDsCascade::m_DxMassUpper

private

Definition at line 63 of file MuPlusDsCascade.h.

◆ m_eventInfo_key

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::MuPlusDsCascade::m_eventInfo_key {this, "EventInfo", "EventInfo", "Input event information"}

private

Definition at line 77 of file MuPlusDsCascade.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_hypoName

std::string DerivationFramework::MuPlusDsCascade::m_hypoName

private

name of the mass hypothesis.

E.g. Jpis, Upsi, etc. Will be used as a prefix for decorations

Definition at line 87 of file MuPlusDsCascade.h.

◆ m_iVertexFitter

ToolHandle< Trk::TrkVKalVrtFitter > DerivationFramework::MuPlusDsCascade::m_iVertexFitter

private

Definition at line 78 of file MuPlusDsCascade.h.

◆ m_MassLower

double DerivationFramework::MuPlusDsCascade::m_MassLower

private

Definition at line 64 of file MuPlusDsCascade.h.

◆ m_MassUpper

double DerivationFramework::MuPlusDsCascade::m_MassUpper

private

Definition at line 65 of file MuPlusDsCascade.h.

◆ m_mcpCuts

bool DerivationFramework::MuPlusDsCascade::m_mcpCuts

private

Definition at line 97 of file MuPlusDsCascade.h.

◆ m_muonCollectionKey

std::string DerivationFramework::MuPlusDsCascade::m_muonCollectionKey

private

Definition at line 94 of file MuPlusDsCascade.h.

◆ m_partPropSvc

ServiceHandle<IPartPropSvc> DerivationFramework::MuPlusDsCascade::m_partPropSvc {this, "PartPropSvc", "PartPropSvc"}

private

Definition at line 83 of file MuPlusDsCascade.h.

◆ m_PV_max

int DerivationFramework::MuPlusDsCascade::m_PV_max

private

Definition at line 90 of file MuPlusDsCascade.h.

◆ m_PV_minNTracks

size_t DerivationFramework::MuPlusDsCascade::m_PV_minNTracks

private

Definition at line 92 of file MuPlusDsCascade.h.

◆ m_pvRefitter

ToolHandle< Analysis::PrimaryVertexRefitter > DerivationFramework::MuPlusDsCascade::m_pvRefitter

private

Definition at line 79 of file MuPlusDsCascade.h.

◆ m_refitPV

bool DerivationFramework::MuPlusDsCascade::m_refitPV

private

Definition at line 85 of file MuPlusDsCascade.h.

◆ m_refPVContainerName

std::string DerivationFramework::MuPlusDsCascade::m_refPVContainerName

private

Definition at line 86 of file MuPlusDsCascade.h.

◆ m_thresholdPt

double DerivationFramework::MuPlusDsCascade::m_thresholdPt

private

Definition at line 96 of file MuPlusDsCascade.h.

◆ m_trackSelectionTools

std::unique_ptr<InDet::InDetTrackSelectionTool> DerivationFramework::MuPlusDsCascade::m_trackSelectionTools

private

Definition at line 82 of file MuPlusDsCascade.h.

◆ m_trkSelector

ToolHandle< Trk::ITrackSelectorTool > DerivationFramework::MuPlusDsCascade::m_trkSelector

private

Definition at line 95 of file MuPlusDsCascade.h.

◆ m_useCombMeasurement

bool DerivationFramework::MuPlusDsCascade::m_useCombMeasurement

private

Definition at line 99 of file MuPlusDsCascade.h.

◆ m_V0Tools

ToolHandle< Trk::V0Tools > DerivationFramework::MuPlusDsCascade::m_V0Tools

private

Definition at line 80 of file MuPlusDsCascade.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexContainerKey

std::string DerivationFramework::MuPlusDsCascade::m_vertexContainerKey

private

Definition at line 56 of file MuPlusDsCascade.h.

◆ m_vertexDxContainerKey

std::string DerivationFramework::MuPlusDsCascade::m_vertexDxContainerKey

private

Definition at line 57 of file MuPlusDsCascade.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vtx0Daug1MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx0Daug1MassHypo

private

Definition at line 68 of file MuPlusDsCascade.h.

◆ m_vtx0MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx0MassHypo

private

Definition at line 66 of file MuPlusDsCascade.h.

◆ m_vtx1Daug1MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx1Daug1MassHypo

private

Definition at line 69 of file MuPlusDsCascade.h.

◆ m_vtx1Daug2MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx1Daug2MassHypo

private

Definition at line 70 of file MuPlusDsCascade.h.

◆ m_vtx1Daug3MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx1Daug3MassHypo

private

Definition at line 71 of file MuPlusDsCascade.h.

◆ m_vtx1MassHypo

double DerivationFramework::MuPlusDsCascade::m_vtx1MassHypo

private

Definition at line 67 of file MuPlusDsCascade.h.

◆ m_VxPrimaryCandidateName

std::string DerivationFramework::MuPlusDsCascade::m_VxPrimaryCandidateName

private

Name of primary vertex container.

Definition at line 60 of file MuPlusDsCascade.h.

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

Public Member Functions

Static Public Member Functions

Protected Member Functions

Private Types

Private Member Functions

Private Attributes

Detailed Description

Member Typedef Documentation

◆ StoreGateSvc_t

Constructor & Destructor Documentation

◆ MuPlusDsCascade()

◆ ~MuPlusDsCascade()

Member Function Documentation

◆ addBranches()

◆ declareGaudiProperty() [1/4]

◆ declareGaudiProperty() [2/4]

◆ declareGaudiProperty() [3/4]

◆ declareGaudiProperty() [4/4]

◆ declareProperty() [1/6]

◆ declareProperty() [2/6]

◆ declareProperty() [3/6]

◆ declareProperty() [4/6]

◆ declareProperty() [5/6]

◆ declareProperty() [6/6]

◆ detStore()

◆ evtStore() [1/2]

◆ evtStore() [2/2]

◆ extraDeps_update_handler()

◆ FindVertexTrack()

◆ initialize()

◆ inputHandles()

◆ interfaceID()

◆ msg() [1/2]

◆ msg() [2/2]

◆ msgLvl()

◆ outputHandles()

◆ performSearch()

◆ renounce()

◆ renounceArray()

◆ sysInitialize()

◆ sysStart()

◆ updateVHKA()

Member Data Documentation

◆ m_cascadeOutputsKeys

◆ m_CascadeTools

◆ m_chi2cut

◆ m_combOnly

◆ m_constrDx

◆ m_detStore

◆ m_DoVertexType

◆ m_Dx_pid

◆ m_DxMassLower

◆ m_DxMassUpper

◆ m_eventInfo_key

◆ m_evtStore

◆ m_hypoName

◆ m_iVertexFitter

◆ m_MassLower

◆ m_MassUpper

◆ m_mcpCuts

◆ m_muonCollectionKey

◆ m_partPropSvc

◆ m_PV_max

◆ m_PV_minNTracks

◆ m_pvRefitter

◆ m_refitPV

◆ m_refPVContainerName

◆ m_thresholdPt

◆ m_trackSelectionTools

◆ m_trkSelector

◆ m_useCombMeasurement

◆ m_V0Tools

◆ m_varHandleArraysDeclared

◆ m_vertexContainerKey

◆ m_vertexDxContainerKey

◆ m_vhka

◆ m_vtx0Daug1MassHypo

◆ m_vtx0MassHypo

◆ m_vtx1Daug1MassHypo

◆ m_vtx1Daug2MassHypo