 |
ATLAS Offline Software
|
#include <QCDTruthJetFilter.h>
|
| QCDTruthJetFilter (const std::string &name, ISvcLocator *pSvcLocator) |
|
StatusCode | filterInitialize () |
|
StatusCode | filterFinalize () |
|
StatusCode | filterEvent () |
|
virtual StatusCode | sysInitialize () override |
| Override sysInitialize. More...
|
|
virtual const DataObjIDColl & | extraOutputDeps () const override |
| Return the list of extra output dependencies. 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 | 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, V, H > &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 |
|
|
CLHEP::HepRandomEngine * | getRandomEngine (const std::string &streamName, const EventContext &ctx) const |
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKey> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyArrayType &) |
| specialization for handling Gaudi::Property<SG::VarHandleKeyArray> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &hndl, const SG::VarHandleType &) |
| specialization for handling Gaudi::Property<SG::VarHandleBase> More...
|
|
Gaudi::Details::PropertyBase & | declareGaudiProperty (Gaudi::Property< T, V, H > &t, const SG::NotHandleType &) |
| specialization for handling everything that's not a Gaudi::Property<SG::VarHandleKey> or a <SG::VarHandleKeyArray> More...
|
|
|
double | m_minPtCut {-1.*Gaudi::Units::GeV} |
|
Gaudi::Property< double > | m_MinPt {this, "MinPt", -1.*Gaudi::Units::GeV} |
| Min pT for the truth jets. More...
|
|
double | m_maxPtCut {7000.*Gaudi::Units::GeV} |
|
Gaudi::Property< double > | m_MaxPt {this, "MaxPt", 7000*Gaudi::Units::GeV} |
| Max pT for the truth jets. More...
|
|
Gaudi::Property< double > | m_MinEta {this, "MinEta", s_startMinEta} |
| Min eta for the truth jets. More...
|
|
double | m_minEtaCut {s_startMinEta} |
|
Gaudi::Property< double > | m_MaxEta {this, "MaxEta", 999.0} |
| Max eta for the truth jets. More...
|
|
double | m_maxEtaCut {999.0} |
|
Gaudi::Property< double > | m_MinPhi {this, "MinPhi", -999.0} |
| Min phi for the lead truth jet. More...
|
|
Gaudi::Property< double > | m_MaxPhi {this, "MaxPhi", 999.0} |
|
Gaudi::Property< bool > | m_SymEta {this, "SymEta", false} |
| Use symmetric cut for min eta? (Default false for p-Pb run filters) More...
|
|
Gaudi::Property< bool > | m_doShape {this, "DoShape", false} |
| Attempt to flatten the pT distribution. More...
|
|
SG::ReadHandleKey< xAOD::JetContainer > | m_TruthJetContainerName {this, "TruthJetContainer", "AntiKt4TruthWZJets"} |
|
ServiceHandle< IAthRNGSvc > | m_rndmSvc {this, "RndmSvc", "AthRNGSvc"} |
|
long | m_total {0} |
| Total number of events tested. More...
|
|
long | m_passed {0} |
| Number of events passing all cuts. More...
|
|
long | m_ptfailed {0} |
| Number of events failing the pT cuts. More...
|
|
double | m_norm {1.} |
| Normalization for weights. More...
|
|
double | m_high {1.} |
| High-side function level. More...
|
|
DataObjIDColl | m_extendedExtraObjects |
|
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 |
|
Definition at line 19 of file QCDTruthJetFilter.h.
◆ StoreGateSvc_t
◆ QCDTruthJetFilter()
QCDTruthJetFilter::QCDTruthJetFilter |
( |
const std::string & |
name, |
|
|
ISvcLocator * |
pSvcLocator |
|
) |
| |
◆ ATLAS_NOT_CONST_THREAD_SAFE() [1/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 76 of file GenBase.h.
78 ATH_MSG_ERROR(
"McEventCollection is empty during first event access");
79 return *(
events()->begin());
◆ ATLAS_NOT_CONST_THREAD_SAFE() [2/2]
◆ cmTomm()
void GenBase::cmTomm |
( |
HepMC::GenEvent * |
evt | ) |
|
|
protectedinherited |
Scale event lengths by x 10.
Definition at line 78 of file GenBase.cxx.
79 for (HepMC::GenEvent::vertex_iterator vtx =
evt->vertices_begin(); vtx !=
evt->vertices_end(); ++vtx) {
80 const HepMC::FourVector fv((*vtx)->position().x() * 10,
81 (*vtx)->position().y() * 10,
82 (*vtx)->position().z() * 10,
83 (*vtx)->position().t() * 10);
84 (*vtx)->set_position(fv);
◆ declareGaudiProperty() [1/4]
specialization for handling Gaudi::Property<SG::VarHandleKeyArray>
Definition at line 170 of file AthCommonDataStore.h.
175 hndl.documentation());
◆ declareGaudiProperty() [2/4]
specialization for handling Gaudi::Property<SG::VarHandleKey>
Definition at line 156 of file AthCommonDataStore.h.
161 hndl.documentation());
◆ declareGaudiProperty() [3/4]
specialization for handling Gaudi::Property<SG::VarHandleBase>
Definition at line 184 of file AthCommonDataStore.h.
189 hndl.documentation());
◆ declareGaudiProperty() [4/4]
◆ declareProperty() [1/6]
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.
250 this->declare(hndl.
vhKey());
251 hndl.
vhKey().setOwner(
this);
253 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [2/6]
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.
229 return PBASE::declareProperty(
name,hndl,
doc);
◆ declareProperty() [3/6]
◆ declareProperty() [4/6]
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.
338 return PBASE::declareProperty(
name, property,
doc);
◆ declareProperty() [5/6]
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.
◆ declareProperty() [6/6]
◆ detStore()
◆ event_const()
const HepMC::GenEvent* GenBase::event_const |
( |
| ) |
const |
|
inlineinherited |
Access the current signal event (const)
Definition at line 83 of file GenBase.h.
85 ATH_MSG_ERROR(
"Const McEventCollection is empty during first event access");
◆ events_const() [1/2]
◆ events_const() [2/2]
◆ evtStore() [1/2]
◆ evtStore() [2/2]
◆ execute()
StatusCode GenFilter::execute |
( |
| ) |
|
|
virtualinherited |
- Todo:
- Probably the filter should only look at the first event... right?
Reimplemented from GenBase.
Definition at line 29 of file GenFilter.cxx.
32 return StatusCode::FAILURE;
35 ATH_MSG_WARNING(
"More than one event in current McEventCollection -- which is valid?");
51 sc = StatusCode::FAILURE;
◆ extraDeps_update_handler()
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 & AthAlgorithm::extraOutputDeps |
( |
| ) |
const |
|
overridevirtualinherited |
Return the list of extra output dependencies.
This list is extended to include symlinks implied by inheritance relations.
Definition at line 50 of file AthAlgorithm.cxx.
57 return Algorithm::extraOutputDeps();
◆ filterEvent()
StatusCode QCDTruthJetFilter::filterEvent |
( |
| ) |
|
|
virtual |
Implements GenFilter.
Definition at line 58 of file QCDTruthJetFilter.cxx.
62 const EventContext& ctx = Gaudi::Hive::currentContext();
65 ATH_MSG_WARNING(
"Failed to retrieve random number engine QCDTruthJetFilter");
66 setFilterPassed(
false);
67 return StatusCode::SUCCESS;
73 if (!truthjetTES.isValid()) {
76 return StatusCode::SUCCESS;
80 double pt_lead = -1, phi_lead = 0;
82 if (!(*it_truth))
continue;
85 if (pt_lead < (*it_truth)->pt()) {
86 pt_lead = (*it_truth)->pt();
87 phi_lead = (*it_truth)->phi();
96 setFilterPassed(
false);
98 return StatusCode::SUCCESS;
103 setFilterPassed(
false);
105 if (phi_lead < m_MinPhi || phi_lead >
m_MaxPhi) {
107 return StatusCode::SUCCESS;
110 setFilterPassed(
true);
118 double rnd = rndm->flat();
120 setFilterPassed(
false);
121 ATH_MSG_DEBUG(
"Event failed weighting cut. Weight is " <<
w <<
" for pt_lead of " << pt_lead <<
" high end is " <<
m_high <<
" rnd is " << rnd);
122 return StatusCode::SUCCESS;
127 setFilterPassed(
true);
132 ATH_MSG_DEBUG(
"Event passed. Will mod event weights by " <<
w*
m_norm <<
" for pt_lead of " << pt_lead <<
" norm " <<
m_norm <<
" w " <<
w <<
" high " <<
m_high <<
" rnd " << rnd);
135 for (
unsigned int i=0;
i<mec->
size();++
i) {
136 if ( !(*mec)[
i] )
continue;
137 orig = (*mec)[
i]->weights().
size()>0?(*mec)[
i]->weights()[0]:1.;
138 if ((*mec)[
i]->weights().size()>0) (*mec)[
i]->weights()[0] = orig*
w*
m_norm;
139 else (*mec)[
i]->weights().push_back(
w*
m_norm*orig );
142 (*mec)[
i]->add_attribute(
"filterWeight", std::make_shared<HepMC3::DoubleAttribute>(
w*
m_norm));
146 return StatusCode::SUCCESS;
◆ filterFinalize()
StatusCode QCDTruthJetFilter::filterFinalize |
( |
| ) |
|
|
virtual |
◆ filterInitialize()
StatusCode QCDTruthJetFilter::filterInitialize |
( |
| ) |
|
|
virtual |
◆ finalize()
StatusCode GenFilter::finalize |
( |
| ) |
|
|
inherited |
Definition at line 56 of file GenFilter.cxx.
59 return StatusCode::SUCCESS;
◆ fitFn()
static double QCDTruthJetFilter::fitFn |
( |
const double |
x | ) |
|
|
inlinestatic |
- Todo:
- Move to an inline in the cc
Definition at line 68 of file QCDTruthJetFilter.h.
69 static const double p[7] =
70 { 1000000, -5.5, 0.0001, 0.00000012, 1./4150., 23, 0.4 };
◆ fitFnR()
static double QCDTruthJetFilter::fitFnR |
( |
const double * |
x, |
|
|
const double * |
|
|
) |
| |
|
inlinestatic |
- Todo:
- Move to an inline in the cc
Definition at line 61 of file QCDTruthJetFilter.h.
62 static const double p[7] =
63 { 1000000, -5.5, 0.0001, 0.00000012, 1./4150., 23, 0.4 };
◆ getRandomEngine()
CLHEP::HepRandomEngine * QCDTruthJetFilter::getRandomEngine |
( |
const std::string & |
streamName, |
|
|
const EventContext & |
ctx |
|
) |
| const |
|
private |
◆ GeVToMeV()
void GenBase::GeVToMeV |
( |
HepMC::GenEvent * |
evt | ) |
|
|
protectedinherited |
Scale event energies/momenta by x 1000.
- Todo:
- Add HepMC units awareness and do it differently when HepMC provides this functionality directly (and reference-based FourVector accessors)
Definition at line 58 of file GenBase.cxx.
59 for (HepMC::GenEvent::particle_iterator
p =
evt->particles_begin();
p !=
evt->particles_end(); ++
p) {
60 const HepMC::FourVector fv((*p)->momentum().px() * 1000,
61 (*p)->momentum().py() * 1000,
62 (*p)->momentum().pz() * 1000,
63 (*p)->momentum().e() * 1000);
64 (*p)->set_momentum(fv);
65 (*p)->set_generated_mass(1000 * (*p)->generated_mass());
◆ initialize()
StatusCode GenFilter::initialize |
( |
| ) |
|
|
virtualinherited |
◆ inputHandles()
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.
◆ MeVToGeV()
void GenBase::MeVToGeV |
( |
HepMC::GenEvent * |
evt | ) |
|
|
protectedinherited |
Scale event energies/momenta by x 1/1000.
Definition at line 68 of file GenBase.cxx.
69 for (HepMC::GenEvent::particle_iterator
p =
evt->particles_begin();
p !=
evt->particles_end(); ++
p) {
70 const HepMC::FourVector fv((*p)->momentum().px() / 1000,
71 (*p)->momentum().py() / 1000,
72 (*p)->momentum().pz() / 1000,
73 (*p)->momentum().e() / 1000);
74 (*p)->set_momentum(fv);
75 (*p)->set_generated_mass((*p)->generated_mass() / 1000);
◆ mmTocm()
void GenBase::mmTocm |
( |
HepMC::GenEvent * |
evt | ) |
|
|
protectedinherited |
Scale event lengths by x 1/10.
Definition at line 87 of file GenBase.cxx.
88 for (HepMC::GenEvent::vertex_iterator vtx =
evt->vertices_begin(); vtx !=
evt->vertices_end(); ++vtx) {
89 const HepMC::FourVector fv((*vtx)->position().x() / 10,
90 (*vtx)->position().y() / 10,
91 (*vtx)->position().z() / 10,
92 (*vtx)->position().t() / 10);
93 (*vtx)->set_position(fv);
◆ msg() [1/2]
◆ msg() [2/2]
◆ msgLvl()
◆ nFailed()
int GenFilter::nFailed |
( |
| ) |
const |
|
inlineinherited |
◆ nNeeded()
int GenFilter::nNeeded |
( |
| ) |
const |
|
inlineinherited |
◆ nPassed()
int GenFilter::nPassed |
( |
| ) |
const |
|
inlineinherited |
◆ outputHandles()
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.
◆ particleData()
const HepPDT::ParticleData* GenBase::particleData |
( |
int |
pid | ) |
const |
|
inlineinherited |
Access an element in the particle data table.
Definition at line 126 of file GenBase.h.
◆ particleTable()
const HepPDT::ParticleDataTable& GenBase::particleTable |
( |
| ) |
const |
|
inlineinherited |
Get a particle data table.
Definition at line 118 of file GenBase.h.
◆ partPropSvc()
Access the particle property service.
Definition at line 113 of file GenBase.h.
◆ pdt()
const HepPDT::ParticleDataTable& GenBase::pdt |
( |
| ) |
const |
|
inlineinherited |
Shorter alias to get a particle data table.
Definition at line 123 of file GenBase.h.
◆ renounce()
◆ renounceArray()
◆ sysInitialize()
StatusCode AthAlgorithm::sysInitialize |
( |
| ) |
|
|
overridevirtualinherited |
◆ sysStart()
Handle START transition.
We override this in order to make sure that conditions handle keys can cache a pointer to the conditions container.
◆ updateVHKA()
◆ m_detStore
◆ m_doShape
Gaudi::Property<bool> QCDTruthJetFilter::m_doShape {this, "DoShape", false} |
|
private |
◆ m_evtStore
◆ m_extendedExtraObjects
DataObjIDColl AthAlgorithm::m_extendedExtraObjects |
|
privateinherited |
◆ m_high
double QCDTruthJetFilter::m_high {1.} |
|
private |
◆ m_MaxEta
Gaudi::Property<double> QCDTruthJetFilter::m_MaxEta {this, "MaxEta", 999.0} |
|
private |
◆ m_maxEtaCut
double QCDTruthJetFilter::m_maxEtaCut {999.0} |
|
private |
◆ m_MaxPhi
Gaudi::Property<double> QCDTruthJetFilter::m_MaxPhi {this, "MaxPhi", 999.0} |
|
private |
◆ m_MaxPt
Gaudi::Property<double> QCDTruthJetFilter::m_MaxPt {this, "MaxPt", 7000*Gaudi::Units::GeV} |
|
private |
◆ m_maxPtCut
◆ m_mcEventKey
std::string GenBase::m_mcEventKey {} |
|
protectedinherited |
StoreGate key for the MC event collection (defaults to GEN_EVENT)
Definition at line 137 of file GenBase.h.
◆ m_mcevents_const
Const handle to the MC event collection.
Definition at line 163 of file GenBase.h.
◆ m_MinEta
Gaudi::Property<double> QCDTruthJetFilter::m_MinEta {this, "MinEta", s_startMinEta} |
|
private |
◆ m_minEtaCut
◆ m_MinPhi
Gaudi::Property<double> QCDTruthJetFilter::m_MinPhi {this, "MinPhi", -999.0} |
|
private |
◆ m_MinPt
Gaudi::Property<double> QCDTruthJetFilter::m_MinPt {this, "MinPt", -1.*Gaudi::Units::GeV} |
|
private |
◆ m_minPtCut
◆ 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 139 of file GenBase.h.
◆ m_nFail
◆ m_nNeeded
◆ m_norm
double QCDTruthJetFilter::m_norm {1.} |
|
private |
◆ m_nPass
◆ m_passed
long QCDTruthJetFilter::m_passed {0} |
|
private |
◆ m_ppSvc
ServiceHandle<IPartPropSvc> GenBase::m_ppSvc {this, "PartPropSvc", "PartPropSvc"} |
|
privateinherited |
Handle on the particle property service.
Definition at line 160 of file GenBase.h.
◆ m_ptfailed
long QCDTruthJetFilter::m_ptfailed {0} |
|
private |
◆ m_rndmSvc
◆ m_SymEta
Gaudi::Property<bool> QCDTruthJetFilter::m_SymEta {this, "SymEta", false} |
|
private |
◆ m_total
long QCDTruthJetFilter::m_total {0} |
|
private |
◆ m_TruthJetContainerName
◆ m_varHandleArraysDeclared
◆ m_vhka
◆ s_startMinEta
constexpr double QCDTruthJetFilter::s_startMinEta {-10} |
|
staticconstexprprivate |
The documentation for this class was generated from the following files:
def retrieve(aClass, aKey=None)
void setSeed(const std::string &algName, const EventContext &ctx)
Set the random seed using a string (e.g.
Gaudi::Property< double > m_MaxPt
Max pT for the truth jets.
Const iterator class for DataVector/DataList.
Gaudi::Property< double > m_MinEta
Min eta for the truth jets.
static double fitFn(const double x)
const_pointer_type cptr()
Dereference the pointer.
const McEventCollection * events_const() const
Access the current event's McEventCollection (const)
std::string m_mcEventKey
StoreGate key for the MC event collection (defaults to GEN_EVENT)
StoreGateSvc_t m_evtStore
Pointer to StoreGate (event store by default)
std::vector< SG::VarHandleKeyArray * > m_vhka
virtual StatusCode filterEvent()=0
double m_high
High-side function level.
const std::string & key() const
Return the StoreGate ID for the referenced object.
ServiceHandle< IAthRNGSvc > m_rndmSvc
virtual StatusCode filterInitialize()
virtual void setOwner(IDataHandleHolder *o)=0
SG::ReadCondHandle< T > makeHandle(const SG::ReadCondHandleKey< T > &key, const EventContext &ctx=Gaudi::Hive::currentContext())
Gaudi::Property< double > m_MinPhi
Min phi for the lead truth jet.
Gaudi::Details::PropertyBase & declareGaudiProperty(Gaudi::Property< T, V, H > &hndl, const SG::VarHandleKeyType &)
specialization for handling Gaudi::Property<SG::VarHandleKey>
ServiceHandle< IPartPropSvc > m_ppSvc
Handle on the particle property service.
ServiceHandle< StoreGateSvc > & evtStore()
The standard StoreGateSvc (event store) Returns (kind of) a pointer to the StoreGateSvc.
Gaudi::Property< double > m_MaxPhi
CLHEP::HepRandomEngine * getRandomEngine(const std::string &streamName, const EventContext &ctx) const
virtual StatusCode sysInitialize() override
Override sysInitialize.
long m_total
Total number of events tested.
virtual std::vector< Gaudi::DataHandle * > outputHandles() const override
Return this algorithm's output handles.
::StatusCode StatusCode
StatusCode definition for legacy code.
Gaudi::Details::PropertyBase & declareProperty(Gaudi::Property< T, V, H > &t)
SG::ReadHandleKey< McEventCollection > m_mcevents_const
Const handle to the MC event collection.
Gaudi::Property< double > m_MaxEta
Max eta for the truth jets.
#define CHECK(...)
Evaluate an expression and check for errors.
StoreGateSvc_t m_detStore
Pointer to StoreGate (detector store by default)
StatusCode initialize(bool used=true)
If this object is used as a property, then this should be called during the initialize phase.
This defines the McEventCollection, which is really just an ObjectVector of McEvent objects.
virtual bool isValid() override final
Can the handle be successfully dereferenced?
virtual void renounce()=0
const HepPDT::ParticleDataTable & particleTable() const
Get a particle data table.
std::conditional< std::is_base_of< SG::VarHandleKeyArray, T >::value, VarHandleKeyArrayType, type2 >::type type
A wrapper class for event-slot-local random engines.
static constexpr double s_startMinEta
virtual StatusCode filterFinalize()
CLHEP::HepRandomEngine * getEngine(const EventContext &ctx) const
Retrieve the random engine corresponding to the provided EventContext.
long m_ptfailed
Number of events failing the pT cuts.
Gaudi::Property< double > m_MinPt
Min pT for the truth jets.
Gaudi::Property< bool > m_doShape
Attempt to flatten the pT distribution.
DataObjIDColl m_extendedExtraObjects
bool empty() const noexcept
long m_passed
Number of events passing all cuts.
double m_norm
Normalization for weights.
#define ATH_MSG_WARNING(x)
const HepPDT::ParticleDataTable & pdt() const
Shorter alias to get a particle data table.
SG::ReadHandleKey< xAOD::JetContainer > m_TruthJetContainerName
SG::VarHandleKey & vhKey()
Return a non-const reference to the HandleKey.
filterPassed
now, add a sequencer: it will only execute seqalg1 and 2 and never 3
GenFilter(const std::string &name, ISvcLocator *pSvcLocator)
constexpr int pow(int base, int exp) noexcept
size_type size() const noexcept
Returns the number of elements in the collection.
virtual StatusCode initialize() override
const_iterator begin() const noexcept
Return a const_iterator pointing at the beginning of the collection.
Gaudi::Property< bool > m_SymEta
Use symmetric cut for min eta? (Default false for p-Pb run filters)