ATLAS Offline Software
Loading...
Searching...
No Matches
xAODLeptonPairFilter Class Reference

#include <xAODLeptonPairFilter.h>

Inheritance diagram for xAODLeptonPairFilter:
Collaboration diagram for xAODLeptonPairFilter:

Public Member Functions

virtual StatusCode filterInitialize () override final
virtual StatusCode filterEvent (const EventContext &ctx) override final
 GenFilter (const std::string &name, ISvcLocator *pSvcLocator)
virtual StatusCode sysInitialize () override
 Override sysInitialize.
virtual bool isClonable () const override
 Specify if the algorithm is clonable.
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
 Get filter decision:
virtual void setFilterPassed (bool state, const EventContext &ctx) const
 Set filter decision:
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
Event loop algorithm methods: not to be overloaded
StatusCode initialize ()
StatusCode execute (const EventContext &ctx)
 Execute method.
StatusCode finalize ()
Gen-specific event loop methods: to be overloaded!
virtual StatusCode filterFinalize ()
Counters
int nPassed () const
int nFailed () const
int nNeeded () const
bool keepAll () const
Event collection accessors (const and non-const)
HepMC::GenEvent *event ATLAS_NOT_CONST_THREAD_SAFE ()
 Access the current signal event (first in the McEventCollection).
McEventCollection *events ATLAS_NOT_CONST_THREAD_SAFE ()
 Access the current event's McEventCollection.
const HepMC::GenEventevent_const (const EventContext &ctx) const
 Access the current signal event (const).
const McEventCollectionevents_const (const EventContext &ctx) const
 Access the current event's McEventCollection (const).

Protected Member Functions

virtual bool isReEntrant () const override final
 Legacy algorithms are not thread-safe.
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.

Protected Attributes

Counters and requirements
int m_nPass
int m_nFail
int m_nNeeded
bool m_keepAll

Private Types

typedef ServiceHandle< StoreGateSvcStoreGateSvc_t

Private Member Functions

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::TruthParticleContainerm_truthPartContKey {this, "TruthParticleContainerKey", "TruthGen"}
Gaudi::Property< int > m_nSFOS_Max {this,"NSFOS_Max",0,"Number of Same Flavor Opposite Sign Max"}
Gaudi::Property< int > m_nSFOS_Min {this,"NSFOS_Min",0,"Number of Same Flavor Opposite Sign Min"}
Gaudi::Property< int > m_nSFSS_Max {this,"NSFSS_Max",-1,"Number of Same Flavor Same Sign Max"}
Gaudi::Property< int > m_nSFSS_Min {this,"NSFSS_Min",-1,"Number of Same Flavor Same Sign Min"}
Gaudi::Property< int > m_nOFOS_Max {this,"NOFOS_Max",-1,"Number of Opposite Flavor Opposite Sign Max"}
Gaudi::Property< int > m_nOFOS_Min {this,"NOFOS_Min",-1,"Number of Opposite Flavor Opposite Sign Min"}
Gaudi::Property< int > m_nOFSS_Max {this,"NOFSS_Max",-1,"Number of Opposite Flavor Same Sign Max"}
Gaudi::Property< int > m_nOFSS_Min {this,"NOFSS_Min",-1,"Number of Opposite Flavor Same Sign Min"}
Gaudi::Property< int > m_nPairSum_Max {this,"NPairSum_Max",-1,"Number of Pair sum Max"}
Gaudi::Property< int > m_nPairSum_Min {this,"NPairSum_Min",-1,"Number of Pair sum Min"}
Gaudi::Property< bool > m_bUseSFOSInSum {this,"UseSFOSInSum",false,"SFOS to be considered in sum"}
Gaudi::Property< bool > m_bUseSFSSInSum {this,"UseSFSSInSum",false,"SFSS to be considered in sum"}
Gaudi::Property< bool > m_bUseOFOSInSum {this,"UseOFOSInSum",false,"OFOS to be considered in sum"}
Gaudi::Property< bool > m_bUseOFSSInSum {this,"UseOFSSInSum",false,"OFSS to be considered in sum"}
Gaudi::Property< bool > m_onlyMassiveParents {this,"OnlyMassiveParents",false,"only count leptons with a parent with M>20 GeV"}
Gaudi::Property< double > m_Ptmin {this,"Ptcut",10000.,"min Pt"}
Gaudi::Property< double > m_EtaRange {this,"Etacut",10.,"Eta range max"}
Gaudi::Property< int > m_nLeptons_Max {this,"NLeptons_Max",-1," max num leptons"}
Gaudi::Property< int > m_nLeptons_Min {this,"NLeptons_Min",0," min num leptons"}
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

Properties

SG::ReadHandleKey< McEventCollectionm_mcevents_const { this, "McEventKey", "GEN_EVENT", "StoreGate key of the MC event collection" }
 Const handle to the MC event collection.
std::string m_mcEventKey {}
 StoreGate key for the MC event collection (defaults to GEN_EVENT).
BooleanProperty m_mkMcEvent {this, "MakeMcEvent", false, "Create a new MC event collection if it doesn't exist"}
 Flag to determine if a new MC event collection should be made if it doesn't exist.

Detailed Description

Definition at line 25 of file xAODLeptonPairFilter.h.

Member Typedef Documentation

◆ StoreGateSvc_t

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

Definition at line 388 of file AthCommonDataStore.h.

Member Function Documentation

◆ ATLAS_NOT_CONST_THREAD_SAFE() [1/2]

McEventCollection *events GenBase::ATLAS_NOT_CONST_THREAD_SAFE ( )
inherited

Access the current event's McEventCollection.

Note
This function will make a new McEventCollection if there is not already a valid one and MakeMcEvent=True.

◆ ATLAS_NOT_CONST_THREAD_SAFE() [2/2]

HepMC::GenEvent *event GenBase::ATLAS_NOT_CONST_THREAD_SAFE ( )
inlineinherited

Access the current signal event (first in the McEventCollection).

Note
This function will make a new McEventCollection if there is not already a valid one and MakeMcEvent=True.

Definition at line 74 of file GenBase.h.

74 {
75 if (events()->empty())
76 ATH_MSG_ERROR("McEventCollection is empty during first event access");
77 return *(events()->begin());
78 }
#define ATH_MSG_ERROR(x)
static const Attributes_t empty

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

158 {
160 hndl.value(),
161 hndl.documentation());
162
163 }
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)

◆ declareProperty()

Gaudi::Details::PropertyBase & AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::declareProperty ( Gaudi::Property< T, V, H > & t)
inlineinherited

Definition at line 145 of file AthCommonDataStore.h.

145 {
146 typedef typename SG::HandleClassifier<T>::type htype;
148 }
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>

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

◆ event_const()

const HepMC::GenEvent * GenBase::event_const ( const EventContext & ctx) const
inlineinherited

Access the current signal event (const).

Definition at line 81 of file GenBase.h.

81 {
82 const McEventCollection* coll = events_const(ctx);
83 if (coll->empty())
84 ATH_MSG_ERROR("Const McEventCollection is empty during first event access");
85 return *(coll->begin());
86 }
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
bool empty() const noexcept
Returns true if the collection is empty.
const McEventCollection * events_const(const EventContext &ctx) const
Access the current event's McEventCollection (const).
Definition GenBase.h:95

◆ events_const()

const McEventCollection * GenBase::events_const ( const EventContext & ctx) const
inlineinherited

Access the current event's McEventCollection (const).

Definition at line 95 of file GenBase.h.

95 {
96 SG::ReadHandle<McEventCollection> ret = SG::makeHandle(m_mcevents_const, ctx);
97 if (!ret.isValid())
98 ATH_MSG_ERROR("No McEventCollection found in StoreGate with key " << m_mcevents_const.key());
99 return ret.cptr();
100 }
SG::ReadHandleKey< McEventCollection > m_mcevents_const
Const handle to the MC event collection.
Definition GenBase.h:119
virtual bool isValid() override final
Can the handle be successfully dereferenced?
const_pointer_type cptr()
Dereference the pointer.
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())

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

◆ execute()

StatusCode GenFilter::execute ( const EventContext & ctx)
virtualinherited

Execute method.

Todo
Probably the filter should only look at the first event... right?

Reimplemented from GenBase.

Definition at line 27 of file GenFilter.cxx.

27 {
28 if (events_const(ctx)->empty()) {
29 ATH_MSG_ERROR("No events found in McEventCollection");
30 return StatusCode::FAILURE;
31 } else if (events_const(ctx)->size() > 1) {
33 ATH_MSG_WARNING("More than one event in current McEventCollection -- which is valid?");
34 }
36 if (filterPassed(ctx) || m_keepAll ) {
37 ATH_MSG_DEBUG("Event passed filter");
38 m_nPass += 1;
39 } else {
40 ATH_MSG_DEBUG("Event failed filter");
41 m_nFail += 1;
42 }
43 // Bail out once we have enough events
44 if (m_nPass >= m_nNeeded && m_nNeeded > 0)
45 sc = StatusCode::FAILURE;
46 return sc;
47}
#define ATH_MSG_WARNING(x)
#define ATH_MSG_DEBUG(x)
static Double_t sc
size_t size() const
Number of registered mappings.
virtual bool filterPassed(const EventContext &ctx) const
bool m_keepAll
Definition GenFilter.h:66
int m_nPass
Definition GenFilter.h:63
int m_nNeeded
Definition GenFilter.h:65
virtual StatusCode filterEvent(const EventContext &ctx)=0
int m_nFail
Definition GenFilter.h:64
::StatusCode StatusCode
StatusCode definition for legacy code.

◆ 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 & AthCommonAlgorithm< 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 89 of file AthCommonAlgorithm.cxx.

54{
55 // If we didn't find any symlinks to add, just return the collection
56 // from the base class. Otherwise, return the extended collection.
57 if (!m_extendedExtraObjects.empty()) {
59 }
61}
Common base class for algorithms.

◆ filterEvent()

StatusCode xAODLeptonPairFilter::filterEvent ( const EventContext & ctx)
finaloverridevirtual

Implements GenFilter.

Definition at line 36 of file xAODLeptonPairFilter.cxx.

36 {
37//---------------------------------------------------------------------------
38
39 // Loop over all events in McEventCollection
40 std::vector<int> vLeptonPDGIDs;
41 std::vector<double> vLeptonPt;
42 std::vector<double> vLeptonEta;
43 std::vector< std::vector < size_t > > vLeptonParentPDGIDs;
44
45 // Retrieve TruthGen container from xAOD Gen slimmer, contains all particles witout barcode_zero and
46 // duplicated barcode ones
47 SG::ReadHandle<xAOD::TruthParticleContainer> xTruthParticleContainer{m_truthPartContKey, ctx};
48 CHECK(xTruthParticleContainer.isValid());
49
50 // Loop over all particles in the event
51 for (const xAOD::TruthParticle* pitr : *xTruthParticleContainer) {
52 if( !MC::isStable(pitr) ) continue;
53 // check stable particles only
54 // We do not place requirements on their origins (updated: optionally rejecting hadron decays)
55 // save pdg ids of found leptons
56 // do not consider taus
57 if( !MC::isElectron(pitr) && !MC::isMuon(pitr) ) continue;
58 //only consider leptons which satisfy pt and eta requirements
59 if( (pitr->pt() < m_Ptmin) || std::abs(pitr->eta()) > m_EtaRange) continue;
61 {
62 auto p = pitr;
63 bool massiveParent = false;
64 while(p)
65 {
66 auto vxp = p->prodVtx();
67 if(!vxp) break;
68 if(vxp->nIncomingParticles()!=1) break;
69 p = vxp->incomingParticle(0);
70 const int pdg = std::abs(p->pdgId());
71 if(!(MC::isSMLepton(pdg) || MC::isPhoton(pdg)))
72 {
73 massiveParent = (p->m()>20000);
74 break;
75 }
76 }
77 if(!massiveParent) continue;
78 }
79
80 vLeptonPDGIDs.push_back(pitr->pdgId());
81 vLeptonPt.push_back(pitr->pt());
82 vLeptonEta.push_back(pitr->eta());
83
84 std::vector<size_t> parentPDG_tmp;
85 for(size_t thisParent_id=0; thisParent_id<pitr->prodVtx()->nIncomingParticles();thisParent_id++) {
86 parentPDG_tmp.push_back(pitr->prodVtx()->incomingParticle(thisParent_id)->pdgId());
87 }
88 vLeptonParentPDGIDs.push_back(std::move(parentPDG_tmp));
89 } //loop over TruthParticles
90
91
92 int nLeptons = vLeptonPDGIDs.size();
93
94 //Filter automatically passes if number of leptons is different than expected
95 if (!( (nLeptons >= m_nLeptons_Min) &&
96 (m_nLeptons_Max < 0 || nLeptons <= m_nLeptons_Max) )) {
97 ATH_MSG_INFO("# Lep = "<<nLeptons << " Pass" );
98 return StatusCode::SUCCESS;
99 }
100
101 int nSFOS = 0;
102 int nSFSS = 0;
103 int nOFOS = 0;
104 int nOFSS = 0;
105
106 int id1 = 0;
107 int id2 = 0;
108 //loop over all possible pair combinations of the leptons
109 //which were found.
110 for (unsigned int i = 0 ; i < vLeptonPDGIDs.size() ; i++){
111 //don't double count
112 for (unsigned int j = i ; j < vLeptonPDGIDs.size() ; j++){
113 //Do not compare the lepton with itself
114 if(i==j) continue;
115 id1 = vLeptonPDGIDs[i];
116 id2 = vLeptonPDGIDs[j];
117 //classify the pair and count it
118 if(std::abs(id1)==std::abs(id2) && id1*id2 < 0) nSFOS+=1;
119 else if(std::abs(id1)==std::abs(id2) && id1*id2 > 0) nSFSS+=1;
120 else if(std::abs(id1)!=std::abs(id2) && id1*id2 < 0) nOFOS+=1;
121 else if(std::abs(id1)!=std::abs(id2) && id1*id2 > 0) nOFSS+=1;
122 else ATH_MSG_ERROR( "Couldn't classify lepton pair" );
123 }
124 }
125
126 //based on configuration
127 //create sum of multiple pair combinations
128 //which can alternatively be cut on instead
129 //of the individual ones
130 int nPairSum = 0;
131 if(m_bUseSFOSInSum) nPairSum+=nSFOS;
132 if(m_bUseSFSSInSum) nPairSum+=nSFSS;
133 if(m_bUseOFOSInSum) nPairSum+=nOFOS;
134 if(m_bUseOFSSInSum) nPairSum+=nOFSS;
135
136
137 bool passSFOS = false;
138 bool passSFSS = false;
139 bool passOFOS = false;
140 bool passOFSS = false;
141 bool passPairSum = false;
142 //Test if number of lepton pairs satisfies requirements
143 //The maximum requirement is ignored if it is negative
144 //No requirement is made if both min and max values are negative
145 ATH_MSG_INFO("# Lep " << vLeptonPDGIDs.size() << ", "<< nSFOS << " SFOS, "<<nSFSS<< " SFSS, " << nOFOS << " OFOS, " << nOFSS << " OFSS pairs ," << nPairSum << "summed pairs" );
146
147 if(nSFOS >= m_nSFOS_Min && (m_nSFOS_Max<0 || nSFOS <= m_nSFOS_Max)) passSFOS=true;
148 if(nSFSS >= m_nSFSS_Min && (m_nSFSS_Max<0 || nSFSS <= m_nSFSS_Max)) passSFSS=true;
149 if(nOFOS >= m_nOFOS_Min && (m_nOFOS_Max<0 || nOFOS <= m_nOFOS_Max)) passOFOS=true;
150 if(nOFSS >= m_nOFSS_Min && (m_nOFSS_Max<0 || nOFSS <= m_nOFSS_Max)) passOFSS=true;
151 if(nPairSum >= m_nPairSum_Min && (m_nPairSum_Max<0 || nPairSum <= m_nPairSum_Max)) passPairSum=true;
152
153 if(passSFOS && passSFSS && passOFOS && passOFSS && passPairSum){
154 ATH_MSG_INFO("Pass" );
155 for (unsigned int i = 0;i<vLeptonPDGIDs.size();i++){
156 int pdg = vLeptonPDGIDs[i];
157 double pt = vLeptonPt[i];
158 double eta = vLeptonEta[i];
159 msg(MSG::DEBUG) << pdg << ": Pt = "<<pt<<", Eta = "<<eta << ", Parent PDG: ";
160 for (unsigned int j=0;j<vLeptonParentPDGIDs[i].size();j++) msg(MSG::DEBUG) << vLeptonParentPDGIDs[i][j];
161 msg(MSG::DEBUG) << endmsg;
162 }
163 return StatusCode::SUCCESS;
164 }
165
166 // if we get here we have failed
167 setFilterPassed(false, ctx);
168 ATH_MSG_INFO("Fail" );
169 return StatusCode::SUCCESS;
170}
Scalar eta() const
pseudorapidity method
#define endmsg
#define ATH_MSG_INFO(x)
#define CHECK(...)
Evaluate an expression and check for errors.
virtual void setFilterPassed(bool state, const EventContext &ctx) const
Gaudi::Property< int > m_nSFSS_Min
Gaudi::Property< int > m_nOFOS_Max
Gaudi::Property< double > m_Ptmin
Gaudi::Property< int > m_nPairSum_Min
SG::ReadHandleKey< xAOD::TruthParticleContainer > m_truthPartContKey
Gaudi::Property< int > m_nOFSS_Min
Gaudi::Property< double > m_EtaRange
Gaudi::Property< bool > m_bUseSFOSInSum
Gaudi::Property< bool > m_onlyMassiveParents
Gaudi::Property< bool > m_bUseOFSSInSum
Gaudi::Property< int > m_nSFOS_Max
Gaudi::Property< int > m_nLeptons_Min
Gaudi::Property< int > m_nLeptons_Max
Gaudi::Property< bool > m_bUseOFOSInSum
Gaudi::Property< int > m_nPairSum_Max
Gaudi::Property< int > m_nSFSS_Max
Gaudi::Property< int > m_nOFOS_Min
Gaudi::Property< int > m_nOFSS_Max
Gaudi::Property< int > m_nSFOS_Min
Gaudi::Property< bool > m_bUseSFSSInSum
bool isSMLepton(const T &p)
APID: the fourth generation leptons are not standard model leptons.
bool isPhoton(const T &p)
bool isElectron(const T &p)
bool isStable(const T &p)
Identify if the particle is stable, i.e. has not decayed.
bool isMuon(const T &p)
float j(const xAOD::IParticle &, const xAOD::TrackMeasurementValidation &hit, const Eigen::Matrix3d &jab_inv)
TruthParticle_v1 TruthParticle
Typedef to implementation.
MsgStream & msg
Definition testRead.cxx:32

◆ filterFinalize()

◆ filterInitialize()

StatusCode xAODLeptonPairFilter::filterInitialize ( )
finaloverridevirtual

Reimplemented from GenFilter.

Definition at line 29 of file xAODLeptonPairFilter.cxx.

29 {
30//---------------------------------------------------------------------------
31 CHECK(m_truthPartContKey.initialize());
32 return StatusCode::SUCCESS;
33}

◆ filterPassed()

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

Get filter decision:

Definition at line 93 of file AthCommonAlgorithm.h.

93 {
94 return execState( ctx ).filterPassed();
95 }

◆ finalize()

StatusCode GenFilter::finalize ( )
inherited

Definition at line 50 of file GenFilter.cxx.

50 {
51 ATH_MSG_INFO("Events passed = " << m_nPass << " Events failed = " << m_nFail);
53 return StatusCode::SUCCESS;
54}
virtual StatusCode filterFinalize()
Definition GenFilter.h:47

◆ GenFilter()

GenFilter::GenFilter ( const std::string & name,
ISvcLocator * pSvcLocator )

Definition at line 33 of file GenFilter.cxx.

9 : GenBase(name, pSvcLocator)
10{
11 declareProperty("TotalPassed", m_nNeeded=-1);
12 declareProperty("KeepAllEvents", m_keepAll=false);
13 m_nPass = 0;
14 m_nFail = 0;
15}
GenBase(const std::string &name, ISvcLocator *pSvcLocator)
Constructor.
Definition GenBase.cxx:11

◆ initialize()

StatusCode GenFilter::initialize ( )
virtualinherited

Reimplemented from GenBase.

Definition at line 18 of file GenFilter.cxx.

18 {
20 m_nPass = 0;
21 m_nFail = 0;
23 return StatusCode::SUCCESS;
24}
virtual StatusCode initialize() override
Definition GenBase.cxx:17
virtual StatusCode filterInitialize()
Definition GenFilter.h:45

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

virtual bool AthCommonAlgorithm< Gaudi::Algorithm >::isClonable ( ) const
inlineoverridevirtualinherited

Specify if the algorithm is clonable.

Only relevant for non-reentrant algorithms. Actual number of clones needs to be set via the "Cardinality" property.

Reimplemented in AFP_DigiTop, AlgB, AlgT, BCM_Digitization, CscDigitBuilder, CscDigitToCscRDO, G4AtlasAlg, G4RunAlg, HGTD_Digitization, HiveAlgBase, InDet::GNNSeedingTrackMaker, InDet::SCT_Clusterization, InDet::SiSPGNNTrackMaker, InDet::SiSPSeededTrackFinder, InDet::SiTrackerSpacePointFinder, ISF::SimKernelMT, ITk::StripDigitization, ITkPixelCablingAlg, ITkStripCablingAlg, LArHitEMapMaker, LArTTL1Maker, LUCID_DigiTop, LVL1::L1TopoSimulation, MergeCalibHits, MergeGenericMuonSimHitColl, MergeHijingPars, MergeMcEventCollection, MergeTrackRecordCollection, MergeTruthJets, MergeTruthParticles, MuonDigitizer, PileUpMTAlg, PixelDigitization, RoIBResultToxAOD, SCT_ByteStreamErrorsTestAlg, SCT_CablingCondAlgFromCoraCool, SCT_CablingCondAlgFromText, SCT_ConditionsParameterTestAlg, SCT_ConditionsSummaryTestAlg, SCT_ConfigurationConditionsTestAlg, SCT_Digitization, 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, SGInputLoader, Simulation::BeamEffectsAlg, TileHitVecToCnt, TileMuonFitter, TilePulseForTileMuonReceiver, TileRawChannelMaker, TRTDigitization, and ZDC_DigiTop.

Definition at line 68 of file AthCommonAlgorithm.h.

68 {
69 return true;
70 }

◆ isReEntrant()

virtual bool AthAlgorithm::isReEntrant ( ) const
inlinefinaloverrideprotectedvirtualinherited

Legacy algorithms are not thread-safe.

Definition at line 47 of file AthAlgorithm.h.

47{ return false; }

◆ keepAll()

bool GenFilter::keepAll ( ) const
inlineinherited

Definition at line 55 of file GenFilter.h.

55{return m_keepAll; }

◆ msg()

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

Definition at line 24 of file AthCommonMsg.h.

24 {
25 return this->msgStream();
26 }

◆ msgLvl()

bool AthCommonMsg< Gaudi::Algorithm >::msgLvl ( const MSG::Level lvl) const
inlineinherited

Definition at line 30 of file AthCommonMsg.h.

30 {
31 return this->msgLevel(lvl);
32 }

◆ nFailed()

int GenFilter::nFailed ( ) const
inlineinherited

Definition at line 53 of file GenFilter.h.

53{ return m_nFail; }

◆ nNeeded()

int GenFilter::nNeeded ( ) const
inlineinherited

Definition at line 54 of file GenFilter.h.

54{ return m_nNeeded; }

◆ nPassed()

int GenFilter::nPassed ( ) const
inlineinherited

Definition at line 52 of file GenFilter.h.

52{ return m_nPass; }

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

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

381 {
382 h.renounce();
384 }
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)

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

364 {
366 }

◆ setFilterPassed()

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

Set filter decision:

Reimplemented in AthFilterAlgorithm.

Definition at line 99 of file AthCommonAlgorithm.h.

99 {
101 }

◆ sysExecute()

StatusCode AthCommonAlgorithm< 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.

Reimplemented in AthAnalysisAlgorithm.

Definition at line 80 of file AthCommonAlgorithm.cxx.

41{
42 return BaseAlg::sysExecute (ctx);
43}

◆ sysInitialize()

StatusCode AthCommonAlgorithm< 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 AthAnalysisAlgorithm, AthFilterAlgorithm, AthHistogramAlgorithm, HypoBase, InputMakerBase, and PyAthena::Alg.

Definition at line 60 of file AthCommonAlgorithm.cxx.

71 {
73
74 if (sc.isFailure()) {
75 return sc;
76 }
77
78 ServiceHandle<ICondSvc> cs("CondSvc",name());
79 for (auto h : outputHandles()) {
80 if (h->isCondition() && h->mode() == Gaudi::DataHandle::Writer) {
81 // do this inside the loop so we don't create the CondSvc until needed
82 if ( cs.retrieve().isFailure() ) {
83 ATH_MSG_WARNING("no CondSvc found: won't autoreg WriteCondHandles");
85 }
86 if (cs->regHandle(this,*h).isFailure()) {
88 ATH_MSG_ERROR("unable to register WriteCondHandle " << h->fullKey()
89 << " with CondSvc");
90 }
91 }
92 }
93 return sc;
94}
virtual StatusCode sysInitialize() override
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override

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

◆ updateVHKA()

void AthCommonDataStore< AthCommonMsg< Gaudi::Algorithm > >::updateVHKA ( Gaudi::Details::PropertyBase & )
inlineinherited

Definition at line 308 of file AthCommonDataStore.h.

308 {
309 // debug() << "updateVHKA for property " << p.name() << " " << p.toString()
310 // << " size: " << m_vhka.size() << endmsg;
311 for (auto &a : m_vhka) {
313 for (auto k : keys) {
314 k->setOwner(this);
315 }
316 }
317 }

Member Data Documentation

◆ m_bUseOFOSInSum

Gaudi::Property<bool> xAODLeptonPairFilter::m_bUseOFOSInSum {this,"UseOFOSInSum",false,"OFOS to be considered in sum"}
private

Definition at line 65 of file xAODLeptonPairFilter.h.

65{this,"UseOFOSInSum",false,"OFOS to be considered in sum"};

◆ m_bUseOFSSInSum

Gaudi::Property<bool> xAODLeptonPairFilter::m_bUseOFSSInSum {this,"UseOFSSInSum",false,"OFSS to be considered in sum"}
private

Definition at line 66 of file xAODLeptonPairFilter.h.

66{this,"UseOFSSInSum",false,"OFSS to be considered in sum"};

◆ m_bUseSFOSInSum

Gaudi::Property<bool> xAODLeptonPairFilter::m_bUseSFOSInSum {this,"UseSFOSInSum",false,"SFOS to be considered in sum"}
private

Definition at line 63 of file xAODLeptonPairFilter.h.

63{this,"UseSFOSInSum",false,"SFOS to be considered in sum"};

◆ m_bUseSFSSInSum

Gaudi::Property<bool> xAODLeptonPairFilter::m_bUseSFSSInSum {this,"UseSFSSInSum",false,"SFSS to be considered in sum"}
private

Definition at line 64 of file xAODLeptonPairFilter.h.

64{this,"UseSFSSInSum",false,"SFSS to be considered in sum"};

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

Gaudi::Property<double> xAODLeptonPairFilter::m_EtaRange {this,"Etacut",10.,"Eta range max"}
private

Definition at line 75 of file xAODLeptonPairFilter.h.

75{this,"Etacut",10.,"Eta range max"};

◆ 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 AthCommonAlgorithm< Gaudi::Algorithm >::m_extendedExtraObjects
privateinherited

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

Empty if no symlinks were found.

Definition at line 108 of file AthCommonAlgorithm.h.

◆ m_keepAll

bool GenFilter::m_keepAll
protectedinherited

Definition at line 66 of file GenFilter.h.

◆ m_mcEventKey

std::string GenBase::m_mcEventKey {}
protectedinherited

StoreGate key for the MC event collection (defaults to GEN_EVENT).

Definition at line 110 of file GenBase.h.

110{};

◆ m_mcevents_const

SG::ReadHandleKey<McEventCollection> GenBase::m_mcevents_const { this, "McEventKey", "GEN_EVENT", "StoreGate key of the MC event collection" }
privateinherited

Const handle to the MC event collection.

Definition at line 119 of file GenBase.h.

119{ this, "McEventKey", "GEN_EVENT", "StoreGate key of the MC event collection" };

◆ m_mkMcEvent

BooleanProperty GenBase::m_mkMcEvent {this, "MakeMcEvent", false, "Create a new MC event collection if it doesn't exist"}
protectedinherited

Flag to determine if a new MC event collection should be made if it doesn't exist.

Definition at line 112 of file GenBase.h.

112{this, "MakeMcEvent", false, "Create a new MC event collection if it doesn't exist"};

◆ m_nFail

int GenFilter::m_nFail
protectedinherited

Definition at line 64 of file GenFilter.h.

◆ m_nLeptons_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nLeptons_Max {this,"NLeptons_Max",-1," max num leptons"}
private

Definition at line 81 of file xAODLeptonPairFilter.h.

81{this,"NLeptons_Max",-1," max num leptons"};

◆ m_nLeptons_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nLeptons_Min {this,"NLeptons_Min",0," min num leptons"}
private

Definition at line 82 of file xAODLeptonPairFilter.h.

82{this,"NLeptons_Min",0," min num leptons"};

◆ m_nNeeded

int GenFilter::m_nNeeded
protectedinherited

Definition at line 65 of file GenFilter.h.

◆ m_nOFOS_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nOFOS_Max {this,"NOFOS_Max",-1,"Number of Opposite Flavor Opposite Sign Max"}
private

Definition at line 49 of file xAODLeptonPairFilter.h.

49{this,"NOFOS_Max",-1,"Number of Opposite Flavor Opposite Sign Max"};

◆ m_nOFOS_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nOFOS_Min {this,"NOFOS_Min",-1,"Number of Opposite Flavor Opposite Sign Min"}
private

Definition at line 50 of file xAODLeptonPairFilter.h.

50{this,"NOFOS_Min",-1,"Number of Opposite Flavor Opposite Sign Min"};

◆ m_nOFSS_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nOFSS_Max {this,"NOFSS_Max",-1,"Number of Opposite Flavor Same Sign Max"}
private

Definition at line 53 of file xAODLeptonPairFilter.h.

53{this,"NOFSS_Max",-1,"Number of Opposite Flavor Same Sign Max"};

◆ m_nOFSS_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nOFSS_Min {this,"NOFSS_Min",-1,"Number of Opposite Flavor Same Sign Min"}
private

Definition at line 54 of file xAODLeptonPairFilter.h.

54{this,"NOFSS_Min",-1,"Number of Opposite Flavor Same Sign Min"};

◆ m_nPairSum_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nPairSum_Max {this,"NPairSum_Max",-1,"Number of Pair sum Max"}
private

Definition at line 58 of file xAODLeptonPairFilter.h.

58{this,"NPairSum_Max",-1,"Number of Pair sum Max"};

◆ m_nPairSum_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nPairSum_Min {this,"NPairSum_Min",-1,"Number of Pair sum Min"}
private

Definition at line 59 of file xAODLeptonPairFilter.h.

59{this,"NPairSum_Min",-1,"Number of Pair sum Min"};

◆ m_nPass

int GenFilter::m_nPass
protectedinherited

Definition at line 63 of file GenFilter.h.

◆ m_nSFOS_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nSFOS_Max {this,"NSFOS_Max",0,"Number of Same Flavor Opposite Sign Max"}
private

Definition at line 41 of file xAODLeptonPairFilter.h.

41{this,"NSFOS_Max",0,"Number of Same Flavor Opposite Sign Max"};

◆ m_nSFOS_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nSFOS_Min {this,"NSFOS_Min",0,"Number of Same Flavor Opposite Sign Min"}
private

Definition at line 42 of file xAODLeptonPairFilter.h.

42{this,"NSFOS_Min",0,"Number of Same Flavor Opposite Sign Min"};

◆ m_nSFSS_Max

Gaudi::Property<int> xAODLeptonPairFilter::m_nSFSS_Max {this,"NSFSS_Max",-1,"Number of Same Flavor Same Sign Max"}
private

Definition at line 45 of file xAODLeptonPairFilter.h.

45{this,"NSFSS_Max",-1,"Number of Same Flavor Same Sign Max"};

◆ m_nSFSS_Min

Gaudi::Property<int> xAODLeptonPairFilter::m_nSFSS_Min {this,"NSFSS_Min",-1,"Number of Same Flavor Same Sign Min"}
private

Definition at line 46 of file xAODLeptonPairFilter.h.

46{this,"NSFSS_Min",-1,"Number of Same Flavor Same Sign Min"};

◆ m_onlyMassiveParents

Gaudi::Property<bool> xAODLeptonPairFilter::m_onlyMassiveParents {this,"OnlyMassiveParents",false,"only count leptons with a parent with M>20 GeV"}
private

Definition at line 70 of file xAODLeptonPairFilter.h.

70{this,"OnlyMassiveParents",false,"only count leptons with a parent with M>20 GeV"};

◆ m_Ptmin

Gaudi::Property<double> xAODLeptonPairFilter::m_Ptmin {this,"Ptcut",10000.,"min Pt"}
private

Definition at line 74 of file xAODLeptonPairFilter.h.

74{this,"Ptcut",10000.,"min Pt"};

◆ m_truthPartContKey

SG::ReadHandleKey<xAOD::TruthParticleContainer> xAODLeptonPairFilter::m_truthPartContKey {this, "TruthParticleContainerKey", "TruthGen"}
private

Definition at line 34 of file xAODLeptonPairFilter.h.

34{this, "TruthParticleContainerKey", "TruthGen"};

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


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