DerivationFramework::FourLeptonVertexingAlgorithm Class Reference

#include <FourLeptonVertexingAlgorithm.h>

Inheritance diagram for DerivationFramework::FourLeptonVertexingAlgorithm:

Collaboration diagram for DerivationFramework::FourLeptonVertexingAlgorithm:

Public Types
using	LeptonQuadruplet = std::array<const xAOD::IParticle*, 4>

Public Member Functions
	FourLeptonVertexingAlgorithm (const std::string &n, ISvcLocator *p)
StatusCode	initialize () override
StatusCode	execute (const EventContext &ctx) const override
StatusCode	finalize () override
virtual StatusCode	sysInitialize () override
	Override sysInitialize.
virtual bool	isClonable () const override
	Specify if the algorithm is clonable.
virtual unsigned int	cardinality () const override
	Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.
virtual StatusCode	sysExecute (const EventContext &ctx) override
	Execute an algorithm.
virtual const DataObjIDColl &	extraOutputDeps () const override
	Return the list of extra output dependencies.
virtual bool	filterPassed (const EventContext &ctx) const
virtual void	setFilterPassed (bool state, const EventContext &ctx) const
ServiceHandle< StoreGateSvc > &	evtStore ()
	The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
const ServiceHandle< StoreGateSvc > &	detStore () const
	The standard StoreGateSvc/DetectorStore Returns (kind of) a pointer to the StoreGateSvc.
virtual StatusCode	sysStart () override
	Handle START transition.
virtual std::vector< Gaudi::DataHandle * >	inputHandles () const override
	Return this algorithm's input handles.
virtual std::vector< Gaudi::DataHandle * >	outputHandles () const override
	Return this algorithm's output handles.
Gaudi::Details::PropertyBase &	declareProperty (Gaudi::Property< T, V, H > &t)
void	updateVHKA (Gaudi::Details::PropertyBase &)
MsgStream &	msg () const
bool	msgLvl (const MSG::Level lvl) const

Protected Member Functions
void	renounceArray (SG::VarHandleKeyArray &handlesArray)
	remove all handles from I/O resolution
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.

Private Types
using	MuonTrk = xAOD::Muon::TrackParticleType
typedef ServiceHandle< StoreGateSvc >	StoreGateSvc_t

Private Member Functions
bool	passSelection (const xAOD::Electron *elec) const
bool	passSelection (const xAOD::Muon *muon) const
bool	passSelection (const LeptonQuadruplet &quad) const
std::vector< LeptonQuadruplet >	buildAllQuadruplets (const EventContext &ctx) const
std::unique_ptr< xAOD::Vertex >	fitQuadruplet (const EventContext &ctx, const LeptonQuadruplet &quad) const
const xAOD::TrackParticle *	trackParticle (const xAOD::IParticle *part) const
int	charge (const xAOD::IParticle *part) const
Gaudi::Details::PropertyBase &	declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
	specialization for handling Gaudi::Property<SG::VarHandleKey>

Private Attributes
SG::ReadHandleKey< xAOD::MuonContainer >	m_muonKey
	Input containers.
SG::ReadHandleKey< xAOD::ElectronContainer >	m_elecKey
SG::ReadHandleKey< xAOD::EventInfo >	m_evtKey {this, "EventInfoKey", "EventInfo"}
ToolHandle< CP::IMuonSelectionTool >	m_muonSelTool
ToolHandle< IAsgElectronLikelihoodTool >	m_elecSelTool
ToolHandle< Trk::IVertexFitter >	m_fitter {this, "VertexFitter", "", "Vertex fitter tool"}
Gaudi::Property< float >	m_minMuonPt
Gaudi::Property< float >	m_minElecPt {this, "MinElecPt", 4.5 * Gaudi::Units::GeV, " Minimum pt cut applied on the electron"}
Gaudi::Property< int >	m_muonTrkProp
MuonTrk	m_muonTrk {MuonTrk::Primary}
Gaudi::Property< bool >	m_elecUseGSF
Gaudi::Property< float >	m_z0Cut
Gaudi::Property< float >	m_lowSFOS_Cut {this, "LowSFOSMass", 3.5 * Gaudi::Units::GeV, "Minimal mass for the lower SFOS pair"}
Gaudi::Property< float >	m_highPair_Cut {this, "HighMassPair", 40.*Gaudi::Units::GeV, "Minimum mass for the largest lepton pair"}
Gaudi::Property< float >	m_LeadPtCut {this, "LeadingLeptonPt", 15. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "}
Gaudi::Property< float >	m_SubLeadPtCut {this, "SubLeadingLeptonPt", 10. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "}
Gaudi::Property< float >	m_VtxChi2Cut {this, "VertexChi2Cut", 20, "Maximal chii n.D.o.F cut on the reconstructed vertex" }
SG::WriteHandleKey< xAOD::VertexContainer >	m_vtxKey {this, "OutVertexContainer", "FourLeptonVertices", "Output vertex container"}
Gaudi::Property< bool >	m_pruneCov {this, "PruneCovariance", true, "Clears the covariance vector"}
Gaudi::Property< bool >	m_pruneWeight {this, "PruneTrackWeights", true, "Clears the track weight vector"}
std::array< std::atomic< unsigned int >, 8 > m_num_lep	ATLAS_THREAD_SAFE {}
	Simple counter to evaluate the number of leptons per event.
std::array< std::atomic< unsigned int >, 8 > m_num_ele	ATLAS_THREAD_SAFE {}
std::array< std::atomic< unsigned int >, 8 > m_num_muo	ATLAS_THREAD_SAFE {}
std::atomic< unsigned int >	m_tried_fits {0}
	How many vertex fits were attempted.
std::atomic< unsigned int >	m_good_fits {0}
	How many vertex fits were actually successful.
DataObjIDColl	m_extendedExtraObjects
	Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.
StoreGateSvc_t	m_evtStore
	Pointer to StoreGate (event store by default)
StoreGateSvc_t	m_detStore
	Pointer to StoreGate (detector store by default)
std::vector< SG::VarHandleKeyArray * >	m_vhka
bool	m_varHandleArraysDeclared

Detailed Description

Definition at line 26 of file FourLeptonVertexingAlgorithm.h.

Member Typedef Documentation

LeptonQuadruplet

using DerivationFramework::FourLeptonVertexingAlgorithm::LeptonQuadruplet = std::array<const xAOD::IParticle*, 4>

Definition at line 33 of file FourLeptonVertexingAlgorithm.h.

MuonTrk

using DerivationFramework::FourLeptonVertexingAlgorithm::MuonTrk = xAOD::Muon::TrackParticleType

private

Definition at line 65 of file FourLeptonVertexingAlgorithm.h.

StoreGateSvc_t

typedef ServiceHandle<StoreGateSvc> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::StoreGateSvc_t

privateinherited

Definition at line 388 of file AthCommonDataStore.h.

Constructor & Destructor Documentation

FourLeptonVertexingAlgorithm()

DerivationFramework::FourLeptonVertexingAlgorithm::FourLeptonVertexingAlgorithm	(	const std::string &	n,
		ISvcLocator *	p )

Definition at line 43 of file FourLeptonVertexingAlgorithm.cxx.

: AthReentrantAlgorithm(n, p) {}

Member Function Documentation

buildAllQuadruplets()

std::vector< LeptonQuadruplet > DerivationFramework::FourLeptonVertexingAlgorithm::buildAllQuadruplets ( const EventContext & ctx ) const

private

Retrieve the container from the store gate

Dump the leptons in a common vector to build all possible quadruplets

Increment the lepton counter

None of the leptons has enough momentum to be selected in the offline analysis

The leading lepton has too little momentum. All subsequent combinations will have even less. Bail out of the loop

Definition at line 95 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                                               {
        std::vector<LeptonQuadruplet> to_ret{};

        SG::ReadHandle<xAOD::MuonContainer> muonContainer{m_muonKey, ctx};
        if (!muonContainer.isValid()) {
            ATH_MSG_FATAL("Failed to retrieve the muon container for input " << m_muonKey.fullKey());
            throw std::runtime_error("Invalid key access");
        }
        SG::ReadHandle<xAOD::ElectronContainer> elecContainer{m_elecKey, ctx};
        if (!elecContainer.isValid()) {
            ATH_MSG_FATAL("Failed to retrieve the electron container for input " << m_elecKey.fullKey());
            throw std::runtime_error("Invalid key access");
        }
        using PartVec = std::vector<const xAOD::IParticle*>;
        PartVec selected_lep{};
        const size_t num_lep = elecContainer->size() + muonContainer->size();
        if (num_lep < 4) {
            ATH_MSG_DEBUG("Less than four objects reconstructed");
            ++m_num_lep[num_lep];
            ++m_num_ele[elecContainer->size()];
            ++m_num_muo[muonContainer->size()];            
            return to_ret;
        }
        selected_lep.reserve(num_lep);
        size_t n_ele{0},n_muo{0};
        for (const xAOD::Muon* muon : *muonContainer) {
            if (passSelection(muon)) {
                ATH_MSG_DEBUG("Add " << muon);
                selected_lep.push_back(muon);
                ++n_ele;
            }
        }
        for (const xAOD::Electron* elec : *elecContainer) {
            if (passSelection(elec)) {
                ATH_MSG_DEBUG("Add " << elec);
                selected_lep.push_back(elec);
                ++n_muo;
            }
        }
        ++m_num_lep[std::min(m_num_lep.size() -1 , n_ele + n_muo)];
        ++m_num_ele[std::min(m_num_ele.size() -1 , n_ele)];
        ++m_num_muo[std::min(m_num_muo.size() -1 , n_muo)];
        
        if (selected_lep.size() < 4) {
            ATH_MSG_DEBUG("Less than four leptons survived.");
            return to_ret;
        }
        std::sort(selected_lep.begin(),selected_lep.end(),
                        [](const xAOD::IParticle* a, const xAOD::IParticle* b) {                            
                            return a->pt() > b->pt();});
        if (selected_lep[0]->pt() < m_LeadPtCut || selected_lep[1]->pt() < m_SubLeadPtCut) return to_ret;
 
        to_ret.reserve(std::pow(selected_lep.size(), 4));
        for (PartVec::const_iterator itr1 = selected_lep.begin() + 3; itr1 != selected_lep.end(); ++itr1) {
            for (PartVec::const_iterator itr2 = selected_lep.begin() + 2; itr2 != itr1; ++itr2) {
                for (PartVec::const_iterator itr3 = selected_lep.begin() + 1; itr3 != itr2; ++itr3) {
                    if ( (*itr3)->pt() < m_SubLeadPtCut) return to_ret;
                    for (PartVec::const_iterator itr4 = selected_lep.begin(); itr4 != itr3; ++itr4) {
                        LeptonQuadruplet quad{*itr4, *itr3, *itr2, *itr1};
                        if (quad[0]->pt() < m_LeadPtCut) return to_ret;
                        ATH_MSG_DEBUG("Create new quaduplet " << quad);
                        to_ret.push_back(std::move(quad));
                    }
                }
            }
        }
        return to_ret;
    }

cardinality()

unsigned int AthCommonReentrantAlgorithm< Gaudi::Algorithm >::cardinality ( ) const

overridevirtualinherited

Cardinality (Maximum number of clones that can exist) special value 0 means that algorithm is reentrant.

Override this to return 0 for reentrant algorithms.

Definition at line 75 of file AthCommonReentrantAlgorithm.cxx.

{
  return 0;
}

charge()

int DerivationFramework::FourLeptonVertexingAlgorithm::charge ( const xAOD::IParticle * part ) const

private

Definition at line 176 of file FourLeptonVertexingAlgorithm.cxx.

                                                                           {
        static const SG::AuxElement::ConstAccessor<float> acc_charge{"charge"};
        if (acc_charge.isAvailable(*part)) return acc_charge(*part);
        return trackParticle(part)->charge();
    }

declareGaudiProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareGaudiProperty ( Gaudi::Property< T, V, H > & 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());
 
  }

declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & 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< Gaudi::Algorithm > >::detStore ( ) const

inlineinherited

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

Definition at line 95 of file AthCommonDataStore.h.

{ return m_detStore; }

evtStore()

ServiceHandle< StoreGateSvc > & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::evtStore ( )

inlineinherited

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

Definition at line 85 of file AthCommonDataStore.h.

{ return m_evtStore; }

execute()

StatusCode DerivationFramework::FourLeptonVertexingAlgorithm::execute ( const EventContext & ctx ) const

override

Setup the output container

Definition at line 74 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                  {
        SG::WriteHandle<xAOD::VertexContainer> vtxContainer{m_vtxKey, ctx};
        ATH_CHECK(vtxContainer.record(std::make_unique<xAOD::VertexContainer>(), std::make_unique<xAOD::VertexAuxContainer>()));
        xAOD::VertexContainer* out_container = vtxContainer.ptr();
 
        std::vector<LeptonQuadruplet> all_quads = buildAllQuadruplets(ctx);
        for (const LeptonQuadruplet& quad : all_quads) {
            std::unique_ptr<xAOD::Vertex> candidate = fitQuadruplet(ctx, quad);
            if (candidate) out_container->push_back(std::move(candidate));
        }
        return StatusCode::SUCCESS;
    }

extraDeps_update_handler()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

extraOutputDeps()

const DataObjIDColl & AthCommonReentrantAlgorithm< Gaudi::Algorithm >::extraOutputDeps ( ) const

overridevirtualinherited

Return the list of extra output dependencies.

This list is extended to include symlinks implied by inheritance relations.

Definition at line 94 of file AthCommonReentrantAlgorithm.cxx.

{
  // If we didn't find any symlinks to add, just return the collection
  // from the base class.  Otherwise, return the extended collection.
  if (!m_extendedExtraObjects.empty()) {
    return m_extendedExtraObjects;
  }
  return BaseAlg::extraOutputDeps();
}

filterPassed()

virtual bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::filterPassed ( const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 96 of file AthCommonReentrantAlgorithm.h.

                                                           {
    return execState( ctx ).filterPassed();
  }

finalize()

StatusCode DerivationFramework::FourLeptonVertexingAlgorithm::finalize ( )

override

Definition at line 62 of file FourLeptonVertexingAlgorithm.cxx.

                                                      {
        unsigned int n_processed = std::accumulate(m_num_lep.begin(), m_num_lep.end(), 0, [](const std::atomic<unsigned int>& l, unsigned int n){
            return n +l;
        });        
        ATH_MSG_INFO("Proccessed in total "<<n_processed<<" events. Lepton multiplicities observed:");
        for (unsigned int l = 0 ; l < m_num_lep.size(); ++l) {
            ATH_MSG_INFO("    ---- Events with "<<l<<" leptons: "<<m_num_lep[l]<<" (elec/muon) "<<m_num_ele[l]<<"/"<<m_num_muo[l]<<" ("<<(100.*m_num_lep[l] / std::max(1u,n_processed))<<"%)" );            
        }
        const unsigned int n_trials = m_tried_fits;
        ATH_MSG_INFO("---> Attempted fits "<<m_tried_fits<<" out of which "<<m_good_fits<<" ("<<(100.* m_good_fits / std::max(1u, n_trials) )<<"%) succeeded." );
        return StatusCode::SUCCESS;
    }

fitQuadruplet()

std::unique_ptr< xAOD::Vertex > DerivationFramework::FourLeptonVertexingAlgorithm::fitQuadruplet ( const EventContext & ctx, const LeptonQuadruplet & quad ) const

private

Use the beam spot as an initial constraint

Definition at line 181 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                                                                     {
        if (!passSelection(quad)) return nullptr;
        ++m_tried_fits;
        std::vector<const xAOD::TrackParticle*> trks{};
        trks.reserve(4);
        for (const xAOD::IParticle* lep : quad) { trks.push_back(trackParticle(lep)); }
        SG::ReadHandle<xAOD::EventInfo> evtInfo{m_evtKey, ctx};
        if (!evtInfo.isValid()) {
            ATH_MSG_FATAL("Failed to retrieve the event info " << m_evtKey.fullKey());
            throw std::runtime_error("Invalid key access");
        }
 
        const Amg::Vector3D beam_spot{evtInfo->beamPosX(), evtInfo->beamPosY(), evtInfo->beamPosZ()};
        std::unique_ptr<xAOD::Vertex> common_vtx = m_fitter->fit(ctx, trks, beam_spot);
        if (!common_vtx) {
            ATH_MSG_DEBUG("Fit from  " << quad << " failed.");
            return nullptr;
        }
        const float redChi2 = (common_vtx->chiSquared() / common_vtx->numberDoF());
        ATH_MSG_DEBUG("Fit from " << quad << " gave a vertex with position at " << common_vtx->position() << " with chi2 "
                                  << redChi2 << " nDoF: " << common_vtx->numberDoF());
        if (redChi2 > m_VtxChi2Cut) {
            ATH_MSG_DEBUG("Reject due to bad chi2");
            return nullptr;
        }
        static const std::vector<float> empty_vec{};
        if (m_pruneWeight) common_vtx->setTrackWeights(empty_vec);
        if (m_pruneCov) common_vtx->setCovariance(empty_vec);
        std::vector<TrkLink> track_links{};
        for (const xAOD::TrackParticle* trk : trks) {
            const xAOD::TrackParticleContainer* cont = dynamic_cast<const xAOD::TrackParticleContainer*>(trk->container());
            TrkLink link{*cont, trk->index(), ctx};
            link.toPersistent();
            track_links.push_back(std::move(link));
        }
        common_vtx->setTrackParticleLinks(track_links);
        ++m_good_fits;
        return common_vtx;
    }

initialize()

StatusCode DerivationFramework::FourLeptonVertexingAlgorithm::initialize ( )

override

Definition at line 45 of file FourLeptonVertexingAlgorithm.cxx.

                                                        {
        ATH_CHECK(m_muonKey.initialize());
        ATH_CHECK(m_elecKey.initialize());
        ATH_CHECK(m_evtKey.initialize());
        ATH_CHECK(m_vtxKey.initialize());
        ATH_CHECK(m_fitter.retrieve());
 
        if (!m_elecSelTool.empty()) ATH_CHECK(m_elecSelTool.retrieve());
        if (!m_muonSelTool.empty()) ATH_CHECK(m_muonSelTool.retrieve());
 
        if (m_muonTrkProp < MuonTrk::Primary || m_muonTrk > MuonTrk::MSOnlyExtrapolatedMuonSpectrometerTrackParticle) {
            ATH_MSG_FATAL("A bogous muon track particle has been picked " << m_muonTrkProp);
            return StatusCode::FAILURE;
        }
        m_muonTrk = static_cast<MuonTrk>(m_muonTrkProp.value());
        return StatusCode::SUCCESS;
    }

inputHandles()

virtual std::vector< Gaudi::DataHandle * > AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

isClonable()

bool AthCommonReentrantAlgorithm< Gaudi::Algorithm >::isClonable ( ) const

overridevirtualinherited

Specify if the algorithm is clonable.

Reentrant algorithms are clonable.

Reimplemented in InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ITkPixelCablingAlg, ITkStripCablingAlg, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_FlaggedConditionTestAlg, SCT_LinkMaskingTestAlg, SCT_MajorityConditionsTestAlg, SCT_ModuleVetoTestAlg, SCT_MonitorConditionsTestAlg, SCT_PrepDataToxAOD, SCT_RawDataToxAOD, SCT_ReadCalibChipDataTestAlg, SCT_ReadCalibDataTestAlg, SCT_RODVetoTestAlg, SCT_SensorsTestAlg, SCT_SiliconConditionsTestAlg, SCT_StripVetoTestAlg, SCT_TdaqEnabledTestAlg, SCT_TestCablingAlg, SCTEventFlagWriter, SCTRawDataProvider, SCTSiLorentzAngleTestAlg, SCTSiPropertiesTestAlg, and Simulation::BeamEffectsAlg.

Definition at line 68 of file AthCommonReentrantAlgorithm.cxx.

{
  // Reentrant algorithms are clonable.
  return true;
}

msg()

MsgStream & AthCommonMsg< Gaudi::Algorithm >::msg ( ) const

inlineinherited

Definition at line 24 of file AthCommonMsg.h.

                                {
    return this->msgStream();
  }

msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::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< Gaudi::Algorithm > >::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.

passSelection() [1/3]

bool DerivationFramework::FourLeptonVertexingAlgorithm::passSelection ( const LeptonQuadruplet & quad ) const

private

Check the di lepton masses in the quad

If the pair is SFOS check the mass

The high mass pair must at least satisfy a certain threshold

Definition at line 221 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                       {
        float max_m{-1.};
        for (unsigned int i = 1; i < quad.size(); ++i) {
            for (unsigned int j = 0; j < i; ++j) {
                const xAOD::TrackParticle* trk_a = trackParticle(quad[i]);
                const xAOD::TrackParticle* trk_b = trackParticle(quad[j]);                
                
                if (std::abs(trk_b->z0() - trk_a->z0()) > m_z0Cut) return false;
                
                const float di_m1 = (quad[i]->p4() + quad[j]->p4()).M();
                const bool isSFOS1 = charge(quad[i])*charge(quad[j]) < 0 && quad[i]->type() == quad[j]->type();
                if (isSFOS1 && di_m1 < m_lowSFOS_Cut) return false;
                max_m = std::max(di_m1, max_m);                
                for (unsigned int k = 1 ; k < quad.size(); ++k) {
                    if (k == i || k == j) continue;
                    for (unsigned int l = 0; l < k ; ++l) {
                        if (l == i || l == j) continue;
                        const float di_m2 =  (quad[k]->p4() + quad[l]->p4()).M();
                        const bool isSFOS2 = charge(quad[k])*charge(quad[l]) < 0 && quad[k]->type() == quad[l]->type();
                        if (isSFOS2 && di_m2 < m_lowSFOS_Cut) return false;
                        max_m = std::max(max_m, di_m2);
                    }
                }          
            }
        }      
        if (m_highPair_Cut > max_m) return false;
 
        return true;
    }

passSelection() [2/3]

bool DerivationFramework::FourLeptonVertexingAlgorithm::passSelection ( const xAOD::Electron * elec ) const

private

Definition at line 88 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                   {
        return elec->pt() >= m_minElecPt && (m_elecSelTool.empty() || m_elecSelTool->accept(elec)) &&
               !xAOD::EgammaHelpers::getTrackParticlesVec(elec, !m_elecUseGSF).empty();
    }

passSelection() [3/3]

bool DerivationFramework::FourLeptonVertexingAlgorithm::passSelection ( const xAOD::Muon * muon ) const

private

Definition at line 92 of file FourLeptonVertexingAlgorithm.cxx.

                                                                               {
        return muon->pt() >= m_minMuonPt && (m_muonSelTool.empty() || m_muonSelTool->accept(*muon)) && muon->trackParticle(m_muonTrk);
    }

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< Gaudi::Algorithm > >::renounce ( T & h )

inlineprotectedinherited

Definition at line 380 of file AthCommonDataStore.h.

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

renounceArray()

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

inlineprotectedinherited

remove all handles from I/O resolution

Definition at line 364 of file AthCommonDataStore.h.

                                                          {
    handlesArray.renounce();
  }

setFilterPassed()

virtual void AthCommonReentrantAlgorithm< Gaudi::Algorithm >::setFilterPassed ( bool state, const EventContext & ctx ) const

inlinevirtualinherited

Definition at line 100 of file AthCommonReentrantAlgorithm.h.

                                                                            {
    execState( ctx ).setFilterPassed( state );
  }

sysExecute()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysExecute ( const EventContext & ctx )

overridevirtualinherited

Execute an algorithm.

We override this in order to work around an issue with the Algorithm base class storing the event context in a member variable that can cause crashes in MT jobs.

Definition at line 85 of file AthCommonReentrantAlgorithm.cxx.

{
  return BaseAlg::sysExecute (ctx);
}

sysInitialize()

StatusCode AthCommonReentrantAlgorithm< Gaudi::Algorithm >::sysInitialize ( )

overridevirtualinherited

Override sysInitialize.

Override sysInitialize from the base class.

Loop through all output handles, and if they're WriteCondHandles, automatically register them and this Algorithm with the CondSvc

Scan through all outputHandles, and if they're WriteCondHandles, register them with the CondSvc

Reimplemented from AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >.

Reimplemented in HypoBase, and InputMakerBase.

Definition at line 61 of file AthCommonReentrantAlgorithm.cxx.

                                                               {
  StatusCode sc=AthCommonDataStore<AthCommonMsg<BaseAlg>>::sysInitialize();
 
  if (sc.isFailure()) {
    return sc;
  }
  
  ServiceHandle<ICondSvc> cs("CondSvc",name());
  for (auto h : outputHandles()) {
    if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
      // do this inside the loop so we don't create the CondSvc until needed
      if ( cs.retrieve().isFailure() ) {
        ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
        return StatusCode::SUCCESS;
      }      
      if (cs->regHandle(this,*h).isFailure()) {
        sc = StatusCode::FAILURE;
        ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
                      << " with CondSvc");
      }
    }
  }
  return sc;  
}

sysStart()

virtual StatusCode AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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.

trackParticle()

const xAOD::TrackParticle * DerivationFramework::FourLeptonVertexingAlgorithm::trackParticle ( const xAOD::IParticle * part ) const

private

Definition at line 169 of file FourLeptonVertexingAlgorithm.cxx.

                                                                                                        {
        if (part->type() == xAOD::Type::Muon)
            return static_cast<const xAOD::Muon*>(part)->trackParticle(m_muonTrk);
        else if (part->type() == xAOD::Type::Electron)
            return xAOD::EgammaHelpers::getTrackParticlesVec(static_cast<const xAOD::Electron*>(part), !m_elecUseGSF)[0];
        return dynamic_cast<const xAOD::TrackParticle*>(part);
    }

updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::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

ATLAS_THREAD_SAFE [1/3]

std::array<std::atomic<unsigned int>, 8> m_num_lep DerivationFramework::FourLeptonVertexingAlgorithm::ATLAS_THREAD_SAFE {}

mutableprivate

Simple counter to evaluate the number of leptons per event.

Definition at line 91 of file FourLeptonVertexingAlgorithm.h.

{};

ATLAS_THREAD_SAFE [2/3]

std::array<std::atomic<unsigned int>, 8> m_num_ele DerivationFramework::FourLeptonVertexingAlgorithm::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 92 of file FourLeptonVertexingAlgorithm.h.

{};

ATLAS_THREAD_SAFE [3/3]

std::array<std::atomic<unsigned int>, 8> m_num_muo DerivationFramework::FourLeptonVertexingAlgorithm::ATLAS_THREAD_SAFE {}

mutableprivate

Definition at line 93 of file FourLeptonVertexingAlgorithm.h.

{};

m_detStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_detStore

privateinherited

Pointer to StoreGate (detector store by default)

Definition at line 393 of file AthCommonDataStore.h.

m_elecKey

SG::ReadHandleKey<xAOD::ElectronContainer> DerivationFramework::FourLeptonVertexingAlgorithm::m_elecKey

private

Initial value:

{this, "ElectronContainer", "Electrons",
                                                             "Location of the electron container in the store gate"}

Definition at line 50 of file FourLeptonVertexingAlgorithm.h.

                                                        {this, "ElectronContainer", "Electrons",
                                                             "Location of the electron container in the store gate"};

m_elecSelTool

ToolHandle<IAsgElectronLikelihoodTool> DerivationFramework::FourLeptonVertexingAlgorithm::m_elecSelTool

private

Initial value:

{this, "ElectronSelectionTool", "",
                                                             "Applied preselection on the electrons. Will be ignored if empty"}

Definition at line 56 of file FourLeptonVertexingAlgorithm.h.

                                                            {this, "ElectronSelectionTool", "",
                                                             "Applied preselection on the electrons. Will be ignored if empty"};

m_elecUseGSF

Gaudi::Property<bool> DerivationFramework::FourLeptonVertexingAlgorithm::m_elecUseGSF

private

Initial value:

{this, "PickElecGsfTrk", true,
                                           "If this property is set to true it will pick up the electron GSF track"}

Definition at line 70 of file FourLeptonVertexingAlgorithm.h.

                                        {this, "PickElecGsfTrk", true,
                                           "If this property is set to true it will pick up the electron GSF track"};

m_evtKey

SG::ReadHandleKey<xAOD::EventInfo> DerivationFramework::FourLeptonVertexingAlgorithm::m_evtKey {this, "EventInfoKey", "EventInfo"}

private

Definition at line 52 of file FourLeptonVertexingAlgorithm.h.

{this, "EventInfoKey", "EventInfo"};

m_evtStore

StoreGateSvc_t AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_evtStore

privateinherited

Pointer to StoreGate (event store by default)

Definition at line 390 of file AthCommonDataStore.h.

m_extendedExtraObjects

DataObjIDColl AthCommonReentrantAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects

privateinherited

Extra output dependency collection, extended by AthAlgorithmDHUpdate to add symlinks.

Empty if no symlinks were found.

Definition at line 114 of file AthCommonReentrantAlgorithm.h.

m_fitter

ToolHandle<Trk::IVertexFitter> DerivationFramework::FourLeptonVertexingAlgorithm::m_fitter {this, "VertexFitter", "", "Vertex fitter tool"}

private

Definition at line 59 of file FourLeptonVertexingAlgorithm.h.

{this, "VertexFitter", "", "Vertex fitter tool"};

m_good_fits

std::atomic<unsigned int> DerivationFramework::FourLeptonVertexingAlgorithm::m_good_fits {0}

mutableprivate

How many vertex fits were actually successful.

Definition at line 98 of file FourLeptonVertexingAlgorithm.h.

{0};

m_highPair_Cut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_highPair_Cut {this, "HighMassPair", 40.*Gaudi::Units::GeV, "Minimum mass for the largest lepton pair"}

private

Definition at line 78 of file FourLeptonVertexingAlgorithm.h.

{this, "HighMassPair", 40.*Gaudi::Units::GeV, "Minimum mass for the largest lepton pair"};

m_LeadPtCut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_LeadPtCut {this, "LeadingLeptonPt", 15. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "}

private

Definition at line 80 of file FourLeptonVertexingAlgorithm.h.

{this, "LeadingLeptonPt", 15. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "};

m_lowSFOS_Cut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_lowSFOS_Cut {this, "LowSFOSMass", 3.5 * Gaudi::Units::GeV, "Minimal mass for the lower SFOS pair"}

private

Definition at line 77 of file FourLeptonVertexingAlgorithm.h.

{this, "LowSFOSMass", 3.5 * Gaudi::Units::GeV, "Minimal mass for the lower SFOS pair"};

m_minElecPt

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_minElecPt {this, "MinElecPt", 4.5 * Gaudi::Units::GeV, " Minimum pt cut applied on the electron"}

private

Definition at line 63 of file FourLeptonVertexingAlgorithm.h.

{this, "MinElecPt", 4.5 * Gaudi::Units::GeV, " Minimum pt cut applied on the electron"};

m_minMuonPt

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_minMuonPt

private

Initial value:

{this, "MinMuonPt", 3. * Gaudi::Units::GeV,
                                           "Minimum pt cut applied on the muon -- Unfortunately not handled by the MST"}

Definition at line 61 of file FourLeptonVertexingAlgorithm.h.

                                        {this, "MinMuonPt", 3. * Gaudi::Units::GeV,
                                           "Minimum pt cut applied on the muon -- Unfortunately not handled by the MST"};

m_muonKey

SG::ReadHandleKey<xAOD::MuonContainer> DerivationFramework::FourLeptonVertexingAlgorithm::m_muonKey

private

Initial value:

{this, "MuonContainer", "Muons",
                                                         "Location of the input muon container in the StoreGate"}

Input containers.

Definition at line 48 of file FourLeptonVertexingAlgorithm.h.

                                                    {this, "MuonContainer", "Muons",
                                                         "Location of the input muon container in the StoreGate"};

m_muonSelTool

ToolHandle<CP::IMuonSelectionTool> DerivationFramework::FourLeptonVertexingAlgorithm::m_muonSelTool

private

Initial value:

{this, "MuonSelectionTool", "",
                                                         "Applied preselection on the muons. Will be ignored if empty"}

Definition at line 54 of file FourLeptonVertexingAlgorithm.h.

                                                      {this, "MuonSelectionTool", "",
                                                         "Applied preselection on the muons. Will be ignored if empty"};

m_muonTrk

MuonTrk DerivationFramework::FourLeptonVertexingAlgorithm::m_muonTrk {MuonTrk::Primary}

private

Definition at line 68 of file FourLeptonVertexingAlgorithm.h.

{MuonTrk::Primary};

m_muonTrkProp

Gaudi::Property<int> DerivationFramework::FourLeptonVertexingAlgorithm::m_muonTrkProp

private

Initial value:

{this, "PickMuonTrk", MuonTrk::InnerDetectorTrackParticle,
                                           "Pick the proper track particle from the muon"}

Definition at line 66 of file FourLeptonVertexingAlgorithm.h.

                                        {this, "PickMuonTrk", MuonTrk::InnerDetectorTrackParticle,
                                           "Pick the proper track particle from the muon"};

m_pruneCov

Gaudi::Property<bool> DerivationFramework::FourLeptonVertexingAlgorithm::m_pruneCov {this, "PruneCovariance", true, "Clears the covariance vector"}

private

Definition at line 87 of file FourLeptonVertexingAlgorithm.h.

{this, "PruneCovariance", true, "Clears the covariance vector"};

m_pruneWeight

Gaudi::Property<bool> DerivationFramework::FourLeptonVertexingAlgorithm::m_pruneWeight {this, "PruneTrackWeights", true, "Clears the track weight vector"}

private

Definition at line 88 of file FourLeptonVertexingAlgorithm.h.

{this, "PruneTrackWeights", true, "Clears the track weight vector"};

m_SubLeadPtCut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_SubLeadPtCut {this, "SubLeadingLeptonPt", 10. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "}

private

Definition at line 81 of file FourLeptonVertexingAlgorithm.h.

{this, "SubLeadingLeptonPt", 10. *Gaudi::Units::GeV,"Minimal momentum of the leading lepton in the quad "};

m_tried_fits

std::atomic<unsigned int> DerivationFramework::FourLeptonVertexingAlgorithm::m_tried_fits {0}

mutableprivate

How many vertex fits were attempted.

Definition at line 96 of file FourLeptonVertexingAlgorithm.h.

{0};

m_varHandleArraysDeclared

bool AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_varHandleArraysDeclared

privateinherited

Definition at line 399 of file AthCommonDataStore.h.

m_vhka

std::vector<SG::VarHandleKeyArray*> AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::m_vhka

privateinherited

Definition at line 398 of file AthCommonDataStore.h.

m_VtxChi2Cut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_VtxChi2Cut {this, "VertexChi2Cut", 20, "Maximal chii n.D.o.F cut on the reconstructed vertex" }

private

Definition at line 83 of file FourLeptonVertexingAlgorithm.h.

{this, "VertexChi2Cut", 20, "Maximal chii n.D.o.F cut on the reconstructed vertex" };

m_vtxKey

SG::WriteHandleKey<xAOD::VertexContainer> DerivationFramework::FourLeptonVertexingAlgorithm::m_vtxKey {this, "OutVertexContainer", "FourLeptonVertices", "Output vertex container"}

private

Definition at line 86 of file FourLeptonVertexingAlgorithm.h.

{this, "OutVertexContainer", "FourLeptonVertices", "Output vertex container"};

m_z0Cut

Gaudi::Property<float> DerivationFramework::FourLeptonVertexingAlgorithm::m_z0Cut

private

Initial value:

{this, "DeltaZ0Cut", 7.5 * Gaudi::Units::mm,
                                       "All leptons in the quadruplet have to be seperated from each other by maximum"}

Definition at line 73 of file FourLeptonVertexingAlgorithm.h.

                                    {this, "DeltaZ0Cut", 7.5 * Gaudi::Units::mm,
                                       "All leptons in the quadruplet have to be seperated from each other by maximum"};

---

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

• FourLeptonVertexingAlgorithm.h
• FourLeptonVertexingAlgorithm.cxx