#include <PsiPlusPsiCascade.h>

Inheritance diagram for DerivationFramework::PsiPlusPsiCascade:

Collaboration diagram for DerivationFramework::PsiPlusPsiCascade:

Public Member Functions
	PsiPlusPsiCascade (const std::string &type, const std::string &name, const IInterface *parent)

virtual	~PsiPlusPsiCascade ()=default

virtual StatusCode	initialize () override

StatusCode	performSearch (std::vector< Trk::VxCascadeInfo * > cascadeinfoContainer, std::vector< Trk::VxCascadeInfo > *cascadeinfoContainer_noConstr) 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 ()

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
SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexPsi1ContainerKey

SG::ReadHandleKey< xAOD::VertexContainer >	m_vertexPsi2ContainerKey

std::vector< std::string >	m_vertexPsi1HypoNames

std::vector< std::string >	m_vertexPsi2HypoNames

SG::WriteHandleKeyArray< xAOD::VertexContainer >	m_cascadeOutputsKeys

SG::ReadHandleKey< xAOD::VertexContainer >	m_VxPrimaryCandidateName
	Name of primary vertex container. More...

SG::ReadHandleKey< xAOD::TrackParticleContainer >	m_trackContainerName

SG::ReadHandleKey< xAOD::EventInfo >	m_eventInfo_key

double	m_jpsi1MassLower

double	m_jpsi1MassUpper

double	m_jpsi2MassLower

double	m_jpsi2MassUpper

double	m_diTrack1MassLower

double	m_diTrack1MassUpper

double	m_diTrack2MassLower

double	m_diTrack2MassUpper

double	m_psi1MassLower

double	m_psi1MassUpper

double	m_psi2MassLower

double	m_psi2MassUpper

double	m_MassLower

double	m_MassUpper

int	m_vtx1Daug_num

double	m_vtx1Daug1MassHypo

double	m_vtx1Daug2MassHypo

double	m_vtx1Daug3MassHypo

double	m_vtx1Daug4MassHypo

int	m_vtx2Daug_num

double	m_vtx2Daug1MassHypo

double	m_vtx2Daug2MassHypo

double	m_vtx2Daug3MassHypo

double	m_vtx2Daug4MassHypo

double	m_massPsi1

double	m_massPsi2

double	m_massJpsi1

double	m_massJpsi2

double	m_massDiTrk1

double	m_massDiTrk2

bool	m_constrPsi1

bool	m_constrPsi2

bool	m_constrJpsi1

bool	m_constrJpsi2

bool	m_constrDiTrk1

bool	m_constrDiTrk2

double	m_chi2cut_Psi1

double	m_chi2cut_Psi2

double	m_chi2cut

bool	m_removeDuplicatePairs

unsigned int	m_maxCandidates

ToolHandle< Trk::TrkVKalVrtFitter >	m_iVertexFitter

ToolHandle< Analysis::PrimaryVertexRefitter >	m_pvRefitter

ToolHandle< Trk::V0Tools >	m_V0Tools

ToolHandle< DerivationFramework::CascadeTools >	m_CascadeTools

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

bool	m_refitPV

SG::WriteHandleKey< xAOD::VertexContainer >	m_refPVContainerName

std::string	m_hypoName

int	m_PV_max

int	m_DoVertexType

size_t	m_PV_minNTracks

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 33 of file PsiPlusPsiCascade.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

◆ PsiPlusPsiCascade()

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

Definition at line 258 of file PsiPlusPsiCascade.cxx.

                                                                                                              : AthAlgTool(type,name,parent),
     m_vertexPsi1ContainerKey(""),
     m_vertexPsi2ContainerKey(""),
     m_cascadeOutputsKeys({"PsiPlusPsiCascadeVtx1", "PsiPlusPsiCascadeVtx2", "PsiPlusPsiCascadeVtx3"}),
     m_VxPrimaryCandidateName("PrimaryVertices"),
     m_trackContainerName("InDetTrackParticles"),
     m_eventInfo_key("EventInfo"),
     m_jpsi1MassLower(0.0),
     m_jpsi1MassUpper(20000.0),
     m_jpsi2MassLower(0.0),
     m_jpsi2MassUpper(20000.0),
     m_diTrack1MassLower(-1.0),
     m_diTrack1MassUpper(-1.0),
     m_diTrack2MassLower(-1.0),
     m_diTrack2MassUpper(-1.0),
     m_psi1MassLower(0.0),
     m_psi1MassUpper(25000.0),
     m_psi2MassLower(0.0),
     m_psi2MassUpper(25000.0),
     m_MassLower(0.0),
     m_MassUpper(31000.0),
     m_vtx1Daug_num(4),
     m_vtx1Daug1MassHypo(-1),
     m_vtx1Daug2MassHypo(-1),
     m_vtx1Daug3MassHypo(-1),
     m_vtx1Daug4MassHypo(-1),
     m_vtx2Daug_num(4),
     m_vtx2Daug1MassHypo(-1),
     m_vtx2Daug2MassHypo(-1),
     m_vtx2Daug3MassHypo(-1),
     m_vtx2Daug4MassHypo(-1),
     m_massPsi1(-1),
     m_massPsi2(-1),
     m_massJpsi1(-1),
     m_massJpsi2(-1),
     m_massDiTrk1(-1),
     m_massDiTrk2(-1),
     m_constrPsi1(false),
     m_constrPsi2(false),
     m_constrJpsi1(false),
     m_constrJpsi2(false),
     m_constrDiTrk1(false),
     m_constrDiTrk2(false),
     m_chi2cut_Psi1(-1.0),
     m_chi2cut_Psi2(-1.0),
     m_chi2cut(-1.0),
     m_removeDuplicatePairs(false),
     m_maxCandidates(0),
     m_iVertexFitter("Trk::TrkVKalVrtFitter"),
     m_pvRefitter("Analysis::PrimaryVertexRefitter", this),
     m_V0Tools("Trk::V0Tools"),
     m_CascadeTools("DerivationFramework::CascadeTools")
   {
     declareProperty("Psi1Vertices",             m_vertexPsi1ContainerKey);
     declareProperty("Psi2Vertices",             m_vertexPsi2ContainerKey);
     declareProperty("Psi1VtxHypoNames",         m_vertexPsi1HypoNames);
     declareProperty("Psi2VtxHypoNames",         m_vertexPsi2HypoNames);
     declareProperty("VxPrimaryCandidateName",   m_VxPrimaryCandidateName);
     declareProperty("TrackContainerName",       m_trackContainerName);
     declareProperty("RefPVContainerName",       m_refPVContainerName = "RefittedPrimaryVertices");
     declareProperty("Jpsi1MassLowerCut",        m_jpsi1MassLower);
     declareProperty("Jpsi1MassUpperCut",        m_jpsi1MassUpper);
     declareProperty("Jpsi2MassLowerCut",        m_jpsi2MassLower);
     declareProperty("Jpsi2MassUpperCut",        m_jpsi2MassUpper);
     declareProperty("DiTrack1MassLower",        m_diTrack1MassLower); // only effective when m_vtx1Daug_num=4
     declareProperty("DiTrack1MassUpper",        m_diTrack1MassUpper); // only effective when m_vtx1Daug_num=4
     declareProperty("DiTrack2MassLower",        m_diTrack2MassLower); // only effective when m_vtx2Daug_num=4
     declareProperty("DiTrack2MassUpper",        m_diTrack2MassUpper); // only effective when m_vtx2Daug_num=4
     declareProperty("Psi1MassLowerCut",         m_psi1MassLower);
     declareProperty("Psi1MassUpperCut",         m_psi1MassUpper);
     declareProperty("Psi2MassLowerCut",         m_psi2MassLower);
     declareProperty("Psi2MassUpperCut",         m_psi2MassUpper);
     declareProperty("MassLowerCut",             m_MassLower);
     declareProperty("MassUpperCut",             m_MassUpper);
     declareProperty("HypothesisName",           m_hypoName = "TQ");
     declareProperty("NumberOfPsi1Daughters",    m_vtx1Daug_num); // 3 or 4 only
     declareProperty("Vtx1Daug1MassHypo",        m_vtx1Daug1MassHypo);
     declareProperty("Vtx1Daug2MassHypo",        m_vtx1Daug2MassHypo);
     declareProperty("Vtx1Daug3MassHypo",        m_vtx1Daug3MassHypo);
     declareProperty("Vtx1Daug4MassHypo",        m_vtx1Daug4MassHypo);
     declareProperty("NumberOfPsi2Daughters",    m_vtx2Daug_num); // 3 or 4 only
     declareProperty("Vtx2Daug1MassHypo",        m_vtx2Daug1MassHypo);
     declareProperty("Vtx2Daug2MassHypo",        m_vtx2Daug2MassHypo);
     declareProperty("Vtx2Daug3MassHypo",        m_vtx2Daug3MassHypo);
     declareProperty("Vtx2Daug4MassHypo",        m_vtx2Daug4MassHypo);
     declareProperty("Jpsi1Mass",                m_massJpsi1);
     declareProperty("Jpsi2Mass",                m_massJpsi2);
     declareProperty("DiTrack1Mass",             m_massDiTrk1);
     declareProperty("DiTrack2Mass",             m_massDiTrk2);
     declareProperty("Psi1Mass",                 m_massPsi1);
     declareProperty("Psi2Mass",                 m_massPsi2);
     declareProperty("ApplyPsi1MassConstraint",  m_constrPsi1);
     declareProperty("ApplyPsi2MassConstraint",  m_constrPsi2);
     declareProperty("ApplyJpsi1MassConstraint", m_constrJpsi1);
     declareProperty("ApplyJpsi2MassConstraint", m_constrJpsi2);
     declareProperty("ApplyDiTrk1MassConstraint",m_constrDiTrk1); // only effective when m_vtx1Daug_num=4
     declareProperty("ApplyDiTrk2MassConstraint",m_constrDiTrk2); // only effective when m_vtx2Daug_num=4
     declareProperty("Chi2CutPsi1",              m_chi2cut_Psi1);
     declareProperty("Chi2CutPsi2",              m_chi2cut_Psi2);
     declareProperty("Chi2Cut",                  m_chi2cut);
     declareProperty("RemoveDuplicatePairs",     m_removeDuplicatePairs); // only effective when m_vertexPsi1ContainerKey == m_vertexPsi2ContainerKey
     declareProperty("MaxCandidates",            m_maxCandidates);
     declareProperty("RefitPV",                  m_refitPV         = true);
     declareProperty("MaxnPV",                   m_PV_max          = 1000);
     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);
   }

◆ ~PsiPlusPsiCascade()

virtual DerivationFramework::PsiPlusPsiCascade::~PsiPlusPsiCascade ( )

virtualdefault

Member Function Documentation

◆ addBranches()

StatusCode DerivationFramework::PsiPlusPsiCascade::addBranches ( ) const

overridevirtual

Pass the thinning service

Implements DerivationFramework::IAugmentationTool.

Definition at line 66 of file PsiPlusPsiCascade.cxx.

                                                   {
     if ((m_vtx1Daug_num != 3 && m_vtx1Daug_num != 4) || (m_vtx2Daug_num != 3 && m_vtx2Daug_num != 4)) {
       ATH_MSG_FATAL("Incorrect number of Psi daughters (should be 3 or 4)");
       return StatusCode::FAILURE;
     }
  
     constexpr int topoN = 3;
     if(m_cascadeOutputsKeys.size() != topoN) {
       ATH_MSG_FATAL("Incorrect number of VtxContainers");
       return StatusCode::FAILURE;
     }
     std::array<SG::WriteHandle<xAOD::VertexContainer>, topoN> VtxWriteHandles; int ikey(0);
     for(const SG::WriteHandleKey<xAOD::VertexContainer>& key : m_cascadeOutputsKeys) {
       VtxWriteHandles[ikey] = SG::WriteHandle<xAOD::VertexContainer>(key);
       ATH_CHECK( VtxWriteHandles[ikey].record(std::make_unique<xAOD::VertexContainer>(), std::make_unique<xAOD::VertexAuxContainer>()) );
       ikey++;
     }
  
     //----------------------------------------------------
     // retrieve primary vertices
     //----------------------------------------------------
     SG::ReadHandle<xAOD::VertexContainer> pvContainer(m_VxPrimaryCandidateName);
     ATH_CHECK( pvContainer.isValid() );
     if (pvContainer.cptr()->size()==0) {
       ATH_MSG_WARNING("You have no primary vertices: " << pvContainer.cptr()->size());
       return StatusCode::RECOVERABLE;
     }
  
     //----------------------------------------------------
     // Record refitted primary vertices
     //----------------------------------------------------
     SG::WriteHandle<xAOD::VertexContainer> refPvContainer;
     if(m_refitPV) {
       refPvContainer = SG::WriteHandle<xAOD::VertexContainer>(m_refPVContainerName);
       ATH_CHECK( refPvContainer.record(std::make_unique<xAOD::VertexContainer>(), std::make_unique<xAOD::VertexAuxContainer>()) );
     }
  
     std::vector<Trk::VxCascadeInfo*> cascadeinfoContainer;
     std::vector<Trk::VxCascadeInfo*> cascadeinfoContainer_noConstr;
     ATH_CHECK(performSearch(&cascadeinfoContainer,&cascadeinfoContainer_noConstr));
  
     SG::ReadHandle<xAOD::EventInfo> evt(m_eventInfo_key);
     ATH_CHECK( evt.isValid() );
     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<VertexLinkVector> Psi1LinksDecor("Psi1VertexLinks");
     SG::AuxElement::Decorator<VertexLinkVector> Psi2LinksDecor("Psi2VertexLinks");
     SG::AuxElement::Decorator<float> chi2_decor("ChiSquared");
     SG::AuxElement::Decorator<int> ndof_decor("nDoF");
     SG::AuxElement::Decorator<float> chi2_nc_decor("ChiSquared_nc");
     SG::AuxElement::Decorator<int> ndof_nc_decor("nDoF_nc");
     SG::AuxElement::Decorator<float> Pt_decor("Pt");
     SG::AuxElement::Decorator<float> PtErr_decor("PtErr");
     SG::AuxElement::Decorator<float> chi2_SV1_decor("ChiSquared_SV1");
     SG::AuxElement::Decorator<float> chi2_nc_SV1_decor("ChiSquared_nc_SV1");
     SG::AuxElement::Decorator<float> chi2_V1_decor("ChiSquared_V1");
     SG::AuxElement::Decorator<int> ndof_V1_decor("nDoF_V1");
     SG::AuxElement::Decorator<float> lxy_SV1_decor("lxy_SV1");
     SG::AuxElement::Decorator<float> lxyErr_SV1_decor("lxyErr_SV1");
     SG::AuxElement::Decorator<float> a0xy_SV1_decor("a0xy_SV1");
     SG::AuxElement::Decorator<float> a0xyErr_SV1_decor("a0xyErr_SV1");
     SG::AuxElement::Decorator<float> a0z_SV1_decor("a0z_SV1");
     SG::AuxElement::Decorator<float> a0zErr_SV1_decor("a0zErr_SV1");
     SG::AuxElement::Decorator<float> chi2_SV2_decor("ChiSquared_SV2");
     SG::AuxElement::Decorator<float> chi2_nc_SV2_decor("ChiSquared_nc_SV2");
     SG::AuxElement::Decorator<float> chi2_V2_decor("ChiSquared_V2");
     SG::AuxElement::Decorator<int> ndof_V2_decor("nDoF_V2");
     SG::AuxElement::Decorator<float> lxy_SV2_decor("lxy_SV2");
     SG::AuxElement::Decorator<float> lxyErr_SV2_decor("lxyErr_SV2");
     SG::AuxElement::Decorator<float> a0xy_SV2_decor("a0xy_SV2");
     SG::AuxElement::Decorator<float> a0xyErr_SV2_decor("a0xyErr_SV2");
     SG::AuxElement::Decorator<float> a0z_SV2_decor("a0z_SV2");
     SG::AuxElement::Decorator<float> a0zErr_SV2_decor("a0zErr_SV2");
  
     // Get the container and identify the input Psi's
     SG::ReadHandle<xAOD::VertexContainer> psi1Container(m_vertexPsi1ContainerKey);
     ATH_CHECK( psi1Container.isValid() );
     SG::ReadHandle<xAOD::VertexContainer> psi2Container(m_vertexPsi2ContainerKey);
     ATH_CHECK( psi2Container.isValid() );
  
     for(size_t ic=0; ic<cascadeinfoContainer.size(); ic++) {
       Trk::VxCascadeInfo* cascade_info = cascadeinfoContainer[ic];
       if(cascade_info==nullptr) {
         ATH_MSG_ERROR("CascadeInfo is null");
         continue;
       }
  
       Trk::VxCascadeInfo* cascade_info_noConstr = cascadeinfoContainer_noConstr[ic];
  
       const std::vector<xAOD::Vertex*> &cascadeVertices = cascade_info->vertices();
       if(cascadeVertices.size() != topoN) ATH_MSG_ERROR("Incorrect number of vertices");
       if(cascadeVertices[0]==nullptr || cascadeVertices[1]==nullptr || cascadeVertices[2]==nullptr) ATH_MSG_ERROR("Error null vertex");
       // Keep vertices
       for(int i=0; i<topoN; i++) VtxWriteHandles[i].ptr()->push_back(cascadeVertices[i]);
  
       cascade_info->setSVOwnership(false); // Prevent Container from deleting vertices
       const auto mainVertex = cascadeVertices[2]; // this is the mother vertex
       const std::vector< std::vector<TLorentzVector> > &moms = cascade_info->getParticleMoms();
  
       // Set links to cascade vertices
       std::vector<VertexLink> precedingVertexLinks;
       VertexLink vertexLink1;
       vertexLink1.setElement(cascadeVertices[0]);
       vertexLink1.setStorableObject(*VtxWriteHandles[0].ptr());
       if( vertexLink1.isValid() ) precedingVertexLinks.push_back( vertexLink1 );
       VertexLink vertexLink2;
       vertexLink2.setElement(cascadeVertices[1]);
       vertexLink2.setStorableObject(*VtxWriteHandles[1].ptr());
       if( vertexLink2.isValid() ) precedingVertexLinks.push_back( vertexLink2 );
       CascadeLinksDecor(*mainVertex) = precedingVertexLinks;
  
       // Identify the input Psi2
       const xAOD::Vertex* psi2Vertex(0);
       if(m_vtx2Daug_num==4) psi2Vertex = BPhysPVCascadeTools::FindVertex<4>(psi2Container.cptr(), cascadeVertices[1]);
       else psi2Vertex = BPhysPVCascadeTools::FindVertex<3>(psi2Container.cptr(), cascadeVertices[1]);
       // Identify the input Psi1
       const xAOD::Vertex* psi1Vertex(0);
       if(m_vtx1Daug_num==4) psi1Vertex = BPhysPVCascadeTools::FindVertex<4>(psi1Container.cptr(), cascadeVertices[0]);
       else psi1Vertex = BPhysPVCascadeTools::FindVertex<3>(psi1Container.cptr(), cascadeVertices[0]);
  
       // Set links to input vertices
       std::vector<const xAOD::Vertex*> psi2VerticestoLink;
       if(psi2Vertex) psi2VerticestoLink.push_back(psi2Vertex);
       else ATH_MSG_WARNING("Could not find linking Jpsi");
       if(!BPhysPVCascadeTools::LinkVertices(Psi2LinksDecor, psi2VerticestoLink, psi2Container.cptr(), mainVertex)) ATH_MSG_ERROR("Error decorating with Psi2 vertex");
  
       std::vector<const xAOD::Vertex*> psi1VerticestoLink;
       if(psi1Vertex) psi1VerticestoLink.push_back(psi1Vertex);
       else ATH_MSG_WARNING("Could not find linking Psi1");
       if(!BPhysPVCascadeTools::LinkVertices(Psi1LinksDecor, psi1VerticestoLink, psi1Container.cptr(), mainVertex)) ATH_MSG_ERROR("Error decorating with Psi1 vertex");
  
       xAOD::BPhysHypoHelper vtx(m_hypoName, mainVertex);
  
       // Get refitted track momenta from all vertices, charged tracks only
       BPhysPVCascadeTools::SetVectorInfo(vtx, cascade_info);
       vtx.setPass(true);
  
       //
       // Decorate main vertex
       //
       // mass, mass error
       // https://gitlab.cern.ch/atlas/athena/-/blob/21.2/Tracking/TrkVertexFitter/TrkVKalVrtFitter/TrkVKalVrtFitter/VxCascadeInfo.h
       BPHYS_CHECK( vtx.setMass(m_CascadeTools->invariantMass(moms[2])) );
       BPHYS_CHECK( vtx.setMassErr(m_CascadeTools->invariantMassError(moms[2],cascade_info->getCovariance()[1])) );
       // pt and pT error (the default pt of mainVertex is != the pt of the full cascade fit!)
       Pt_decor(*mainVertex)       = m_CascadeTools->pT(moms[2]);
       PtErr_decor(*mainVertex)    = m_CascadeTools->pTError(moms[2],cascade_info->getCovariance()[1]);
       // chi2 and ndof (the default chi2 of mainVertex is != the chi2 of the full cascade fit!)
       chi2_decor(*mainVertex)     = cascade_info->fitChi2();
       ndof_decor(*mainVertex)     = cascade_info->nDoF();
       chi2_nc_decor(*mainVertex)  = cascade_info_noConstr ? cascade_info_noConstr->fitChi2() : -999999.;
       ndof_nc_decor(*mainVertex)  = cascade_info_noConstr ? cascade_info_noConstr->nDoF() : -1;
  
       // decorate the Psi1 vertex
       chi2_SV1_decor(*cascadeVertices[0])    = m_V0Tools->chisq(cascadeVertices[0]);
       chi2_nc_SV1_decor(*cascadeVertices[0]) = cascade_info_noConstr ? m_V0Tools->chisq(cascade_info_noConstr->vertices()[0]) : -999999.;
       chi2_V1_decor(*cascadeVertices[0])     = m_V0Tools->chisq(psi1Vertex);
       ndof_V1_decor(*cascadeVertices[0])     = m_V0Tools->ndof(psi1Vertex);
       lxy_SV1_decor(*cascadeVertices[0])     = m_CascadeTools->lxy(moms[0],cascadeVertices[0],mainVertex);
       lxyErr_SV1_decor(*cascadeVertices[0])  = m_CascadeTools->lxyError(moms[0],cascade_info->getCovariance()[0],cascadeVertices[0],mainVertex);
       a0z_SV1_decor(*cascadeVertices[0])     = m_CascadeTools->a0z(moms[0],cascadeVertices[0],mainVertex);
       a0zErr_SV1_decor(*cascadeVertices[0])  = m_CascadeTools->a0zError(moms[0],cascade_info->getCovariance()[0],cascadeVertices[0],mainVertex);
       a0xy_SV1_decor(*cascadeVertices[0])    = m_CascadeTools->a0xy(moms[0],cascadeVertices[0],mainVertex);
       a0xyErr_SV1_decor(*cascadeVertices[0]) = m_CascadeTools->a0xyError(moms[0],cascade_info->getCovariance()[0],cascadeVertices[0],mainVertex);
  
       // decorate the Psi2 vertex
       chi2_SV2_decor(*cascadeVertices[1])    = m_V0Tools->chisq(cascadeVertices[1]);
       chi2_nc_SV2_decor(*cascadeVertices[1]) = cascade_info_noConstr ? m_V0Tools->chisq(cascade_info_noConstr->vertices()[1]) : -999999.;
       chi2_V2_decor(*cascadeVertices[1])     = m_V0Tools->chisq(psi2Vertex);
       ndof_V2_decor(*cascadeVertices[1])     = m_V0Tools->ndof(psi2Vertex);
       lxy_SV2_decor(*cascadeVertices[1])     = m_CascadeTools->lxy(moms[1],cascadeVertices[1],mainVertex);
       lxyErr_SV2_decor(*cascadeVertices[1])  = m_CascadeTools->lxyError(moms[1],cascade_info->getCovariance()[1],cascadeVertices[1],mainVertex);
       a0z_SV2_decor(*cascadeVertices[1])     = m_CascadeTools->a0z(moms[1],cascadeVertices[1],mainVertex);
       a0zErr_SV2_decor(*cascadeVertices[1])  = m_CascadeTools->a0zError(moms[1],cascade_info->getCovariance()[1],cascadeVertices[1],mainVertex);
       a0xy_SV2_decor(*cascadeVertices[1])    = m_CascadeTools->a0xy(moms[1],cascadeVertices[1],mainVertex);
       a0xyErr_SV2_decor(*cascadeVertices[1]) = m_CascadeTools->a0xyError(moms[1],cascade_info->getCovariance()[1],cascadeVertices[1],mainVertex);
  
       double Mass_Moth = m_CascadeTools->invariantMass(moms[2]); // size=2
       ATH_CHECK(helper.FillCandwithRefittedVertices(m_refitPV, pvContainer.cptr(), m_refitPV ? refPvContainer.ptr() : 0, &(*m_pvRefitter), m_PV_max, m_DoVertexType, cascade_info, 2, Mass_Moth, vtx));
     } // 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 cascade_info : cascadeinfoContainer) delete cascade_info;
     for (auto cascade_info_noConstr : cascadeinfoContainer_noConstr) delete cascade_info_noConstr;
  
     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

◆ initialize()

StatusCode DerivationFramework::PsiPlusPsiCascade::initialize ( )

overridevirtual

Definition at line 23 of file PsiPlusPsiCascade.cxx.

                                            {
     // retrieving vertex Fitter
     ATH_CHECK( m_iVertexFitter.retrieve() );
  
     // retrieve PV refitter
     ATH_CHECK( m_pvRefitter.retrieve() );
  
     // retrieving the V0 tools
     ATH_CHECK( m_V0Tools.retrieve() );
  
     // retrieving the Cascade tools
     ATH_CHECK( m_CascadeTools.retrieve() );
  
     ATH_CHECK( m_vertexPsi1ContainerKey.initialize() );
     ATH_CHECK( m_vertexPsi2ContainerKey.initialize() );
     ATH_CHECK( m_VxPrimaryCandidateName.initialize() );
     ATH_CHECK( m_trackContainerName.initialize() );
     ATH_CHECK( m_refPVContainerName.initialize() );
     ATH_CHECK( m_cascadeOutputsKeys.initialize() );
     ATH_CHECK( m_eventInfo_key.initialize() );
  
     ATH_CHECK( m_partPropSvc.retrieve() );
     auto pdt = m_partPropSvc->PDT();
  
     // retrieve particle masses
     // https://gitlab.cern.ch/atlas/athena/-/blob/main/Generators/TruthUtils/TruthUtils/AtlasPID.h
     if(m_massJpsi1 < 0.) m_massJpsi1 = BPhysPVCascadeTools::getParticleMass(pdt, MC::JPSI);
     if(m_massJpsi2 < 0.) m_massJpsi2 = BPhysPVCascadeTools::getParticleMass(pdt, MC::JPSI);
     if(m_massPsi1 < 0.) m_massPsi1 = BPhysPVCascadeTools::getParticleMass(pdt, MC::PSI2S);
     if(m_massPsi2 < 0.) m_massPsi2 = BPhysPVCascadeTools::getParticleMass(pdt, MC::PSI2S);
  
     if(m_vtx1Daug1MassHypo < 0.) m_vtx1Daug1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::MUON);
     if(m_vtx1Daug2MassHypo < 0.) m_vtx1Daug2MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::MUON);
     if(m_vtx1Daug3MassHypo < 0.) m_vtx1Daug3MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS);
     if(m_vtx1Daug4MassHypo < 0.) m_vtx1Daug4MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS);
     if(m_vtx2Daug1MassHypo < 0.) m_vtx2Daug1MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::MUON);
     if(m_vtx2Daug2MassHypo < 0.) m_vtx2Daug2MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::MUON);
     if(m_vtx2Daug3MassHypo < 0.) m_vtx2Daug3MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS);
     if(m_vtx2Daug4MassHypo < 0.) m_vtx2Daug4MassHypo = BPhysPVCascadeTools::getParticleMass(pdt, MC::PIPLUS);
  
     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::PsiPlusPsiCascade::interfaceID ( )

inlinestatic

Definition at line 36 of file PsiPlusPsiCascade.h.

36 { return IID_PsiPlusPsiCascade;}

◆ 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::PsiPlusPsiCascade::performSearch	(	std::vector< Trk::VxCascadeInfo * > *	cascadeinfoContainer,
		std::vector< Trk::VxCascadeInfo * > *	cascadeinfoContainer_noConstr
	)		const

Definition at line 371 of file PsiPlusPsiCascade.cxx.

                                                                                                                                                                    {
     ATH_MSG_DEBUG( "PsiPlusPsiCascade::performSearch" );
     assert(cascadeinfoContainer!=nullptr && cascadeinfoContainer_noConstr!=nullptr);
  
     // Get TrackParticle container (for setting links to the original tracks)
     SG::ReadHandle<xAOD::TrackParticleContainer> trackContainer(m_trackContainerName);
     ATH_CHECK( trackContainer.isValid() );
  
     std::vector<const xAOD::TrackParticle*> tracksJpsi1;
     std::vector<const xAOD::TrackParticle*> tracksJpsi2;
     std::vector<const xAOD::TrackParticle*> tracksDiTrk1;
     std::vector<const xAOD::TrackParticle*> tracksDiTrk2;
     std::vector<const xAOD::TrackParticle*> tracksPsi1;
     std::vector<const xAOD::TrackParticle*> tracksPsi2;
     std::vector<double> massesPsi1;
     massesPsi1.push_back(m_vtx1Daug1MassHypo);
     massesPsi1.push_back(m_vtx1Daug2MassHypo);
     massesPsi1.push_back(m_vtx1Daug3MassHypo);
     if(m_vtx1Daug_num==4) massesPsi1.push_back(m_vtx1Daug4MassHypo);
     std::vector<double> massesPsi2;
     massesPsi2.push_back(m_vtx2Daug1MassHypo);
     massesPsi2.push_back(m_vtx2Daug2MassHypo);
     massesPsi2.push_back(m_vtx2Daug3MassHypo);
     if(m_vtx2Daug_num==4) massesPsi2.push_back(m_vtx2Daug4MassHypo);
  
     // Get Psi1 container
     SG::ReadHandle<xAOD::VertexContainer> psi1Container(m_vertexPsi1ContainerKey);
     ATH_CHECK( psi1Container.isValid() );
  
     // Get Psi2 container
     SG::ReadHandle<xAOD::VertexContainer> psi2Container(m_vertexPsi2ContainerKey);
     ATH_CHECK( psi2Container.isValid() );
  
     // Select the Psi2 candidates before calling cascade fit
     std::vector<const xAOD::Vertex*> selectedPsi2Candidates;
     for(auto vxcItr=psi2Container.cptr()->cbegin(); vxcItr!=psi2Container.cptr()->cend(); ++vxcItr) {
       // Check the passed flag first
       const xAOD::Vertex* vtx = *vxcItr;
       bool passed = false;
       for(size_t i=0; i<m_vertexPsi2HypoNames.size(); i++) {
     SG::AuxElement::Accessor<Char_t> flagAcc("passed_"+m_vertexPsi2HypoNames[i]);
     if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
       passed |= 1;
     }
       }
       if(m_vertexPsi2HypoNames.size() && !passed) continue;
  
       // Check Psi2 candidate invariant mass and skip if need be
       double mass_psi2 = m_V0Tools->invariantMass(*vxcItr, massesPsi2);
       if (mass_psi2 < m_psi2MassLower || mass_psi2 > m_psi2MassUpper) continue;
  
       // Add loose cut on Jpsi2 mass from Psi2 -> Jpsi2 pi+ pi-, or on phi mass from Ds+ -> phi pi+
       TLorentzVector p4_mu1, p4_mu2;
       p4_mu1.SetPtEtaPhiM( (*vxcItr)->trackParticle(0)->pt(), 
                (*vxcItr)->trackParticle(0)->eta(),
                (*vxcItr)->trackParticle(0)->phi(), m_vtx2Daug1MassHypo);
       p4_mu2.SetPtEtaPhiM( (*vxcItr)->trackParticle(1)->pt(), 
                (*vxcItr)->trackParticle(1)->eta(),
                (*vxcItr)->trackParticle(1)->phi(), m_vtx2Daug2MassHypo);
       double mass_jpsi2 = (p4_mu1 + p4_mu2).M();
       if (mass_jpsi2 < m_jpsi2MassLower || mass_jpsi2 > m_jpsi2MassUpper) continue;
  
       if(m_vtx2Daug_num==4 && m_diTrack2MassLower>=0 && m_diTrack2MassUpper>=0) {
     TLorentzVector p4_trk1, p4_trk2;
     p4_trk1.SetPtEtaPhiM( (*vxcItr)->trackParticle(2)->pt(),
                   (*vxcItr)->trackParticle(2)->eta(),
                   (*vxcItr)->trackParticle(2)->phi(), m_vtx2Daug3MassHypo);
     p4_trk2.SetPtEtaPhiM( (*vxcItr)->trackParticle(3)->pt(),
                   (*vxcItr)->trackParticle(3)->eta(),
                   (*vxcItr)->trackParticle(3)->phi(), m_vtx2Daug4MassHypo);
     double mass_diTrk2 = (p4_trk1 + p4_trk2).M();
     if (mass_diTrk2 < m_diTrack2MassLower || mass_diTrk2 > m_diTrack2MassUpper) continue;
       }
  
       double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
       if(m_chi2cut_Psi2>0 && chi2DOF>m_chi2cut_Psi2) continue;
  
       selectedPsi2Candidates.push_back(*vxcItr);
     }
     if(selectedPsi2Candidates.size()==0) return StatusCode::SUCCESS;
  
     // Select the Psi1 candidates before calling cascade fit
     std::vector<const xAOD::Vertex*> selectedPsi1Candidates;
     for(auto vxcItr=psi1Container.cptr()->cbegin(); vxcItr!=psi1Container.cptr()->cend(); ++vxcItr) {
       // Check the passed flag first
       const xAOD::Vertex* vtx = *vxcItr;
       bool passed = false;
       for(size_t i=0; i<m_vertexPsi1HypoNames.size(); i++) {
     SG::AuxElement::Accessor<Char_t> flagAcc("passed_"+m_vertexPsi1HypoNames[i]);
     if(flagAcc.isAvailable(*vtx) && flagAcc(*vtx)) {
       passed |= 1;
     }
       }
       if(m_vertexPsi1HypoNames.size() && !passed) continue;
  
       // Check Psi candidate invariant mass and skip if need be
       double mass_psi1 = m_V0Tools->invariantMass(*vxcItr,massesPsi1);
       if(mass_psi1 < m_psi1MassLower || mass_psi1 > m_psi1MassUpper) continue;
  
       // Add loose cut on Jpsi1 mass from Psi1 -> Jpsi1 pi+ pi-, or on phi mass from Ds+ -> phi pi+
       TLorentzVector p4_mu1, p4_mu2;
       p4_mu1.SetPtEtaPhiM( (*vxcItr)->trackParticle(0)->pt(), 
                (*vxcItr)->trackParticle(0)->eta(),
                (*vxcItr)->trackParticle(0)->phi(), m_vtx1Daug1MassHypo);
       p4_mu2.SetPtEtaPhiM( (*vxcItr)->trackParticle(1)->pt(), 
                (*vxcItr)->trackParticle(1)->eta(),
                (*vxcItr)->trackParticle(1)->phi(), m_vtx1Daug2MassHypo);
       double mass_jpsi1 = (p4_mu1 + p4_mu2).M();
       if (mass_jpsi1 < m_jpsi1MassLower || mass_jpsi1 > m_jpsi1MassUpper) continue;
  
       if(m_vtx1Daug_num==4 && m_diTrack1MassLower>=0 && m_diTrack1MassUpper>=0) {
     TLorentzVector p4_trk1, p4_trk2;
     p4_trk1.SetPtEtaPhiM( (*vxcItr)->trackParticle(2)->pt(),
                   (*vxcItr)->trackParticle(2)->eta(),
                   (*vxcItr)->trackParticle(2)->phi(), m_vtx1Daug3MassHypo);
     p4_trk2.SetPtEtaPhiM( (*vxcItr)->trackParticle(3)->pt(),
                   (*vxcItr)->trackParticle(3)->eta(),
                   (*vxcItr)->trackParticle(3)->phi(), m_vtx1Daug4MassHypo);
     double mass_diTrk1 = (p4_trk1 + p4_trk2).M();
     if (mass_diTrk1 < m_diTrack1MassLower || mass_diTrk1 > m_diTrack1MassUpper) continue;
       }
  
       double chi2DOF = (*vxcItr)->chiSquared()/(*vxcItr)->numberDoF();
       if(m_chi2cut_Psi1>0 && chi2DOF>m_chi2cut_Psi1) continue;
  
       selectedPsi1Candidates.push_back(*vxcItr);
     }
     if(selectedPsi1Candidates.size()==0) return StatusCode::SUCCESS;
  
     std::vector<std::pair<const xAOD::Vertex*, const xAOD::Vertex*> > candidatePairs;
     for(auto psi1Itr=selectedPsi1Candidates.cbegin(); psi1Itr!=selectedPsi1Candidates.cend(); ++psi1Itr) {
       tracksPsi1.clear();
       for(size_t i=0; i<(*psi1Itr)->nTrackParticles(); i++) tracksPsi1.push_back((*psi1Itr)->trackParticle(i));
       for(auto psi2Itr=selectedPsi2Candidates.cbegin(); psi2Itr!=selectedPsi2Candidates.cend(); ++psi2Itr) {
     bool skip = false;
     for(size_t j=0; j<(*psi2Itr)->nTrackParticles(); j++) {
       if(std::find(tracksPsi1.cbegin(), tracksPsi1.cend(), (*psi2Itr)->trackParticle(j)) != tracksPsi1.cend()) { skip = true; break; }
     }
     if(skip) continue;
     if(m_vertexPsi1ContainerKey.key() == m_vertexPsi2ContainerKey.key()) {
       for(size_t ic=0; ic<candidatePairs.size(); ic++) {
         const xAOD::Vertex* psi1Vertex = candidatePairs[ic].first;
         const xAOD::Vertex* psi2Vertex = candidatePairs[ic].second;
         if((psi1Vertex == *psi1Itr && psi2Vertex == *psi2Itr) || (psi1Vertex == *psi2Itr && psi2Vertex == *psi1Itr)) { skip = true; break; }
       }
     }
     if(skip) continue;
     candidatePairs.push_back(std::pair<const xAOD::Vertex*, const xAOD::Vertex*>(*psi1Itr,*psi2Itr));
       }
     }
  
     std::sort( candidatePairs.begin(), candidatePairs.end(), [](std::pair<const xAOD::Vertex*, const xAOD::Vertex*> a, std::pair<const xAOD::Vertex*, const xAOD::Vertex*> b) { return a.first->chiSquared()/a.first->numberDoF()+a.second->chiSquared()/a.second->numberDoF() < b.first->chiSquared()/b.first->numberDoF()+b.second->chiSquared()/b.second->numberDoF(); } );
     if(m_maxCandidates>0 && candidatePairs.size()>m_maxCandidates) {
       candidatePairs.erase(candidatePairs.begin()+m_maxCandidates, candidatePairs.end());
     }
  
     for(size_t ic=0; ic<candidatePairs.size(); ic++) {
       const xAOD::Vertex* psi1Vertex = candidatePairs[ic].first;
       const xAOD::Vertex* psi2Vertex = candidatePairs[ic].second;
  
       tracksPsi1.clear();
       for(size_t it=0; it<psi1Vertex->nTrackParticles(); it++) tracksPsi1.push_back(psi1Vertex->trackParticle(it));
       if (tracksPsi1.size() != massesPsi1.size()) {
     ATH_MSG_ERROR("Problems with Psi1 input: number of tracks or track mass inputs is not correct!");
       }
       tracksPsi2.clear();
       for(size_t it=0; it<psi2Vertex->nTrackParticles(); it++) tracksPsi2.push_back(psi2Vertex->trackParticle(it));
       if (tracksPsi2.size() != massesPsi2.size()) {
     ATH_MSG_ERROR("Problems with Psi2 input: number of tracks or track mass inputs is not correct!");
       }
  
       tracksJpsi1.clear();
       tracksJpsi1.push_back(psi1Vertex->trackParticle(0));
       tracksJpsi1.push_back(psi1Vertex->trackParticle(1));
       tracksDiTrk1.clear();
       if(m_vtx1Daug_num==4) {
     tracksDiTrk1.push_back(psi1Vertex->trackParticle(2));
     tracksDiTrk1.push_back(psi1Vertex->trackParticle(3));
       }
       tracksJpsi2.clear();
       tracksJpsi2.push_back(psi2Vertex->trackParticle(0));
       tracksJpsi2.push_back(psi2Vertex->trackParticle(1));
       tracksDiTrk2.clear();
       if(m_vtx2Daug_num==4) {
     tracksDiTrk2.push_back(psi2Vertex->trackParticle(2));
     tracksDiTrk2.push_back(psi2Vertex->trackParticle(3));
       }
  
       TLorentzVector p4_moth;
       TLorentzVector tmp;
       for(size_t it=0; it<psi1Vertex->nTrackParticles(); it++) {
     tmp.SetPtEtaPhiM(psi1Vertex->trackParticle(it)->pt(),psi1Vertex->trackParticle(it)->eta(),psi1Vertex->trackParticle(it)->phi(),massesPsi1[it]);
     p4_moth += tmp;
       }
       for(size_t it=0; it<psi2Vertex->nTrackParticles(); it++) {
     tmp.SetPtEtaPhiM(psi2Vertex->trackParticle(it)->pt(),psi2Vertex->trackParticle(it)->eta(),psi2Vertex->trackParticle(it)->phi(),massesPsi2[it]);
     p4_moth += tmp;
       }
       if (p4_moth.M() < m_MassLower || p4_moth.M() > m_MassUpper) continue;
     
       // Apply the user's settings to the fitter
       std::unique_ptr<Trk::IVKalState> state = m_iVertexFitter->makeState();
       // Robustness: http://cdsweb.cern.ch/record/685551
       int robustness = 0;
       m_iVertexFitter->setRobustness(robustness, *state);
       // Build up the topology
       // Vertex list
       std::vector<Trk::VertexID> vrtList;
       // Psi1 vertex
       Trk::VertexID vID1;
       // https://gitlab.cern.ch/atlas/athena/-/blob/21.2/Tracking/TrkVertexFitter/TrkVKalVrtFitter/TrkVKalVrtFitter/IVertexCascadeFitter.h
       if (m_constrPsi1) {
     vID1 = m_iVertexFitter->startVertex(tracksPsi1,massesPsi1,*state,m_massPsi1);
       } else {
     vID1 = m_iVertexFitter->startVertex(tracksPsi1,massesPsi1,*state);
       }
       vrtList.push_back(vID1);
       // Psi2 vertex
       Trk::VertexID vID2;
       if (m_constrPsi2) {
     vID2 = m_iVertexFitter->nextVertex(tracksPsi2,massesPsi2,*state,m_massPsi2);
       } else {
     vID2 = m_iVertexFitter->nextVertex(tracksPsi2,massesPsi2,*state);
       }
       vrtList.push_back(vID2);
       // Mother vertex including Psi1 and Psi2
       std::vector<const xAOD::TrackParticle*> tp; tp.clear();
       std::vector<double> tp_masses; tp_masses.clear();
       m_iVertexFitter->nextVertex(tp,tp_masses,vrtList,*state);
       if (m_constrJpsi1) {
     std::vector<Trk::VertexID> cnstV; cnstV.clear();
     if ( !m_iVertexFitter->addMassConstraint(vID1,tracksJpsi1,cnstV,*state,m_massJpsi1).isSuccess() ) {
       ATH_MSG_WARNING("addMassConstraint for Jpsi1 failed");
     }
       }
       if (m_constrDiTrk1 && m_vtx1Daug_num==4 && m_massDiTrk1>0) {
     std::vector<Trk::VertexID> cnstV; cnstV.clear();
     if ( !m_iVertexFitter->addMassConstraint(vID1,tracksDiTrk1,cnstV,*state,m_massDiTrk1).isSuccess() ) {
       ATH_MSG_WARNING("addMassConstraint for DiTrk1 failed");
     }
       }
       if (m_constrJpsi2) {
     std::vector<Trk::VertexID> cnstV; cnstV.clear();
     if ( !m_iVertexFitter->addMassConstraint(vID2,tracksJpsi2,cnstV,*state,m_massJpsi2).isSuccess() ) {
       ATH_MSG_WARNING("addMassConstraint for Jpsi2 failed");
     }
       }
       if (m_constrDiTrk2 && m_vtx2Daug_num==4 && m_massDiTrk2>0) {
     std::vector<Trk::VertexID> cnstV; cnstV.clear();
     if ( !m_iVertexFitter->addMassConstraint(vID2,tracksDiTrk2,cnstV,*state,m_massDiTrk2).isSuccess() ) {
       ATH_MSG_WARNING("addMassConstraint for DiTrk2 failed");
     }
       }
       // Do the work
       std::unique_ptr<Trk::VxCascadeInfo> result(m_iVertexFitter->fitCascade(*state));
  
       bool pass = false;
       if (result != nullptr) {
     for(auto v : result->vertices()) {
       if(v->nTrackParticles()==0) {
         std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
         v->setTrackParticleLinks(nullLinkVector);
       }
     }
     // reset links to original tracks
     BPhysPVCascadeTools::PrepareVertexLinks(result.get(), trackContainer.cptr());
  
     // necessary to prevent memory leak
     result->setSVOwnership(true);
  
     // Chi2/DOF cut
     double chi2DOF = result->fitChi2()/result->nDoF();
     bool chi2CutPassed = (m_chi2cut <= 0.0 || chi2DOF < m_chi2cut);
  
     if(chi2CutPassed) {
       cascadeinfoContainer->push_back(result.release());
       pass = true;
     }
       }
  
       // do cascade fit again without any mass constraints
       if(pass) {
     if(m_constrJpsi1 || m_constrPsi1 || m_constrJpsi2 || m_constrPsi2 || (m_constrDiTrk1 && m_vtx1Daug_num==4 && m_massDiTrk1>0) || (m_constrDiTrk2 && m_vtx2Daug_num==4 && m_massDiTrk2>0)) {
       std::unique_ptr<Trk::IVKalState> state (m_iVertexFitter->makeState());
       m_iVertexFitter->setRobustness(robustness, *state);
       std::vector<Trk::VertexID> vrtList_nc;
       // Psi1 vertex
       Trk::VertexID vID1_nc = m_iVertexFitter->startVertex(tracksPsi1,massesPsi1,*state);
       vrtList_nc.push_back(vID1_nc);
       // Psi2 vertex
       Trk::VertexID vID2_nc = m_iVertexFitter->nextVertex(tracksPsi2,massesPsi2,*state);
       vrtList_nc.push_back(vID2_nc);
       // Mother vertex including Psi1 and Psi2
       std::vector<const xAOD::TrackParticle*> tp; tp.clear();
       std::vector<double> tp_masses; tp_masses.clear();
       m_iVertexFitter->nextVertex(tp,tp_masses,vrtList_nc,*state);
       // Do the work
       std::unique_ptr<Trk::VxCascadeInfo> result_nc(m_iVertexFitter->fitCascade(*state));
  
       if (result_nc != nullptr) {
         for(auto v : result_nc->vertices()) {
           if(v->nTrackParticles()==0) {
         std::vector<ElementLink<xAOD::TrackParticleContainer> > nullLinkVector;
         v->setTrackParticleLinks(nullLinkVector);
           }
         }
         // reset links to original tracks
         BPhysPVCascadeTools::PrepareVertexLinks(result_nc.get(), trackContainer.cptr());
  
         // necessary to prevent memory leak
         result_nc->setSVOwnership(true);
         cascadeinfoContainer_noConstr->push_back(result_nc.release());
       }
       else cascadeinfoContainer_noConstr->push_back(0);
     }
     else cascadeinfoContainer_noConstr->push_back(0);
       }
     } //Iterate over candidatePairs
  
     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

SG::WriteHandleKeyArray<xAOD::VertexContainer> DerivationFramework::PsiPlusPsiCascade::m_cascadeOutputsKeys

private

Definition at line 48 of file PsiPlusPsiCascade.h.

◆ m_CascadeTools

ToolHandle< DerivationFramework::CascadeTools > DerivationFramework::PsiPlusPsiCascade::m_CascadeTools

private

Definition at line 99 of file PsiPlusPsiCascade.h.

◆ m_chi2cut

double DerivationFramework::PsiPlusPsiCascade::m_chi2cut

private

Definition at line 92 of file PsiPlusPsiCascade.h.

◆ m_chi2cut_Psi1

double DerivationFramework::PsiPlusPsiCascade::m_chi2cut_Psi1

private

Definition at line 90 of file PsiPlusPsiCascade.h.

◆ m_chi2cut_Psi2

double DerivationFramework::PsiPlusPsiCascade::m_chi2cut_Psi2

private

Definition at line 91 of file PsiPlusPsiCascade.h.

◆ m_constrDiTrk1

bool DerivationFramework::PsiPlusPsiCascade::m_constrDiTrk1

private

Definition at line 88 of file PsiPlusPsiCascade.h.

◆ m_constrDiTrk2

bool DerivationFramework::PsiPlusPsiCascade::m_constrDiTrk2

private

Definition at line 89 of file PsiPlusPsiCascade.h.

◆ m_constrJpsi1

bool DerivationFramework::PsiPlusPsiCascade::m_constrJpsi1

private

Definition at line 86 of file PsiPlusPsiCascade.h.

◆ m_constrJpsi2

bool DerivationFramework::PsiPlusPsiCascade::m_constrJpsi2

private

Definition at line 87 of file PsiPlusPsiCascade.h.

◆ m_constrPsi1

bool DerivationFramework::PsiPlusPsiCascade::m_constrPsi1

private

Definition at line 84 of file PsiPlusPsiCascade.h.

◆ m_constrPsi2

bool DerivationFramework::PsiPlusPsiCascade::m_constrPsi2

private

Definition at line 85 of file PsiPlusPsiCascade.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_diTrack1MassLower

double DerivationFramework::PsiPlusPsiCascade::m_diTrack1MassLower

private

Definition at line 57 of file PsiPlusPsiCascade.h.

◆ m_diTrack1MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_diTrack1MassUpper

private

Definition at line 58 of file PsiPlusPsiCascade.h.

◆ m_diTrack2MassLower

double DerivationFramework::PsiPlusPsiCascade::m_diTrack2MassLower

private

Definition at line 59 of file PsiPlusPsiCascade.h.

◆ m_diTrack2MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_diTrack2MassUpper

private

Definition at line 60 of file PsiPlusPsiCascade.h.

◆ m_DoVertexType

int DerivationFramework::PsiPlusPsiCascade::m_DoVertexType

private

Definition at line 106 of file PsiPlusPsiCascade.h.

◆ m_eventInfo_key

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::PsiPlusPsiCascade::m_eventInfo_key

private

Definition at line 51 of file PsiPlusPsiCascade.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::PsiPlusPsiCascade::m_hypoName

private

Definition at line 104 of file PsiPlusPsiCascade.h.

◆ m_iVertexFitter

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

private

Definition at line 96 of file PsiPlusPsiCascade.h.

◆ m_jpsi1MassLower

double DerivationFramework::PsiPlusPsiCascade::m_jpsi1MassLower

private

Definition at line 53 of file PsiPlusPsiCascade.h.

◆ m_jpsi1MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_jpsi1MassUpper

private

Definition at line 54 of file PsiPlusPsiCascade.h.

◆ m_jpsi2MassLower

double DerivationFramework::PsiPlusPsiCascade::m_jpsi2MassLower

private

Definition at line 55 of file PsiPlusPsiCascade.h.

◆ m_jpsi2MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_jpsi2MassUpper

private

Definition at line 56 of file PsiPlusPsiCascade.h.

◆ m_massDiTrk1

double DerivationFramework::PsiPlusPsiCascade::m_massDiTrk1

private

Definition at line 82 of file PsiPlusPsiCascade.h.

◆ m_massDiTrk2

double DerivationFramework::PsiPlusPsiCascade::m_massDiTrk2

private

Definition at line 83 of file PsiPlusPsiCascade.h.

◆ m_massJpsi1

double DerivationFramework::PsiPlusPsiCascade::m_massJpsi1

private

Definition at line 80 of file PsiPlusPsiCascade.h.

◆ m_massJpsi2

double DerivationFramework::PsiPlusPsiCascade::m_massJpsi2

private

Definition at line 81 of file PsiPlusPsiCascade.h.

◆ m_MassLower

double DerivationFramework::PsiPlusPsiCascade::m_MassLower

private

Definition at line 65 of file PsiPlusPsiCascade.h.

◆ m_massPsi1

double DerivationFramework::PsiPlusPsiCascade::m_massPsi1

private

Definition at line 78 of file PsiPlusPsiCascade.h.

◆ m_massPsi2

double DerivationFramework::PsiPlusPsiCascade::m_massPsi2

private

Definition at line 79 of file PsiPlusPsiCascade.h.

◆ m_MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_MassUpper

private

Definition at line 66 of file PsiPlusPsiCascade.h.

◆ m_maxCandidates

unsigned int DerivationFramework::PsiPlusPsiCascade::m_maxCandidates

private

Definition at line 94 of file PsiPlusPsiCascade.h.

◆ m_partPropSvc

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

private

Definition at line 100 of file PsiPlusPsiCascade.h.

◆ m_psi1MassLower

double DerivationFramework::PsiPlusPsiCascade::m_psi1MassLower

private

Definition at line 61 of file PsiPlusPsiCascade.h.

◆ m_psi1MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_psi1MassUpper

private

Definition at line 62 of file PsiPlusPsiCascade.h.

◆ m_psi2MassLower

double DerivationFramework::PsiPlusPsiCascade::m_psi2MassLower

private

Definition at line 63 of file PsiPlusPsiCascade.h.

◆ m_psi2MassUpper

double DerivationFramework::PsiPlusPsiCascade::m_psi2MassUpper

private

Definition at line 64 of file PsiPlusPsiCascade.h.

◆ m_PV_max

int DerivationFramework::PsiPlusPsiCascade::m_PV_max

private

Definition at line 105 of file PsiPlusPsiCascade.h.

◆ m_PV_minNTracks

size_t DerivationFramework::PsiPlusPsiCascade::m_PV_minNTracks

private

Definition at line 107 of file PsiPlusPsiCascade.h.

◆ m_pvRefitter

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

private

Definition at line 97 of file PsiPlusPsiCascade.h.

◆ m_refitPV

bool DerivationFramework::PsiPlusPsiCascade::m_refitPV

private

Definition at line 102 of file PsiPlusPsiCascade.h.

◆ m_refPVContainerName

SG::WriteHandleKey<xAOD::VertexContainer> DerivationFramework::PsiPlusPsiCascade::m_refPVContainerName

private

Definition at line 103 of file PsiPlusPsiCascade.h.

◆ m_removeDuplicatePairs

bool DerivationFramework::PsiPlusPsiCascade::m_removeDuplicatePairs

private

Definition at line 93 of file PsiPlusPsiCascade.h.

◆ m_trackContainerName

SG::ReadHandleKey<xAOD::TrackParticleContainer> DerivationFramework::PsiPlusPsiCascade::m_trackContainerName

private

Definition at line 50 of file PsiPlusPsiCascade.h.

◆ m_V0Tools

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

private

Definition at line 98 of file PsiPlusPsiCascade.h.

◆ m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< AlgTool > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

◆ m_vertexPsi1ContainerKey

SG::ReadHandleKey<xAOD::VertexContainer> DerivationFramework::PsiPlusPsiCascade::m_vertexPsi1ContainerKey

private

Definition at line 44 of file PsiPlusPsiCascade.h.

◆ m_vertexPsi1HypoNames

std::vector<std::string> DerivationFramework::PsiPlusPsiCascade::m_vertexPsi1HypoNames

private

Definition at line 46 of file PsiPlusPsiCascade.h.

◆ m_vertexPsi2ContainerKey

SG::ReadHandleKey<xAOD::VertexContainer> DerivationFramework::PsiPlusPsiCascade::m_vertexPsi2ContainerKey

private

Definition at line 45 of file PsiPlusPsiCascade.h.

◆ m_vertexPsi2HypoNames

std::vector<std::string> DerivationFramework::PsiPlusPsiCascade::m_vertexPsi2HypoNames

private

Definition at line 47 of file PsiPlusPsiCascade.h.

◆ m_vhka

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

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

◆ m_vtx1Daug1MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx1Daug1MassHypo

private

Definition at line 68 of file PsiPlusPsiCascade.h.

◆ m_vtx1Daug2MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx1Daug2MassHypo

private

Definition at line 69 of file PsiPlusPsiCascade.h.

◆ m_vtx1Daug3MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx1Daug3MassHypo

private

Definition at line 70 of file PsiPlusPsiCascade.h.

◆ m_vtx1Daug4MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx1Daug4MassHypo

private

Definition at line 71 of file PsiPlusPsiCascade.h.

◆ m_vtx1Daug_num

int DerivationFramework::PsiPlusPsiCascade::m_vtx1Daug_num

private

Definition at line 67 of file PsiPlusPsiCascade.h.

◆ m_vtx2Daug1MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx2Daug1MassHypo

private

Definition at line 73 of file PsiPlusPsiCascade.h.

◆ m_vtx2Daug2MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx2Daug2MassHypo

private

Definition at line 74 of file PsiPlusPsiCascade.h.

◆ m_vtx2Daug3MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx2Daug3MassHypo

private

Definition at line 75 of file PsiPlusPsiCascade.h.

◆ m_vtx2Daug4MassHypo

double DerivationFramework::PsiPlusPsiCascade::m_vtx2Daug4MassHypo

private

Definition at line 76 of file PsiPlusPsiCascade.h.

◆ m_vtx2Daug_num

int DerivationFramework::PsiPlusPsiCascade::m_vtx2Daug_num

private

Definition at line 72 of file PsiPlusPsiCascade.h.

◆ m_VxPrimaryCandidateName

SG::ReadHandleKey<xAOD::VertexContainer> DerivationFramework::PsiPlusPsiCascade::m_VxPrimaryCandidateName

private

Name of primary vertex container.

Definition at line 49 of file PsiPlusPsiCascade.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

◆ PsiPlusPsiCascade()

◆ ~PsiPlusPsiCascade()

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()

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

◆ m_chi2cut_Psi2

◆ m_constrDiTrk1

◆ m_constrDiTrk2

◆ m_constrJpsi1

◆ m_constrJpsi2

◆ m_constrPsi1

◆ m_constrPsi2

◆ m_detStore

◆ m_diTrack1MassLower

◆ m_diTrack1MassUpper

◆ m_diTrack2MassLower

◆ m_diTrack2MassUpper

◆ m_DoVertexType

◆ m_eventInfo_key

◆ m_evtStore

◆ m_hypoName

◆ m_iVertexFitter

◆ m_jpsi1MassLower

◆ m_jpsi1MassUpper

◆ m_jpsi2MassLower

◆ m_jpsi2MassUpper

◆ m_massDiTrk1

◆ m_massDiTrk2

◆ m_massJpsi1

◆ m_massJpsi2

◆ m_MassLower

◆ m_massPsi1

◆ m_massPsi2

◆ m_MassUpper

◆ m_maxCandidates

◆ m_partPropSvc

◆ m_psi1MassLower

◆ m_psi1MassUpper

◆ m_psi2MassLower